基于MATLAB的DPSK通信系統(tǒng)仿真

#目錄TOC\o"1-5"\h\z1設(shè)計(jì)題目32設(shè)計(jì)要求33設(shè)計(jì)原理3DPSK調(diào)制?DPSK解調(diào)?4設(shè)計(jì)流程55代碼實(shí)現(xiàn)5基帶信號(hào)的獲取55.2差分編碼6DPSK調(diào)制7AWGN信道8限噪處理85.6差分相干解調(diào)9AWGN信道115.7功率譜密度圖的繪制6測(cè)試結(jié)果126.1調(diào)制解調(diào)正確性的測(cè)試12系統(tǒng)抗噪性能的測(cè)試147總結(jié)178心得體會(huì)17附:完整代碼181設(shè)計(jì)題目用MATLAB仿真DPSK通信系統(tǒng)2設(shè)計(jì)要求基本參數(shù)如下：1）Fc=1800Hz2）Rb=1200bps3）考慮加入AWGN信道，Eb/N0=20dB4）自定義一小段信息序列結(jié)果要求：1）繪制各階段信號(hào)波形2）繪制信號(hào)的PSD3設(shè)計(jì)原理二進(jìn)制差分相移鍵控常簡稱為二相相對(duì)調(diào)相，記為2DPSK。它不是利用載波相位的絕對(duì)數(shù)值傳送數(shù)字信息，而是用前后碼元的相對(duì)載波相位值傳送數(shù)字信息。所謂相對(duì)載波相位是只本碼元初相與前一碼元初相之差。DPSK調(diào)制DPSK方式即是利用前后相鄰碼元的相對(duì)載波相位值去表示數(shù)字信息的一種方式。例如，假設(shè)相位值用相位偏移△屮表示（△屮定義為本碼元初相與前一碼兀初相只差），并設(shè)△屮=nf數(shù)字信息1△屮=0—數(shù)字信息0則數(shù)字信息序列與2DPSK信號(hào)的碼元相位關(guān)系可舉例表示如如下：數(shù)字信息：00111001012DPSK信號(hào)相位：000n0nnn00n或nnn0n000nn0畫出的2PSK及DPSK信號(hào)的波形如圖3-1所示。

碼變?yōu)橄鄬?duì)碼。2DPSK產(chǎn)生的原理框圖如圖2-2所示⑴，圖（a）為模擬調(diào)制法,圖（b）為鍵控法。（a）模擬調(diào)制法（b）鍵控法圖3-22DPSK信號(hào)的調(diào)制原理圖從上面分析可見，無論接收信號(hào)是2DPSK還是2PSK信號(hào)，單從接收端看是區(qū)分不開的。因此2DPSK信號(hào)的功率譜密度和2PSK信號(hào)的功率譜密度是完全一樣的。DPSK解調(diào)DPSK信號(hào)可以采用相干解調(diào)法（極性比較法）和差分相干解調(diào)法（相位比較法）。本次設(shè)計(jì)采用差分相干解調(diào)法，圖3-3為差分相干解調(diào)法原理圖。差分相干解調(diào)法的思想是直接比較前后碼元的相位差。由于此時(shí)的解調(diào)已同時(shí)完成了碼變換作業(yè)，故無需另加碼變換器。由于這種解調(diào)方法又無需專門的相干載波，故事一種比較實(shí)用的方法。

