北郵通原軟件實(shí)驗(yàn)報(bào)告.

1、信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告i信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告2013年通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告學(xué) 院：信息與通信工程學(xué)院班 級(jí)：2011211104姓 名：學(xué) 號(hào)： 班內(nèi)序號(hào):#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告目錄實(shí)驗(yàn)一調(diào)幅信號(hào)波形頻譜仿真4一. 實(shí)驗(yàn)題目4二. 基本原理41. AM調(diào)制原理 42. DSB調(diào)制原理43.SSB調(diào)制原理4三

2、. 仿真思路4四. 程序框圖5五. 仿真源代碼5六. 仿真結(jié)果及分析 8實(shí)驗(yàn)二調(diào)頻信號(hào)頻譜仿真10一. 實(shí)驗(yàn)題目10二. 基本原理10三. 仿真思路10四. 程序框圖11五. 仿真源代碼11六. 仿真結(jié)果分析 12實(shí)驗(yàn)三 單雙極性波形及功率譜仿真 13一. 實(shí)驗(yàn)題目13二. 基本原理131. 單極性歸零碼 132. 雙極性歸零碼 133. 各種碼的比較13三. 仿真思路131. 產(chǎn)生RZ碼132. 仿真功率譜密度 13四. 程序框圖14五. 仿真源代碼 14六. 仿真及結(jié)果分析 17實(shí)驗(yàn)四根升余弦滾降功率譜密度及眼圖仿真 25一. 實(shí)驗(yàn)題目25二. 基本原理251. 升余弦滾降252. 眼圖2

