1、- QPSK通信系統(tǒng)的 Monte Carlo 仿真 一、實(shí)驗(yàn)?zāi)康?、提高獨(dú)立學(xué)習(xí)的能力；2、培養(yǎng)發(fā)現(xiàn)問題、解決問題和分析問題的能力；3、學(xué)習(xí)Matlab 的使用；4、掌握4PSK通信系統(tǒng)的Monte Carlo仿真方法；5、掌握4PSK通信系統(tǒng)的的組成原理；6、比較編碼信號(hào)和為編碼信號(hào)在隨機(jī)信道中的傳輸，加深對(duì)糾錯(cuò)編碼的理解；二、系統(tǒng)框圖及編程原理實(shí)驗(yàn)原理 PSK是利用載波的不同相位表示相應(yīng)的數(shù)字信息。對(duì)于二進(jìn)制相位調(diào)制（M=2)來說，兩個(gè)載波相位是0和。對(duì)于M相相位調(diào)制來說M=2k，這里k是每個(gè)傳輸符號(hào)的信息比特?cái)?shù)。4PSK是M=4的載波相位調(diào)制。這里，將理論差錯(cuò)概率與仿真的差錯(cuò)概率比較，

2、進(jìn)一步觀察仿真與理論值之間的差別。同時(shí)，用不同的判決準(zhǔn)則對(duì)接受信號(hào)進(jìn)行判決。并比較兩種判別方法的差別。一QPSK調(diào)制原理1信號(hào)能量分析一組M載波相位調(diào)制信號(hào)波形的一般表示式為 m=0，M-1式中是發(fā)送濾波器的脈沖形狀，它決定了傳輸信號(hào)的頻譜特性，A是信號(hào)的幅度。注意到，PSK信號(hào)對(duì)所有m都具有相等的能量，即代表每個(gè)傳輸符號(hào)的能量。2噪聲分析傳輸信號(hào)的信道假設(shè)被加性噪聲n(t)所污損，這樣信號(hào)在接收端將產(chǎn)生誤碼。因?yàn)閚(t)是功率譜為的白高斯過程的一個(gè)樣本函數(shù)，所以噪聲分量就是零均值高斯型的，即3信號(hào)判決分析最佳檢測器將接收信號(hào)向量r投射到M個(gè)可能的傳輸信號(hào)向量之一上去，并選取對(duì)應(yīng)于最大投影的向

3、量。據(jù)此，得到相關(guān)準(zhǔn)則為 m=0,M-1檢測器觀察到接收信號(hào)向量，并計(jì)算r在4種可能的信號(hào)向量上的投影。根據(jù)選取對(duì)應(yīng)于最大投影的信號(hào)點(diǎn)作為判決，從而判決出信號(hào)。同時(shí)，檢測器的判決準(zhǔn)則也可采用最小距離法，即利用星座圖上符號(hào)間的距離進(jìn)行判決，從而得到判決結(jié)果。二Monte Carlo 仿真過程仿真框圖如圖（一）圖（一） 用于 Monte Carlo仿真的4PSK系統(tǒng)的方框圖如圖所示，利用一個(gè)隨機(jī)數(shù)發(fā)生器，產(chǎn)生（0，1）范圍內(nèi)的隨機(jī)數(shù)。再將這個(gè)范圍分成四個(gè)相等的區(qū)間（0，0.25），（0.25，0.5），（0.5，0.75），（0.75，1.0），這些子區(qū)間分別對(duì)應(yīng)于00，01，11，10信息比特對(duì)

