接下来我会发一段代码，请增加代码注释并返回代码 %% BPSK调制 fs = 2000; ts = 0:0.00001:3.5-0.00001; % 为了使信号看起来更光滑，作图时采样频率为100kHz s_b = rectpulse(s, 1000); % 将冲激信号补成矩形信号 s_bpsk = (1-2.*s_b).*cos(2*pi*fs*ts); % 扩频后信号BPSK调制时域波形，(1-2.*s_b)是1，-1序列 %% 画出扩频前和扩频后BPSK信号时域波形 figure(3); subplot(2, 1, 2); plot(ts, s_bpsk); xlabel('s'); axis([0.055, 0.085, -1.2, 1.2]) title('扩频后bpsk信号时域波形'); subplot(2, 1, 1); s_bb = rectpulse(x, 7000); s_bpskb = (1-2.*s_bb).*cos(2*pi*fs*ts); % 无扩频信号BPSK调制时域波形 plot(ts, s_bpskb); xlabel('s'); axis([0.055, 0.085, -1.2, 1.2]); title('扩频前bpsk信号时域波形'); %% 画出扩频前和扩频后BPSK调制信号的频谱图 figure(4); N = 400000; ybb = fft(s_bpskb, N); % 无扩频信号BPSK调制频谱 magb = abs(ybb); fbb = (1:N/2) * 100000 / N; subplot(2, 1, 1); plot(fbb, magb(1:N/2) * 2 / N); axis([1700, 2300, 0, 0.8]); title('扩频前调制信号频谱图'); xlabel('Hz'); yb = fft(s_bpsk, N); % 扩频信号BPSK调制频谱 mag = abs(yb); fb = (1:N/2) * 100000 / N; subplot(2, 1, 2); plot(fb, mag(1:N/2) * 2 / N); axis([1700, 2300, 0, 0.8]); title('扩频后调制信号频谱图'); xlabel('Hz');

%% BPSK调制
fs = 2000; % 采样频率
ts = 0:0.00001:3.5-0.00001; % 时间序列

% 为了使信号看起来更光滑，作图时采样频率为100kHz
s_b = rectpulse(s, 1000); % 将冲激信号补成矩形信号

s_bpsk = (1-2.*s_b).*cos(2*pi*fs*ts); % 扩频后信号BPSK调制时域波形，(1-2.*s_b)是1，-1序列

%% 画出扩频前和扩频后BPSK信号时域波形
figure(3);
subplot(2, 1, 2);
plot(ts, s_bpsk);
xlabel('时间 (s)');
axis([0.055, 0.085, -1.2, 1.2]);
title('扩频后BPSK信号时域波形');

subplot(2, 1, 1);
s_bb = rectpulse(x, 7000); % 将冲激信号补成矩形信号
s_bpskb = (1-2.*s_bb).*cos(2*pi*fs*ts); % 无扩频信号BPSK调制时域波形
plot(ts, s_bpskb);
xlabel('时间 (s)');
axis([0.055, 0.085, -1.2, 1.2]);
title('扩频前BPSK信号时域波形');

%% 画出扩频前和扩频后BPSK调制信号的频谱图
figure(4);
N = 400000; % 傅里叶变换点数
ybb = fft(s_bpskb, N); % 无扩频信号BPSK调制频谱
magb = abs(ybb); % 幅度谱
fbb = (1:N/2) * 100000 / N; % 频率序列
subplot(2, 1, 1);
plot(fbb, magb(1:N/2) * 2 / N);
axis([1700, 2300, 0, 0.8]);
title('扩频前调制信号频谱图');
xlabel('频率 (Hz)');

yb = fft(s_bpsk, N); % 扩频信号BPSK调制频谱
mag = abs(yb); % 幅度谱
fb = (1:N/2) * 100000 / N; % 频率序列
subplot(2, 1, 2);
plot(fb, mag(1:N/2) * 2 / N);
axis([1700, 2300, 0, 0.8]);
title('扩频后调制信号频谱图');
xlabel('频率 (Hz)');