FM調(diào)制解調(diào)系統(tǒng)設(shè)計(jì)與仿真

1、電子信息系統(tǒng)仿真課程設(shè)計(jì)09級(jí)電子信息工程專業(yè)083班級(jí)題目FM調(diào)制解調(diào)系統(tǒng)設(shè)計(jì)與仿真內(nèi)容摘要頻率調(diào)制（FM在常應(yīng)用通信系統(tǒng)中。FM廣泛應(yīng)用于電視信號(hào)的傳輸、衛(wèi)星和電話系統(tǒng)等。FM調(diào)制解調(diào)系統(tǒng)設(shè)計(jì)主要是通過對(duì)模擬通信系統(tǒng)主要原理和技術(shù)進(jìn)行研究，理解FM調(diào)制原理和FM系統(tǒng)調(diào)制解調(diào)的基本過程，學(xué)會(huì)建立FM調(diào)制模型并利用集成環(huán)境下的M文件，對(duì)FM調(diào)制解調(diào)系統(tǒng)進(jìn)行設(shè)計(jì)和仿真，并分別繪制出基帶信號(hào)，載波信號(hào)，已調(diào)信號(hào)的時(shí)域波形；再進(jìn)一步分別繪制出對(duì)已調(diào)信號(hào)疊加噪聲后信號(hào)，相干解調(diào)后信號(hào)和解調(diào)基帶信號(hào)的時(shí)域波形；最后繪出FM基帶信號(hào)通過上述信道和調(diào)制和解調(diào)系統(tǒng)后的誤碼率與信噪比的關(guān)系，并通過與理論結(jié)果波形