4、，再用這些比特對(duì)來選擇信號(hào)相位向量。加性噪聲的同相分量和正交分量 ，在上面討論過，即為零均值，方差為 的統(tǒng)計(jì)獨(dú)立的高斯隨機(jī)變量。在檢測器觀察到的接收信號(hào)向量 ，利用上面討論的兩種檢測方法，得到判決結(jié)果，并與傳輸符號(hào)作比較，最后對(duì)符號(hào)差錯(cuò)和比特差錯(cuò)計(jì)數(shù)三、實(shí)驗(yàn)內(nèi)容及程序分析（以下程序皆以N=1000為例）%映射比較子函數(shù)%函數(shù)分為四步第一步產(chǎn)生隨機(jī)序列，進(jìn)行4PSK映射。%第二步：調(diào)用高斯高斯噪聲子函數(shù)，產(chǎn)生正交兩路高斯噪聲，與輸出符號(hào)序列相加%第三步：檢測接受信號(hào)%第四步：計(jì)算誤碼率和誤比特率N=1000; %符號(hào)長度E=1; %計(jì)算噪聲方差 sgma=input(方差=);sgma=sqrt

5、(sgma);% 4PSK比特映射 s00=1 0; s01=0 1;s11=-1 0;s10=0 -1; %第一步產(chǎn)生隨機(jī)序列，進(jìn)行4PSK映射。生成隨機(jī)信源for i=1:N, %生成隨機(jī)信源 temp=rand; if (temp0.25), % With probability 1/4, source output is 00. dsource1(i)=0; dsource2(i)=0; elseif (temp0.5), % With probability 1/4, source output is 01. dsource1(i)=0; dsource2(i)=1; elseif

