【正文】 title(39。基帶信號 p39。)%抽取 IQ 信號 QPSK 兩個基帶符號為一組，映射到一對 IQpQI=reshape(p,2,5)。pQ=pQI(2,:)。pI=pQI(1,:)。d2=[0:1:t(end)1]。%IQ 信號波形%dQ=[d2。pQ]39。%dI=[d2。pI]39。%mI=pulstran(,dI,39。rectpuls39。,1)。%mQ=pulstran(,dQ,39。rectpuls39。,1)。 %IQ 值映射 “0”映射成1， “1”映射成 1I(pI==1)=1/sqrt(2)。I(pI==0)=1/sqrt(2)。Q(pQ==1)=1/sqrt(2)。Q(pQ==0)=1/sqrt(2)。dImap=[d2。I]39。dQmap=[d2。Q]39。mImap=pulstran(,dImap,39。rectpuls39。,1)。mQmap=pulstran(,dQmap,39。rectpuls39。,1)。 subplot(4,1,2)plot(t,mImap)axis([t(1) t(end) ])grid ontitle(39。QPSK 映射后的 I 路信號39。)subplot(4,1,3)plot(t,mQmap)axis([t(1) t(end) ])grid ontitle(39。QPSK 映射后的 Q 路信號39。)%載波信號fc=10 %載波頻率c1=cos(2*pi*fc*t)。c2=sin(2*pi*fc*t)。%已調(diào)信號s1=mImap.*c1。s2=mQmap.*c2。s=s1s2。subplot(4,1,4)plot(t,s)axis([t(1) t(end) ])title(39。QPSK 進入信道的信號 s=s1s239。)%QDSK 信號波形clc。clear all。close all。%全局變量Ts=。 %時間分辨率(時點間間隔） ，抽樣時間間隔t=0:Ts:5。 %時間坐標 sN=length(t)。fs=1/Ts。df=fs/(N1)。 %df 為頻率分辨率（頻率間隔）n=(N1)/2:(N1)/2。f=n*df。 %頻率坐標%基帶信號figure。p=[1 0 1 1 0 1 0 0 1 1]。 %待傳送的編碼串d1=[0::t(end)]。 %每個編碼的發(fā)送延遲時間p1=[p, zeros(length(d1)length(p))]。 %在發(fā)送串后補零d=[d1。p1]39。 %產(chǎn)生 d 矩陣m=pulstran(,d,39。rectpuls39。,)。subplot (4,1,1)。plot(t,m)axis([t(1) t(end) ])grid ontitle(39?；鶐盘?p39。)%抽取 IQ 信號 QDPSK 兩個基帶符號為一組 ，映射到一對 IQpQI=reshape(p,2,5)。pQ=pQI(2,:)。pI=pQI(1,:)。d2=[0:1:t(end)1]。%IQ 信號波形%dQ=[d2。pQ]39。%dI=[d2。pI]39。%mI=pulstran(,dI,39。rectpuls39。,1)。%mQ=pulstran(,dQ,39。rectpuls39。,1)。 %IQ 值映射 “0”映射成+1， “1”映射成1I(pI==1)=1/sqrt(2)。I(pI==0)=+1/sqrt(2)。Q(pQ==1)=1/sqrt(2)。Q(pQ==0)=+1/sqrt(2)。dImap=[d2。I]39。dQmap=[d2。Q]39。mImap=pulstran(,dImap,39。rectpuls39。,1)。mQmap=pulstran(,dQmap,39。rectpuls39。,1)。 subplot(4,1,2)plot(t,mImap)axis([t(1) t(end) ])grid ontitle(39。QDPSK 映射后的 I 路信號39。)subplot(4,1,3)plot(t,mQmap)axis([t(1) t(end) ])grid ontitle(39。QDPSK 映射后的 Q 路信號39。)%載波信號fc=10 %載波頻率c1=cos(2*pi*fc*t+pi/2)。c2=sin(2*pi*fc*tpi/2)。%已調(diào)信號s1=mImap.*c1。s2=mQmap.*c2。s=s1+s2。subplot(4,1,4)plot(t,s)axis([t(1) t(end) ])title(39。QDPSK 進入信道的信號 s=s1+s239。)%QPSK 功率譜密度m=1。p1=10。Fs=1000。nfft=1024。 cs_n=xcorr(s,39。unbiased39。)。 %計算序列的自相關函數(shù) CSk=fft(cs_n,nfft)。 %傅里葉變換Pxx=abs(CSk)。 index=0:round(nfft/21)。 k=index*Fs/nfft。 plot_Pxx=10*log10(Pxx(index+1))。 figureplot(k,plot_Pxx)。 title(39。QPSK 功率譜密度39。)%QDPSK 功率譜密度m=。p1=10。Fs=1000。nfft=1024。 cs_n=xcorr(s,39。unbiased39。)。 %計算序列的自相關函數(shù) CSk=fft(cs_n,nfft)。 %傅里葉變換Pxx=abs(CSk)。 index=0:round(nfft/21)。 k=index*Fs/nfft。 plot_Pxx=10*log10(Pxx(index+1))。 figureplot(k,plot_Pxx)。 title(39。QDPSK 功率譜密度 39。)%QPSK 誤碼率clear all。clc。M = 4。 x = randint(100000,1,M)。 % 隨機序列s= pskmod(x,M)。 % 調(diào)制SNR = 0:2:30。 % 信噪比范圍for n = 1:length(SNR)r = awgn(s,SNR(n),39。measured39。)。 % 加高斯噪聲r1 = pskdemod(r,M)。 % 解調(diào)[nErrors, BER(n)] = biterr(x(1:end),r1(1:end))。%計算仿真誤比特率endfor i=1:length(SNR) SNR1(i)=10^(SNR(i)/10)。 % 信噪比 BERtheory(i)=1(*erfc(sqrt(SNR1(i)/2)))^2。 BERtheory(i)=BERtheory(i)/2。endfiguresemilogy(SNR,BERtheory,39。b39。)。 %不同顏色畫出誤比特率曲線對比hold onsemilogy(SNR,BER,39。r39。)。 %不同顏色畫出誤比特率曲線對比legend(39。 理論 BER39。,39。仿真 BER39。)。 %曲線名稱grid onxlabel(39。SNR (dB)39。)。 ylabel(39。BER39。)。 %橫縱坐標title(39。QPSK 通過 AWGN 信道39。)。%QDPSK 誤碼率clear all。close all。clc。%%參數(shù)設置Fc = 5e6。data_len=10000。 %%隨機產(chǎn)生傳輸信號rand_data_len =randn(1,data_len)。for i=1:data_len if rand_data_len(i)= input(i)=1。 else input(i)=0。 endend%%串并轉(zhuǎn)換for i=1:data_len if rem(i,2)==1 if input(i)==1 I(i)=1。 I(i+1)=1。 else I(i)=1。 I(i+1)=1。 end else if input(i)==1 Q(i1)=1。 Q(i)=1。 else Q(i1)=1。 Q(i)=1。 end endend%% 碼變換I_diff=zeros(1,data_len)。Q_diff=zeros(1,data_len)。I_diff(1)=I(1)。Q_diff(1)=Q(1)。for i=2:data_len if I(i)==1amp。amp。Q(i)==1 % 0+1 11 if I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 else I_diff(i)=1。 Q_diff(i)=1。 end elseif I(i)==1amp。amp。Q(i)==1 if I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 else I_diff(i)=1。 Q_diff(i)=1。 end elseif I(i)==1amp。amp。Q(i)==1 if I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 else I_diff(i)=1。 Q_diff(i)=1。 end else if I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 elseif I_diff(i1)==1amp。amp。Q_diff(i1)==1 I_diff(i)=1。 Q_diff(i)=1。 else I_diff(i)=1。 Q_diff(i)=1。 end endend%%零內(nèi)插zero=5。