通信原理i課程設(shè)計(jì)報(bào)告

1、<p><b>　　北京信息科技大學(xué)</b></p><p><b>　　通信原理I課程設(shè)計(jì)</b></p><p>　　題 目: 模擬通信系統(tǒng)（PM調(diào)制）Matlab仿真平臺的</p><p><b>　　設(shè)計(jì)和實(shí)現(xiàn)</b></p><p>　　學(xué) 院

2、： 信息與通信工程學(xué)院 </p><p>　　專 業(yè)： 通信工程 </p><p>　　學(xué)生姓名： 班級/學(xué)號 </p><p>　　指導(dǎo)老師：

3、 </p><p>　　起止時間： 2015年1月 </p><p>　　《通信原理I課程設(shè)計(jì)》任務(wù)書</p><p><b>　　目錄</b></p><p>　　概述--------------------------------------------

4、-----4</p><p>　　設(shè)計(jì)目的與要求---------------------------------------4</p><p>　　設(shè)計(jì)內(nèi)容---------------------------------------------4</p><p>　　設(shè)計(jì)原理--------------------------------------------

5、-5</p><p>　　實(shí)驗(yàn)結(jié)果---------------------------------------------9</p><p>　　心得與體會------------------------------------------11</p><p>　　參考文獻(xiàn)--------------------------------------------12

6、</p><p>　　附錄：源代碼----------------------------------------12</p><p><b>　　概述</b></p><p>　　本次課程設(shè)計(jì)是對通信原理所學(xué)的最主要知識，即模擬或數(shù)字通信系統(tǒng)的各種調(diào)制方式的設(shè)計(jì)仿真實(shí)現(xiàn)。</p><p><b>　　二．設(shè)計(jì)

7、目的和要求</b></p><p>　　通信原理I課程設(shè)計(jì)的目的是為了使學(xué)生加深對所學(xué)的通信原理知識的理解，比較扎實(shí)地掌握通信原理的基礎(chǔ)知識和基本理論，增強(qiáng)分析問題和解決問題的能力，培養(yǎng)學(xué)生的專業(yè)素質(zhì)，提高其利用通信原理知識處理通信系統(tǒng)問題的能力，為今后專業(yè)課程的學(xué)習(xí)、畢業(yè)設(shè)計(jì)和工作打下良好的基礎(chǔ)。</p><p>　　本課設(shè)要求，學(xué)生根據(jù)所學(xué)知識獨(dú)立完成基本設(shè)計(jì)任務(wù)；經(jīng)老師審

8、核同意并在條件允許的情況下，可以自行命題。本課程設(shè)計(jì)以上機(jī)為主，大部分時間由學(xué)生上機(jī)操作，必要時配合少量的理論講授。</p><p><b>　　設(shè)計(jì)內(nèi)容</b></p><p>　　輸入：輸入模擬信號，例如正弦型單音頻信號等，至少選擇兩種模擬信號。給出其時域波形和功率譜密度。將此模擬信號數(shù)字化（PCM編碼），得到數(shù)字信號，接著進(jìn)行PCM解碼，將編碼前和解碼后的信號波

9、形進(jìn)行比較，驗(yàn)證所作PCM編碼的正確性。</p><p>　　調(diào)制：對輸入的模擬信號進(jìn)行DSB、SSB、PM調(diào)制，給出調(diào)制后信號的時域波形和功率譜密度。</p><p>　　信道：假定信道屬于加性高斯信道，或自行設(shè)計(jì)。自行設(shè)計(jì)可以加分。</p><p>　　解調(diào)： DSB、SSB、PM解調(diào)，仿真獲得該系統(tǒng)的輸出波形，并得到該模擬傳輸系統(tǒng)的性能指標(biāo)，即該系統(tǒng)的輸出信噪

10、比隨輸入信噪比的變化曲線。</p><p><b>　　設(shè)計(jì)原理</b></p><p><b>　　PCM編碼</b></p><p>　　PCM就是把一個時間連續(xù),取值連續(xù)的模擬信號變換成時間離散,取值離散的數(shù)字信號后在信道中傳輸.脈沖編碼調(diào)制就是對模擬信號先抽樣,再對樣值幅度量化,編碼的過程.</p>

