FM調(diào)制和相干解調(diào),非相干解調(diào)

粹1通信原理課程設(shè)計(jì)報(bào)告題目：模擬通信系統(tǒng)頻率調(diào)制的買現(xiàn)專業(yè)班級(jí)：13電信科技2班姓名：學(xué)號(hào)：指導(dǎo)教師： 彭文娟【摘要】本設(shè)計(jì)主要是用MATLAB軟件仿真實(shí)現(xiàn)語音信號(hào)作為輸入信號(hào)的FM調(diào)制與解調(diào)，并繪制出原調(diào)制信號(hào)、已調(diào)信號(hào)、解調(diào)信號(hào)的時(shí)域和頻域圖；再進(jìn)一步繪制出已調(diào)信號(hào)通過高斯白噪聲信道的時(shí)域圖，最后繪出FM基帶信號(hào)通過上述信道后的解調(diào)波形和時(shí)域圖。最后通過理論與波形對(duì)比分析不同信噪比對(duì)信號(hào)解調(diào)的分析。關(guān)鍵詞：FM；調(diào)制；解調(diào)；MATLAB；高斯白噪聲目錄TOC\o"1-5"\h\z1弓丨言 11.1課程設(shè)計(jì)題目 11.2課程設(shè)計(jì)要求 1\o"CurrentDocument"2FM調(diào)制解調(diào)系統(tǒng)建立 12.1調(diào)制模型 1\o"CurrentDocument"FM解調(diào) 3相干解調(diào) 3非相干解調(diào) 5\o"CurrentDocument"3系統(tǒng)抗噪聲性能分析 7\o"CurrentDocument"3.1高斯白噪聲特性 73.2不同信噪比對(duì)信號(hào)傳輸?shù)挠绊?8\o"CurrentDocument"4部分程序代碼 9\o"CurrentDocument"5心得體會(huì) 12\o"CurrentDocument"6參考文獻(xiàn) 14第第#頁信道為高斯白噪聲信道，其功率為a2$cope4$cope4調(diào)制信號(hào)的時(shí)域圖Q密?1氓*飛 a5cope3-口觀tfl,. , , -| i i rU"?????■?,?-???■?????????????????T??????'?--*'***???-?T f?1■ 1???^Bi1|?“.??????? ??I，■????■?**I???■?■?** ? ?????? ??■??■■■?■??■ ?，??。?**? ??■????* ??????■■???■■■??????，?■??Qo????■????????■??????！?。閒刪刪『刪刪刪輸刪I刪刪訛呱刪?刪刪1曲1刪酬刪l刪??■?■：??■??????■■????????■??“?■??_??? ???■????????■???《??? ￡■???■????■??????????■?■??????■???????????????????-??????■TOC\o"1-5"\h\zI■ W*???| f??■■■■ ?■ 1?? I?|*，??丨$■?■ (??■ 丨■?: ?■? "?<tt?????|?^tI■???{ ??? 丨.■>| ???<| ?1t99?I| ??I*1??■■ ?■?? ?? ?■-€■■4 P * k > r- :?0-.-? .?已調(diào)的時(shí)域圖FM解調(diào)調(diào)制信號(hào)的解調(diào)有相干解調(diào)和非相干解調(diào)兩種。相干解調(diào)即同步解調(diào)，僅僅適用于窄帶調(diào)頻信號(hào)，且需要同步信號(hào)。非相干解調(diào)不需要同步信號(hào)，應(yīng)用較為廣泛。相干解調(diào)相干解調(diào)模型

