資源描述:
《基帶傳輸實(shí)驗(yàn)圖形》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、%每種信噪比下的發(fā)送符號數(shù)%信噪比,E/N0%消息數(shù)據(jù)%比特0在發(fā)送消息中的位置%比特0影射為發(fā)送波形sO%比特1在發(fā)送消息中的位置%比特1映射為發(fā)送波形si%總的發(fā)送波形%數(shù)據(jù)轉(zhuǎn)罝,方便接收端處理%通過AWGN信道%與s0相關(guān)%與si相關(guān)%判決正交信號在AWGN信道下的傳輸性能代碼:clearallnsamp=10;sO=ones(1,nsamp);si=[ones(l,nsamp/2)-ones(1,nsamp/2)];nsymbol=100000;EbN0=0:12;msg=randint(l,nsymbol);sOO=zeros(nsymb
2、ol,1);si1=zeros(nsymbol,1);indx=find(msg==0);s00(indx)=l;s00=s00*s0;indxl=find(msg==1);si1(indxl)=1;sil=sll*sl;s=s00+sll;s=$?;forindx=1:length(EbN0)decmsg=zeros(1,nsymbol);r=awgn(s,EbN0(indx)-7);r00=s0*r;rll=sl*r;indxl=fmd(rl1>=r00);decmsg(indxl)=l;[err,ber(indx)]=biterr(msg,d
3、ecmsg);endsemilogy(EbNO,ber;-koEbNO,qfunc(sqrt(l0.八(EbNO/10))));title(’二進(jìn)制正交信號在AWGN信道下的誤比特率性能’)xlabelCEb/N(T);ylabel(’誤比特率Pe')legend仿真結(jié)果7理論結(jié)果閣形Figure1口D運(yùn)S曇Dp_h0.二進(jìn)制正交信號在AWGN信道下的誤比特率性能仿真結(jié)果理論結(jié)果鑭O?2O-3O-4O-5468Eb/NO1012改善代碼可以不考慮發(fā)送信號波形的影響:%發(fā)送符號數(shù)%信噪比%消息數(shù)據(jù)clearallnsymbol=100000;EbN
4、0=0:12;msg=randint(1,nsymbol);E=l;rO=zeros(l,nsymbol);rl=zeros(l,nsymbol);indx=find(msg==O);rO(indx)=E;indxl=find(msg==1);rl(indxl)=E;forindx=l:length(EbN0)%dB轉(zhuǎn)換為線性值%噪聲方差%相關(guān)器的輸出dec=zeros(1,length(msg));snr=10.A(EbN0(indx)/10);sigma=l/(2*snr);rOO=rO+sqrt(sigma)*randn(1,length(m
5、sg));%判決rll=r1+$qrt(sigma)*randn(1,length(m$g));indxl=find(rll>=rOO);dec(indxl)=l;[err,ber(indx)]=biterr(msg,dec);endfiguresemilogy(EbN0,ber,'-ko’,EbN0,qfunc(sqrt(10.A(EbN0/10))));title(’二進(jìn)制正交信號在AWGN信道下的誤比特率性能’)xlabel(.Eb/NO.);ylabel(’誤比特率Pe')legendC仿真結(jié)果理論結(jié)果)圖形??WFigureFileEdi
6、tViewInsertToolsDesktopWindowHelp□GI□圓■0這種修改可以使仿真的運(yùn)行時間變短,效率商用simulink仿真上邊的代碼:仿真吋間設(shè)為100000雙極性信號在AWGN信道下的傳輸性能程序代碼如下clearalln$amp=10;%每個脈沖信號的抽樣點(diǎn)數(shù)sl=-$O;nsymbol=100000;EbN0=0:10;msg=randint(l,nsymbol);sOO=zeros(nsymbol,1);si1=zeros(nsymbol,1);indx=find(msg==0);s00(indx)=l;s00=s00*
7、s0;indxl=find(msg==l);si1(indxl)=1;sil=sll*sl;s=s00+sll;s=s.’;%每種信噪比下的發(fā)送符號數(shù)%信噪比,E/N0%消息數(shù)據(jù)%比特0在發(fā)送消息中的位置%比特0影射為發(fā)送波形sO%比特1在發(fā)送消息中的位置%比特1映射為發(fā)送波形si%總的發(fā)送波形%數(shù)裾轉(zhuǎn)賈,方便接收端處理forindx=1:length(EbN0)decmsg=zeros(1,nsymbol);r=awgn(s,EbN0(indx)-7);r00=s0*r;indxl=find(r00<0);decmsg(indxl)=l;[err
8、,ber(indx)]=biterr(msg,decmsg);end%通過AWGN信道%與s0相關(guān)%判決semilogy(