fsk_pe.m

一个基于matlab的fsk仿真

for(snr=0:10)
    m=0;
    for(l=1:100000)
        a=rand(1,1);
        a=2*a-1;
        a=ceil(a);
        if(a==0)
            b=zeros(1,1001);
        else
            b=ones(1,1001);
        end;
        c=b;
        n=[0:0.001:1];
        d=cos((4+2*c)*pi.*n);
        y=awgn(d,snr-27,'measured');
        z1=y.*cos(6*pi*n)-y.*cos(4*pi*n);
        z=sum(z1)*0.001;
        if(z<0)
            o=0;
        else
            o=1;
        end;
        if(o~=a)
            m=m+1;
        end;
    end;
    pe(snr+1)=m/100000;
end;
snr=[0:10];
plot(snr,pe,'V-')
hold on;
plot(snr,0.5*erfc(sqrt((10.^(snr/10))/2)),'.r-')