4設(shè)計(jì)流程設(shè)計(jì)流程如下框圖所示，箭頭上的信號(hào)表示實(shí)際傳遞的信號(hào)名稱。bn差分編碼［dn=>DPSKbn差分編碼［dn=>DPSK調(diào)制modzmodAWGN信道=>BPF對(duì)比Bnmtdemod延時(shí)相乘Bnmtdemod延時(shí)相乘差分相干解調(diào)5代碼實(shí)現(xiàn)5.1基帶信號(hào)的獲取利用MATLAB自帶函數(shù)randint產(chǎn)生定長度的二進(jìn)制隨機(jī)序列sign，為了將sign轉(zhuǎn)換為單極性NRZ信號(hào)，要對(duì)sign的每個(gè)隨機(jī)碼進(jìn)行T點(diǎn)增樣，最后得到基帶數(shù)字信號(hào)bn代碼實(shí)現(xiàn)如下:%%%%%%%%%%%%fc=1800;Rb=1200;Tb=1/Rb;T=50;M=8;N=9;初始參量定義%%%%%%%%%%%載波頻率%基帶數(shù)據(jù)傳輸速率%每個(gè)Tb內(nèi)進(jìn)行的增樣點(diǎn)數(shù)%原始碼字序列的長度%差分碼字序列的長度sign=randint(1,M);%產(chǎn)生長度為M的隨機(jī)二進(jìn)制序列SNR=20;%信噪比%%%%%%%%%%%%bn=[];forn=1:M對(duì)sign進(jìn)行T點(diǎn)增樣后的基帶數(shù)字信號(hào)%%%%%%%%%%ifsign(n)==O;F=zeros(1,T);elseF=ones(1,T);endbn=[bn,F];end繪圖基帶數(shù)字信號(hào)囪2i1111111111111111TOC\o"1-5"\h\z1；t……3IIIIIIIIIIIIIIIIJIIII-I111IIIIIIIIIIIII■2!!!!'01002003004005.2差分編碼設(shè)差分碼cfm(二進(jìn)制序列)，差分編碼原理為cfm[n]=b[n]xorcfm[n-l]例如，如果sign=[10010011],那么cfm=[100011101];同時(shí)，還要對(duì)得到的差分碼cfm進(jìn)行T點(diǎn)增樣，得到差分信號(hào)dn(單極性NRZ信號(hào))。代碼實(shí)現(xiàn)如下：%%%%%%%%%%%%%差分編碼得到差分碼cfm%%%%%%%%%%cfm=[];cfm(1)=1;%設(shè)。壬皿的初值為1fori=1:Mcfm(i+1)=xor(sign(i),cfm(i));end%%%%%%%%%%%%%對(duì)cfm進(jìn)行T點(diǎn)增樣%%%%%%%%%%%%%%dn=[];%差分信號(hào)forn=1:N;ifcfm(n)==0;A=zeros(1,T);%每個(gè)symbol進(jìn)行T點(diǎn)增樣elseA=ones(1,T);enddn=[dnA];end繪圖DPSK調(diào)制DPSK調(diào)制由于PSK調(diào)制的主要原因是DPSK調(diào)制解決了“不確定性反相”問題，即DPSK信號(hào)通過相鄰時(shí)隙載波相位的變化與否來“攜帶”信息。在代碼實(shí)現(xiàn)這一步的過程中，我假設(shè)了一個(gè)相鄰時(shí)隙相位差變量delta，易知,當(dāng)sign的一個(gè)碼字為'1'時(shí)，delta=pi，否則delta=O.通過delta就能反映前后時(shí)隙的載波相位變化，進(jìn)而實(shí)現(xiàn)DPSK調(diào)制。代碼實(shí)現(xiàn)如下：%%%%%%%%%%%%%%%%dpsK調(diào)制%%%%%%%%%%%%%%%t=Tb/T:Tb/T:Tb；%對(duì)Tb=1/Rb間隔進(jìn)行T點(diǎn)劃分mod=cos(2*fc*pi*t);%調(diào)制后的信號(hào)delta=0;%相位差forn=1:Mifsign(n)==1;delta=delta+pi;elsedelta=delta;endc2=cos(2*fc*pi*t+delta);mod=[modc2];end繪圖AWGN信道利用MATLAB自帶函數(shù)awgn實(shí)現(xiàn)對(duì)已調(diào)信號(hào)mod進(jìn)行加噪處理。代碼實(shí)現(xiàn)如下：%%%%%%%%%%%調(diào)制信號(hào)過高斯白噪聲的信道%%%%%%%modz=awgn(mod,SNR);%^AWGN信道繪圖5.5限噪處理利用MATLAB自帶函數(shù)firl構(gòu)造FIRI型帶通濾波器，使其對(duì)準(zhǔn)信號(hào)頻帶,讓信號(hào)幾乎無失真通過，同時(shí)抑制帶外噪聲，此處取BPF的通帶寬度為2*Rb，中心頻率為fc。(由于設(shè)計(jì)的是數(shù)字濾波器，所以不要忘記對(duì)頻率進(jìn)行歸一化處理)值得注意的是，由于設(shè)計(jì)的濾波器是100階，所以通過BPF后的信號(hào)會(huì)有出現(xiàn)接近于一個(gè)symbol寬度的時(shí)延。所以在信號(hào)modz接入BPF前，要對(duì)其進(jìn)行一個(gè)拓寬處理，對(duì)modz最后一個(gè)symbol進(jìn)行復(fù)制，得到modz_broad，再接入至【」BPF上。同時(shí)，不要忘記對(duì)從BPF出來的信號(hào)r0進(jìn)行截取，因?yàn)閞0的第一個(gè)symbol是濾波器時(shí)延帶來的。代碼實(shí)現(xiàn)如下：%%%%%%%%%%%BPF%%%%%%%%%%%%%modz_broad=[modzmodz((N-1)*T+1:N*T)];fs=Rb*T;hBPF=fir1(100,2*pi*[fc-Rbfc+Rb]/fs);rO=filter(hBPF,1,modz_broad);r=r0(T+1:(N+1)*T);%對(duì)r0進(jìn)行截取繪圖5.6差分相干解調(diào)差分相干解調(diào)的核心思想是比較兩個(gè)相鄰時(shí)隙上信號(hào)的相位，從而直接還原出信號(hào)Bn，分為三步：(1)延時(shí)相乘代碼實(shí)現(xiàn)如下：%%%%%%%%%%%延時(shí)相乘%%%%%%%r_shift=r(T+1:N*T);demod=[];fori=1:M*T;a=r_shift(i).*r(i);demod=[demoda];end繪圖⑵LPF利用MATLAB自帶函數(shù)firl構(gòu)造FIRI型低通濾波器，此處取LPF的通帶寬度為Rb。同BPF一樣，由于設(shè)計(jì)的濾波器是100階，所以通過LPF后的信號(hào)會(huì)有出現(xiàn)接近于一個(gè)symbol寬度的時(shí)延。所以要對(duì)信號(hào)demod接入LPF前，進(jìn)行拓寬處理，得到demod_broad，再接入到LPF上。同樣，不要忘記對(duì)從LPF出來的信號(hào)mt0進(jìn)行截取，因?yàn)閙t0的第一個(gè)symbol是濾波器時(shí)延帶來的。代碼實(shí)現(xiàn)如下：%%%%%%%%%%%%%LPF%%%%%%%%%%demod_broad=[demoddemod((M-1)*T+1:M*T)]hLPF=fir1(100,2*pi*[Rb]/fs);mtO=filter(hLPF,1,demod_broad);mt=mt0(T+1:(M+1)*T);%^mt0進(jìn)行截取繪圖⑶抽樣判決實(shí)現(xiàn)思想：對(duì)mt—個(gè)symbol的T個(gè)采樣點(diǎn)進(jìn)行求和，如果sum大于0,則Bn對(duì)應(yīng)這個(gè)symbol內(nèi)的所有點(diǎn)的值為1,反之為0.代碼實(shí)現(xiàn)如下：%%%%%%%%%%%%%抽樣判決%%%%%%%%%%Bn=zeros(1,M*T);fori=1:M;ifsum(mt((i-1)*T+1:i*T))>0;Bn((i-1)*T+1:i*T)=0;elseBn((i-1)*T+1:i*T)=1;endend繪圖從上圖可知，最終經(jīng)DPSK通信系統(tǒng)調(diào)制解調(diào)后得到的信號(hào)Bn與一開始的基帶數(shù)字信號(hào)bn完全一致。5.[?功率譜密度圖的繪制采用Welch法繪制DPSK調(diào)制信號(hào)mod的功率譜密度圖(Powerspectrummagnitude)。選擇適當(dāng)?shù)拇昂瘮?shù)w(n)，并再周期圖計(jì)算前直接加進(jìn)去，加窗的優(yōu)點(diǎn)是無論什么樣的窗函數(shù)均可使譜估計(jì)非負(fù)。二是在分段時(shí)，可使各段之間有重疊，這樣會(huì)使方差減小。本設(shè)計(jì)中嘗試了三種窗，分別是矩形窗、海明窗和Blacman窗。代碼實(shí)現(xiàn)如下：%%%%%%%%%%%%%%%PSD%%%%%%%%%%%%%%%window=boxcar(100);%矩形窗window1=hamming(100);%海明窗window2=blackman(100);%Blackman窗Nfft=N*T*Rb;noverlap=0;%數(shù)據(jù)無重疊range='onesided';%計(jì)算單邊帶功率譜密度[Pxx,f]=pwelch(mod,window,noverlap,Nfft,fs,range);[Pxx1,f]=pwelch(mod,window1,noverlap,Nfft,fs,range);[Pxx2,f]=pwelch(mod,window2,noverlap,Nfft,fs,range);figure(3)subplot(311)plot(f,Pxx);axis([0500000.0008]);title('rectangularwindowPSD');xlabel('f/Hz');subplot(312)plot(f,Pxx1);axis([0500000.0006]);title('HammingPSD');xlabel('f/Hz');subplot(313)plot(f,Pxx2);axis([0500000.0006]);title('BlackmanPSD');xlabel('f/Hz');