11、<p>　　所謂抽樣,就是對模擬信號進(jìn)行周期性掃描,把時間上連續(xù)的信號變成時間上離散的信號.該模擬信號經(jīng)過抽樣后還應(yīng)當(dāng)包含原信號中所有信息,也就是說能無失真的恢復(fù)原模擬信號.它的抽樣速率的下限是由抽樣定理確定的.設(shè)一個頻帶限制的（0，fH）Hz內(nèi)的時間連續(xù)信號m（t）如果它不少于2fH次每秒的速率進(jìn)行抽樣，則m(t)可以由抽樣值完全確定。</p><p>　　抽樣定理指出，由樣值序列無失真恢復(fù)原信號的條

12、件是f S≥2 f h ，為了滿足抽樣定理，要求模擬信號的頻譜限制在0～f h之內(nèi)（fh為模擬信號的最高頻率）。為此，在抽樣之前，先設(shè)置一個前置低通濾波器，將模擬信號的帶寬限制在fh以下，如果前置低通濾波器特性不良或者抽樣頻率過低都會產(chǎn)生折疊噪聲。抽樣頻率小于2倍頻譜最高頻率時，信號的頻譜有混疊。抽樣頻率大于2倍頻譜最高頻率時，信號的頻譜無混疊。另外要注意的是，采樣間隔的 周期要足夠的小，采樣率要做夠的大，要不然會出現(xiàn)如下圖所示的混疊現(xiàn)

13、象，一般情況下TsWs=2π,Wn>2Wm。在該實(shí)驗(yàn)中,抽樣速率采用8Kbit/s. </p><p>　　所謂量化,就是把經(jīng)過抽樣得到的瞬時值將其幅度離散,即用一組規(guī)定的電平,把瞬時抽樣值用最接近的電平值來表示. 一個模擬信號經(jīng)過抽樣量化后,得到已量化的脈沖幅度調(diào)制信號,它僅為有限個數(shù)值. </p><p>　　所謂編碼,就是用一組二進(jìn)制碼組來表示每一個有固定電平的量化值.然而,實(shí)

14、際上量化是在編碼過程中同時完成的,故編碼過程也稱為模/數(shù)變換,可記作A/D.</p><p>　　話音信號先經(jīng)防混疊低通濾波器，進(jìn)行脈沖抽樣，變成8KHz重復(fù)頻率的抽樣信號（即離散的脈沖調(diào)幅PAM信號），然后將幅度連續(xù)的PAM信號用“四舍五入”辦法量化為有限個幅度取值的信號，再經(jīng)編碼后轉(zhuǎn)換成二進(jìn)制碼。對于電話，CCITT規(guī)定抽樣率為8KHz，每抽樣值編8位碼，即共有2∧8=256個量化值，因而每話路PCM編碼后的

15、標(biāo)準(zhǔn)數(shù)碼率是64kb/s。為解決均勻量化時小信號量化誤差大，音質(zhì)差的問題，在實(shí)際中采用不均勻選取量化間隔的非線性量化方法，即量化特性在小信號時分層密，量化間隔小，而在大信號時分層疏，量化間隔大。 　　</p><p>　　在實(shí)際中使用的是兩種對數(shù)形式的壓縮特性：A律和U律，A律編碼主要用于30/32路一次群系統(tǒng)，U律編碼主要用于24路一次群系統(tǒng)。A律PCM用于歐洲和中國，U律PCM用于北美和日本。</p&g

16、t;<p>　　PCM幀結(jié)構(gòu)，一個復(fù)幀包括16個幀，一個幀為125μs，分為32個時隙，其中偶幀的零時隙傳輸同步信息碼0011011，奇幀的零時隙傳輸對告碼，16時隙傳輸信令信息，其它各時隙傳輸數(shù)據(jù)，每個時隙傳輸8比特?cái)?shù)據(jù)。 　　可采用u率或者是A率進(jìn)行編碼。我國采用的是A率13折線編碼。</p><p><b>　　2.PM調(diào)制與解調(diào)</b></p><p