6、(temp %j化星座圖函數(shù)sgma=input(方差=);%輸入方差sgma=sqrt(sgma);%求sgmaN=1000;%輸入數(shù)據(jù)點(diǎn)數(shù)E=1;% 單個(gè)符號(hào)能量%4PSK映射規(guī)則%s00=1 0;s01=0 1;s11=-1 0;s10=0 -1;% 化星座圖 %產(chǎn)生信源序列% for i=1:N,% a uniform random variable between 0 and 1 temp=rand;%隨機(jī)序列if (temp0.25), % With probability 1/4, source output is 00. dsource1(i)=0; dsource2(i)=0

7、; elseif (temp0.5), % With probability 1/4, source output is 01. dsource1(i)=0; dsource2(i)=1; elseif (temp0.5) source(i)=1; else source(i)=0; endend%(7,4)漢明編碼%第二步：實(shí)現(xiàn)4碼到7碼的變換for k=1:N/4 out(7*k-6)=source(4*k-3); out(7*k-5)=source(4*k-2); out(7*k-4)=source(4*k-1); out(7*k-3)=source(4*k); out(7*k-2)=x

8、or(xor(out(7*k-6),out(7*k-5),out(7*k-4); out(7*k-1)=xor(xor(out(7*k-6),out(7*k-5),out(7*k-3); out(7*k)=xor(xor(out(7*k-6),out(7*k-4),out(7*k-3);end%第三步：進(jìn)行4PSK映射。s00=1 0;%s00=1 0s01=0 1;%s01=0 1s11=-1 0;%s11=-1 0s10=0 -1;%s10=0 -1%進(jìn)行4PSK映射%for i=1:N/2 if(source(2*i-1)=0&source(2*i)=0) s=s00; elseif(s

9、ource(2*i-1)=0&source(2*i)=1) s=s01; elseif(source(2*i-1)=1&source(2*i)=0) s=s10; elseif(source(2*i-1)=1&source(2*i)=1) s=s11; end%第四步：調(diào)用高斯高斯噪聲子函數(shù)，產(chǎn)生正交兩路高斯噪聲，與輸出符號(hào)序列相加n(1),n(2)=gnguass(0,sgma);%調(diào)用gnguass函數(shù) r=s+n;%(7,4)漢明解碼%第五步：7碼到4碼的解碼for k=1:N/4 jiema(4*k-3)=out(7*k-6); jiema(4*k-2)=out(7*k-5); jie

10、ma(4*k-1)=out(7*k-4); jiema(4*k)=out(7*k-3);end%第六步：進(jìn)行碼元檢測%最大投影點(diǎn)準(zhǔn)則%c00=dot(r,s00);%取r在s00上的投影c01=dot(r,s01);%取r在s01上的投影c10=dot(r,s10);%取r在s10上的投影c11=dot(r,s11);%取r在s11上的投影c_max=max(c00 c01 c10 c11);%取c00，c01，c10，c11中的最大值 if (c00=c_max), decis1=0;decis2=0; elseif(c01=c_max), decis1=0;decis2=1; elseif

11、(c10=c_max), decis1=1;decis2=0; else decis1=1;decis2=1; end out(2*i-1)=decis1; out(2*i)=decis2;end%第七步：計(jì)算誤碼率和誤比特率%計(jì)算符號(hào)差錯(cuò)概率%symbolerror=0;for i=1:N/2 if(out(2*i-1)=source(2*i-1)|out(2*i)=source(2*i) symbolerror=symbolerror+1; endendps=2*symbolerror/N;%計(jì)算比特差錯(cuò)概率 biterror=0;for i=1:N if(out(i)=source(i)

12、 biterror=biterror+1; endendpb=biterror/N; 1、高斯白噪聲子函數(shù)functiongsrv1,gsrv2=gnguass(m,sgma)if nargin=0 m=0; sgma=1;elseif nargin=1 sgma=m; m=0;endu=rand;z=sgma*sqrt(2*log10(1/(1-u);u=rand;gsrv1=m+z*cos(2*pi*u);gsrv2=m+z*sin(2*pi*u);函數(shù)評(píng)注：編程產(chǎn)生正交兩路高斯白噪聲2、映射檢測比較子函數(shù)（1）、采用最大投影準(zhǔn)則function pb,ps=cm_sm32(snr_in_

13、dB)% pb,ps=cm_sm32(snr_in_dB)%pb誤符號(hào)率% ps誤比特率N=1000; %符號(hào)長度E=1; %計(jì)算噪聲方差 snr=10(snr_in_dB/10); sgma=sqrt(E/snr/2); s00=1 0; % 比特映射 s01=0 1;s11=-1 0;s10=0 -1; for i=1:N, %生成隨機(jī)信源 temp=rand; if (temp0.25), % With probability 1/4, source output is 00. dsource1(i)=0; dsource2(i)=0; elseif (temp0.5), % With

14、probability 1/4, source output is 01. dsource1(i)=0; dsource2(i)=1; elseif (temp0.75), % With probability 1/4, source output is 10. dsource1(i)=1; dsource2(i)=0; else % With probability 1/4, source output is 11. dsource1(i)=1; dsource2(i)=1; end;end;numofsymbolerror=0; %檢測錯(cuò)誤并計(jì)算錯(cuò)誤率 numofbiterror=0;fo

15、r i=1:N, gsrv1,gsrv2=gnguass(0,sgma); %調(diào)用高斯噪聲子函數(shù) n(1)=gsrv1; n(2)=gsrv2; if (dsource1(i)=0) & (dsource2(i)=0), r=s00+n; elseif (dsource1(i)=0) & (dsource2(i)=1), r=s01+n; elseif (dsource1(i)=1) & (dsource2(i)=0), r=s10+n; else r=s11+n; end; c00=dot(r,s00); %最大投影點(diǎn)準(zhǔn)則，即向量點(diǎn)乘積 c01=dot(r,s01); c10=dot(r,s

16、10); c11=dot(r,s11); c_max=max(c00 c01 c10 c11); if (c00=c_max), decis1=0; decis2=0; elseif (c01=c_max), decis1=0; decis2=1; elseif (c10=c_max), decis1=1; decis2=0; else decis1=1; decis2=1; end; symbolerror=0; %設(shè)置符號(hào)錯(cuò)誤標(biāo)志，以統(tǒng)計(jì)錯(cuò)誤個(gè)數(shù) if (decis1=dsource1(i), numofbiterror=numofbiterror+1; symbolerror=1; end; if (decis2=dsource2(i), numofbiterror=numofbiterror+1; symbolerror=1; end; if (symbolerror=1), numofsymbolerror = numofsymbolerror+1; end;end;ps=numofsymbolerror/N; % 錯(cuò)誤率計(jì)算pb=numofbiterror/(2*N); （2）最小距離判別法（只需將上面紅色部分改為下面程序即可）c00=sqrt(r(1)-s00(1)2+(r(2)-s00(2)2); %最小距離判決 c01=sqrt(r(1)-s01(1)2