帶通濾波器輸入為通過噪聲信道的已調(diào)頻信號(hào)，帶通濾波器是為了濾除頻帶外信道加性噪聲。乘法器的作用是將信號(hào)與相干載波相乘。低通濾波器可以濾除高頻分量。微分器是為了去掉m(t)外的高頻分量。本設(shè)計(jì)的相干解調(diào)是基于matlab中simulink仿真電路實(shí)現(xiàn)。Simulink仿真電路□Scope2Simulink仿真電路□Scope2調(diào)制、已調(diào)、解調(diào)信號(hào)的時(shí)域圖非相干解調(diào)鑒頻器非相干解調(diào)器由限幅器、鑒頻器和低通濾波器等組成，其流程如圖所示。限幅器輸入為已調(diào)頻信號(hào)和噪聲，限幅器是為了消除接收信號(hào)在幅度上可能出現(xiàn)的畸變；帶通濾波器的作用是用來限制帶外噪聲，使調(diào)頻信號(hào)順利通過。鑒頻器中的微分器把調(diào)頻信號(hào)變成調(diào)幅調(diào)頻波，然后由包絡(luò)檢波器檢出包絡(luò)，最后通過低通濾波器取出調(diào)制信號(hào)。(引用百度文庫)輸入調(diào)頻信號(hào)tt=FM= 0 [fTTS(t)S(t)Acos(tICm()d微分器輸岀信號(hào)為c/、dSi(t)dSFMt()Sd(t)=一二 L7dtdtt=-bc+Kfm(t)!in(閒ct+Kf \包絡(luò)檢波的作用是從輸出信號(hào)的幅度變化中檢出調(diào)制信號(hào)。 包絡(luò)檢波器輸出為So(t)KdCKfm()t=Kdfe+KdKi(汕尹0+Kd稱為鑒頻靈敏度，是已調(diào)信號(hào)單位頻偏對(duì)應(yīng)的調(diào)制信號(hào)的幅度，隔除無用的直流，得mo(t)=KdKim()t

本過程基于mdtbb代碼實(shí)現(xiàn):程序流程圖微分器通過程序?qū)崿F(xiàn)，代碼如下:fori=l:length(t)-l %接受信號(hào)通過微分器處理diffLnsfm(i)=(nsfm(i4-l)-nsfm(i))./dt;enddifILnsfrnn=abs(hilbert(diffLnsfm));%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度%(包絡(luò)檢波)(引用)通過代碼實(shí)現(xiàn)解調(diào)波形如下：|戶?回fat鍛制信號(hào)的時(shí)3或團(tuán)061■I-icn0&時(shí)間t|戶?回fat鍛制信號(hào)的時(shí)3或團(tuán)061■I-icn0&時(shí)間t無嗓蘆茶件下己調(diào)常號(hào)的時(shí)域圖時(shí)間I無噪聲采件下解関信號(hào)的時(shí)域圖5?時(shí)間IBlF蚣ute3I>■■ —i文件迫翔馳JBEfitV)駆入壯iTfl(D￡DO?3CW)二i凸uj飛.亍、心⑩啜比-a口0：Tb3系統(tǒng)的抗噪聲性能分析3.1高斯白噪聲特性MATLAB自身帶有標(biāo)準(zhǔn)高斯分布的內(nèi)部函數(shù) 。函數(shù)產(chǎn)生的隨機(jī)randnrandn序列服從均值為m方，方差◎J］的高斯分布設(shè)信號(hào)通過高斯白噪聲信道后的信號(hào)變?yōu)?/p>

r(t)=AcocJSc(6 +)nt(A2其有用信號(hào)功率為 S=——A2T-2*1010A2T-2*1010噪聲功率為信噪比為可得3.2不同信噪比對(duì)信號(hào)對(duì)信號(hào)傳輸影響本報(bào)告重點(diǎn)討論非相干解調(diào)時(shí)，大小信噪比對(duì)信號(hào)傳輸?shù)挠绊?。S調(diào)頻信號(hào)(t)FM二=Acosqt+:cIZtf[m(叫T輸入功率為Si—A2，噪聲功率為NinoBFM2SiA2加以—=-Ni2Bfm在大信噪比條件下，信號(hào)和噪聲的相互作用可以忽略，這時(shí)可以把信號(hào)和噪。3A2K2H12(t)-5—= ￡ 聲分開來算，這里，我們可以得到解調(diào)器的輸岀信噪比N。8^2nofm3上面的式子中，A上面的式子中，A是載波的振幅,Kf是調(diào)頻器靈敏度，fm是調(diào)制信號(hào)m(t)的GFMFlofin最高頻率，n。是噪聲單邊功率譜密度。GFMFlofin得到解調(diào)器的制度增益為考慮在寬帶調(diào)頻時(shí)，信號(hào)帶寬為Bfm=2(mf+1)f釧公(f+fin)可以得到2GFM"3mf(mr1)(引用百度文庫)可以看出，大信噪比時(shí)寬帶調(diào)頻系統(tǒng)的信噪比增益是很高的， 它與調(diào)頻指數(shù)的立方成正比。可見，加大調(diào)頻指數(shù) ，可使調(diào)頻系統(tǒng)的抗噪聲性能迅速改善。在仿真過程中采用 10db在仿真過程中采用 10db和30db兩種不同信噪比以示區(qū)別，其時(shí)域圖如下MFigure4 … I.口回！叢?HFigure5 —— 一- [oi回tBSWl4部分程序代碼解調(diào)fori=l:length(t)-l4部分程序代碼解調(diào)fori=l:length(t)-ldiftLnsfm(i)=(nslTn(i4-l)-nsfm(i))./dt;%接受信號(hào)通過微分器處理enddiftlnsfmn=abs(hilbert(diflLnsfiTi));%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(包絡(luò)檢波)zero=(max(diftLnslmn)-min(diftLnsfmn))/2;diftLnsfmnl=difLnslTnn-zero;g、[?1*■匕?匕?1*■”enddiftlnsfmn=abs(hilbert(diflLnsfiTi));%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(包絡(luò)檢波)zero=(max(diftLnslmn)-min(diftLnsfmn))/2;diftLnsfmnl=difLnslTnn-zero;g、[?1*■匕?匕?1*■”?上■匕^1^?丄小^2^^1^?*?*?*?*?*?*?*?*?*?*?*?*?*?*?%****材**材**材ts=0.001;fs=l/ts;d?=0.25;I甘域轉(zhuǎn)換為頻域**材**材**材**%*****對(duì)調(diào)制信號(hào)m(t)求傅里葉變換*****m=am*cos(2*pi*fim*t);fc=l/ts;ifnargin==2nl=0;elsenl=fc/df;endn2=length(m);n=2A(max(nextpow2(nl),nextpow2(n2)));M=fft(m,n);m=[m,zeros(l,n?n2)];dfl=fs/n;M=M/fs;f=[0:dfl:dfl*(length(m)-1)]-fs/2;%抽樣間隔%抽樣頻率%所需的頻率分辨率，用在求傅里葉變換%時(shí)，它表示FFT的最小頻率間隔%調(diào)制信號(hào)%以上程序是對(duì)調(diào)制信號(hào)m(t)求傅里變換%縮放，便于在頻譜圖上觀察整體%時(shí)間向量對(duì)應(yīng)的頻率向量%以上是對(duì)已調(diào)信號(hào)u求傅里變換%縮放%接受信號(hào)通過微分器處理%************ 對(duì)已調(diào)信號(hào)u求傅里變換顯示頻域**********fs=l/ts;ifnargin==2nl=0;elsenl=fe/df;endn2=length(sfm);n=2A(max(nextpow2(nl),nextpow2(n2)));U=fft(sfm,n);u=[sfm,zeros(l,n?n2)];dfl=fs/n;U=U/fs;%■* ■匕 .匕■卜■匕■p<Tw*TwrT.<Tw <T* *T*<T*nsfm=sfiTi;fori=l:length(t)-l

diftLnsfm(i)=(nslrn(i4-1)-nsfm(i))./dt;enddiftLnsfmn=abs(hilbert(difLnsfm)); %hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(包絡(luò)檢波)zero=(max(diiT_nsltnn)-min(diftLnsfmn))/2;diiT_nsfmnl=difLnsfmn-zero;subplot(3,l,3);%繪制無噪聲條件下解調(diào)信號(hào)的時(shí)域圖plot((l:length(difLnsfmn1))71000,difLnsfmnl.MOO/r');xlabelf時(shí)間t);mieC無噪聲條件下解調(diào)信號(hào)的時(shí)域圖 ');.匕?匕dispf按任意鍵可以看到原調(diào)制信號(hào)、小信噪比高斯白噪聲條件下已調(diào)信號(hào)和解調(diào)信號(hào)已調(diào)信號(hào)的曲線*)pause**scope(4)scope(4)subplot(3,l,l);plot(t,mt);xlabel(*時(shí)間tscope(4)subplot(3,l,l);plot(t,mt);xlabel(*時(shí)間t1);title(*調(diào)制信號(hào)的時(shí)域圖*);db1=amA2/(2*(l0A(snl/10)));n1=sqrt(db1)*randn(size(t));nsfrnl=nl+sfm;%繪制調(diào)制信號(hào)的時(shí)域圖fori=l:length(t)-ldiftLnsfm1(i)=(nsfml(i+1)-nsfm1(i))./dt;end%計(jì)算對(duì)應(yīng)的小信噪比高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號(hào)(信號(hào)通%過信道傳輸)%接受信號(hào)通過微分器處理diflLnsfmnl=abs(hilbert(diff_nsfml));%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(包絡(luò)檢波)zero=(max(dilT_nslTnn)-min(dilT_nsfmn))/2;diflLnsfmn1=difLnslrnn1?zero;subplot(3,l,2);plot(l:length(diif_nsfm),difLnsfm); %繪制含小信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖xlabelf時(shí)間t);title(*含小信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖 *);subplot⑶1,3); %繪制含小信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖plot((l:length(difLnsfmn1))71000,diflLnsfmnl.MOO/r*);xlabelf時(shí)間t);titleC含小信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖 ’)；f■上?土dispf按任意鍵可以看到原調(diào)制信號(hào)、大信噪比高斯白噪聲條件下已調(diào)信號(hào)和解調(diào)信號(hào)已調(diào)信號(hào)的波形*)pause%****材**材**5<°^cope(5)材**********材**材%繪制調(diào)制信號(hào)的時(shí)域圖scope(5)%繪制調(diào)制信號(hào)的時(shí)域圖subplot(3,1,1);plot(t,mt);xlabel('時(shí)間t);

山1乂調(diào)制信號(hào)的時(shí)域圖?);db山1乂調(diào)制信號(hào)的時(shí)域圖?);db1=amA2/(2*(10A(sn2/10)));n1=sqrt(db1)*randn(size(t));nsltnl=nl+sltn;過信道傳輸)fori=l:length(t)-ldiff_nsfm1(i)=(nsfml(i+1)-nsfm1(i))./dt;enddilT_nsfmnl=abs(hilbert(dilf_nsfml));zero=(max(diiT_nslTnn)-min(diftLnsfmn))/2;dilT_nsfmn1=difLnsfmn1-zero;subplot(3,l,2);plot(l:length(difLnsfml),difLnsfml);xlabelf時(shí)間tj;titleC含大信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖subplot(3,l,3)；%計(jì)算對(duì)應(yīng)的大信噪比高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號(hào)（信號(hào)通%接受信號(hào)通過微分器處理%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度（包%絡(luò)檢波）%繪制含大信噪比高斯白噪聲已調(diào)信號(hào)%的時(shí)域圖%繪制含大信噪比高斯白噪聲解調(diào)信號(hào)%的時(shí)域圖plot((l:length(diflLnsfmn1))71000,diflLnsfmnl7400,Y);xlabelf時(shí)間t);title(?含大信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖 ’)；5心得體會(huì)本次課程設(shè)計(jì)中，我主要了解了FM調(diào)制和解調(diào)的原理，對(duì)信號(hào)的傳輸有了更深的認(rèn)識(shí)。此次課程設(shè)計(jì)也花費(fèi)了我大量的時(shí)間，主要用于學(xué)習(xí)motbb的simulink和代碼的學(xué)習(xí)。不得不說matbb功能的強(qiáng)大和復(fù)雜。設(shè)計(jì)之初我們小組查花費(fèi)了大量時(shí)間查閱相關(guān)資料弄懂了FM調(diào)制和解調(diào)的過程以及基本的matlab電路搭建，但是當(dāng)我們搭建好FM調(diào)制和解調(diào)的電路時(shí)，每次運(yùn)行都會(huì)出錯(cuò),得不到想要的結(jié)果，通過對(duì)錯(cuò)誤的反復(fù)排查及分析終于找到錯(cuò)誤的原因是simuKnk中集成模塊的內(nèi)部電路參數(shù)未設(shè)置正確。還有對(duì)matlab代碼實(shí)現(xiàn)非相干解調(diào)部分,代碼的編寫需要非常細(xì)心，一個(gè)標(biāo)點(diǎn)的錯(cuò)誤都會(huì)導(dǎo)致運(yùn)行結(jié)果的失敗，由于對(duì)matlab的不熟悉，代碼的編寫只能參照產(chǎn)考文獻(xiàn)和網(wǎng)上資料進(jìn)行借用和更該。但是經(jīng)過小組成員的不懈努力，最終實(shí)現(xiàn)了課程設(shè)計(jì)的相關(guān)要求。這次課程設(shè)計(jì)讓我深刻的認(rèn)識(shí)到的了科學(xué)的嚴(yán)謹(jǐn)性，任何錯(cuò)誤都是不允許的。還有學(xué)習(xí)交流的重要性，設(shè)計(jì)過程中很多困難都是在大家的討論和請(qǐng)教中解決的。而且在這次設(shè)計(jì)也讓我發(fā)現(xiàn)了自己很多不足之出，多所學(xué)知識(shí)的理解膚淺，學(xué)習(xí)能力較差，這也是我以后要改進(jìn)的地方，理論和實(shí)踐必須有機(jī)的結(jié)合起來才能學(xué)的更好?！緟⒖嘉墨I(xiàn)】【1】樊昌信，通信原理，電子工業(yè)出版社，2008【2】邵玉斌，Matlab/Simulink通信原理建模與仿真實(shí)例分析，清華大學(xué)出版社，2008【3】張靜，MATLAB在控制系統(tǒng)中的應(yīng)用，電子工業(yè)出版社【4】田敏、李江全、劉軍,案例解說MATLAB典型控制應(yīng)用，電子工業(yè)出版社【5】吳偉鈴，龐沁華，通信原理，北京郵電大學(xué)出版社，2005【6】劉浩、韓晶，matlab2012a―本通[7]曹志剛，錢亞生，現(xiàn)代通信原理，清華大學(xué)出版社， 1992

