%參數(shù)設(shè)定
Fs?=?2000;????????????%?采樣頻率
T?=?1/Fs;?????????????%?采樣周期
L?=1000;?????????????%?信號(hào)長(zhǎng)度
t?=?（0:L-1）*T;????????%?時(shí)間向量


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


%產(chǎn)生一個(gè)余弦載波信號(hào)，頻率為fc
fc=500;
u=4*sin（2*pi*fc*t）;
subplot（312）;
plot（u）;

%BPSK調(diào)制
f1=baseband.*u;
subplot（313）;
plot（f1）;


%基帶信號(hào)頻譜
baseband_spectrum?=?fft（baseband）;
%將基帶信號(hào)的頻譜的橫軸換算成頻率
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）;