figure(3) subplot(2,1,1) plot(t, modulating) xlim([0 0.02]) title('原始基带信号') xlabel('时间') ylabel('幅度') subplot(2,1,2) plot(t, demodulated) title('解调后信号') xlabel('时间') ylabel('幅度') xlim([0 0.02])

这段代码是用 MATLAB 绘制信号调制和解调的波形图。subplot 函数用于将绘图窗口分成多个区域，其中第一个参数表示总的行数，第二个参数表示总的列数，第三个参数表示当前绘图的区域编号。plot 函数用于绘制图像，第一个参数表示横轴数据，第二个参数表示纵轴数据。xlim 函数用于设置横轴的范围，title 函数用于设置图像的标题，xlabel 和 ylabel 函数用于设置横轴和纵轴的标签。

相关问题

%% 产生基带信号 i=10;%基带信号码元数 j=5000; t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即将[0,5]分成5000份 f1=10;%载波1频率 f2=5;%载波2频率 fm=i/5;%基带信号频率 码元数是10，而时域长度是5，也就是一个单位2个码元 a=round(rand(1,i));%产生随机序列 st1=t; for n=1:10 if a(n)<1 for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end figure(1); subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]); %% FSK调制 %基带信号求反 st2=t; for n=1:j if st1(n)==1 st2(n)=0; else st2(n)=1; end end subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]); %载波信号 s1=cos(2*pi*f1*t); s2=cos(2*pi*f2*t); subplot(413),plot(s1); title('载波信号s1'); subplot(414),plot(s2); title('载波信号s2'); %调制 F1=st1.*s1;%加入载波1 F2=st2.*s2;%加入载波2 figure(2); subplot(411); plot(t,F1); title('F1=s1*st1'); subplot(412); plot(t,F2); title('F2=s2*st2'); e_fsk=F1+F2; subplot(413); plot(t,e_fsk); title('2FSK信号');%键控法产生的信号在相邻码元之间相位不一定连续 f = linspace(-1000, 1000, length(e_fsk)); Y = fftshift(abs(fft(e_fsk))); subplot(414); plot(f, Y); title('2FSK频域图'); %% 加入噪声 snr = 10 % 信噪比 fsk_noisy = awgn(e_fsk, snr, 'measured'); % 加入噪声 figure(3) plot(t, fsk_noisy); % 绘制加入噪声后的时域波形 title('受干扰时域波形');

这段代码存在一些小问题，主要是缺少一些变量的定义，例如 `awgn` 函数中的 `snr` 变量需要定义。另外，还需要添加一些注释来解释代码的作用。

以下是添加了注释并修正了 `snr` 变量定义的代码：

%% 产生基带信号
i=10; % 基带信号码元数
j=5000;
t=linspace(0,5,j); % 0-5之间产生5000个点行矢量，即将[0,5]分成5000份
f1=10; % 载波1频率
f2=5; % 载波2频率
fm=i/5; % 基带信号频率，码元数是10，而时域长度是5，也就是一个单位2个码元
a=round(rand(1,i)); % 产生随机序列
st1=t;
for n=1:10
    if a(n)<1
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=0;
        end
    else
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=1;
        end
    end
end

figure(1); 
subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]);

%% FSK调制
% 基带信号求反
st2=t;
for n=1:j
    if st1(n)==1
        st2(n)=0;
    else
        st2(n)=1;
    end
end
subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]);

% 载波信号
s1=cos(2*pi*f1*t);
s2=cos(2*pi*f2*t);
subplot(413),plot(s1); title('载波信号s1');
subplot(414),plot(s2); title('载波信号s2');

% 调制
F1=st1.*s1; % 加入载波1
F2=st2.*s2; % 加入载波2
figure(2); 
subplot(411); plot(t,F1); title('F1=s1*st1');
subplot(412); plot(t,F2); title('F2=s2*st2');
e_fsk=F1+F2;
subplot(413); plot(t,e_fsk); title('2FSK信号'); % 键控法产生的信号在相邻码元之间相位不一定连续
f = linspace(-1000, 1000, length(e_fsk));
Y = fftshift(abs(fft(e_fsk)));
subplot(414); plot(f, Y); title('2FSK频域图');

%% 加入噪声
snr = 10; % 信噪比
fsk_noisy = awgn(e_fsk, snr, 'measured'); % 加入噪声
figure(3); 
plot(t, fsk_noisy); % 绘制加入噪声后的时域波形
title('受干扰时域波形');