17、><b>　　調(diào)相信號</b></p><p>　　在模擬調(diào)制中，一個連續(xù)波有三個參數(shù)可以用來攜帶信息而構(gòu)成已調(diào)信號。當(dāng)幅度和頻率保持不變時，改變載波的相位使之隨未調(diào)信號的大小而改變，這就是調(diào)相的概念。</p><p>　　角度調(diào)制信號的一般表示形式為：</p><p>　　S (t)=Acos[ωt+φ(t)]</p>&

18、lt;p>　　式中，A是載波的恒定振幅；[ωt+φ(t)]是信號的瞬時相位，而φ(t)稱為瞬時相位偏移；d[ωt+φ(t)]/dt為信號的瞬時頻率，而dφ(t)/dt稱為瞬時頻率偏移，即相對于ω的瞬時頻率偏移。</p><p>　　設(shè)高頻載波為u=Ucosωt，調(diào)制信號為UΩ(t)，則調(diào)相信號的瞬時相位</p><p>　　φ(t)=ω+KUΩ(t)</p><p

19、>　　瞬時角頻率 ω(t)==ω+K</p><p>　　調(diào)相信號 u=Ucos［ωt+KuΩ(t)］ </p><p>　　將信號的信息加在載波的相位上則形成調(diào)相信號，調(diào)相的表達(dá)式為：</p><p>　　S(t)=Acos[ωt+Kf(t)+φ]</p><p>　　這里K稱為相移指數(shù)，這種調(diào)制方式，載波的幅度和角頻率不變，

20、而瞬時相位偏移是調(diào)制信號f(t)的線性函數(shù)，稱為相位調(diào)制。</p><p>　　調(diào)相與調(diào)頻有著相當(dāng)密切的關(guān)系，我們知道相位與頻率有如下關(guān)系式：</p><p>　　ω==ω+Kf(t)</p><p><b>　　φ(t)=ωt+K</b></p><p>　　所以在調(diào)相時可以先將調(diào)制信號進(jìn)行微分后在進(jìn)行頻率調(diào)制，這樣等

21、效于調(diào)相，此方法稱為間接調(diào)相，與此相對應(yīng)，上述方法稱為直接調(diào)相。調(diào)相信號的產(chǎn)生如圖所示：</p><p><b>　　PM調(diào)相信號的產(chǎn)生</b></p><p><b>　　調(diào)制原理</b></p><p>　　實(shí)現(xiàn)相位調(diào)制的基本原理是使角頻率為ω的高頻載波u(t)通過一個可控相移網(wǎng)絡(luò), 此網(wǎng)絡(luò)產(chǎn)生的相移Δφ受調(diào)制電壓uΩ

22、(t)控制, 滿足Δφ=KuΩ(t)的關(guān)系, 所以網(wǎng)絡(luò)輸出就是調(diào)相信號，可控相移網(wǎng)絡(luò)調(diào)相原理圖如圖所示：</p><p>　　可控相移網(wǎng)絡(luò)調(diào)相原理圖</p><p>　　調(diào)相信號的解調(diào)——相干解調(diào)</p><p>　　由于調(diào)相信號可以分解成同相分量與正交分量之和，因而可以采用線性調(diào)制中的相干解調(diào)法來進(jìn)行解調(diào)，如圖：</p><p><b

23、>　　PM信號相干解調(diào)</b></p><p>　　根據(jù)公式可以設(shè)調(diào)相信號 </p><p>　　并設(shè)相干載波 </p><p>　　則相乘器的輸出為 </p><p>　　經(jīng)過低通濾波器取出其低頻分量</p><p>　　再經(jīng)微分器，即得解調(diào)輸出&l

