資源描述:
《Matlab仿真DF和AF源程序》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、版權(quán)所有,如轉(zhuǎn)載,請(qǐng)注明出處。說(shuō)明:?????該程序是使用Matlab仿真協(xié)同通信的DF(解碼轉(zhuǎn)發(fā))基本性能,基本模型S-R-D三個(gè)節(jié)點(diǎn),接收端使用MRC(最大比合并)。固定DF由于SR信道的錯(cuò)誤,其性能受到限制,在大SNR情況下不能獲得有效分集。而當(dāng)SR距離足夠近的時(shí)候,可能獲得一定的分集增益。????該程序也可以作為學(xué)習(xí)Matlab通信仿真的入門實(shí)例。使用的基本函數(shù):???隨機(jī)數(shù)函數(shù):randrandn??文件操作:fopenfprintffclose??時(shí)間函數(shù):datestr(now)????畫圖函數(shù)
2、:semilogy??legend??axis?????通信基本概念和函數(shù):??瑞利信道及其產(chǎn)生??高斯噪聲的產(chǎn)生???蒙特卡羅仿真???MRC最大比合并?????程序一:DF_SNR_SR_distance.m主程序%writtenbyxiyong826025/08/2009%thiscodeisforfixedDFperformanceevaluation,comparethediffirentSRdistance%該仿真用來(lái)評(píng)估不同的SR距離下的固定DF轉(zhuǎn)發(fā)性能%forS-R-Dthreenodes%t
3、heSDandRDchannelisRayleighfading,andtheaverageSNRisvaried.clearall;%thetimesofMonteCarlo?%蒙特卡羅仿真的次數(shù)Monte_MAX=10^(7);%thelengthofdataframeM2=10^(0);MIN_SNR=0;MAX_SNR=30;INTERVAL=2;%功率劃分因子,在協(xié)同情況下,為了保證總的功率一定,每個(gè)節(jié)點(diǎn)使用1/2的功率發(fā)送POW_DIV=1/2;%openanewrecordfile%打開一個(gè)新的
4、文件,如果文件已經(jīng)存在,那么清除以前的內(nèi)容fid=fopen('record.txt','w');fprintf(fid,'%%%s',datestr(now));fclose(fid);%settheSDdistanceandRDdistance%設(shè)置SD距離和RD的距離,歸一化為1sd_distance=1;rd_distance=1;%變化不同的SR距離sr_distances=[0.1,0.5,0.7,0.9];%使用循環(huán)的方式仿真不同的SR距離性能forsr_INDEX=1:2???????%通
5、過(guò)sr_INDEX變量獲得真正的SR距離???sr_distance=sr_distances(sr_INDEX);????????%snrcount這個(gè)變量用來(lái)對(duì)下面的SNR循環(huán)進(jìn)行計(jì)數(shù)???snrcount=0;???????%變化SNR循環(huán)???forSNR=MIN_SNR:INTERVAL:MAX_SNR???????sig=10^(SNR/10);???????????????%justdisplaysomethingonthescreen.????????%僅僅是在計(jì)算機(jī)上顯示一些信息,因?yàn)榉抡嫱?/p>
6、往很漫長(zhǎng),只是想知道仿真到哪里了^.^???????datestr(now)???????M2???????SNR??????????????snrcount=snrcount+1;???????????????%用來(lái)對(duì)錯(cuò)誤的bit計(jì)數(shù)???????err_num_dsd=0;???????err_num_coop=0;??????????????fortries=0:Monte_MAX???????????????????????%flagindicateifthecooperateisused,if1,u
7、secooperation,elsedirecttranisperferred.???????????%一個(gè)標(biāo)志,表明是否轉(zhuǎn)發(fā)???????????tx_coop=1;???????????????????????%generatethesourceBPSKsignal,0,1????????????%產(chǎn)生一個(gè)長(zhǎng)度為M2的0,1隨機(jī)數(shù)???????????X1=rand(1,M2)>0.5;????????????%generatethesourceBPSKsignal,+1,-1????????????%產(chǎn)
8、生一個(gè)長(zhǎng)度為M2的BPSK隨機(jī)信號(hào)???????????Xs=X1*2-1;????????????%generatethesource-relaychannel,thechannelkeepconstantinaframe????????????%產(chǎn)生SR信道,假設(shè)信道在一幀內(nèi)保持不變???????????CH_sr=xy_RayleighCH(1)/(sr_distance)^2;???????