2、對(duì)比來分析該仿真調(diào)制與解調(diào)系統(tǒng)的正確性及噪聲對(duì)信號(hào)解調(diào)的影響。在課程設(shè)計(jì)中，系統(tǒng)開發(fā)平臺(tái)為WindowsXP,使用工具軟件為7.0。在該平臺(tái)運(yùn)行程序完成了對(duì)FM調(diào)制和解調(diào)以及對(duì)疊加噪聲后解調(diào)結(jié)果的觀察。通過該課程設(shè)計(jì)，達(dá)到了實(shí)現(xiàn)FM信號(hào)通過噪聲信道，調(diào)制和解調(diào)系統(tǒng)的仿真目的。從而了解FM調(diào)制解調(diào)系統(tǒng)的優(yōu)點(diǎn)和缺點(diǎn)，有利于以后設(shè)計(jì)應(yīng)用。關(guān)鍵詞FM;調(diào)制；解調(diào)；MATLAB仿真；信噪比、MATLAB軟件簡(jiǎn)介MATLAB是由美國(guó)mathworks公司發(fā)布的主要面對(duì)科學(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在

3、一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語(yǔ)言（如GFortran）的編輯模式，代表了當(dāng)今國(guó)際科學(xué)計(jì)算軟件的先進(jìn)水平。MATLAB和mathematicsmaple并稱為三大數(shù)學(xué)軟件。它以矩陣為基本數(shù)據(jù)單位，在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算方面首屈一指。MATLAB可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語(yǔ)言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號(hào)處理與通訊、圖像處理、信號(hào)檢測(cè)、金融建模設(shè)計(jì)與分析等領(lǐng)域。二、理論分析2.1一般通信系統(tǒng)通信的目的是傳輸信息

4、。一般通信系統(tǒng)的作用就是將信息從信息源發(fā)送到一個(gè)或多個(gè)目的地。對(duì)于任何一個(gè)通信系統(tǒng)，均可視為由發(fā)送端、信道和接收端三大部分組成(如圖1所示)圖1通信系統(tǒng)一般模型2.2FM調(diào)制原理調(diào)制在通信系統(tǒng)中具有十分重要的作用。一方面，通過調(diào)制可以把基帶信號(hào)的頻譜搬移到所希望的位置上去，從而將調(diào)制信號(hào)轉(zhuǎn)換成適合于信道傳輸或便于信道多路復(fù)用的已調(diào)信號(hào)。另一方面，通過調(diào)制可以提高信號(hào)通過信道傳輸時(shí)的抗干擾能力，同時(shí)，它還和傳輸效率有關(guān)。具體地講，不同的調(diào)制方式產(chǎn)生的已調(diào)信號(hào)的帶寬不同，因此調(diào)制影響傳輸帶寬的利用率?？梢姡{(diào)制方式往往決定一個(gè)通信系統(tǒng)的性能。在本仿真的過程中我們選擇用調(diào)頻調(diào)制方法進(jìn)行調(diào)制。調(diào)制過程

5、是一個(gè)頻譜搬移的過程，它是將低頻信號(hào)的頻譜搬移到載頻位置。而解調(diào)是將位于載頻的信號(hào)頻譜再搬回來，并且不失真地恢復(fù)出原始基帶信號(hào)。在本仿真的過程中我們選擇用非相干解調(diào)方法進(jìn)行解調(diào)。而頻率(FM調(diào)制的名稱源于m(t)與已調(diào)信號(hào)的頻率呈線性關(guān)系。FM調(diào)制就是將調(diào)制信號(hào)的變化映射到已調(diào)信號(hào)的頻率大小。設(shè)調(diào)制信號(hào)為m(t),調(diào)頻信號(hào)的數(shù)學(xué)表達(dá)式為Sfm(t)=AcosRfct2kFMm(t)dt例如：m(t)的時(shí)域波形為m(t)100.51t-1LFM調(diào)頻波如下：FM信號(hào)00512.3FM調(diào)制模型的建立現(xiàn)(f)FM圖2FM調(diào)制模型其中，m(t)為基帶調(diào)制信號(hào)，設(shè)調(diào)制信號(hào)為m(t)=Acos(2fmt)設(shè)

6、正弦載波為c(t)=cos(2fct)可得到已調(diào)調(diào)頻信號(hào)為Sfm(t廣Accos2二fct2kFMm(t)dt2假設(shè)信號(hào)傳輸信道為局斯白噪聲信道，其功率為。2.2調(diào)制過程分析在調(diào)制時(shí)，調(diào)制信號(hào)的頻率去控制載波的頻率的變化，載波的瞬時(shí)頻偏隨調(diào)制信號(hào)m(t)成正比例變化，即d(t)=Kfm(t)dtrad式中，Kf為調(diào)頻靈敏度(7(s，V)。這時(shí)相位偏移為(t)=Kfm()d則可得到調(diào)頻信號(hào)為sfm(t)=Acosct+Kfm()d"'%設(shè)定時(shí)間步長(zhǎng)%產(chǎn)生時(shí)間向量%設(shè)定調(diào)制信號(hào)幅度?？筛恼{(diào)制信號(hào)產(chǎn)生的M文件：dt=0.001;t=0:dt:1.5;am=15;fm=15;mt=

7、am*cos(2*pi*fm*t);fc=50;ct=cos(2*pi*fc*t);kf=10;int_mt(1)=0;fori=1:length(t)-1int_mt(i+1)=int_mt(i)+mt(i)*dt;end%設(shè)JE調(diào)制信號(hào)頻率?？筛?生成調(diào)制信號(hào)%設(shè)定載波頻率可更改%生成載波%設(shè)定調(diào)頻指數(shù)%對(duì)mt進(jìn)行積分%調(diào)制，產(chǎn)生已調(diào)信號(hào)調(diào)制信號(hào)的時(shí)域圖時(shí)間t載波的時(shí)域圖時(shí)間t已調(diào)信號(hào)的時(shí)域圖10-10000.5時(shí)間t11.5sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);圖3FM調(diào)制2.3FM解調(diào)模型的建立調(diào)制信號(hào)的解調(diào)分為相干解調(diào)和非相干解調(diào)兩種。相干解調(diào)僅

8、僅適用于窄帶調(diào)頻信號(hào)，且需同步信號(hào)，故應(yīng)用范圍受限；而非相干解調(diào)不需同步信號(hào)，且對(duì)于NBFM言號(hào)和WBF暗號(hào)均適用，因此是FM系統(tǒng)的主要解調(diào)方式。在本仿真的過程中我們選擇用非相干解調(diào)方法進(jìn)行解調(diào)屆II11鑒頻器Sd(t)1111曲1(/)限幅器11if包絡(luò)嬴(f)低通及帶通I1瞰jjr-fifr檢波11波波器F"圖4FM解調(diào)模型非相干解調(diào)器由限幅器、鑒頻器和低通濾波器等組成，其方框圖如圖5所示。限幅器輸入為已調(diào)頻信號(hào)和噪聲，限幅器是為了消除接收信號(hào)在幅度上可能出現(xiàn)的畸變；帶通濾波器的作用是用來限制帶外噪聲，使調(diào)頻信號(hào)順利通過。鑒頻器中的微分器把調(diào)頻信號(hào)變成調(diào)幅調(diào)頻波，然后由包絡(luò)檢波器

9、檢出包絡(luò)，最后通過低通濾波器取出調(diào)制信號(hào)。2.4解調(diào)過程分析設(shè)輸入調(diào)頻信號(hào)為tS(t廣SFM(t廣Acos(ctKfm()d)cO微分器的作用是把調(diào)頻信號(hào)變成調(diào)幅調(diào)頻波。微分器輸出為dS(t)dtdSFM(t)dt-1cKfm(t)】sin(ctKf_m()d)包絡(luò)檢波的作用是從輸出信號(hào)的幅度變化中檢出調(diào)制信號(hào)。包絡(luò)檢波器輸出為So(t)=KdIcKm(t)1=KdcKdKm(t)心稱為鑒頻靈敏度(%z),是已調(diào)信號(hào)單位頻偏對(duì)應(yīng)的調(diào)制信號(hào)的幅度，經(jīng)低通濾波器后加隔直流電容，隔除無用的直流，得mo(t)=KdKfm(t)微分器通過程序?qū)崿F(xiàn)，代碼如下：fori=1:length(t)-1%接受信號(hào)

10、通過微分器處理diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn=abs(hilbert(diff_nsfm);%hilbert變換,求絕對(duì)值得到瞬時(shí)幅度一一(包絡(luò)檢波)通過M文件繪制出兩種不同信噪比解調(diào)的輸出波形如下：調(diào)制信號(hào)的時(shí)域圖51:-51100.511.5時(shí)間t時(shí)間t無噪聲條件下解調(diào)信號(hào)的時(shí)域圖00.51時(shí)間t1.5圖5FM解調(diào)2.5高斯白噪聲信道特性設(shè)正弦波通過加性高斯白噪聲信道后的信號(hào)為r(t)=Acos(ct口)n(t)其中，白噪聲n的取值的概率分布服從高斯分布。MATLAB身白帶了標(biāo)準(zhǔn)高斯分布的內(nèi)部函數(shù)randn。randn函數(shù)

11、產(chǎn)生的隨機(jī)序列服從均值為m=0,方差§2=1的高斯分布。正弦波通過加性高斯白噪聲信道后的信號(hào)為r(t)=Acos(ct口)n(t)故其有用信號(hào)功率為A2S2噪聲功率為N二，2信噪比%滿足公式B=10logio*)則可得到公式A2B21010我們可以通過這個(gè)公式方便的設(shè)置高斯白噪聲的方差。在本仿真過程中，我們選擇了10db和30db兩種不同信噪比以示區(qū)別，其時(shí)域圖如圖7和圖8無噪聲條件下解調(diào)信號(hào)的時(shí)域圖511k,/XXXXX/fif/ffifIfX/0-/l-51100.511.5時(shí)間t圖6無噪聲條件下已調(diào)信號(hào)的時(shí)域圖調(diào)制信號(hào)的時(shí)域圖00.511.5時(shí)間t含小信噪比高斯白噪聲已調(diào)信號(hào)的

12、時(shí)域圖時(shí)間t含小信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖-200.51.5圖7含小信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖調(diào)制信號(hào)的時(shí)域圖時(shí)間t含大信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖時(shí)間t含大信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖圖8含大信噪比高斯白噪聲已調(diào)信號(hào)的時(shí)域圖2.6調(diào)頻系統(tǒng)的抗噪聲性能分析從前面的分析可知，調(diào)頻信號(hào)的解調(diào)有相干解調(diào)和非相干解調(diào)兩種。相干解調(diào)僅適用于窄帶調(diào)頻信號(hào)，且需同步信號(hào)；而非相干解調(diào)適用于窄帶和寬帶調(diào)頻信號(hào)，而且不需同步信號(hào)，因而是FM系統(tǒng)的主要解調(diào)方式，所以這里僅僅討論非相干解調(diào)系統(tǒng)的抗噪聲性能，其分析模型如圖9所示。1解調(diào)器圖9調(diào)頻系統(tǒng)抗噪聲性能分析模型圖中帶通濾波器的作用是抑制信號(hào)

13、帶寬以外的噪聲。n是均值為零，單邊功率譜密度為n。的高斯白噪聲，經(jīng)過帶通濾波器后變?yōu)檎瓗Ц咚乖肼昻i(t)。限幅器是為了消除接收信號(hào)在幅度上可能出現(xiàn)的畸變?cè)O(shè)調(diào)頻信號(hào)為t&M(t)=Acos(ctKfm()d)-joO故其輸入功率為Ni=noBFM輸入噪聲功率為因此輸入信噪比為SA2Ni2Bfm在大信噪比條件下,信號(hào)和噪聲的相互作用可以忽略，這時(shí)可以把信號(hào)和噪聲分開來算，這里，我們可以得到解調(diào)器的輸出信噪比So3A2K；m2(t)No82nofm3上式中，A為載波的振幅，Kf為調(diào)頻器靈敏度，fm為調(diào)制信號(hào)m的最高頻率，no為噪聲單邊功率譜密度。我們?nèi)缛艨紤]m（t）為單一頻率余弦波時(shí)的情況

14、，可得到解調(diào)器的制度增益為Gfm80No3奇七Gfm-mf8Ni2nofm考慮在寬帶調(diào)頻時(shí)，信號(hào)帶寬為Bm=2（%m+1）B則可以得到Gm=3言帝（t）/m（t）max可以看出，大信噪比時(shí)寬帶調(diào)頻系統(tǒng)的信噪比增益是很高的，它一P一，、與調(diào)頻指數(shù)的立方成正比。可見，加大調(diào)頻指數(shù)fm,可使調(diào)頻系統(tǒng)的抗噪聲性能迅速改善。但加大調(diào)頻指數(shù)pfm,又會(huì)增大信號(hào)占用帶寬Bfm，所以說fm調(diào)制的大信噪比是用帶寬來?yè)Q取的。各種調(diào)制解調(diào)方式的傳輸增益和所占帶寬各有不同，我們可根據(jù)實(shí)際情況選用比較合適的調(diào)制解調(diào)方式，下表列出了幾種方式的各項(xiàng)指標(biāo)（表2.1）:調(diào)制方式傳輸帶寬GsYS=（N）O/（S）/baseban