%FM調(diào)制解調(diào)系統(tǒng).m%頻率調(diào)制與解調(diào)的Matlab演不源程序%可以任意改原調(diào)制信號(hào)函數(shù) m⑴%信息工程 陳麗丹07323202?*?*?*?*?*?*?*?*?*?*?*?*?*?*??*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?*?初始化echooffcloseallclearallclc?*?*?*?*?*?*?*?*?*?*?*?*?*?*?(yc材******材***^**fm調(diào)制dt=0.001;t=0:dt:1.5;am=5;flTl=5；mt=am叱os(2巾i*fhi*t);250;cAcosQ^i^fc吐)；kf=10;int_mt(l)=0;fori=l:length(t)-lint_mt(i+l)=int_mt(i)+mt(i)5Mt;endsfm=am叱os(2尊^枕*t+2巾詳1<1巧n(_mt);%■* .匕 .匕■”.匕?b.匕■卜■匕%設(shè)定時(shí)間步長%產(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)制信號(hào)?*?*?*?*?*?*?*?*?*?*?*?*?*?*?%?*■丄■ ?丄snl=10;sn2=30;snl=10;sn2=30;sn=0;db=amA2/(2*(10A(sn/10)));n=sqrt(db)*randn(size(t));nsfm=n+sfm;%設(shè)定信躁比（小信噪比）%設(shè)定信躁比（大信噪比）%設(shè)定信躁比（無信噪比）%計(jì)算對(duì)應(yīng)的高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號(hào)（信號(hào)通%過信道傳輸）%料**材**材**材**卩\|解調(diào)********林*材**材**%接受信號(hào)通過微分器處理fori=l:length(t)-l%接受信號(hào)通過微分器處理ditT_nsfm(i)=(nslTn(i4-l)-nsfm(i))./dt;end

diftLnsfmn=abs(hilbert(diftLnsfm)); %hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(包絡(luò)檢波)zero=(max(difLnsftnn)-min(dilT_nsfmn))/2;difiLnsfmnl=diftlnsltnn-zero;*If *?*?*?*?*?*?*?*?*?*?*?*?*?*?2^*1* ■丄*?*?*?*?*?*?*?*?*?*?*?*?*?*??*?*?*?*?*?*?*?*?*?*?*?*?*?*?時(shí)域至0頻域轉(zhuǎn)換******材**祕(mì)**ts=0.001;fc=l/ts;dM).25;%材***對(duì)調(diào)制信號(hào)m(f)求傅里葉變換*****m=am叱os(2*pi*fm*t);fs=l/ts;ifnargin==2nl=0;elsenl=fs/df;endn2=length(m);n=2A(max(nextpow2(n1),nextpow2(n2)));M=fft(m,n);m=[m,zeros(l,n-n2)];dfl=fs/n;換M=M/fs;f=[O:dfl:dfl*(length(m)-1)]-fs/2;%抽樣間隔%抽樣頻率%所需的頻率分辨率，用在求傅里葉變換%時(shí)，它表示FFT的最小頻率間隔%原調(diào)信號(hào)%抽樣間隔%抽樣頻率%所需的頻率分辨率，用在求傅里葉變換%時(shí)，它表示FFT的最小頻率間隔%原調(diào)信號(hào)%以上程序是對(duì)調(diào)制后的信號(hào) u求傅里變%縮放，便于在頻鋪圖上整體觀察%時(shí)間向量對(duì)應(yīng)的頻率向量fs=l/ts;ifnargin==2nl=0;elsenl=fe/df;endn2=length(sfm);n=2A(max(nextpow2(nl),nextpow2(n2)));u=[sfm,zeros(l,n-n2)];dfl=fs/n;U=U/fs;%■* ■匕 .匕■卜■匕■p *TwrT. *T*<T*C*T*?7* *I**T**I**7*r*Iw*1**7***TW*iW*?**T**TW*iwr rTw<Tw<?w?Tw*Tw?Tw?Tw<^rTw