24、t;/p><p>　　可見，相干解調(diào)可以恢復(fù)原調(diào)制信號。這種解調(diào)方法需要本地載波與調(diào)制載波同步，否則將使解調(diào)信號失真。</p><p><b>　　實(shí)驗(yàn)結(jié)果</b></p><p>　　1.PCM編碼與解碼：</p><p>　　對原始的模擬基帶信號進(jìn)行PM調(diào)制</p><p><b>　　

25、對調(diào)相信號進(jìn)行解調(diào)</b></p><p><b>　　六．心得與體會</b></p><p><b>　　七．參考文獻(xiàn)</b></p><p>　　[1] 《通信原理》 ，周炯槃等，2005年11月 北京郵電大學(xué)出版社[3] 《通信原理》 ，樊昌信等，國防工業(yè)出版社，2006第6版</p>

26、<p>　　八．附錄：實(shí)驗(yàn)源代碼</p><p><b>　　close all</b></p><p>　　clear all;</p><p><b>　　clc;</b></p><p>　　Fs=100000; %抽取電壓值樣本的頻率</p><p>　　t

27、=[0:1/Fs:0.00999]; %每隔1/Fs秒抽取一個樣本,共抽取從0～0.001時間段內(nèi)的樣本點(diǎn),用來組成"模擬信號"</p><p>　　f=300; %該信號的頻率</p><p>　　m=cos(2*pi*f*t); %調(diào)制信號</p><p><b>　　figure(1)</b></p>

28、<p>　　subplot(3,1,1);</p><p>　　plot(t,m);grid on;</p><p>　　title('模擬信號的波形')</p><p>　　axis([0 1/f*2 -2 2]);</p><p>　　for i=1:length(m);</p><p>

29、　　%設(shè)計(jì)輸入范圍是-1.2～1.2V,對模擬信號抽取1000個樣本</p><p>　　%循環(huán)對每個樣本進(jìn)行歸一化并計(jì)算出量化單位x</p><p>　　%設(shè)量化器最大分層電平是2048個量化單位.</p><p>　　x(i)=m(i)/1.2*2048;</p><p><b>　　end</b></p>

30、;<p>　　%對每個樣本進(jìn)行編碼</p><p>　　for i=1:length(m);%length(y)=n1</p><p>　　if x(i)>0 %x(i)=A(i)</p><p>　　out(1)=1;%抽樣值大于0，極性碼為1 %code(i,1)=out()</p><p><b>　　e

31、lse</b></p><p>　　out(1)=0;%抽樣值小于0，極性碼為0</p><p><b>　　end</b></p><p><b>　　%確定段落碼</b></p><p>　　if abs(x(i))>=0 & abs(x(i))<32</p

32、><p>　　out(2)=0;out(3)=0;out(4)=0;step=1;st=0;%x的絕對值在0到32之間，則落在第一段內(nèi)，段落碼為001，段落起始電平st為16，量化間隔step為1</p><p>　　elseif 32<=abs(x(i)) & abs(x(i))<64</p><p>　　out(2)=0;out(3)=0;out

33、(4)=1;step=2;st=32;%x的絕對值在32到64之間，則落在第二段內(nèi)，段落碼為010，段落起始電平st為32，量化間隔step為2</p><p>　　elseif 64<=abs(x(i)) & abs(x(i))<128</p><p>　　out(2)=0;out(3)=1;out(4)=0;step=4;st=64;%x的絕對值在64到128之間，

34、則落在第三段內(nèi)，段落碼為011，段落起始電平st為64，量化間隔step為4</p><p>　　elseif 128<=abs(x(i)) & abs(x(i))<256</p><p>　　out(2)=0;out(3)=1;out(4)=1;step=8;st=128;%x的絕對值在128到256之間，則落在第四段內(nèi)，段落碼為100，段落起始電平st為128，量化

35、間隔step為8</p><p>　　elseif 256<=abs(x(i)) & abs(x(i))<512</p><p>　　out(2)=1;out(3)=0;out(4)=0;step=16;st=256;%x的絕對值在256到512之間，則落在第五段內(nèi)，段落碼為101，段落起始電平st為256，量化間隔step為16</p><p>

