【導讀】求通信距離遠、通信容量大、傳輸質量好。作為其關鍵技術之一的調制解調技。術一直是人們研究的一個重要方向。從模擬調制到數(shù)字調制，從二進制發(fā)展到。一個系統(tǒng)的通信質量，很大程度上依賴于所采用的調制方式。信號進行調制解調，從而達到傳遞數(shù)字信息的目的。目前，相移鍵控一般采用四相和多相制。字調制方案，傳遞數(shù)據(jù)，改變或者調制階段的載波信號。用PSK有限數(shù)目的階段，每個階。段分配一個唯一的二進制數(shù)字模式。一般情況下，各相位編碼的比特數(shù)相等。每一個模式的位元符號的形式是特定的階段。這種解調器，其目的是特別為使。制解調的輸出數(shù)據(jù)，從而恢復到原始數(shù)據(jù)。這就要求接收器能夠比較接收信號。的參考信號——這樣一個系統(tǒng)被稱為相干。本文就四相制移相。鍵控技術的概念，應用和原理方面進行了詳細介紹。行了系統(tǒng)仿真和實現(xiàn)方面也做了詳盡的闡述。隨著經濟?；牟粩喟l(fā)展，人們對通信的要求也越來越高。據(jù)格式、通信協(xié)議等。使用SUMLINK對QPSK調制、解調進行模擬。

　　

【正文】 單位寬度的高電平） Ray_demodata1=zeros(1,nb/delta_T)。 for q=1:nb Ray_demodata1((q1)/delta_T+1:q/delta_T)=Ray_demodata(q)。 end % 累計誤碼數(shù) % abs(demodatadata)求接收端和發(fā)射端 % 數(shù)據(jù)差的絕對值，累計之后就是誤碼個數(shù) Ray_num_BER=sum(abs(Ray_demodatadata)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 誤碼率計算 %% 調用了 cm_sm32()。和 cm_sm33(）函數(shù) %%聲明：函數(shù)聲明在另外倆個 M文件中 %%作用： cm_sm32()用于瑞利信道誤碼率的計算 %% cm_sm33()用于高斯信道誤碼率的計算 %% ecoh on/off 作用在于決定是否顯示指令內容 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% SNRindB1=0:1:6。 SNRindB2=0::6。 % 瑞利衰落信道 for i=1:length(SNRindB1), [pb,ps]=cm_sm32(SNRindB1(i))。 % 比特誤碼率 smld_bit_ray_err_prb(i)=pb。 34 smld_symbol_ray_err_prb(i)=ps。 disp([ps,pb])。 echo off。 end。 % 高斯信道 echo on。 for i=1:length(SNRindB1), [pb1,ps1]=cm_sm33(SNRindB1(i))。 smld_bit_awgn_err_prb(i)=pb1。 smld_symbol_awgn_err_prb(i)=ps1。 disp([ps1,pb1])。 echo off。 end。 % 理論曲線 echo on。 for i=1:length(SNRindB2), SNR=exp(SNRindB2(i)*log(10)/10)。 theo_err_awgn_prb(i)=*erfc(sqrt(SNR))。 theo_err_ray_prb(i)=*(11/sqrt(1+1/SNR))。 echo off。 end。 h = 。 % 輸出顯示部分 % 第一部分（理想） figure(1) subplot(3,2,1)。 plot(data0),title(39?；鶐盘?39。)。 axis([0 20210 2 2])。 subplot(3,2,2)。 psd(h,data1,39。fs39。,fs),title(39?；鶐盘柟β首V密度 39。)。 subplot(3,2,3)。 plot(s),title(39。調制信號 39。)。 axis([0 500 3 3])。 35 subplot(3,2,4)。 psd(h,s,39。fs39。,fs),title(39。調制信號功率譜密度 39。)。 subplot(3,2,5)。 plot(demodata1),title(39。解調輸出 39。)。 axis([0 20210 2 2])。 subplot(3,2,6)。 psd(h,demodata1,39。fs39。,fs),title(39。解調輸出功率譜密度 39。)。 % 通過高斯信道 figure(2) subplot(2,2,1)。 plot(s1),title(39。調制信號 (Awgn)39。)。 axis([0 500 5 5])。 subplot(2,2,2)。 psd(h,s1,39。fs39。,fs),title(39。調制信號功率譜密度 (Awgn)39。)。 subplot(2,2,3)。 plot(s111),title(39。高斯噪聲曲線 39。)。 axis([0 2021 5 5])。 figure(3) for i=1:nb/2 plot(idata(i),qdata(i),39。r+39。),title(39。QPSK信號星座圖（ Awgn） 39。)。hold on。 axis([2 2 2 2])。 plot(Awgn_ichsum(i),Awgn_qchsum(i),39。*39。)。hold on。 legend(39。理論值（發(fā)射端） 39。,39。實際值（接收端） 39。)。 end %通過高斯信道再通過瑞利衰落信道 figure(4) subplot(2,2,1) plot(Ray_s),title(39。調制信號 (Ray+Awgn)39。)。 axis([0 500 5 5])。 subplot(2,2,2)。 psd(h,Ray_s,39。fs39。,fs),title(39。調制信號功率譜密度 (Ray)39。)。 figure(5) for i=1:nb/2 plot(idata(i),qdata(i),39。r+39。),title(39。QPSK信號星座圖（ Awgn+Ray） 39。)。hold on。 axis([2 2 2 2])。 36 plot(Ray_ichsum(i),Ray_qchsum(i),39。*39。)。hold on。 legend(39。理論值（發(fā)射端） 39。,39。實際值（ 接收端） 39。)。 end figure(6) semilogy(SNRindB2,theo_err_awgn_prb,39。r39。),title(39。誤碼率曲線 39。)。hold on。 semilogy(SNRindB1,smld_bit_awgn_err_prb,39。r*39。)。hold on。 semilogy(SNRindB2,theo_err_ray_prb)。hold on。 semilogy(SNRindB1,smld_bit_ray_err_prb,39。+39。)。 xlabel(39。Eb/No39。)。ylabel(39。BER39。)。 legend(39。理論 AWGN39。,39。仿真 AWGN39。,39。理論 Rayleigh39。,39。仿真 Rayleigh39。)。 cm_sm32()用于瑞利信道誤碼率的計算 ： function [pb,ps]=cm_sm32(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 numofsymbolerror=0。 numofbiterror=0。 counter=0。 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 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 37 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, ray=raylrnd()。 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。 38 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 %