echo on
Lc=20;   %每比特碼片數目
A1=3;     %第一個正弦干擾信號的幅度
A2=7;%第二個正弦干擾信號的幅度
A3=12;%第三個正弦干擾信號的幅度
A4=0;%第四種情況：無干擾
w0=1;%以弧度表達的正弦干擾信號頻率
SNRindB=0:2:30;
for i=1:length(SNRindB),
    %計算誤碼率
    smld_err_prb1(i)=ss_Pe94(SNRindB(i),Lc,A1,w0);
    smld_err_prb2(i)=ss_Pe94(SNRindB(i),Lc,A2,w0);
    smld_err_prb3(i)=ss_Pe94(SNRindB(i),Lc,A3,w0);
end;
SNRindB4=0:1:8;
for i=1:length(SNRindB4),
    %計算無干擾的情況下的誤碼率
    smld_err_prb4(i)=ss_Pe94(SNRindB(i),Lc,A4,w0);
end;
semilogy (SNRindB,smld_err_prb1,SNRindB,smld_err_prb2,SNRindB,smld_err_prb3);
ylabel('誤碼率Pe');
xlabel('信噪比S/N');
text(15,10^-3,'A=3');
text(20,8*(10^-3),'A=7');
text(25,10^-2,'A=12');