36、;　　elseif 512<=abs(x(i)) & abs(x(i))<1024</p><p>　　out(2)=1;out(3)=0;out(4)=1;step=32;st=512;%x的絕對值在512到1024之間，則落在第六段內(nèi)，段落碼為110，段落起始電平st為512，量化間隔step為32</p><p>　　elseif 1024<=abs(x(i

37、)) & abs(x(i))<2048</p><p>　　out(2)=1;out(3)=1;out(4)=0;step=64;st=1024;%x的絕對值在1024到2048之間，則落在第七段內(nèi)，段落碼為111，段落起始電平st為1024，量化間隔step為64</p><p><b>　　else</b></p><p>　

38、　out(2)=1;out(3)=1;out(4)=1;step=64;st=1024;%x的絕對值超出了2048，則段落碼為111，段落起始電平st為1024，量化間隔step為64</p><p><b>　　end</b></p><p><b>　　%確定段內(nèi)碼</b></p><p>　　if (abs(x(i)

39、)>=2048)</p><p>　　out(2:8)=[ 1 1 1 1 1 1 1];%x的絕對值超出了2048，段落碼和段內(nèi)碼均為1111111</p><p><b>　　else</b></p><p>　　tmp=floor((abs(x(i))-st)/step);%確定x落在某段的第tmp級內(nèi)</p><

40、;p>　　t=dec2bin(tmp,4)-48;%將tmp轉(zhuǎn)化為四位二進(jìn)制碼</p><p>　　out(5:8)=t(1:4);%段內(nèi)碼為t</p><p><b>　　end</b></p><p>　　a(i*8-7:i*8)=out(1:8);</p><p><b>　　end</b&g

41、t;</p><p><b>　　%輸出PCM編碼</b></p><p>　　t=0:0.001:799.999;</p><p>　　y=a(ceil(10*t+0.01));</p><p><b>　　figure(1)</b></p><p>　　subplot(3

42、,1,2)</p><p>　　plot(t,y);</p><p>　　axis([2 4 -0.2 1.2]);</p><p><b>　　grid on;</b></p><p>　　title('PCM編碼后的信號(局部)');</p><p>　　%對該P(yáng)CM編碼進(jìn)行反

43、PCM編碼</p><p>　　%%抽取其中的PCM編碼序列</p><p>　　for i=1:length(m)*8</p><p>　　pcm(i)=y(i*100-50);</p><p><b>　　end</b></p><p>　　%循環(huán)將PCM編碼變成量化電平數(shù)</p>

44、<p>　　for i=1:length(m);</p><p><b>　　%初始化計(jì)算參數(shù)</b></p><p><b>　　dlm=0;</b></p><p><b>　　dp=0;</b></p><p>　　%將2-4位段落碼換算成起始電平</

45、p><p>　　if pcm(i*8-6)==1;</p><p>　　dlm=dlm+4;</p><p><b>　　end</b></p><p>　　if pcm(i*8-5)==1;</p><p>　　dlm=dlm+2;</p><p><b>　　en

46、d</b></p><p>　　if pcm(i*8-4)==1;</p><p>　　dlm=dlm+1;</p><p><b>　　end</b></p><p>　　%得到由段落碼換算來的起始電平</p><p>　　dp=2^(dlm+4);</p><p

47、>　　%先加上最小量化電平的一半</p><p>　　dp=dp+2^(dlm)/2;</p><p>　　%將5-8位段內(nèi)碼結(jié)合段落碼一起轉(zhuǎn)換成電平數(shù)</p><p><b>　　%第5位</b></p><p>　　if pcm(i*8-3)==1;</p><p>　　dp=dp+2^

48、(dlm+3);</p><p><b>　　end</b></p><p><b>　　%第6位</b></p><p>　　if pcm(i*8-2)==1;</p><p>　　dp=dp+2^(dlm+2);</p><p><b>　　end</b&

