【文章內容簡介】 仿真流程 仿真中用到的 MATLAB 函數(shù)介紹 Randsrc： x1=randsrc(1,1000,[1,0。,])。%產(chǎn)生隨機信源，信源符號 1和 0的概率分別為 ， Symerr : [wnum(1,SNR+11),pe(1,SNR+11)]=symerr(x1(1:length(R_X)),R_X)。%計算誤碼率 Xor : b(n+1)=xor(x1(n),b(n))。 %差分碼編碼規(guī)則 b(n+1)是 x1(n)與b(n)相異 或 O_a=[]。O_a(1)=xor(R_X(1),1)。%碼反變換 Rcosfir: Num=rcosfir(Alfa,N_Filter,Rate,Tc,39。sqrt39。)。%根升余弦濾波器單位沖激序列 Filter : x=filter(Num,[1],c)。 %濾波 Awgn : y=awgn(x,snr)。 %向功率為 0dB 的信號 x 添加高斯白噪聲，輸出信號 y的信噪比 SNR為參數(shù) snr (單位 dB)。如果信號 x為復信號，該函數(shù)添加復高斯白噪聲 。 R_Y=awgn(y,SNR,39。measured39。)。%加上高斯白噪聲 誤碼率與信噪比的關系結果分析 與仿真系統(tǒng)參數(shù)設計 分析結果 —— 誤碼率與信噪比的關系 圖 51 誤碼率與信噪比的關系 由 結果圖可知隨信噪比增加誤碼率降低 。 PSK 比 DPSK 抗噪性能好，但由于 2PSK系統(tǒng)存在 “ 反向工作 ” 問題，而 2DPSK 系統(tǒng)不存在 “ 反向工作 ” 問題，因此在實際應用中，真正作為傳輸用的數(shù)字調相信號幾乎都是 DPSK 信號。 2PSK 信號的帶寬、頻帶利用率和 2DPSK 信號的相同。 仿真參數(shù)設計 抽樣速率 fs=256Hz 載波速率 fc=64Hz 碼元速率 sigv=16B 碼元周期 Tc=1/sigv。 升余弦濾波器參數(shù) Alfa= 采用了符號 1和 0概率分別是 ， 的隨機信源，并且信噪比分別設置在10dB— 10dB 變化范圍內 附錄、相移鍵控（ PSK）和差分相移鍵控（ DPSK）的 程序 設計 function pskvsdpsk clc。clear all。close all wnum=zeros(2,21)。 %用于存儲不同信噪比不同調制方式下的錯誤符號個數(shù) pe=zeros(2,21)。 %存儲誤碼率 i=1。 for SNR=10:10 %信噪比變化范圍 x1=randsrc(1,1000,[1,0。,])。%產(chǎn)生隨機信源，信源符號 1和 0的概率分別為 ， if SNR==0 i=0。 figure,subplot(211),plot(x1),title(39。原信號波形 39。) subplot(212),plot(abs(fftshift(x1))),title(39。原信號波形頻譜 39。) else i=1。 end NSig=length(x1)。%碼元個數(shù) % %基于 2PSK 的調制與解調 for n=1:NSig %2PSK 信 號產(chǎn)生 if x1(n)==1 Sig(n)=1。%2PSK 信號 else Sig(n)=1。 end end R_X=runup(Sig,SNR,i)。%調用子函數(shù)，完成調制與解調，返回解調信號 %判決 for n=1:length(R_X) if R_X(n)0 R_X(n)=1。 else R_X(n)=0。 end end if SNR==0 figure,subplot(411),plot(R_X),title(39。2PSK 抽樣判決后信號波形 ,SNR=039。) axis([0,100,]) subplot(412),plot(abs(fft(R_X))),title(39。2PSK 抽樣判決后信號頻譜 ,SNR=039。) axis([0,100,inf,inf]) subplot(413),plot(x1),title(39。原信號波形 39。) axis([0,100,]) subplot(414),plot(abs(fft(x1))),title(39。原信號頻譜 39。) axis([0,100,inf,inf]) else i=1。 end [wnum(1,SNR+11),pe(1,SNR+11)]=symerr(x1(1:length(R_X)),R_X)。%計算誤碼率 % %基于 2DPSK 的調制與解調 if SNR==0 i=1。 end b=[]。b(1)=1。 %b 為差分碼，設第一個符號為 1 for n=1:NSig % b(n+1)=xor(x1(n),b(n))。 %差分碼編碼規(guī)則 b(n+1)是 x1(n)與 b(n)相異或 if b(n+1)==1 Sig(n)=1。 %2DPSK 信號 else Sig(n)=1。 end end R_X=runup(Sig,SNR,i)。%調用子函數(shù)，完成調制與解調，返回解調信號 %判決 for n=1:length(R_X) if R_X(n)0 R_X(n)=1。 else R_X(n)=0。 end end O_a=[]。O_a(1)=xor(R_X(1),1)。%碼反變換 if O_a(1)==1 O_Y(1)=1。 else O_Y(1)=0。 end for n=2:length(R_X) O_a(n)=xor(R_X(n),R_X(n1))。 if O_a(n)==1 O_Y(n)=1。 %輸出信號 else O_Y(n)=0。 end end if SNR==0 figure,subplot(411),plot(O_Y),title(39。2DPSK 抽樣判決后信號波形 ,SNR=039。) axis([0,100,]) subplot(412),plot(abs(fft(O_Y))),title(39。2DPSK 抽樣判決后信號頻譜 ,SNR=039。) axis([0,100,inf,inf]) subplot(413),plot(x1),title(39。原信號波形 39。) axis([0,100,]) subplot(414),plot(abs(fft(x1))),title(39。原信號頻譜 39。) axis([0,100,inf,inf]) else end [wnum(2,SNR+11),pe(2,SNR+11)]=symerr(x1(1:length(O_Y)),O_Y)。%計算誤碼率 end wnum,pe figure t=10:10。 plot(t,pe(1,:),39。*39。)。 %畫出誤碼率曲線 gtext(39。2psk39。) hold on plot(t,pe(2,:),39。.39。)。 %畫出誤碼率曲線 gtext(39。2dpsk39。) xlabel(39。信噪比