rectangularwindowPSDx10rectangularwindowPSDx10LLLLLLLLLr■-r-rrrr―r—rr0500100015002000250030003500400045005000f/Hzx10-4HammingPSD-LLLLLI-LLL-rrrrrrr0500100015002000250030003500400045005000f/Hzx10-4BlackmanPSDLLLLL1-LLL■--r-一二rf~?一一rr"r~rr—F0500100015002000250030003500400045005000505050f/Hz繪圖繪圖從第一幅圖即rectangularPSD中可以看出，單邊帶功率譜密度圖的中心頻率大致在1800Hz即載波頻率處，第一零點(diǎn)大致在600Hz和3000Hz,即fc-Rb和fc+Rb處。6測(cè)試結(jié)果6.1調(diào)制解調(diào)正確性的測(cè)試在前面我取用了一個(gè)特例sign=[10010011]來說明每一步的正確性。下面利用randint產(chǎn)生定長度的二進(jìn)制隨機(jī)序列sign，多次運(yùn)行結(jié)果都顯示最終差分相干解調(diào)得出的信號(hào)Bn與基帶數(shù)字信號(hào)bn一致，說明此次基于MATLAB設(shè)計(jì)的DPSK通信系統(tǒng)能夠成功實(shí)現(xiàn)DPSK信號(hào)的調(diào)制與解調(diào)。某幾次運(yùn)行結(jié)果的截圖如下:基帶數(shù)字信號(hào)bn差分信號(hào)dn21』J|il|1'l''m\0DPSK調(diào)制信號(hào)mod-20200400過AWGN信道后的modz信號(hào)過BPF后的信號(hào)r0200400J1f-/W71V120-20200400-10Luv-1111-￥ii.4■-■r,0200400延時(shí)相乘后的信號(hào)demod1過LPF后的信號(hào)mt0200400經(jīng)抽樣判決后信號(hào)Bn基帶數(shù)字信號(hào)bn2|r-0-2L0200400差分信號(hào)dnI'j■■\l9n\卩I-■I)DPSK調(diào)制信號(hào)mod20-2020040012「-10JULKfl02004000-20200400基帶數(shù)字信號(hào)bn20差分信號(hào)dnDPSK調(diào)制信號(hào)mod0200400202004002「0-102004000-1-20200400基帶數(shù)字信號(hào)bn2V差分信號(hào)基帶數(shù)字信號(hào)bn2V差分信號(hào)dn2211^1[lil[iblin0:川兒』1mmv.-202004000-2*-0200400DPSK調(diào)制信號(hào)mod6.2系統(tǒng)抗噪性能的測(cè)試在前面調(diào)制解調(diào)正確性的測(cè)試中，設(shè)定信噪比SNR=20?下面通過檢驗(yàn)在不同信噪比值下，系統(tǒng)調(diào)制解調(diào)的正確性，來簡單評(píng)定系統(tǒng)的抗噪性能。SNR=10時(shí)，正確解調(diào)0"-2020040012-10■V20040000-2200400SNR=5時(shí)，正確解調(diào)-202差分信號(hào)dn02004000訓(xùn)加DPSK調(diào)制信號(hào)mod2I11—12-20200400-10■P十、■J卜02004000-2200400SNR=2時(shí)，正確解調(diào)20基帶數(shù)字信號(hào)bn差分信號(hào)dn20DPSK調(diào)制信號(hào)mod2Y??i\;■Ji'iliJ.j-r-,c*人〕；V￥U￥'H'f'l■r,UU.if!i[\0-2L-2-200200400過LPF后的信號(hào)mt2004000200400經(jīng)抽樣判決后信號(hào)Bn21-10uf45八02004000-2200400SNR=1.5時(shí)，有誤碼出現(xiàn)基帶數(shù)字信號(hào)bn20■111U20差分信號(hào)dn21,1>J\，!|1ARr■rifi?■h11J.1.11UUmI)00200400DPSK調(diào)制信號(hào)mod-220200400-202004001經(jīng)抽樣判決后信號(hào)Bn2nr0-1200400延時(shí)相乘后的信號(hào)demod10-1r'f一"尸.Aftl1JIp過LPF后的信號(hào)mt02004000——-20200400所以，可以定性判斷該系統(tǒng)的抗噪性能良好。7總結(jié)經(jīng)過兩周的MATLAB課程設(shè)計(jì)，我收獲多多?？偨Y(jié)后主要有一下幾點(diǎn)可談1）良好編程習(xí)慣的培養(yǎng)我們編程時(shí)要有良好的風(fēng)格，源代碼的邏輯簡明清晰，易讀易懂是好程序的重要標(biāo)準(zhǔn)。我的幾點(diǎn)小收獲是：標(biāo)識(shí)符命名盡量含義鮮明，適當(dāng)?shù)淖⒔?，注意程序的視覺組織。2）信息篩查能力的訓(xùn)練拿到一個(gè)課程設(shè)計(jì)題目后，應(yīng)該從何下手，如何獲取有效資源、篩查信息，是對(duì)個(gè)人能力的一種考驗(yàn)。經(jīng)過這次試驗(yàn)后，我認(rèn)為首先需要把理論知識(shí)弄透徹（尤其是課本上的），再閱讀別人的設(shè)計(jì)報(bào)告，重復(fù)別人所做的工作是一種不錯(cuò)學(xué)習(xí)方法，但是需要懂得去粗取精，提取有效信息，否則會(huì)浪費(fèi)大量時(shí)間。3）調(diào)試思想的建立這次試驗(yàn)中我被很多問題難住，卻常常不知如何下手找出錯(cuò)誤，經(jīng)老師指點(diǎn)后才懂得要按一定方法進(jìn)行調(diào)試糾錯(cuò)，結(jié)果這些錯(cuò)誤常常是一些拼寫錯(cuò)誤或符號(hào)不對(duì)稱等小錯(cuò)。明確整體的架構(gòu)，有針對(duì)地進(jìn)行測(cè)試。調(diào)試方法多種多樣，但調(diào)試思想的建立卻是最為基礎(chǔ)、根本。8心得體會(huì)課程實(shí)驗(yàn)是培養(yǎng)學(xué)生綜合運(yùn)用所學(xué)知識(shí)，發(fā)現(xiàn)，提出,分析和解決實(shí)際問題,鍛煉實(shí)踐能力的重要環(huán)節(jié),是對(duì)學(xué)生實(shí)際工作能力的具體訓(xùn)練和考察過程。在這兩周的時(shí)間里我切切實(shí)實(shí)的學(xué)到了許多知識(shí)，尤其是提高了在MATLAB中

