資源描述：
《QPSK通過Rayleigh信道多徑衰落的Matlab仿真.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、參照《通信系統(tǒng)仿真原理與無線應(yīng)用》351頁(yè)例14-1在這個(gè)例子里，我們對(duì)有3條固定路徑的AWGN多徑信道中的QPSK系統(tǒng)進(jìn)行BER性能仿真，并與在理想的AWGN信道（沒有多徑）中同樣系統(tǒng)地BER性能進(jìn)行比較……書上有比較詳細(xì)的數(shù)學(xué)推導(dǎo)，不抄了。這個(gè)例子似乎沒有考慮多普勒頻移。待我繼續(xù)學(xué)習(xí)下一個(gè)例子，這個(gè)也沒太看懂。下面是該例子的源程序，P0、P1、P2分別是LOS路徑和兩條延遲瑞利分量的相對(duì)功率級(jí)。當(dāng)p0=0且delay！=0時(shí)為瑞利頻率選擇性衰落，delay==0時(shí)為瑞利平坦衰落。主程序scr
2、iptfile：%兩徑瑞利衰落信道仿真%設(shè)定默認(rèn)參數(shù)NN=256;%傳輸符號(hào)個(gè)數(shù)tb=0.5;%一比特時(shí)間fs=10;%每符號(hào)采樣數(shù)ebn0db=[1:2:15];%設(shè)定Eb/N0%建立QPSK信號(hào)x=random_binary(NN,fs)+i*random_binary(NN,fs);%x為QPSK信號(hào)%輸入功率和延遲p0=0;%視距LOS分量p1=20;%第一路徑分量p2=1;%第二路徑分量delay=1;%按照每符號(hào)采樣數(shù)決定的延遲delay0=0;delay1=0;delay2=del
3、ay;%設(shè)定復(fù)高斯（瑞利）衰減gain1=sqrt(p1)*abs(randn(1,NN)+i*randn(1,NN));gain2=sqrt(p2)*abs(randn(1,NN)+i*randn(1,NN));fork=1:NNforkk=1:fsindex=(k-1)*fs+kk;ggain1(1,index)=gain1(1,k);ggain2(1,index)=gain2(1,k);endendy1=x;fork=1:delay2y2(1,k)=y1(1,k)*sqrt(p0);end
4、fork=(delay2+1):(NN*fs)y2(1,k)=y1(1,k)*sqrt(p0)+y1(1,k-delay1)*ggain1(1,k)+y1(1,k-delay2)*ggain2(1,k);end%匹配濾波器b=-ones(1,fs);b=b/fs;a=1;y=filter(b,a,y2);%仿真結(jié)束%UsethesemianalyticBERestimator.Thefollowingsetsupthesemi%analyticestimator.Findthemaximunma
5、gnitudeofthecrosscorrelation%andthecorrespondinglag.[corlags]=vxcorr(x,y);cmax=max(max(abs(cor)));nmax=find(abs(cor)==cmax);timelag=lags(nmax);corrmag=cmax;theta=angle(cor(nmax));y=y*exp(-i*theta);%derotate%NoiseBWcalibrationhh=impz(b,a);ts=1/16;nbw=
6、(fs/2)*sum(hh.^2);%Delaytheinput,anddoBERestimationonthelast128bits.Usemiddle%sample.Makesuretheindexdoesnotexceednumberofinputpoints.Eb%shouldbecomputedatthereceiverinput.index=(10*fs+8:fs:(NN-10)*fs+8);xx=x(index);yy=y(index-timelag+1);[n1n2]=size(
7、y2);ny2=n1*n2;eb=tb*sum(sum(abs(y2).^2))/ny2;eb=eb/2;[peideal,pesystem]=qpsk_berest(xx,yy,ebn0db,eb,tb,nbw);figuresemilogy(ebn0db,peideal,'b*-',ebn0db,pesystem,'r+-')xlabel('Eb/N0(db)');ylabel('ProbabilityofError');gridonaxis([01410^(-10)1]);%Endofsc
8、riptfile.相關(guān)的一些調(diào)用程序（4個(gè)）：[1]vxcorr.mfunction[c,lags]=vxcorr(a,b)%Thisfunctioncalculatestheunscaledcross-correlationof2vectorsof%thesamelength.Theoutputlength(c)islength(a)+length(b)-1.Itisa%simplifiedfunctionofxcorrfunctioninmatlabR12usingthedefinition
