b=1; else b=-1; end
h1=randn; h2=randn; h3=randn; h4=randn; n1=randn; n2=randn; n3=randn; n4=randn;
s1=snr*abs(h1)*b+n1*h1/abs(h1); s2=snr*abs(h2)*b+n2*h2/abs(h2); s3=snr*abs(h3)*b+n3*h3/abs(h3); s4=snr*abs(h4)*b+n4*h4/abs(h4); S=s1+s2+s3+s4; if S>0; z=1; else z=-1; end if z~=b
err=err+1; end end
ber=err/10000
我们把A依次输入1到15，可以得到15个运行结果，利用以下程序可以画出曲线图： x=1:1:15;
y=[0.0743,0.0675,0.0600,0.0484,0.0463,0.0380,0.0345,0.0292,0.0239,0.0190,0.0183,0.0145,0.0116,0.0104,0.0083]; semilogy(x,y)
6根接收天线时
A=1; err=0;
snr=sqrt(10^(A/20)); for k=1:10000 s=rand;
29
if s>0.5 b=1; else b=-1; end
h1=randn; h2=randn; h3=randn; h4=randn; h5=randn; h6=randn; n1=randn; n2=randn; n3=randn; n4=randn; n5=randn; n6=randn;
s1=snr*abs(h1)*b+n1*h1/abs(h1); s2=snr*abs(h2)*b+n2*h2/abs(h2); s3=snr*abs(h3)*b+n3*h3/abs(h3); s4=snr*abs(h4)*b+n4*h4/abs(h4); s5=snr*abs(h5)*b+n5*h5/abs(h5); s6=snr*abs(h6)*b+n6*h6/abs(h6); S=s1+s2+s3+s4+s5+s6; if S>0; z=1; else z=-1; end if z~=b
err=err+1; end end
ber=err/10000
我们把A依次输入1到15，可以得到15个运行结果，利用以下程序可以画出曲线图: x=1:1:15;
y=[0.0360,0.0343,0.0281,0.0219,0.0191,0.0144,0.0140,0.0101,0.0084,0.0068,0.0056,0.0043,0.0030,0.0028,0.0016]; semilogy(x,y)
30