的編程能力?；仡櫰鸫舜瓮ㄐ旁碚n程實(shí)驗(yàn)，至今我仍感慨頗多，在整整兩周的時(shí)間里，我的心情和課程設(shè)計(jì)的進(jìn)展起起伏伏，一度懷疑自己的能力開始否定自己，直到最后把每處都弄通了、達(dá)到完美了，心情才真正舒暢開來。但是，我確實(shí)學(xué)到很多很多的的東西，同時(shí)鞏固了所學(xué)知識(shí)。整個(gè)設(shè)計(jì)的過程，可以說得是困難重重，這畢竟算是自己第一次獨(dú)立完成，在設(shè)計(jì)的過程中發(fā)現(xiàn)了自己的不足之處，對(duì)軟件不熟悉，動(dòng)手能力比較欠缺。設(shè)計(jì)的過程很糾結(jié)，調(diào)試的過程也很痛苦，但是知識(shí)就是在這個(gè)過程中學(xué)習(xí)的，經(jīng)驗(yàn)也是在這個(gè)過程中積累的。經(jīng)歷過困難才會(huì)有更多的收獲。掌握正確的方法，保持認(rèn)真的作風(fēng)和平和的心態(tài)，這樣才能真正的學(xué)到知識(shí)。初始參量定義這次課程實(shí)驗(yàn)終于完成了，在設(shè)計(jì)中遇到了很多問題，多謝李曉峰老師的耐心講解，從老師的身上我學(xué)到了很多實(shí)用的知識(shí)，也由衷地慶幸自己能遇到這樣一位好老師，不僅專業(yè)知識(shí)扎實(shí)，而且品格高尚、謙和有禮，這就是我理想中的好老師。再次感謝李老師。同時(shí)，對(duì)給予我?guī)椭耐瑢W(xué)們也表示忠心的感謝。初始參量定義附：完整代碼%%%%%%%%%%%%DPSK調(diào)制解調(diào)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fc=1800;%%%%%%%%%%%%fc=1800;Rb=1200;Tb=1/Rb;T=50;M=8;N=9;sign=zeros(1,M);sign=randint(1,N);SNR=1.5;%載波頻率%基帶數(shù)據(jù)傳輸速率%每個(gè)Tb內(nèi)進(jìn)行的增樣點(diǎn)數(shù)%原始碼字序列的長度%差分碼字序列的長度%產(chǎn)生長度為M的隨機(jī)二進(jìn)制序列%信噪比%%%%%%%%%%%%對(duì)sign進(jìn)行T點(diǎn)增樣后的基帶數(shù)字信號(hào)%%%%%%%%%%bn=[];forn=1:Mifsign(n)==0;F=zeros(1,T);elseF=ones(1,T);end

