%參數設定
Fs?=?2000;????????????%?采樣頻率
T?=?1/Fs;?????????????%?采樣周期
L?=1000;?????????????%?信號長度
t?=?（0:L-1）*T;????????%?時間向量


%產生一個方波，頻率為16，占空比為50?
f=4;
baseband=?square（2*pi*f*t50）;
figure;
subplot（311）;
plot（baseband）;


%產生一個余弦載波信號，頻率為fc
fc=500;
u=4*sin（2*pi*fc*t）;
subplot（312）;
plot（u）;

%BPSK調制
f1=baseband.*u;
subplot（313）;
plot（f1）;


%基帶信號頻譜
baseband_spectrum?=?fft（baseband）;
%將基帶信號的頻譜的橫軸換算成頻率
P_b?=?abs（baseband_spectrum/L）;
P_B?=?P_b（1:L/2+1）;
P_B（2:end-1）?=?2*P_B（2:end-1）;
f?=?Fs*（0:（L/2））/L;
figure;
subplot（211）;