shiyan1.m

用MATLAB语言编程设计 本文介绍了一个用MATLAB语言编写程序

clear all

T=2;
j=1;
delt_t=0.1;
Wc=4*pi;
sigma=1;
t=0:delt_t:1;
nCount1=zeros(11,1);
nCount2=0;
sum1=0;
sum2=0;
for SNR=10.^(0:0.1:1)
    A=sqrt(4*SNR*delt_t/T);
    s0=A*sin(Wc*t);
    s1=A*sin(Wc*t+pi);
    while nCount1(j)<100
        m=rand(1);
        if m<0.5
            r0=s0+randn(1,11);
            sum1=sum(r0.*s0*delt_t);
            sum2=sum(r0.*s1*delt_t);
            if sum1<sum2
                nCount1(j)=nCount1(j)+1;
            end
        else
            r1=s1+randn(1,11);
            sum1=sum(r1.*s0*delt_t);
            sum2=sum(r1.*s1*delt_t);
            if sum1>sum2
                nCount1(j)=nCount1(j)+1;
            end
        end
        nCount2=nCount2+1;  
    end
     nCount1(j)=nCount1(j)/nCount2;
     nCount2=0;
     j=j+1;
end
  delt=1:-0.1:0.1;
  SNR=0:10
  semilogy(SNR,nCount1,'r',SNR,0.5*erfc(sqrt(10.^(SNR/10)))),grid