ss_pe.m

一个函数,可用于构建CDMA信源,有数据产生的PN码的形成,以及噪声设置

function[p]=ss_Pe(snr_in_dB,Lc,A,w0)
%运算得出的误码率
snr=10^(snr_in_dB/10);
sgma=1; %噪声的标准方差设置为固定值
Eb=2*sgma^2*snr; %达到设定信噪比所需要的信号幅度
E_chip=Eb/Lc;%每码片的能量
N=10000;%传送的比特数目
%为减少该程序的运算时间,数据的产生 噪声 干扰 译码和差错计算都一
%起执行,这样做有助于超长运算矢量的计算
num_of_err=0;
for i=1:N   %产生下一数据比特
    temp=rand;
    if(temp<0.5),data=-1;
    else data=1;
    end;
    for j=1:Lc  %将其重复Lc次
        repeated_data(j)=data;
    end;
    for j=1:Lc  %产生比特传输使用的PN序列
        temp=rand;
        if(temp<0.5)pn_seq(j)=-1;
        else pn_seq(j)=1;
        end
    end;
    trans_sig=sqrt(E_chip)*repeated_data.*pn_seq;  %发送信号
    noise=sgma*randn(1,Lc); %方差为Sgma^2的高斯白噪声
    n=(i-1)*Lc+1:i*Lc; %干扰
    interference=A*sin(w0*n);
        rec_sig=trans_sig+noise+interference;  %接收信号
    temp=rec_sig.*pn_seq; %从接收信号中产生判决变量
    decision_variable=sum(temp);
        if (decision_variable<0) decision=-1; %进行判决
        else decision=1;
        end;
        if(decision~=data)
            num_of_err=num_of_err+1;
        end;
    end;
    p=num_of_err/N;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%