temp.m

非常好用

clear,clc
for SNR=0:2:30
    sdata=(floor(rand(1,50)*2)-1);
    a=find(sdata==0);
    sdata(a)=1;
    sdata
    [zz,DopCom]=ChannelAeroT1(j*sdata);
    rdata=imag(zz)>0;
    a=find(sdata==-1);
    sdata(a)=0;
    error=sum(xor(rdata,sdata));
    Ber(SNR/2+1)=error/50;
end
plot((0:2:30)/2+1,Ber);