BPSK调制及功率谱，频谱。

%参数设定
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）;