ask_pe.m

来自「一个用matlab实现的ask仿真程序」· M 代码 · 共 33 行

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=c.*cos(4*pi*n);
        d1=cos(4*pi*n);
        y=awgn(d1,snr-24,'measured')-d1+d;
        z1=y.*cos(4*pi*n);
        z=sum(z1)*0.001;
        if(z<0.25)
            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.-')