【導(dǎo)讀】本次方向設(shè)計根據(jù)當(dāng)今現(xiàn)代通信技術(shù)的發(fā)展，對QPSK信號的工作原理進(jìn)行了仿真分析。近年來，軟件無線電作為解決通信體制兼容性問題的重要方法受到各方面的注意。其中軟件的增加、高頻譜效率的調(diào)制解調(diào)模塊是移動通信系統(tǒng)的關(guān)鍵技術(shù)，它的軟件。了解QPSK的實現(xiàn)方法及數(shù)學(xué)原理。個概念有個整體的理解，學(xué)習(xí)數(shù)字調(diào)制中誤碼率測試的標(biāo)準(zhǔn)及計算方法。道下傳輸，以及該方式的誤碼率測試。MATLAB仿真系統(tǒng)的注意事項，并鍛煉自己的編程能力，通過編程完成QPSK調(diào)制解調(diào)系統(tǒng)的仿真，以及誤碼率測試，并得出響應(yīng)波形。在完成要求任務(wù)的條件下，嘗試優(yōu)化程序。前已經(jīng)廣泛應(yīng)用于無線通信中，成為現(xiàn)代通信中一種十分重要的調(diào)制解調(diào)方式。同的相位來分別表示它們。由于每一種載波相位代表2個比特信息，故每個四進(jìn)制碼元。QPSK是在M=4時的調(diào)相技術(shù)，它規(guī)定了四種載波相位，分別為45°，135°，225°，成一組，共有四種組合，即00，01，10，11，其中每一組稱為雙比特碼元。

　　

【正文】 n=sgma*randn(1,2)。 % 2 normal distributed with 0, variance sgma if ((dsource1(i)==0) amp。 (dsource2(i)==0)), r=ray*s00+n。 elseif ((dsource1(i)==0) amp。 (dsource2(i)==1)), r=ray*s01+n。 elseif ((dsource1(i)==1) amp。 (dsource2(i)==0)), r=s10*ray+n。 else r=s11*ray+n。 end。 % The correlation metrics are puted below c00=dot(r,s00)。 c01=dot(r,s01)。 c10=dot(r,s10)。 c11=dot(r,s11)。 % The decision on the ith symbol is made next 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。 % Increment the error counter, if the decision is not correct symbolerror=0。 if (decis1~=dsource1(i)), numofbiterror=numofbiterror+1。 symbolerror=1。 end。 if (decis2~=dsource2(i)), numofbiterror=numofbiterror+1。 symbolerror=1。 21 end。 if (symbolerror==1), numofsymbolerror=numofsymbolerror+1。 end。 end counter=counter+1。 end ps=numofsymbolerror/(N*counter)。 % since there are totally N symbols pb=numofbiterror/(2*N*counter)。 % since 2N bits are transmitted cm_sm33()用于高斯信道誤碼率的計算 ： function [pb1,ps1]=cm_sm33(snr_in_dB) % [pb,ps]=cm_sm32(snr_in_dB) % CM_SM3 finds the probability of bit error and symbol error for % the given value of snr_in_dB, signal to noise ratio in dB. N=100。 E=1。 % energy per symbol snr=10^(snr_in_dB/10)。 % signal to noise ratio sgma=sqrt(E/snr)/2。 % noise variance s00=[1 0]。 s01=[0 1]。 s11=[1 0]。 s10=[0 1]。 % signal mapping % generation of the data source numofsymbolerror=0。 numofbiterror=0。 counter=0。 while(numofbiterror100) for i=1:N, temp=rand。 % a uniform random variable between 0 and 1 if (temp), % with probability 1/4, source output is 00 dsource1(i)=0。 dsource2(i)=0。 elseif (temp), % with probability 1/4, source output is 01 dsource1(i)=0。 dsource2(i)=1。 elseif (temp), % 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。 % detection and the probability of error calculation for i=1:N, % the received signal at the detection, for the ith symbol,is: n=sgma*randn(1,2)。 % 2 normal distributed with 0, variance 22 sgma if ((dsource1(i)==0) amp。 (dsource2(i)==0)), r=s00+n。 elseif ((dsource1(i)==0) amp。 (dsource2(i)==1)), r=s01+n。 elseif ((dsource1(i)==1) amp。 (dsource2(i)==0)), r=s10+n。 else r=s11+n。 end。 % The correlation metrics are puted below c00=dot(r,s00)。 c01=dot(r,s01)。 c10=dot(r,s10)。 c11=dot(r,s11)。 % The decision on the ith symbol is made next 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。 % Increment the error counter, if the decision is not correct symbolerror=0。 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 counter=counter+1。 end ps1=numofsymbolerror/(N*counter)。 % since there are totally N symbols pb1=numofbiterror/(2*N*counter)。 % since 2N bits are transmitted