基于matlab的AMFMPM調(diào)制

1、目錄第1章前言 2第2章 AM , FM,PMM制原理 3AM調(diào)制原理 3FM調(diào)制原理 3PM調(diào)制原理 5第3章 幾種調(diào)制方式的比較 7PM與FM的比較 7幾種不同的模擬調(diào)制方式 7幾種模擬調(diào)制的性能比較 9幾種模擬調(diào)制的特點(diǎn)及應(yīng)用 10第4章AM,FM,PM勺調(diào)制仿真 10AM的調(diào)制仿真 11理想狀態(tài)下的AM調(diào)制仿真 11含噪聲情況下的AM調(diào)制仿真 12FM的調(diào)制仿真 14理想狀態(tài)下的FM調(diào)制仿真 14含噪聲情況下的FM調(diào)制仿真 15PM的調(diào)制仿真 16心得體會(huì) 17參考文獻(xiàn) 18附錄 19AM FM PM實(shí)現(xiàn)及性能比較第1章前言通信系統(tǒng)是為了有效可靠的傳輸信息，信息由信源發(fā)出，以語言、圖像

2、、數(shù)據(jù)為 媒體，通過電(光)信號(hào)將信息傳輸，由信宿接收。通信系統(tǒng)又可分為數(shù)字通信與模擬 通信。基于課程設(shè)計(jì)的要求，下面簡(jiǎn)要介紹模擬通信系統(tǒng)。信源是模擬信號(hào)，信道中傳輸?shù)囊彩悄M信號(hào)的系統(tǒng)為模擬通信。 模擬通信系統(tǒng) 的模型如圖1所示。圖i模擬通信系統(tǒng)模型調(diào)制器：使信號(hào)與信道相匹配，便于頻分復(fù)用等。發(fā)濾波器：濾除調(diào)制器輸出的 無用信號(hào)。收濾波器：濾除信號(hào)頻帶以外的噪聲，一般設(shè)N(t)為高斯白噪聲，則Ni(t) 為窄帶白噪聲。第2章AM, FM,P網(wǎng)制原理AM調(diào)制原理幅度調(diào)制是用調(diào)制信號(hào)去控制高頻正弦載波的幅度，使其按調(diào)制信號(hào)的規(guī)律變化 的過程。幅度調(diào)制器的一般模型如圖所示。哂廣、-岫 rT cos

3、curf圖幅度調(diào)制模型在圖中，若假設(shè)濾波器為全通網(wǎng)絡(luò)（H（ ）=1）,調(diào)制信號(hào)mt疊加直流A后再與載波相乘，則輸出的信號(hào)就是常規(guī)雙邊帶（AM調(diào)幅.AM調(diào)制器模型如圖所示：圖AM調(diào)制模型AM信號(hào)波形的包絡(luò)與輸入基帶信號(hào) mt成正比，故用包絡(luò)檢波的方法很容易恢 復(fù)原始調(diào)制信號(hào)。 但為了保證包絡(luò)檢波時(shí)不發(fā)生失真，必須滿足 A 1mtimax,否則 將出現(xiàn)過調(diào)幅現(xiàn)象而帶來失真。AM信號(hào)的頻譜是由載頻分量和上、下兩個(gè)邊帶組成 （通常稱頻譜中畫斜線的部分為上邊帶，不畫斜線的部分為下邊帶）。上邊帶的頻譜與原調(diào)制信號(hào)的頻譜結(jié)構(gòu)相同，下邊帶是上邊帶的鏡像。顯然，無論是上邊帶還是下 邊帶，都含有原調(diào)制信號(hào)的完整信

