修改代码 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` 函数绘制了误码率曲线。

希望这可以帮助您！