bn=[bn,F];end%%%%%%%%%%%%%差分編碼得到差分碼cfm%%%%%%%%%%cfm=[];cfm(1)=1;%設(shè)。壬皿的初值為1fori=1:Mcfm(i+1)=xor(sign(i),cfm(i));end%%%%%%%%%%%%%對(duì)cfm進(jìn)行T點(diǎn)增樣%%%%%%%%%%%%%%dn=[];forn=1:N;ifcfm(n)==0;A=zeros(1,T);else%差分信號(hào)%每個(gè)symbol進(jìn)行T點(diǎn)增樣A=ones(1,T);enddn=[dnA];end%%%%%%%%%%%%%%%%t=Tb/T:Tb/T:Tb;dpsK調(diào)制%%%%%%%%%%%%%%%%對(duì)Tb=1/Rb間隔進(jìn)行T點(diǎn)劃分mod=cos(2*fc*pi*t);delta=0;%調(diào)制后的信號(hào)%相位差forn=1:Mifsign(n)==1;delta=delta+pi;elsedelta=delta;endc2=cos(2*fc*pi*t+delta);mod=[modc2];end%%%%%%%%%%%調(diào)制信號(hào)過高斯白噪聲的信道%%%%%%%modz=awgn(mod,SNR);%^AWGN信道%%%%%%%%%%%BPF%%%%%%%%%%%%%modz_broad=[modzmodz((N-1)*T+1:N*T)];fs=Rb*T;hBPF=fir1(100,2*pi*[fc-Rbfc+Rb]/fs);rO=filter(hBPF,1,modz_broad);r=rO(T+1:(N+1)*T);%對(duì)r0進(jìn)行截取%%%%%%%%%%%延時(shí)相乘%%%%%%%r_shift=r(T+1:N*T);demod=[];fori=1:M*T;a=r_shift(i).*r(i);demod=[demoda];end%%%%%%%%%%%%%LPF%%%%%%%%%%demod_broad=[demoddemod((M-1)*T+1:M*T)]hLPF=fir1(100,2*pi*[Rb]/fs);mtO=filter(hLPF,1,demod_broad);mt=mt0(T+1:(M+1)*T);%^mt0進(jìn)行截取%%%%%%%%%%%%%抽樣判決%%%%%%%%%%Bn=zeros(1,M*T);fori=1:M;ifsum(mt((i-1)*T+1:i*T))>0;Bn((i-1)*T+1:i*T)=0;elseBn((i-1)*T+1:i*T)=1;endend%%%%%%%%%%%%畫圖%%%%%%%%%%%%%%%%%%%%%%%figured)%%%%%%%%%%figure(1);subplot(3,3,1);plot(bn);gridon;axis([0T*N-22]);title('基帶數(shù)字信號(hào)bn');figure(1);subplot(3,3,2);plot(dn);gridon;axis([0T*N-22]);title('差分信號(hào)dn');figure(1);subplot(3,3,3);plot(mod);gridon;axis([0T*N-22]);title('DPSK調(diào)制信號(hào)mod');figure(1);subplot(3,3,4);plot(modz);gridonaxis([0T*N-22]);title('過AWGN信道后的modz信號(hào)');figure(1);subplot(3,3,5);plot(r);gridonaxis([0T*N-22]);title('過BPF后的信號(hào)r')figure(1);subplot(3,3,6);plot(demod);gridonaxis([0T*N-11]);title('延時(shí)相乘后的信號(hào)demod')figure(1);subplot(3,3,7);plot(mt);gridonaxis([0T*N-11]);title('過LPF后的信號(hào)mt');figure(1);subplot(3,3,8);plot(Bn);gridonaxis([0T*N-22]);title('經(jīng)抽樣判決后信號(hào)Bn')%%%%%%%%%%%%figure(2)%%%%%%%%%%figure(2);subplot(3,3,1);plot(abs(fft(bn)));axis([0T*N050]);title('原始信號(hào)bn頻譜');figure(2);subplot(3,3,2);plot(abs(fft(mod)));axis([OT*N050]);title('DPSK信號(hào)mod頻譜');figure(2);subplot(3,3,3);