49、gt;</p><p><b>　　%第7位</b></p><p>　　if pcm(i*8-1)==1;</p><p>　　dp=dp+2^(dlm+1);</p><p><b>　　end</b></p><p><b>　　%第8位</b>&

50、lt;/p><p>　　if pcm(i*8)==1;</p><p>　　dp=dp+2^(dlm);</p><p><b>　　end</b></p><p>　　%最后判斷第1位碼,決定正負(fù)</p><p>　　if pcm(i*8-7)==0;</p><p><

51、;b>　　dp=-dp;</b></p><p><b>　　end</b></p><p>　　jiedp(i)=dp;</p><p><b>　　end</b></p><p>　　%將量化電平數(shù)變成電壓</p><p>　　for i=1:lengt

52、h(m);</p><p>　　dianya(i)=jiedp(i)/2048*1.2;</p><p><b>　　end</b></p><p>　　Fs=100000; %抽取電壓值樣本的頻率</p><p>　　t3=[0:1/Fs:0.00999]; %每隔1/Fs秒抽取一個樣本,共抽取從0～0.01時間段內(nèi)

53、的樣本點(diǎn),用來組成"模擬信號"</p><p>　　f=300; %該信號的頻率</p><p>　　figure(1);subplot(3,1,3);plot(t3,dianya);grid on</p><p>　　axis([0 1/f*2 -2 2]);title('PCM解碼后的模擬信號')</p><

54、;p>　　Fs=100000; %抽取電壓值樣本的頻率</p><p>　　t=[0:1/Fs:0.00999]; %每隔1/Fs秒抽取一個樣本,共抽取從0～0.001時間段內(nèi)的樣本點(diǎn),用來組成"模擬信號"</p><p><b>　　y=1;</b></p><p><b>　　M=fft(m);<

55、/b></p><p>　　M=abs(M(1:length(M)/2+1));%調(diào)制信號頻譜</p><p>　　frqM=[0:length(M)-1]*Fs/length(M)/2;</p><p><b>　　Fc=30000;</b></p><p>　　m1=amod(y*m,Fc,Fs,'pm

56、'); %pm調(diào)制</p><p>　　m11=awgn(m1,4);%加噪聲</p><p>　　M1=fft(m1);</p><p>　　M1=abs(M1(1:length(M1)/2+1));</p><p>　　frqM1=[0:length(M1)-1]*Fs/length(M1)/2;</p><

57、p>　　M11=fft(m11);</p><p>　　M11=abs(M11(1:length(M11)/2+1));</p><p>　　frqM11=[0:length(M11)-1]*Fs/length(M11)/2;</p><p>　　u1=ademod(m1,Fc,Fs,'pm');%未加高斯白噪聲的解調(diào)</p>&

58、lt;p>　　u11=ademod(m11,Fc,Fs,'pm');%加高斯白噪聲解調(diào)</p><p>　　U1=fft(u1);</p><p>　　U1=abs(U1(1:length(U1)/2+1));</p><p>　　frqU1=[0:length(U1)-1]*Fs/length(U1)/2;</p><p&

59、gt;　　U11=fft(u11);</p><p>　　U11=abs(U11(1:length(U11)/2+1));</p><p>　　frqU11=[0:length(U11)-1]*Fs/length(U11)/2;</p><p><b>　　%功率譜</b></p><p><b>　　fs2=

60、800;</b></p><p>　　ts2=1/fs2;</p><p>　　t2=0:ts2:2;</p><p><b>　　nfft=64;</b></p><p>　　power=(norm(m)^2)/length(m+1);</p><p>　　spow=abs(fft(

61、m,nfft).^2);</p><p>　　f2=(0:nfft-1)/ts2/nfft;%設(shè)置頻率范圍</p><p>　　f2=f2-fs2/2;</p><p>　　power1=(norm(m1)^2)/length(m1+1);</p><p>　　spow1=abs(fft(m1,nfft).^2);</p>&l

62、t;p>　　power2=(norm(m11)^2)/length(m11+1);</p><p>　　spow2=abs(fft(m11,nfft).^2);</p><p>　　power3=(norm(u1)^2)/length(u1+1);</p><p>　　spow3=abs(fft(u1,nfft).^2);</p><p&g

63、t;　　power4=(norm(u11)^2)/length(u11+1);</p><p>　　spow4=abs(fft(u11,nfft).^2);</p><p>　　figure(2);</p><p>　　subplot(3,3,1) %子圖形式顯示結(jié)果</p><p>　　plot(t,m)

64、 %現(xiàn)在的m信號是重新構(gòu)建的信號， %因?yàn)樵趯求傅里葉變換時m=[m,zeros(1,n-n2)]</p><p><b>　　%定義兩軸的刻度</b></p><p>　　xlabel('時間t') </p><p>　　title('原調(diào)制信號的時域圖')</