%以上是對(duì)已調(diào)信號(hào)u求傅里變換%縮放顯示程序?丄dispC按任意鍵可以看到原調(diào)制信號(hào)、載波信號(hào)和已調(diào)信號(hào)的曲線 ?)pause%材水*****5^**5^0^。^!*。(1)*材**材******材***figure(l)subplot(3,l,l);plot(t,mt);xlabel(*時(shí)間t1);titleC調(diào)制信號(hào)的時(shí)域圖?);%繪制調(diào)制信號(hào)的時(shí)域圖subplot(3,l,2);plot(t,ct);xlabelf時(shí)間f);titleC載波的時(shí)域圖?);subplot(3,l,3)；plot(t,slrn);xlabelC時(shí)間f);titleC已調(diào)信號(hào)的時(shí)域圖?);%繪制載波的時(shí)域圖%繪制已調(diào)信號(hào)的時(shí)域圖dispC按任意鍵可以看到原調(diào)制信號(hào)和已調(diào)信號(hào)在頻域內(nèi)的圖形?)pause*figure(2)figure(2)subplot(2,l,l)plot(f,abs(ffishift(M)))xlabelf頻率f)titleC原調(diào)制信號(hào)的頻譜圖?)subplot(2,l,2)plot(f,abs(fttshift(U)))xlabelC頻率f)title(f已調(diào)信號(hào)的頻譜圖?)/?匕■匕.b.匕*1**1**1* *Xj**m*^^*i**1*%ffishift:將FFT中的DC分量移到頻譜中心dispf按任意鍵可以看到原調(diào)制信號(hào)、無噪聲條件下已調(diào)信號(hào)和解調(diào)信號(hào)的曲線 ')pause%****材**材**^?1^0口2(3)*材**材***********figure(3)subplot(3,l,l);plot(t,mt);xlabef時(shí)間t1);5bC調(diào)制信號(hào)的時(shí)域圖?);%繪制調(diào)制信號(hào)的時(shí)域圖subplot(3,l,2);plot(t,sfm);xlabelf時(shí)間tj;titleC無噪聲條件下已調(diào)信號(hào)的時(shí)域圖 ?);%繪制已調(diào)信號(hào)的時(shí)域圖nsfin=sfiTi;fori=l:length(t)-ldiiT_nsfrn(i)=(nsfim(i+1)-nsfm(i))./dt;%接受信號(hào)通過微分器處理enddiiT_nstrnn=abs(hilbert(diftlnsfiTi));%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(包絡(luò)檢波)zero=(max(diff_nsiTnn)-min(diiT_nsfmn))/2;dil1Lnsfmnl=difLnslrnn-zero;subplot⑶1,3); %繪制無噪聲條件下解調(diào)信號(hào)的時(shí)域圖plot((l:length(difLnsfmn1))71000,difLnsfmnl7400,'r');xlabelf時(shí)間t*);titleC無噪聲條件下解調(diào)信號(hào)的時(shí)域圖 ?);dispf按任意鍵可以看到原調(diào)制信號(hào)、小信噪比高斯白噪聲條件下已調(diào)信號(hào)和解調(diào)信號(hào)已調(diào)信號(hào)的曲線*）pause**figure(4)figure(4)subplot(3,l,l);plot(t,mt);xlabel(*時(shí)間t1);titleCfigure(4)subplot(3,l,l);plot(t,mt);xlabel(*時(shí)間t1);titleC調(diào)制信號(hào)的時(shí)域圖*);db1=amA2/(2*(l0A(snl/10)));%繪制調(diào)制信號(hào)的時(shí)域圖n1=sqrt(db1)*randn(size(t));nsfml=nl+slTn;fori=l:length(t)-ldiftLnsfm1(i)=(nsfml(i+1)-nsfm1(i))./dt;end%計(jì)算對(duì)應(yīng)的小信噪比高斯白躁聲的方差%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信號(hào)（信號(hào)通%過信道傳輸）%接受信號(hào)通過微分器處理diflLnsfmnl=abs(hilbert(diftlnsfml));%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(包絡(luò)檢波)zero=(max(dilT_nslTnn)-min(dilT_ns