希望这可以帮到你！

修改代码 i=10;%基带信号码元数 j=5000; t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即将[0,5]分成5000份 f1=10;%载波1频率 f2=5;%载波2频率 fm=i/5;%基带信号频率 码元数是10，而时域长度是5，也就是一个单位2个码元 a=round(rand(1,i));%产生随机序列 %产生基带信号 st1=t; for n=1:10 if a(n)<1 for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end figure(1); subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]); %基带信号求反 st2=t; for n=1:j if st1(n)==1 st2(n)=0; else st2(n)=1; end end subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]); %载波信号 s1=cos(2*pi*f1*t); s2=cos(2*pi*f2*t); subplot(413),plot(s1); title('载波信号s1'); subplot(414),plot(s2); title('载波信号s2'); %调制 F1=st1.*s1;%加入载波1 F2=st2.*s2;%加入载波2 figure(2); subplot(411); plot(t,F1); title('F1=s1*st1'); subplot(412); plot(t,F2); title('F2=s2*st2'); e_fsk=F1+F2; subplot(413); plot(t,e_fsk); title('2FSK信号');%键控法产生的信号在相邻码元之间相位不一定连续 %% %加噪 SNR = 0:0.1:10; % 信噪比范围 error = zeros(1,length(SNR)); % 误码率数组 for i = 1:length(SNR) fsk_noisy = awgn(e_fsk, snr, 'measured'); % 加入噪声 subplot(414); plot(t, fsk_noisy); % 绘制加入噪声后的时域波形 title('受干扰时域波形');

好的，以下是修改后的代码，主要是添加了计算误码率的部分：

i=10;%基带信号码元数
j=5000;
t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即将[0,5]分成5000份
f1=10;%载波1频率
f2=5;%载波2频率
fm=i/5;%基带信号频率    码元数是10，而时域长度是5，也就是一个单位2个码元
a=round(rand(1,i));%产生随机序列

%产生基带信号
st1=t;
for n=1:10
    if a(n)<1
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=0;
        end
    else
        for m=j/i*(n-1)+1:j/i*n
            st1(m)=1;
        end
    end
end

figure(1);
subplot(411); plot(t,st1); title('基带信号st1'); axis([0,5,-1,2]);

%基带信号求反
st2=t;
for n=1:j
    if st1(n)==1
        st2(n)=0;
    else
        st2(n)=1;
    end
end
subplot(412); plot(t,st2); title('基带信号反码st2'); axis([0,5,-1,2]);

%载波信号
s1=cos(2*pi*f1*t);
s2=cos(2*pi*f2*t);
subplot(413),plot(s1); title('载波信号s1');
subplot(414),plot(s2); title('载波信号s2');

%调制
F1=st1.*s1;%加入载波1
F2=st2.*s2;%加入载波2
figure(2);
subplot(411); plot(t,F1); title('F1=s1*st1');
subplot(412); plot(t,F2); title('F2=s2*st2');
e_fsk=F1+F2;
subplot(413); plot(t,e_fsk); title('2FSK信号');%键控法产生的信号在相邻码元之间相位不一定连续

%加噪
SNR = 0:0.1:10; % 信噪比范围
error = zeros(1,length(SNR)); % 误码率数组
for i = 1:length(SNR)
    fsk_noisy = awgn(e_fsk, SNR(i), 'measured');  % 加入噪声
    subplot(414); plot(t, fsk_noisy);  % 绘制加入噪声后的时域波形
    title('受干扰时域波形');
    %解调
    demod1=fsk_noisy.*s1;
    demod2=fsk_noisy.*s2;
    [b,a]=butter(6,2*fm/j,'low');
    df1=filter(b,a,demod1);%低通滤波
    df2=filter(b,a,demod2);%低通滤波
    z1=df1(j/i/2:j/i:j);
    z2=df2(j/i/2:j/i:j);

    %判决
    for k=1:i
        if z1(k)>z2(k)
            c(k)=1;
        else
            c(k)=0;
        end
    end

    %比较
    error(i) = sum(c~=a); % 统计误码数
    BER(i) = error(i)/i; % 计算误码率
end

%绘制误码率曲线
figure(3);
semilogy(SNR, BER);
xlabel('信噪比(dB)');
ylabel('误码率');
title('2FSK误码率曲线');

其中，我们添加了一个 `SNR` 的范围，然后在循环中计算加噪后的误码率。最后，我们使用 `semilogy` 函数绘制了误码率曲线。

希望这可以帮助您！