【正文】 chsum(i),39。subplot(2,2,4)。,fs),title(39。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 通過高斯信道figure(2)subplot(2,2,1)。)。psd(h,s,39。基帶信號功率譜密度39。plot(data0),title(39。% 理論曲線echo on。 end。 % 將極性碼變成對應(yīng)的波形信號end % 累計誤碼數(shù)% abs(demodatadata)求接收端和發(fā)射端% 數(shù)據(jù)差的絕對值，累計之后就是誤碼個數(shù)Ray_num_BER=sum(abs(Ray_demodatadata)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 誤碼率計算%% 調(diào)用了cm_sm32()。 end Ray_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。% else Ray_idata3(n)=0。 % 存放奇數(shù)位demodata(2:ml:nb)=qdata3。for n=1:nb/2 Awgn_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。% 抽樣判決的過程，與0作比較，data=0,則置1，否則置0for n=1:nb/2% A1(n)=sum(idata2((n1)/delta_T+1:n/delta_T))。%%%%%%%%%%%%%%%%%%%%%%%%%%% 高斯信道 s1=awgn(s,SNR)。endqdata1=qch.*b。 % 將奇偶位分開，因此間隔m1為2 qdata=datanrz(2:ml:nb)。 % 采樣數(shù) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 調(diào)制部分%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 基帶信號的產(chǎn)生data=randn(1,nb)。 % 基帶信號寬度，也就是頻率fc=10/T。最后顯示的是高斯信道和瑞利衰落信道的誤碼率對比。在AWGN信道中，判決結(jié)果具有最小的負(fù)號差錯概率。最后，將這兩個二進(jìn)制PSK信號相加，從而得期望的QPSK。下面介紹QPSK信號的產(chǎn)生和檢測。足不出戶便可以做實驗。了解QPSK的實現(xiàn)方法及數(shù)學(xué)原理。這類數(shù)字調(diào)制技術(shù)的優(yōu)點(diǎn)是已調(diào)信號具有相對窄的功率譜和對放大設(shè)備沒有線性要求,不足之處是其頻譜利用率低于線性調(diào)制技術(shù)。在完成要求任務(wù)的條件下，嘗試優(yōu)化程序。 其他其中，i＝1，2，2，4；E為發(fā)射信號的每個符號的能量，T為符號持續(xù)時間，載波頻率f等于nc/T，nc為固定整數(shù)。這兩個二進(jìn)制波形a1（t），和a2（t）被用來調(diào)制一對正交載波或者說正交基本函數(shù)：1（t）＝，2（t）＝。；類似地。星座圖反映可接收信號早高斯噪聲的影響下發(fā)生了誤碼，但是大部分還是保持了原來的特性。希望以后再接再厲，努力學(xué)習(xí)。 % 采樣頻率SNR=0。 % 創(chuàng)建一個1*nb/delta_T的零矩陣for q=1:nb data0((q1)/delta_T+1:q/delta_T)=data(q)。 % 奇數(shù)位數(shù)據(jù)與余弦函數(shù)相乘，得到一路的調(diào)制信號qch=zeros(1,nb/2/delta_T)。Ray_idata=idata1.*ray_ich。 % 這里面其實隱藏了一個串并轉(zhuǎn)換的過程qdata2=s1.*b。 else qdata3(n)=0。 else qdata4(n)=0。 % 對應(yīng)的信號與正余弦信號相乘% Ray_idata3=zeros(1,nb/2)。for n=1:nb/2 Ray_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。 % 存放奇數(shù)位Ray_demodata(2:ml:nb)=Ray_qdata4。 smld_symbol_ray_err_prb(i)=ps。 disp([ps1,pb1])。end。psd(h,data1,39。)。subplot(3,2,5)。,fs),title(39。subplot(2,2,2)。高斯噪聲曲線39。)。理論值（發(fā)射端）39。subplot(2,2,2)。r+39。)。r39。hold on。)。,39。numofbiterror=0。 % signal mapping% generation of the data sourcewhile(numofbiterror100)for i=1:N, temp=rand。 dsource2(i)=1。 (dsource2(i)==1)), r=ray*s01+n。 c11=dot(r,s11)。 else decis1=1。 symbolerror=1。E=1。numofbiterror=0。 dsource2(i)=0。 elseif ((dsource1(i)==0) amp。 c10=dot(r,s10)。 decis2=0。 if (decis2~=dsource2(i)), numofbiterror=numofbiterror+1。 % since 2N bits are transmitted。 end。 decis2=1。 % The decision on the ith symbol is made next c_max=max([c00,c01,c10,c11])。 elseif ((dsource1(i)==1) amp。 dsource2(i)=1。while(numofbiterror100)for i=1:N, temp=rand。 % signal to noise ratiosgma=sqrt(E/snr)/2。 if (symbolerror==1), numofsymbolerror=numofsymbolerror+1。 end。 if (c00==c_max), decis1=0。 (dsource2(i)==0)), r=s10*ray+n。end。 dsource2(i)=0。snr=10^(snr_in_dB/10)。,39。BER39。hold on。誤碼率曲線39。legend(39。QPSK信號星座圖（Awgn+Ray）39。fs39。實際值（接收端）39。axis([2 2 2 2])。axis([0 2000 5 5])。fs39。)。解調(diào)輸出39。subplot(3,2,4)。,fs),title(39。 % 類似于C語言的宏定義，方便以下的調(diào)用 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%