4、息。故 AM信號(hào)是帶有載波的雙邊帶信號(hào)，它的帶 寬信號(hào)帶寬的兩倍。FM調(diào)制原理頻率調(diào)制的一般表達(dá)式為:FMffi PWE常相似，如果預(yù)先不知道調(diào)制信號(hào)的具體形式，則無法判斷已調(diào)信號(hào) 是調(diào)頻信號(hào)還是調(diào)相信號(hào)。圖 直接調(diào)頻法m(t) | dt|PMl Sfm (t)圖間接調(diào)頻法圖所示的產(chǎn)生調(diào)頻信號(hào)的方法稱為直接調(diào)頻法，圖所示的產(chǎn)生調(diào)頻信號(hào)的方法稱 為間接調(diào)頻法4。由于實(shí)際相位調(diào)制器的調(diào)節(jié)范圍不可能超出-柏封，因而間接調(diào)頻 的方法僅適用于相位偏移和頻率偏移不大的窄帶調(diào)制情形，而直接調(diào)頻則適用于寬帶 調(diào)制情形。根據(jù)調(diào)制后載波瞬時(shí)相位偏移的大小，可將頻率調(diào)制分為寬帶調(diào)頻(WBF M與窄 帶調(diào)頻(NBFM

5、Io寬帶與窄帶調(diào)制的區(qū)分并無嚴(yán)格的界限，但通常認(rèn)為由調(diào)頻所引起 的最大瞬時(shí)相位偏移遠(yuǎn)小于300時(shí)，GA"避稱為窄帶調(diào)頻。否則，稱為寬帶調(diào)頻。為方便起見，無妨假設(shè)正弦載波的振幅式，得)d tttSFM(t) cos ct Kf m(t)d =cos ctcosKFm( )d sin csinKFm(2-3)通過化解，利用傅立葉變化公式可得 NBFMJ號(hào)的頻域表達(dá)式:(2-4)在NBFMK由于下邊頻為負(fù)，因而合成矢量不與載波同相，而是存在相位偏移 當(dāng)最大相位偏移滿足式(2-2)時(shí)，合成矢量的幅度基本不變，這樣就形成了 FM信號(hào)。圖NBFM信號(hào)頻譜PM調(diào)制原理在模擬調(diào)制中，一個(gè)連續(xù)波有三個(gè)

6、參數(shù)可以用來攜帶信息而構(gòu)成已調(diào)信號(hào)。當(dāng)幅度和頻率保持不變時(shí)，改變載波的相位使之隨未調(diào)信號(hào)的大小而改變， 這就是調(diào)相的 概念。角度調(diào)制信號(hào)的一般表示形式為：Sm(t)=AcosCt+6(t)(2-5)式中，A是載波的恒定振幅；ct+4(t)是信號(hào)的瞬時(shí)相位，而 小稱為瞬 時(shí)相位偏移；dct+小(t)/dt 為信號(hào)的瞬時(shí)頻率，而d(I)(t)/dt 稱為瞬時(shí)頻率偏 移，即相對(duì)于cdc的瞬時(shí)頻率偏移。設(shè)高頻載波為uc=Ucmcosct,調(diào)制信號(hào)為UQ(t),則調(diào)相信號(hào)的瞬時(shí)相位Ht)= ct+KpUQ (t)瞬時(shí)角頻率 (t)= ddt = 3 c+KP duQdt調(diào)相信號(hào) u PM =Ucm co

7、s ct+KPuQ(t)將信號(hào)的信息加在載波的相位上則形成調(diào)相信號(hào)，調(diào)相的表達(dá)式為：Spm (t)=Acos c t+K pm f(t)+ 小 0(2-6)這里Kpm稱為相移指數(shù)，這種調(diào)制方式，載波的幅度和角頻率不變，而瞬時(shí)相位偏移是調(diào)制信號(hào)f(t)的線性函數(shù)，稱為相位調(diào)制。調(diào)相與調(diào)頻有著相當(dāng)密切的關(guān)系，我們知道相位與頻率有如下關(guān)系式:="%t = C +KPM f(t)(|)(t)= 出 CD c t+K pm f(t)dt所以在調(diào)相時(shí)可以先將調(diào)制信號(hào)進(jìn)行微分后在進(jìn)行頻率調(diào)制，這樣等效于調(diào)相，此方法稱為間接調(diào)相，與此相對(duì)應(yīng)，上述方法稱為直接調(diào)相。調(diào)相信號(hào)的產(chǎn)生如圖所 示：RD d相

8、位調(diào)制器 仔.代)，J微分口-g頻率調(diào)制器 獸4)(0直接調(diào)相(坊間接調(diào)相圖PM調(diào)相信號(hào)的產(chǎn)生實(shí)現(xiàn)相位調(diào)制的基本原理是使角頻率為CDc的高頻載波uc(t)通過一個(gè)可控相移網(wǎng)絡(luò)，此網(wǎng)絡(luò)產(chǎn)生的相移 A小受調(diào)制電壓uQ(t)控制，滿足A(|)=KpuQ(t)的關(guān)系， 所以網(wǎng)絡(luò)輸出就是調(diào)相信號(hào)，可控相移網(wǎng)絡(luò)調(diào)相原理圖如圖所示：圖可控相移網(wǎng)絡(luò)調(diào)相原理圖第3章幾種調(diào)制方式的比較PM與FM的比較表 PM與FM的比較FM(1)瞬時(shí)頻率： okpu (t)(2)瞬時(shí)相位：t(t) ot kF 0u (t)dt o(3)最大頻偏：mkF |U |max kfU(4)取人相位：.,tUmfmkF | u (t)dt

9、 |kF 0max表達(dá)式:ufm (t) U cos (t)tU cos ot kf 0u (t)dt okfUU cos otsin t oU cos ot mf sin t oPMdu (t)(t)okpdt(t)ot kpu (t)odu (t)mkp | max k pUdtmpm kp |u (t)|max kpUupm (t) U cos (t)U cosotkpu(t)oU cosotkpUcostoU cos ot mpcos t o幾種不同的模擬調(diào)制方式假定所有調(diào)制系統(tǒng)在接收機(jī)輸入端具有相等的信號(hào)功率，且加性噪聲都是均值為0、雙邊功率譜密度為/2的高斯白噪聲，基帶信號(hào)帶寬為，

10、在所有系統(tǒng)都滿足稗=0 « 3 （r） = 卜嘰g=1例如，加（力為正弦型信號(hào)。綜合前面的分析，可總結(jié)各種模擬調(diào)制方式的信號(hào)帶寬、 制度增益、輸出信噪比、設(shè)備（調(diào)制與解調(diào)）復(fù)雜程度、主要應(yīng)用等如表所示。表中 還進(jìn)一步假設(shè)了 AM為100%（制。表調(diào)制方工（信號(hào)帶寬制度增益設(shè)備復(fù)雜度主要應(yīng)用DSB%2.%。巾等要求相干 解調(diào)，常與DSB 信號(hào)一起傳輸一 個(gè)小導(dǎo)頻點(diǎn)對(duì)點(diǎn)的專用 通信，低帶寬 信號(hào)多路復(fù)用 系統(tǒng)SSB1顯 絢力較大：要耕目干 解調(diào)，調(diào)制器也 較復(fù)雜短波無線電廣 播，話音頻分 多路通信VSB略人于九近似£犯近似4SB較大；要求相干 解調(diào)，調(diào)制器需 要時(shí)稱浦彼數(shù)據(jù)傳輸

11、F商 用電視廣播AM2 A2I1- -e*3%較小調(diào)制與解 調(diào)包絡(luò)檢波）簡(jiǎn) 單中短波無線電 廣播FM2 + DZw3曷（超廣 + 士）3.國(guó) 一德：2 '公/福中等調(diào)制器有 點(diǎn)復(fù)雜，解調(diào)器 較簡(jiǎn)單微波巾繼、超 短波小功率電 臺(tái)（窄帶）$ 衛(wèi)星逋信、調(diào) 顓立體聲廣播（寬帶）幾種模擬調(diào)制的性能比較就抗噪性能而言，WBFMDSB SSB VS以之,AMR差。NBFMf AMg近。示出了各種模擬調(diào)制系統(tǒng)的性能曲線， 圖中的圓點(diǎn)表示門限點(diǎn)。門限點(diǎn)以下，曲線迅速下跌；門限點(diǎn)以上，DSB SSB的信噪比比AM高以上，而 FM（=6）的信噪比比AM® 22dB0就頻帶利用率而言，SSB最女？

12、，VSBt SSB接近，DSB AM NBFMC之，WBFlMt幾種模擬調(diào)制的特點(diǎn)及應(yīng)用AM 調(diào)制的優(yōu)點(diǎn)是接收設(shè)備簡(jiǎn)單；缺點(diǎn)是功率利用率低，抗干擾能力差，信號(hào)帶 寬較寬，頻帶利用率不高。因此，AM制式用于通信質(zhì)量要求不高的場(chǎng)合，目前主要 用在中波和短波的調(diào)幅廣播中。DSB調(diào)制的優(yōu)點(diǎn)是功率利用率高，但帶寬與 AM相同，頻帶利用率不高，接收要 求同步解調(diào)，設(shè)備較復(fù)雜。只用于點(diǎn)對(duì)點(diǎn)的專用通信及低帶寬信號(hào)多路復(fù)用系統(tǒng)。SSB調(diào)制的優(yōu)點(diǎn)是功率利用率和頻帶利用率都較高，抗干擾能力和抗選擇性衰 落能力均優(yōu)于AM而帶寬只有AM的一半；缺點(diǎn)是發(fā)送和接收設(shè)備都復(fù)雜。 SSB制式 普遍用在頻帶比較擁擠的場(chǎng)合，如短波

13、波段的無線電廣播和頻分多路復(fù)用系統(tǒng)中。VSB調(diào)制性能與SSB相當(dāng)，原則上也需要同步解調(diào)，但在某些 VSB系統(tǒng)中，附 加一個(gè)足夠大的載波，形成（VSB+C合成信號(hào)，就可以用包絡(luò)檢波法進(jìn)行解調(diào)。這 種（VSB+C方式綜合了 AM SS刑DS屋者的優(yōu)點(diǎn)。所以VS*E數(shù)據(jù)傳輸、商用電 視廣播等領(lǐng)域得到廣泛使用。FM波的幅度恒定不變，這使得它對(duì)非線性器件不甚敏感，給FM帶來了抗快衰落能力。利用自動(dòng)增益控制和帶通限幅還可以消除快衰落造成的幅度變化效應(yīng)。這些特點(diǎn)使得NBFM寸微波中繼系統(tǒng)頗具吸引力。WBFMJ抗干擾能力強(qiáng)，可以實(shí)現(xiàn)帶寬與信 噪比的互換，因而WBFMT泛應(yīng)用于長(zhǎng)距離高質(zhì)量的通信系統(tǒng)中，如空間和

14、衛(wèi)星通信、 調(diào)頻立體聲廣播、短波電臺(tái)等。 WBFM勺缺點(diǎn)是頻帶利用率低，存在門限效應(yīng)，因此 在接收信號(hào)弱、干擾大的情況下宜采用NBFM這就是小型通信機(jī)常采用 NBFM勺原因。第4章AM,FM,PM的調(diào)制仿真AM的調(diào)制仿真理想狀態(tài)下的AM調(diào)制仿真圖調(diào)制信號(hào)的時(shí)域波形圖載械此時(shí)城沌甲3避肯號(hào)的時(shí)時(shí)域波10圖 4.3理想狀態(tài)下的已調(diào)信號(hào)的時(shí)域和頻譜圖含噪聲情況下的AM調(diào)制仿真4疊加后的調(diào)制唱號(hào)的域波形4 NIIIIII00020.040.0G0,080.10.12圖疊加噪聲以后的已調(diào)信號(hào)的的時(shí)域和頻譜圖理想狀態(tài)下的FM調(diào)制仿真國(guó)制恰勺的日a國(guó)It nr圖 調(diào)制信號(hào)、載波以及已調(diào)信號(hào)的時(shí)域波形圖圖調(diào)制

15、信號(hào)和已調(diào)信號(hào)的頻譜圖調(diào)制信號(hào)的射揖浮時(shí)機(jī)無唯小條件號(hào)的時(shí)“國(guó)Hint 元喝聲條件下己詞信號(hào)的時(shí)埔圖圖調(diào)制信號(hào)、無噪聲情況下的已調(diào)信號(hào)和解調(diào)信號(hào)時(shí)域圖含噪聲情況下的FM調(diào)制仿真小信噪比為10,大信口比為30周制哨弓r：河域困廠舊1七問圖 調(diào)制信號(hào)、含大信噪比高斯白噪聲情況下的已調(diào)信號(hào)和解調(diào)信號(hào)時(shí)域圖PM的調(diào)制仿真寫圖制信號(hào)的時(shí)域里已調(diào)信號(hào)的幀譜圖時(shí)0044-0.1-0 0500.050 1時(shí)間t題率f第調(diào)后信號(hào)的時(shí)域被而時(shí)間t第5章心得體會(huì)通過這一次課程設(shè)計(jì)，我了解很多關(guān)于專業(yè)的知識(shí)，以前每次學(xué)這些知識(shí)時(shí)，總 是不知道這些東西具體拿來有什么用， 現(xiàn)在才知道，幾個(gè)短短輸入信號(hào)，在有了一個(gè) 簡(jiǎn)單的

16、電路流程后，就能仿真成我們生活中很多常見的東西。 總的來說，這次課程設(shè) 計(jì)過程還是比較愉快輕松的，雖然中間有過一些困難，但是在老師與同學(xué)的指點(diǎn)下我 還是渡過了，在這里我要謝謝幫助我的老師和同學(xué)。參考文獻(xiàn)1 樊昌信.通信原理(第6版).國(guó)防工業(yè)出版社，2006, 092 黎洪松.數(shù)字通彳t原理.西安電子系科技大學(xué)出版社，2005, 073 任嘉偉.數(shù)字頻帶通信系統(tǒng)計(jì)算機(jī)仿真J.電腦知識(shí)與技術(shù)，2008, 074 呂躍廣通信系統(tǒng)仿真.電子工業(yè)出版社，5 席在芳等 基于SIMULINK的現(xiàn)代通信系統(tǒng)仿真分析J.系統(tǒng)仿真學(xué)報(bào)2006,18(10)AM :fm=100;fc=500;fs=5000;Am

17、=1;A=2;N=512;K=N-1;n=0:N-1;t=(0:1/fs:K/fs);yt=Am*cos(2*pi*fm*t);figure(1)subplot(1,1,1),plot(t,yt),title('調(diào)制信號(hào) f1 的時(shí)時(shí)域波');y0=A+yt ;y2=y0.*cos(2*pi*fc*n/fs);y3=fft(y2,N);% fft 變換q1=(0:N/2-1)*fs/N;mx1=abs(y3(1:N/2);figure(2)subplot(2,1,1);plot(t,y2);title('已調(diào)信號(hào)的時(shí)時(shí)域波);subplot(2,1,2);plot(q1

18、,mx1);title('f1已調(diào)信號(hào)的頻譜');yc=cos(2*pi*fc*t);figure(3)subplot(2,1,1),plot(t,yc),title('N=512;n=0:N-1;yc1=Am*cos(2*pi*fc*n/fs);y3=fft(yc1,N);q=(0:N/2-1)*fs/N;mx=abs(y3(1N2);figure(3)subplot(2,1,2),plot(q,mx),title('% 繪圖載波fc時(shí)域波形')載波fc頻譜')y4=*randn(1,length(t);% 用 RAND臚生高斯分布序列w=y4

19、.A2;%噪聲功率figure(4)subplot(2,1,1);plot(t,y4);title('高斯白噪聲時(shí)域波形')y5=fft(y4,N);q2=(0:N/2-1)*fs/N;mx2=abs(y5(1:N/2);subplot(2,1,2),plot(q2,mx2),title('高斯白噪聲頻域波形')y6=y2+y4;figure(5)subplot(2,1,1),plot(t,y6),title('疊加后的調(diào)制信號(hào)時(shí)域波形')q3=q1;mx3=mx1+mx2;subplot(2,1,2),plot(q3,mx3),title(&#

20、39;疊加后的調(diào)制信號(hào)頻譜波形')%調(diào)制yv=y6.*yc; %乘以載波進(jìn)行解調(diào)Ws=yv.A2;p1=fc-fm;k,Wn,beta,ftype=kaiserord(p1 fc,1 0, ,fs); %Fir數(shù)字低通濾波window=kaiser(k+1,beta); % 使用 kaiser 窗函數(shù) b=fir1(k,Wn,ftype,window,'noscale'); %使用標(biāo)準(zhǔn)頻率響應(yīng)的加窗設(shè)計(jì)函數(shù)yt=filter(b,1,yv);yssdb=yt.*2-2;figure(6)subplot(2,1,1),plot(t,yssdb),title('y9

21、=fft(yssdb,N);q=(0:N/2-1)*fs/N;mx=abs(y9(1:N/2);subplot(2,1,2),plot(q,mx),title('ro=y9-yt;W=(yt.A2).*(1/2);R=W/wr=W/roG=r/Rfm：dt=;%t=0:dt:;am=5;%fm=5;%mt=am*cos(2*pi*fm*t);fc=50;%經(jīng)過低通已調(diào)信號(hào)的時(shí)域波形采樣')經(jīng)過低通已調(diào)信號(hào)頻域波形')%解調(diào)設(shè)定時(shí)間步長(zhǎng)% 產(chǎn)生時(shí)間向量設(shè)定調(diào)制信號(hào)幅度設(shè)定調(diào)制信號(hào)頻率%生成調(diào)制信號(hào)設(shè)定載波頻率生成載波設(shè)定調(diào)頻指數(shù)求信號(hào)m（t）的積分調(diào)制，產(chǎn)生已調(diào)信號(hào)調(diào)制信

22、號(hào)設(shè)定信躁比（小信噪比）設(shè)定信躁比（大信噪比）設(shè)定信躁比（無信噪比）計(jì)算對(duì)應(yīng)的高斯白躁聲的方差生成高斯白躁聲生成含高斯白躁聲的已調(diào)信號(hào)（信接受信號(hào)通過微分器處理變換，求絕對(duì)值得到瞬時(shí)ct=cos(2*pi*fc*t);%kf=10;%int_mt=0;for i=1:length(t)-1int_mt(i+1)=int_mt(i)+mt(i)*dt;%end%sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt); %0甘*添力口高斯白噪聲 *sn1=10;%sn2=30;%sn=0;%db=amA2/(2*(10A(sn/10);%n=sqrt(db)*randn(size

23、(t);%nsfm=n+sfm;%號(hào)通%過信道傳輸)%*%*FM 解調(diào) *for i=1:length(t)-1%diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn = abs(hilbert(diff_nsfm); %hilbert幅度(包絡(luò)檢波)zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diffnsfmn1=diffnsfmn-zero;*時(shí)域到頻域轉(zhuǎn)換*ts=;%fs=1/ts;%df=;%換%寸，它表示FFT的最小頻率間隔%*對(duì)調(diào)制信號(hào)m(t)求傅里葉變換*m=am*cos(2*pi*fm*t);%f

24、s=1/ts;if nargin=2抽樣間隔抽樣頻率所需的頻率分辨率，用在求傅里葉變?cè){(diào)信號(hào)n1=0;elsen1=fs/df;endn2=length(m);n=2A(max(nextpow2(n1),nextpow2(n2);M=fft(m,n);m=m,zeros(1,n-n2);df1=fs/n;%傅里變換M=M/fs;%f=0:df1:df1*(length(m)-1)-fs/2;%以上程序是對(duì)調(diào)制后的信號(hào)u求縮放，便于在頻鋪圖上整體觀察時(shí)間向量對(duì)應(yīng)的頻率向量%* 對(duì)已調(diào)信號(hào)u求傅里變換*fs=1/ts;if nargin=2n1=0;elsen1=fs/df;endn2=lengt

25、h(sfm);n=2A(max(nextpow2(n1),nextpow2(n2);U=fft(sfm,n);u=sfm,zeros(1,n-n2);df1=fs/n;%以上是對(duì)已調(diào)信號(hào)u求傅里變換U=U/fs;%縮放%*%*disp('按任意鍵可以看到原調(diào)制信號(hào)、載波信號(hào)和已調(diào)信號(hào)的曲線,)pause*figure(1)*figure(1)subplot(3,1,1);plot(t,mt);xlabel('時(shí)間 t');title(' 調(diào)制信號(hào)的時(shí)域圖');subplot(3,1,2);plot(t,ct);xlabel('時(shí)間 t')

26、;title(' 載波的時(shí)域圖');subplot(3,1,3);plot(t,sfm);xlabel('時(shí)間 t');title(' 已調(diào)信號(hào)的時(shí)域圖');%繪制調(diào)制信號(hào)的時(shí)域圖%繪制載波的時(shí)域圖%繪制已調(diào)信號(hào)的時(shí)域圖%*disp（'按任意鍵可以看到原調(diào)制信號(hào)和已調(diào)信號(hào)在頻域內(nèi)的圖形'）pause*figure(2)*figure(2)將FFT中的DC分subplot(2,1,1)plot(f,abs(fftshift(M)%fftshift:量移到頻譜中心xlabel(' 頻率 f)title(' 原調(diào)制信號(hào)的

27、頻譜圖')subplot(2,1,2)plot(f,abs(fftshift(U)xlabel('頻率 f)title('已調(diào)信號(hào)的頻譜圖')%*disp('按任意鍵可以看到原調(diào)制信號(hào)、無噪聲條件下已調(diào)信號(hào)和解調(diào)信號(hào)的曲線') pause %*figure(3)*繪制調(diào)制信號(hào)的時(shí)域圖figure(3) subplot(3,1,1);plot(t,mt);xlabel('時(shí)間 t');title(' 調(diào)制信號(hào)的時(shí)域圖');subplot(3,1,2);plot(t,sfm);繪制已調(diào)信號(hào)的時(shí)域圖xlabel('

28、;時(shí)間 t');title('無噪聲條件下已調(diào)信號(hào)的時(shí)域圖');接受信號(hào)通過微分器處理nsfm=sfm;for i=1:length(t)-1diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;diff_nsfmn = abs(hilbert(diff_nsfm); %hilbert幅度(包絡(luò)檢波)變換，求絕對(duì)值得到瞬時(shí)zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;subplot(3,1,3);%繪制無噪聲條件下解調(diào)信號(hào)的時(shí)域圖plot(1:length(dif

29、f_nsfmn1)./1000,diff_nsfmn1./400,r,);xlabel('時(shí)間 t');title（,無噪聲條件下解調(diào)信號(hào)的時(shí)域圖,）;0/嚴(yán) *disp（'按任意鍵可以看到原調(diào)制信號(hào)、小信噪比高斯白噪聲條件下已調(diào)信號(hào)和解調(diào)信號(hào)已調(diào)信號(hào)的曲線,）*figure(4)*figure(4)subplot(3,1,1);plot(t,mt);xlabel('時(shí)間 t');title(,調(diào)制信號(hào)的時(shí)域圖,);db1=amA2/(2*(10A(sn1/10);%聲的方差n1=sqrt(db1)*randn(size(t);%nsfm1=n1+sfm

30、;%號(hào)通%過信道傳輸)for i=1:length(t)-1%繪制調(diào)制信號(hào)的時(shí)域圖計(jì)算對(duì)應(yīng)的小信噪比高斯白躁生成高斯白躁聲生成含高斯白躁聲的已調(diào)信號(hào)（信接受信號(hào)通過微分器處理pausediff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i)./dt;變換，求絕對(duì)值得到瞬時(shí)繪制含小信噪比高斯白噪聲已調(diào)信diff_nsfmn1 = abs(hilbert(diff_nsfm1); %hilbert幅度(包絡(luò)檢波)zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot

31、(1:length(diff_nsfm),diff_nsfm); %號(hào)的時(shí)域圖xlabel('時(shí)間 t');title（'含小信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖'）;subplot（3,1,3）;%繪制含小信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,'r');xlabel('時(shí)間 t');title('含小信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖');0/嚴(yán) *disp（'按任意鍵可以看到原調(diào)制信號(hào)、大信噪比高斯白噪聲條件下已調(diào)信

32、號(hào)和解調(diào)信號(hào)已調(diào)信號(hào)的曲線'）pause*figure(5)*figure(5)subplot(3,1,1);plot(t,mt);xlabel('時(shí)間 t');title('調(diào)制信號(hào)的時(shí)域圖');db1=amA2/(2*(10A(sn2/10);%聲的方差n1=sqrt(db1)*randn(size(t);%nsfm1=n1+sfm;%繪制調(diào)制信號(hào)的時(shí)域圖計(jì)算對(duì)應(yīng)的大信噪比高斯白躁生成高斯白躁聲生成含高斯白躁聲的已調(diào)信號(hào)（信號(hào)通過信道傳輸）接受信號(hào)通過微分器處理for i=1:length(t)-1diff_nsfm1(i)=(nsfm1(i+1)-

33、nsfm1(i)./dt;enddiff_nsfmn1 = abs(hilbert(diff_nsfm1);%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度（包zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);繪制含大信噪比高斯白plot(1:length(diff_nsfm1),diff_nsfm1);噪聲已調(diào)信號(hào)%勺時(shí)域圖xlabel(' 時(shí)間 t');title（'含大信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖'）;subplot(3,1,3);繪制含大信

34、噪比高斯白噪聲解調(diào)信號(hào)%勺時(shí)域圖plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,'r');xlabel('時(shí)間 t');title（'含大信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖'）;PM function v,phi=env_phas(x,ts,f0)if nargout=2%nargout為輸出變數(shù)的個(gè)數(shù)z=loweq(x,ts,f0);產(chǎn)生調(diào)制信號(hào)的正交分量phi=angle(z);%angle是對(duì)一個(gè)復(fù)數(shù)求相角的函數(shù)function M,m,df=fftseq(m,ts,df)fs=1/ts;i

35、f nargin=2 n1=0;%nargin為輸入?yún)⒘康膫€(gè)數(shù)else n1=fs/df;endn2=length(m);n=2A(max(nextpow2(n1),nextpow2(n2);%nextpow2(n)取n最接近的較大2次幕M=fft(m,n);%M為信號(hào)m勺傅里葉變換，n為快速傅里葉變換的點(diǎn)數(shù)，及基n-FFT變換m=m,zeros(1,n-n2);構(gòu)建新的mW號(hào)df=fs/n;重新定義頻率分辨率function x1=loweq(x,ts,f0)t=0:ts:ts*(length(x)-1);z=hilbert(x);希爾伯特變換對(duì)的利用-通過實(shí)部來求虛部產(chǎn)生信號(hào)z的正交分量,%t程序t0=;%ts=;%fs=1/ts;%fc=300;%t=-t0/2:ts:t0/2;kf=100;%df=;%x1=z.*exp(-j*2*pi*f0*t);%并將z信號(hào)與它的正交分量加在一起信號(hào)的持續(xù)時(shí)間，用來定義時(shí)間向量抽樣間隔抽樣頻率載波頻率,fc可以任意改變時(shí)間向