65、p><p>　　subplot(3,3,2)</p><p>　　plot(frqM,M) %fftshift：將FFT中的DC分量移到頻譜中心</p><p>　　axis([-1000 50000 0 550])</p><p>　　xlabel('頻率f')</p><p>　　title('

66、原調(diào)制信號的頻譜圖')</p><p>　　subplot(3,3,3)</p><p>　　plot(f2,fftshift(spow));</p><p>　　xlabel('頻率f');title('原調(diào)制信號的功率譜');</p><p>　　disp(['power=',num

67、2str(power),'.']);</p><p>　　subplot(3,3,4)</p><p>　　plot(t,m1)</p><p>　　xlabel('時間t')</p><p>　　title('已調(diào)信號的時域圖')</p><p>　　subplot(3

68、,3,5)</p><p>　　plot(frqM1,M1)</p><p>　　axis([0 50000 0 550])</p><p>　　xlabel('頻率f')</p><p>　　title('已調(diào)信號的頻譜圖')</p><p>　　subplot(3,3,6)</

69、p><p>　　plot(f2,fftshift(spow1));</p><p>　　xlabel('頻率f');title('已調(diào)信號的功率譜');</p><p>　　disp(['power1=',num2str(power1),'.']);</p><p>　　subplo

70、t(3,3,7) </p><p>　　plot(t,m11) </p><p>　　xlabel('時間t') </p><p>　　title('加入高斯噪聲后信號的時域圖')</p><p&g

71、t;　　subplot(3,3,8)</p><p>　　plot(frqM11,M11) </p><p>　　axis([0 50000 0 550])</p><p>　　xlabel('頻率f')</p><p>　　title('信號的頻譜圖')</p><p>　　subpl

72、ot(3,3,9)</p><p>　　plot(f2,fftshift(spow2));</p><p>　　xlabel('頻率f');title('信號的功率譜');</p><p>　　disp(['power2=',num2str(power2),'.']);</p><p

73、>　　figure(3);</p><p>　　subplot(2,3,1)</p><p>　　plot(t,u1)</p><p>　　xlabel('時間t')</p><p>　　title('解調(diào)信號的時域圖（未加噪聲）')</p><p>　　subplot(2,3,2

74、)</p><p>　　plot(frqU1,U1)</p><p>　　axis([-2000 50000 0 2])</p><p>　　xlabel('頻率f')</p><p>　　title('解調(diào)信號的頻譜圖')</p><p>　　subplot(2,3,3)</p&

75、gt;<p>　　plot(f2,fftshift(spow3));</p><p>　　xlabel('頻率f');title('解調(diào)信號的功率譜');</p><p>　　disp(['power3=',num2str(power3),'.']);</p><p>　　subplot(

76、2,3,4)</p><p>　　plot(t,u11)</p><p>　　xlabel('時間t')</p><p>　　title('解調(diào)信號的時域圖（加噪聲）')</p><p>　　subplot(2,3,5)</p><p>　　plot(frqU11,U11)</p&

77、gt;<p>　　axis([-2000 50000 0 2])</p><p>　　xlabel('頻率f')</p><p>　　title('解調(diào)信號的頻譜圖')</p><p>　　subplot(2,3,6)</p><p>　　plot(f2,fftshift(spow4));<