15、d實(shí)現(xiàn)的難易線性調(diào)制AM2B2/、m"m最簡(jiǎn)單DSB-SC2B1復(fù)雜SSBB1最復(fù)雜VSB&B1中等非線性調(diào)制FM2（Efm+1）B3瞌m"t）/m（t）|max簡(jiǎn)單PM2（%m+1）BMm"血匚復(fù)雜表2.1模擬調(diào)制系統(tǒng)的帶寬與噪聲性能三. 仿真實(shí)現(xiàn)圖10程序流程圖3.1MATLAB源代碼%FM調(diào)制解調(diào)系統(tǒng).m%頻率調(diào)制與解調(diào)的Matlab演示源程序%可以任意改原調(diào)制信號(hào)函數(shù)m(t)%電子信息工程杜懷超091308305%*%*初始化*echooffcloseallclearallclc%*FM%*調(diào)制*dt=0.001;t=0:dt:1.5;%設(shè)定時(shí)間步

16、長(zhǎng)%產(chǎn)生時(shí)間向量%設(shè)定調(diào)制信號(hào)幅度%設(shè)定調(diào)制信號(hào)頻率%生成調(diào)制信號(hào)%設(shè)定載波頻率%生成載波%設(shè)定調(diào)頻指數(shù)am=5;fm=5;mt=am*cos(2*pi*fm*t);fc=50;ct=cos(2*pi*fc*t);kf=10;int_mt(1)=0;fori=1:length(t)-1%求信號(hào)m（t）的積分%調(diào)制，產(chǎn)生已調(diào)信號(hào)%調(diào)制信號(hào)int_mt(i+1)=int_mt(i)+mt(i)*dt;endsfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);%*%*%*添加高斯白噪聲*%設(shè)定信躁比（小信噪比）%設(shè)定信躁比（大信噪比）%設(shè)定信躁比（無信噪比）%計(jì)算對(duì)應(yīng)的高斯白躁

