function [p]=ss_Pc94(snr_in_dB,Lc,A,w0)
%計(jì)算在當(dāng)前信噪比,碼片數(shù),干擾正弦幅度,頻率下的誤碼率,不用產(chǎn)生pn序列
% snr_in_dB=10;
% Lc=20;%碼片數(shù)
% A=3;
% w0=1;
snr=10^(snr_in_dB/10);
sgma=1;   %噪聲的標(biāo)準(zhǔn)差
Eb=2*sgma^2*snr;%根據(jù)信噪比,噪聲能量,求出信號(hào)能量,且 N0/2=sgma^2
E_chip=Eb/Lc;%每個(gè)碼元能量

N=10000;
num_of_err=0;%錯(cuò)誤數(shù)計(jì)數(shù)器
for i=1:N
    %生成下一個(gè)bit
    temp=rand;
    if (temp<.5)
        data=-1;
    else
        data=1;
    end
    
    %重復(fù)Lc次,
    for j=1:Lc
        repeated_data(1,j)=data;
    end;
    %生成pn序列
    for j=1:Lc
        temp=rand;
        if(temp<0.5)
            pn_seq(j)=-1;
        else
            pn_seq(j)=1;
        end
    end;
    %得到傳輸信號(hào)
    trans_sig=sqrt(E_chip)*repeated_data.*pn_seq;%假設(shè)碼片間隔Tc=1
    %AWGN的方差為sgma^2
    noise=sgma*randn(1,Lc);
    %干擾
    n=(i-1)*Lc+1:i*Lc;
    interference=A*sin(w0*n);%CW干擾為1*Lc的向量
    %接收信號(hào)
    rec_sig=trans_sig+noise+interference;
    %從接收信號(hào)中判斷信號(hào)
    temp=rec_sig.*pn_seq;
    decision_variable=sum(temp);
    %作出判決,采用硬判決
    if(decision_variable<0);
        decision=-1;
    else decision=1;
    end
    %如果判決錯(cuò)誤者計(jì)數(shù)器加1
    if (decision~=data)
        num_of_err=num_of_err+1;
    end;
end
%得到錯(cuò)誤概率
p=num_of_err/N;