3、5三. 仿真思路26四. 程序框圖26五仿真源代碼26六.仿真結(jié)果及分析27心得體會(huì)294信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)一調(diào)幅信號(hào)波形頻譜仿真一. 實(shí)驗(yàn)題目假設(shè)基帶信號(hào)為m（t）用n（200二t） 2cos（100二t），載波頻率為20kHz，請仿真出AM DSB-SC SSB信號(hào)，觀察已調(diào)信號(hào)的波形及頻譜。二. 基本原理1.AM調(diào)制原理對(duì)于單音頻信號(hào)m(t) = Am sin(2二 fmt)進(jìn)行AM調(diào)制的結(jié)果為Sam 二Ac(AAm si n(2 二fmt)si n 2-f二AcA(1as in (2 二fmt)si n 2 二fct其中調(diào)幅系數(shù)A，要求a乞1以免過調(diào)引起包絡(luò)失真

4、。5信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告由Amax和Amin分別表示AM信號(hào)波形包絡(luò)最大值和最小值，則AM信號(hào)的調(diào)幅系數(shù)為Amax Amin a 二AmaxAmin#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告2.DSB-SC調(diào)制原理DSB信號(hào)的時(shí)域表達(dá)式為頻域表達(dá)式為#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告3.SSB調(diào)制原理SSB信號(hào)只發(fā)送單邊帶，比 DSB節(jié)省一半帶寬，其表達(dá)式為：認(rèn)(/) = &加(/) cos InfJ + Acm (r)sin耳期(/) = Acm (/)cos Infj - Acmt)sin

5、 Ufct三. 仿真思路定義時(shí)域采樣率、截?cái)鄷r(shí)間和采樣點(diǎn)數(shù)，可得到載波和調(diào)制信號(hào)，容易根據(jù)調(diào)制原理寫出各調(diào)制信號(hào)表達(dá)式，由此可以畫出時(shí)域波形圖。另外，對(duì)時(shí)域信號(hào)進(jìn)行FFT變換，此處使用預(yù)先定義的t2f.m函數(shù)替代，進(jìn)行傅立葉變換，得到頻譜，在頻域作圖即可。四. 程序框圖fTh確定載波fA確定f飛確定基本AM,DSB,FFT變換采樣參數(shù)1和基帶待SSB 的一得到信號(hào)作圖調(diào)制信號(hào)條信號(hào)頻譜五. 仿真源代碼1.傅里正變換函數(shù)代碼：fun ctio nS =t2f( s,fs )%傅里葉正變換,s表示輸入信號(hào)，S代表s的頻譜，fs是采樣率N=le ngth(s);%總樣點(diǎn)數(shù)T=1/fs*N;%觀察時(shí)間

6、f=-N/2:(N/2-1)/T;%頻率采樣點(diǎn)tmp1=fft(s)/fs;tmp2=N*ifft(s)/fs;S(1:N/2)=tmp2(N/2+1:-1:2);S(N/2+1:N)=tmp1(1:N/2);S=S.*exp(j*pi*f*T);End 2.傅里葉反變換函數(shù)代碼:7信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告8信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告fun ctio ns = f2t( S,fs )%傅里葉反變換N=le ngth(S);T=N/fs;t=-(T/2):1/fs:(T/2-1/fs); tmp1=fft(S)/T;tmp2=N*ifft(S)/T; s(1:N/2)=

7、tmp1(N/2+1:-1:2); s(N/2+1:N)=tmp2(1:N/2);%時(shí)域采樣點(diǎn)s=s.*exp(-j*pi*t*fs);End3.AM信號(hào)仿真代碼：clear all ;close allfs=800;T=200;%確定基本的采樣參數(shù)N=T*fs;dt=1/fs;df=1/T;t=-T/2:dt:T/2-dt;f=-fs/2:df:fs/2-df;%kHzfm1=1; fm2=0.5;fc=20; a=1/3;%設(shè)定該題目參量%調(diào)幅指數(shù)m=si n(2*pi*fm1*t)+2*cos(2*pi*fm2*t);M=t2f(m,fs);c=cos(2*pi*fc*t);s=(1+a

8、*m).*cos(2*pi*fc*t);S=t2f(s,fs);%待調(diào)制基帶信號(hào)%待調(diào)制基帶信號(hào)傅里葉變換%載波信號(hào)%AM言號(hào)%AMt號(hào)傅里葉變換figure(1)subplot(2,2,1)%做基帶待調(diào)制信號(hào)波形plot(t,m),axis(-2,2,-max(abs(m)+0.5),max(abs(m)+0.5)xlabel( t(s) ),ylabel( m(t)(V) ),title( 基帶待調(diào)制信號(hào)) subplot(2,2,2)%做載波信號(hào)波形plot(t,c),axis(-0.1,0.1,-max(abs(c)+0.5),max(abs(c)+0.5)xlabel( t(s) )

9、,ylabel( c(t)(V) ),title( subplot(2,2,3)%做 AM信號(hào)波形plot(t,s),axis(-2.5,1.5,-3,3)xlabel( t(ms) ),ylabel( s(t)(V) ),title( subplot(2,2,4)%做AM信號(hào)幅頻特性plot(f,abs(S),axis(18,22,0,max(abs(S)+10)xlabel( f(kHz) ),ylabel( |AM(f)|(V/Hz)載波信號(hào))AM時(shí)域波形),title( AM 幅頻特性)9信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告10信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告4.DSB信號(hào)仿真

10、代碼：clear all ;close all ;fs=800;T=200;%設(shè)置采樣基本參數(shù)N=T*fs;dt=1/fs;df=1/T;t=-T/2:dt:T/2-dt;f=-fs/2:df:fs/2-df;%kHzfm仁1;fm2=0.5;fc=20;%設(shè)置該題目參數(shù)m=si n(2*pi*fm1*t)+2*cos(2*pi*fm2*t); M=t2f(m,fs);%基帶待調(diào)制信號(hào)%待調(diào)制信號(hào)頻譜s=m.*cos(2*pi*fc*t);%DSB 言號(hào)S=t2f(s,fs);%DSB言號(hào)頻譜figure(1) subplot(2,1,1) plot(t,s),axis(-2.5,1.5,-4

11、,4) xlabel( t(ms) ),ylabel( s(t)(V) title( DSB-SC 時(shí)域波形) subplot(2,1,2)%故出DSB信號(hào)波形)%故出DSB信號(hào)頻譜plot(f,abs(S),axis(18,22,0,max(abs(S)+10)xlabel( f(kHz) ),ylabel( |DSB-SC(f)|(V/Hz)title( DSB-SC 幅頻特性)5.SSB信號(hào)仿真代碼:clear all ;close allfs=800;T=200;%設(shè)置基本采樣參數(shù)N=T*fs;dt=1/fs;df=1/T;t=-T/2:dt:T/2-dt;f=-fs/2:df:fs/

12、2-df;%kHzfm仁1;fm2=0.5;fc=20;%設(shè)置本題目參數(shù)11信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告m=si n(2*pi*fm1*t)+2*cos(2*pi*fm2*t);M=t2f(m,fs);MH=-j*sig n(f).*M;mh=real(f2t(MH,fs);%基帶待調(diào)制信號(hào)%基帶待調(diào)制信號(hào)傅里葉變換%基帶待調(diào)制信號(hào)的希爾伯特變換#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告%SSB信號(hào)%ssbB號(hào)的傅里葉變換s=m.*cos(2*pi*fc*t)-mh.*si n( 2*pi*fc*t);S=

13、t2f(s,fs);#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告subplot(2,1,2) plot(f,abs(S),axis(18,23,0,max(abs(S)+10) xlabel( f(kHz) ),ylabel( title( SSB幅頻特性)|SSB(f)|(V/Hz)figure(1)subplot(2,1,1) plot(t,s),axis(-3,3,-4,4) xlabel( t(ms) ),ylabel( title( SSB時(shí)域波形)%畫出SSB信號(hào)的時(shí)域波形s(t)(V)%畫出SSB信號(hào)的頻譜12信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)

14、報(bào)告六. 仿真結(jié)果及分析1.AM信號(hào)仿真結(jié)果#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告1S -1-0.50 站 11.51.載液信號(hào)-Q.0S -OH04-0.0200.02孤04 0K .080.1tfrns)#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告AM時(shí)域液冊5Qd-.2iB1營中目 a ns519S20 2052121 622側(cè)朗特性=1III100 -上圖為AM信號(hào)的仿真結(jié)果，在本題目仿真中，取調(diào)幅系數(shù)為1/3,得到的AM信號(hào)可以正確的顯示原信號(hào)的包絡(luò)特性，可以通過AM信號(hào)通過包絡(luò)檢波恢復(fù)原信號(hào)。從幅頻特性上可以看出，AM信號(hào)不僅含有基帶信號(hào)的邊頻

15、分量，還含有離散大載波分量，且載波分 量的幅值很大，說明調(diào)制效率很低。#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告2.DSB信號(hào)仿真結(jié)果13信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告DSSSC時(shí)捕遊搭t(in&562D.59.丄1919篩HWIIJO誥哲一#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告14信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告上圖為位DSB信號(hào)的仿真結(jié)果，其時(shí)域波形無法顯示出原信號(hào)的包絡(luò)特性，但是,從頻譜可以看出，DSB信號(hào)仍攜帶有原信號(hào)的頻率分量，且不含有離散大載波分量，提 高了調(diào)制效率，沒有改變傳輸帶寬，且不能用包絡(luò)檢波解調(diào)。3.SSB信號(hào)仿真結(jié)果

16、SSBAFt域戒邢-U估50上圖為SSB信號(hào)，可以看出SSB信號(hào)的包絡(luò)也無法顯示原信號(hào)的包絡(luò)特性，不能用 包絡(luò)檢波的方式回復(fù)原信號(hào)，但是，從頻譜圖可以看出，SSB信號(hào)中不含有離散大載波分量，相較AM提高了調(diào)頻效率，同時(shí)，這里為上邊帶信號(hào)，相較DSB信號(hào)又節(jié)省了帶寬。綜合以上三種調(diào)幅信號(hào)可知，AM信號(hào)可以利用包絡(luò)展示原信號(hào)特性，可以利用包絡(luò)檢波的方式解調(diào)輸出，但是由于具有離散大載波分量，調(diào)制效率低；DSB信號(hào)不能用包絡(luò)檢波的方式解調(diào)，但是由于去掉了離散大載波分量，因此，調(diào)制效率提高；SSB信號(hào)也不能用包絡(luò)檢波的方式解調(diào)，但是由于去掉了離散大載波且只取了一個(gè)邊帶，因此， 調(diào)制效率高且節(jié)省啦帶寬。實(shí)

17、驗(yàn)二調(diào)頻信號(hào)波形頻譜仿真一. 實(shí)驗(yàn)題目假設(shè)基帶信號(hào) m（t） =sin（200 二t） 2cos（100二t） 4sin（50 二t 二 /3），載波頻率為40kHz，仿真產(chǎn)生FM信號(hào)，觀察波形與頻譜，并與卡松公式做對(duì)照。FM的頻率偏移常數(shù)為 5kHz/V。二. 基本原理單音頻信號(hào)m(t)二 acos(2二 fmt)經(jīng)FM調(diào)制后的表達(dá)式為SFM(t) *cos2二以：(t)其中aKf門:(t)-sin 2二 fmt = : sin 2二 fmtfm調(diào)制指數(shù)fm。由卡松公式可知 FM信號(hào)的帶寬為B 2C1)fm三. 仿真思路同實(shí)驗(yàn)一中相仿，定義必要的仿真參數(shù)，在此基礎(chǔ)上可得到載波信號(hào)和調(diào)制信號(hào)。

18、根t據(jù).；m(t)dt可得到頻偏，由此可寫出最終的fm信號(hào)的表達(dá)式進(jìn)行仿真計(jì)算。對(duì)FM信號(hào)進(jìn)行傅里葉變換可得頻譜特性，變換依舊使用實(shí)驗(yàn)一中給出的t2f.m函數(shù)。四. 程序框圖f確定基本確定載波確定FMFFT變換采樣參數(shù)和基帶待信號(hào)得到信號(hào)作圖調(diào)制信號(hào)頻譜五. 仿真源代碼all ;% 設(shè)置基本采樣參數(shù)%設(shè)置本題目參數(shù)clear all ;closefs=800;T=16;N=T*fs;dt=1/fs; t=-T/2:dt:T/2-dt;df=1/T; f=-fs/2:df:fs/2-df;%kHzfm1= 1;fm2=0.5;fm3=0.25;fc=40;Kf=5;thea=pi/3;%基帶%

19、計(jì)算調(diào)頻產(chǎn)m=si n(2*pi*fm1*t)+2*cos(2*pi*fm2*t)+4*si n( 2*pi*fm3*t+thea);待調(diào)制信號(hào)phi=2*pi*Kf*cumsum(m)*dt;生的相位偏移s=cos(2*pi*fc*t+phi);%FM言號(hào)S=t2f(s,fs);%FM言號(hào)的傅里葉變換17信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告figure(1)subplot(2,1,1)%畫出時(shí)域的調(diào)頻信號(hào)plot(t,s),axis(0,5,-1.2,1.2)subplot(2,1,2)%畫出調(diào)頻信號(hào)的頻譜plot(f,abs(S),axis(0,9

20、0,0,max(abs(S)+0.2)六. 仿真結(jié)果及分析上圖為FM信號(hào)的仿真結(jié)果圖，從時(shí)域波形圖中可以看出，F(xiàn)M信號(hào)時(shí)域波形的頻率疏 密變化。從FM的頻譜圖可以看出，F(xiàn)M信號(hào)的帶寬約為：18信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告W 78kHz-8kHz=70kHz。又由題可知，理論計(jì)算為：fm 4kHz，頻偏為-:fmax 6V 5kHz/V 30kHz。利用卡松公示進(jìn)行B 2(1)fm = 2( fmax fm)2(5 6 4) = 68kHz可以看出理論計(jì)算值可實(shí)際仿真結(jié)果基本相同，驗(yàn)證了卡松公式的有效性對(duì)比試驗(yàn)一，二可以看出，F(xiàn)M以犧牲帶框的代價(jià)換取了高的調(diào)制效率。實(shí)驗(yàn)三單雙極性歸零碼

21、波形及功率譜仿真.實(shí)驗(yàn)題目通過仿真測量占空比為25%、50%、75%以及100%的單雙極性歸零碼波形及其功率譜。二. 基本原理1.單極性歸零碼當(dāng)發(fā)1碼時(shí)，發(fā)出正電流，但持續(xù)時(shí)間短于一個(gè)碼元的時(shí)間寬度，即發(fā)出一個(gè)窄脈沖;當(dāng)發(fā)0碼時(shí)，仍然不發(fā)送電流。單極性歸零碼在符號(hào)等概出現(xiàn)且互不相關(guān)的情況下，功率譜主瓣寬度為2尺，其頻譜含有連續(xù)譜、直流分量、離散始終分量及其奇次諧波分量。2.雙極性歸零碼其中1碼發(fā)正的窄脈沖，0碼發(fā)負(fù)的窄脈沖，兩個(gè)碼元的時(shí)間間隔可以大于每一個(gè) 窄脈沖的寬度，取樣時(shí)間是對(duì)準(zhǔn)脈沖的中心。雙極性歸零碼在符號(hào)等概且不相關(guān)的情況下，功率譜僅含有連續(xù)譜，其主瓣寬度為2RbO20信息與通信工程

22、學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告3. 各種碼的比較不歸零碼(None Return Zero Code)在傳輸中難以確定一位的結(jié)束和另一位的開始，需要 用某種方法使發(fā)送器和接收器之間進(jìn)行定時(shí)或同步。歸零碼(None Return Zero Code)的脈沖較窄，根據(jù)脈沖寬度與傳輸頻帶寬度成反比的關(guān) 系，因而歸零碼在信道上占用的頻帶較寬。單極性碼會(huì)積累直流分量；雙極性碼的直流分量大大減少，這對(duì)數(shù)據(jù)傳輸是很有利的。三. 仿真思路1.產(chǎn)生RZ碼采用歸零矩形脈沖波形的數(shù)字信號(hào)，可以用以下方法產(chǎn)生信號(hào)矢量s。設(shè)a是碼元矢量，n是總?cè)狱c(diǎn)數(shù)，M是總碼元數(shù)，L是每個(gè)碼元內(nèi)的點(diǎn)數(shù)，Rt是要求的占空比，dt是仿真系統(tǒng)的

23、時(shí)域采樣間隔，則a=(ra ndn (1,M)0);% tmp=zeros(L,M); % L1=L*tao; %taoRZ信號(hào)的產(chǎn)生方法是：先產(chǎn)生隨機(jī)的0,1序列將原本不歸零的信號(hào)按照設(shè)定占空比歸零，形成的歸零序列 為占空比tmp(1:L1,:)=o nes(L1,1)*a;2、仿真功率譜密度任意信號(hào)s(t)的功率譜的定義是PS(f)二 limTT：t其中(f)是s(t)截短后的傅氏變換,|ST(f)|2鳥/是的能量譜，t是信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告X是對(duì)應(yīng)的傅氏變換，那么X的功率譜便為&在截短時(shí)間內(nèi)的功率譜。對(duì)于仿真系統(tǒng)，若x是時(shí)域取樣值矢

24、量,P二X.*conj X /T。針對(duì)隨機(jī)過程X(t)，其平均功率譜密度定義為各樣本功率譜密度的數(shù)學(xué)期望PX(f)=+K Pk(f)K km四、程序框圖單極性歸零：確定基本產(chǎn)生0,1按照占空 比產(chǎn)生對(duì)累積求和jr%采樣參數(shù)隨機(jī)碼序應(yīng)的歸零得到功率作圖II列碼序列譜密度L丿14信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告雙極性歸零:確定基本 采樣參數(shù)產(chǎn)生0,1隨機(jī)碼序 列按照占空 比產(chǎn)生對(duì) 應(yīng)的歸零 碼序列作圖22信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告五. 仿真源代碼1單極性歸零碼仿真代碼：clear all ;close all ;L=64; %每隔數(shù)碼元間隔

25、內(nèi)的采樣點(diǎn)數(shù)N=4096; %總采樣點(diǎn)數(shù)M=N/L; %總碼元數(shù)Rs=5;Ts=1/Rs;T=M*Ts;fs=N/T;%設(shè)置基本采樣參數(shù)t=-(T/2):1/fs:(T/2-1/fs);f=-(fs/2):1/T:(fs/2-1/T);tao=1;%占空比為100%，將其改成0.25,0.5,0.75,可以求出對(duì)應(yīng)的仿真結(jié)果EP=zeros(1,N);for loop=1:10000%累積求和產(chǎn)生功率譜密度a=(ra ndn (1,M)0);tmp=zeros(L,M);L1=L*tao;tmp(1:L1,:)=o nes(L1,1)*a;s=tmp(:);S=t2f(s,fs);%樣本信號(hào)的

26、功率譜密度P=abs(SF2/T;EP=EP*(1-1/loop)+P/loop;endfigure(1)plot(t,s),axis(-5,5,-0.2,1.2)xlabel( t(ms) ),ylabel( s(t)(V) ),title( 單極性歸零碼波形(占空比為 100%)figure(2)subplot(2,1,1)plot(f,EP),axis(-30,30,0,max(EP)+0.02)xlabel( f(kHz) ),ylabel( 功率譜(W/kHz) ),title( 單極性歸零碼功率譜(占空比為100%)subplot(2,1,2)plot(f,10*log10(EP+

27、eps),axis(-20,20,-200,0)xlabel( f(kHz) ),ylabel( 功率譜(db) ),title( 單極性歸零碼功率譜(db)( 占空比為100%)2.雙極性歸零碼仿真代碼：clear all ;close all ;%設(shè)置基本采樣參數(shù)L=64; %每隔數(shù)碼元間隔內(nèi)的采樣點(diǎn)數(shù)N=2048; %總采樣點(diǎn)數(shù)M=N/L; %總碼元數(shù)Rs=5;Ts=1/Rs;T=M*Ts;fs=N/T;t=-(T/2):1/fs:(T/2-1/fs);f=-(fs/2):1/T:(fs/2-1/T);%累積求和得到功率譜tao=1; %占空比EP=zeros(1,N);for loop

28、=1:10000a=(ra nd(1,M)0.5);aa=1-2*a;tmp=zeros(L,M);L1=L*tao;tmp(1:L1,:)=o nes(L1,1)*aa;s=tmp(:);S=t2f(s,fs);%樣本信號(hào)的功率譜密度P=abs(S).A2/T;EP=EP*(1-1/loop)+P/loop;endfigure(1)plot(t,s),axis(-3,3,-1.2,1.2)%雙極性歸零碼波形xlabel( t(ms) ),ylabel( s(t)(V)gridtitle(雙極性歸零碼波形)figure(2)%雙極性歸零碼功率譜subplot(2,1,1)plot(f,EP),

29、axis(-50,50,0,max(EP)+0.02)xlabel( f(kHz) ),ylabel( 功率譜(W/kHz)title(雙極性歸零碼功率譜(占空比100%) subplot(2,1,2)plot(f,10*log10(EP+eps),axis(-20,20,-200,0)xlabel( f(kHz) ),ylabel( 功率譜(db)title( 雙極性歸零碼功率譜(db)(占空比100%)實(shí)驗(yàn)結(jié)果及分析25信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告單極性歸零碼波形占空比為1.單極性歸零碼(1 )占空比為100%時(shí)：t(ms)26信息與通信工

30、程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告單極性歸零碼功率譜占空比為100%)27信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告1.510.50-30-20-1001020RkHz)單按性歸霧碼功率譜占空比100%)30#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告可以看出占空比為 100%時(shí)，主瓣帶寬為 5kHz(2 )占空比為75%時(shí)：單極性歸零碼浪形(占空比為了5%#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告t(ms)28信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告29信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告單極性歸零碼功率ig(db)占空比為羽)f

31、(kHz)可以看出占空比為 75%時(shí)，主瓣帶寬為7kHz(3 )占空比為50%時(shí)：單極性歸零碼波形(占空比為別t(ms)30信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告31信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告單扱性歸零碼功率譜(占空比為0%)O4 3 2 1 0*3O,00.OM玉畧B掛IIRiI BIIIAi-20-100102030f(kHz)單極性歸零碼功率譜占空比為50%)可以看出占空比為 50%時(shí)，主瓣帶寬為10kHz(4 )占空比為25%時(shí):單極性歸零碼波形(占空比為茁)t(msl32信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告單極性歸零碼功率譜（占空比為2%）33信息與通信工程學(xué)院通信原理

32、軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告0-30#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告單極性歸零碼功率譜但臥占空比25%）#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告(qp) se#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告可以看出占空比為 50%時(shí)，主瓣帶寬為20kHz2雙極性歸零碼（1 ）占空比為100% :汝極性歸零職波形（占空比10 80.60402 o-0.24 4-0 64 8I- - rIII1T *11iliti|i N k hr - 111_ _ L r L - 11i1 Tl1iiT _ HNh*-3-2-10123t(ms)O

33、2 5 15 0-5O-1-O-o O肯工WQS5掛曰雙極性歸霧碼功率譜（占空比100%）40-30-20-1001020304060f(kHz)雙極性歸零碼功率譜db)(占空比100%)(qp)薛悴O0-220-16 J0 巧05101520ffkHz)35信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告36信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告可以看出當(dāng)雙極性歸零碼占空比為100%時(shí)，主瓣帶寬為 5kHz（2 ）占空比為75%#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告-202t(ms)歡極性歸零碼功率譜（占空比西鳴N工壬AA)ssaeO5 OO40-30-20-1001020304050f(kHz)雙

34、極性歸零碼功率譜(db)(占空比75%)0-2活05ffkHz)10152037信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告可以看出雙極性歸零碼占空比為75%時(shí)，主瓣帶寬約為 7kHz占空比為50%#信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告t(ms)04O02雙扱性歸零碼功率譜(占空比0%500-5040 導(dǎo) 0-20-10 D 10203040fi(kHz)雙極性歸零碼功率譜tdb) C占空比60%)043mJ-協(xié)-100150-20 d5 -10-505101520ffkHz)可以看出雙極性歸零碼占空比為75%時(shí)，主瓣帶寬約為10kHz(4)占空比為25% :38信息與通信工程學(xué)院通信原理軟件實(shí)驗(yàn)報(bào)告02雙極性歸零碼功率譜（占空比25%）懇工乂多塑(1-alpha)/(2*Ts )& abs(f)=(1+alpha)/(2*Ts);Hcos(ii)=Ts/2*(1+cos(pi*Ts/alpha*(abs(f(ii)-(1-alpha)/(2*Ts); ii=fi nd(abs(f)=(1-alpha)/(2*Ts);Hcos(ii)=Ts;%根升余弦特性Hrc