17、聲的方sn1=10;sn2=30;sn=0;n=sqrt(db)*randn(size(t);nsfm=n+sfm;db=amA2/(2*(10A(sn/10);%生成局斯白躁聲%生成含高斯白躁聲的已調(diào)信號(hào)（信號(hào)通%過信道傳輸）%*%.*.%*fm解調(diào)*fori=1:length(t)-1diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn=abs(hilbert(diff_nsfm);(布絡(luò)檢波)zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;%*%.*%*時(shí)

18、域到頻域轉(zhuǎn)換*ts=0.001;fs=1/ts;df=0.25;傅里葉變換%接受信號(hào)通過微分器處理%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度%抽樣間隔%抽樣頻率%所需的頻率分辨率，用在求%時(shí)，它表示FFT的最小頻率間隔%*對(duì)調(diào)制信號(hào)m（t）求傅里葉變換*%原調(diào)信號(hào)m=am*cos（2*pi*fm*t）;fs=1/ts;ifnargin=2n1=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;u求傅里變換M=M/fs;察f=0:df1:d

19、f1*(length(m)-1)-fs/2;%以上程序是對(duì)調(diào)制后的信號(hào)%縮放，便于在頻鋪圖上整體觀%時(shí)間向量對(duì)應(yīng)的頻率向量%*對(duì)已調(diào)信號(hào)u求傅里變換*fs=1/ts;ifnargin=2n1=0;elsen1=fs/df;endn2=length(sfm);n=2A(max(nextpow2(n1),nextpow2(n2);U=fft(sfm,n);u=sfm,zeros(1,n-n2);df1=fs/n;變換U=U/fs;%*%*%*%*顯示程序*disp('按任意鍵可以看到原調(diào)制信號(hào)、載波信號(hào)和已調(diào)信號(hào)的曲線pause%*回佗(1)*figure(1)subplot(3,1,1)

20、;plot(t,mt);xlabel('時(shí)間t');title('調(diào)制信號(hào)的時(shí)域圖');subplot(3,1,2);plot(t,ct);xlabel('時(shí)間t');title('載波的時(shí)域圖');%以上是對(duì)已調(diào)信號(hào)u求傅里%縮放')%繪制調(diào)制信號(hào)的時(shí)域圖%繪制載波的時(shí)域圖subplot（3,1,3）;plot（t,sfm）;%繪制已調(diào)信號(hào)的時(shí)域圖xlabel（'時(shí)間t'）;title（'已調(diào)信號(hào)的時(shí)域圖'）;%*disp（'按任意鍵可以看到原調(diào)制信號(hào)和已調(diào)信號(hào)在頻域內(nèi)的圖形

21、9;）pause*figure(2)*figure(2)subplot(2,1,1)%fftshift:將FFT中的DC分量移plot(f,abs(fftshift(M)到頻譜中心xlabel('頻率f)title('原調(diào)制信號(hào)的頻譜圖')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)*figure(3)subplot(3,1,1);

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

23、度zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;%繪制無噪聲條件下解調(diào)信號(hào)的diff_nsfmn1=diff_nsfmn-zero;subplot(3,1,3);時(shí)域圖plot（1:length（diff_nsfmn1）./1000,diff_nsfmn1./400,'r'）;xlabel（'時(shí)間t'）;title（'無噪聲條件下解調(diào)信號(hào)的時(shí)域圖'）;%*disp（'按任意鍵可以看到原調(diào)制信號(hào)、小信噪比高斯白噪聲條件下已調(diào)信號(hào)和解調(diào)信號(hào)已調(diào)信號(hào)的曲線')pause*figure(4)*figure

24、(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;號(hào)(信號(hào)通fori=1:length(t)-1diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i)./dt;end%繪制調(diào)制信號(hào)的時(shí)域圖%計(jì)算對(duì)應(yīng)的小信噪比高斯白%生成高斯白躁聲%生成含高斯白躁聲的已調(diào)信%過信道傳輸)%接受信號(hào)通過微分器處理diff_nsfmn1=abs(hilbert(

25、diff_nsfm1);%hilbert變換，求絕對(duì)值得到瞬時(shí)幅度(初絡(luò)檢波)-zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot(1:length(diff_nsfm),diff_nsfm);%繪制含小信噪比高斯白噪聲已調(diào)信號(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

26、)./1000,diff_nsfmn1./400,'r');xlabel('時(shí)間t');title('含小信噪比高斯白噪聲解調(diào)信號(hào)的時(shí)域圖');%*disp('按任意鍵可以看到原調(diào)制信號(hào)、大信噪比高斯白噪聲條件下已調(diào)信號(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)*ra

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