基于MATLAB的AM信號的調(diào)制與解調(diào)最終完美版

1、課程設(shè)計任務(wù)書學生姓名：xxx 專業(yè)班級：xxxxxx指導(dǎo)教師:xxx 工作單位：信息工程學院題目：am信號的仿真分析初始條件：1） matlab 軟件2） 通信原理基礎(chǔ)知識設(shè)計內(nèi)容和要求：調(diào)制信號：分別為300hz正弦信號和矩形信號；載波頻率：30khz；解調(diào)方式：分別為包絡(luò)解調(diào)和同步解調(diào)；要求：畫出以下三種情況下調(diào)制信號、已調(diào)信號、解調(diào)信號的波形、頻譜以及解調(diào)器輸入輸出信噪比的關(guān)系曲線；1）調(diào)制信號幅度=0.8×載波幅度；2）調(diào)制信號幅度=載波幅度；3）調(diào)制信號幅度=1.5×載波幅度；指導(dǎo)教師簽名： 年 月 日系主任（或責任教師）簽名： 年 月 日目錄摘要iabstra

2、ctii1 課程設(shè)計內(nèi)容與要求11.1 課程設(shè)計內(nèi)容11.2設(shè)計要求12 調(diào)制和解調(diào)原理22.1調(diào)制與解調(diào)22.2 幅度調(diào)制的原理23 同步解調(diào)與包絡(luò)解調(diào)33.1同步解調(diào)33.2 包絡(luò)解調(diào)44 噪聲與抗噪聲性能模型分析54.1 噪聲的分類54.2 本文噪聲模型54.3 抗噪聲性能的分析模型64.4相干解調(diào)的抗噪聲性能65 am基于matlab的調(diào)制與解調(diào)85.1載波頻率為30kz的載波信號85.2頻率為300hz的正弦波調(diào)制95.2.1調(diào)制信號幅度=0.8×載波幅度95.2.2調(diào)制信號幅度=載波幅度時的調(diào)制信號105.2.3調(diào)制信號幅度=1.5×載波幅度時的調(diào)制信號115.

3、3頻率為300hz的正弦波解調(diào)與輸入輸出信噪比曲線115.3.1 調(diào)制信號幅度=0.8×載波信號幅度115.3.2調(diào)制信號幅度=載波信號幅度145.3.3調(diào)制信號幅度=1.5×載波信號幅度155.4調(diào)制信號為300hz矩形脈沖波175.4.1調(diào)制信號幅度=0.8×載波信號幅度175.4.2調(diào)制信號幅度=載波信號幅度185.4.3調(diào)制信號幅度=1.5×載波信號幅度195.5 頻率為300hz的矩形波解調(diào) 及輸入輸出信噪比的關(guān)系曲線195.5.1 調(diào)制信號幅度=0.8×載波信號幅度195.5.2調(diào)制信號幅度=載波信號幅度225.5.3調(diào)制信號幅度=

4、1.5×載波信號幅度236 小結(jié)與體會257 參考文獻26摘要調(diào)制在通信系統(tǒng)中有十分重要的作用。通過調(diào)制，不僅可以進行頻譜搬移，把調(diào)制信號的頻譜搬移到所希望的位置上，從而將調(diào)制信號轉(zhuǎn)換成適合于傳播的已調(diào)信號，而且它對系統(tǒng)的傳輸有效性和傳輸?shù)目煽啃杂兄艽蟮挠绊?，調(diào)制方式往往決定了一個通信系統(tǒng)的性能。matlab軟件廣泛用于數(shù)字信號分析，系統(tǒng)識別，時序分析與建模，神經(jīng)網(wǎng)絡(luò)、動態(tài)仿真等方面有著廣泛的應(yīng)用。本課題利用matlab軟件對am調(diào)制解調(diào)系統(tǒng)進行模擬仿真，分別利用300hz正弦波和矩形波，對30khz正弦波進行調(diào)制，觀察調(diào)制信號、已調(diào)信號和解調(diào)信號的波形和頻譜分布，并在解調(diào)時引入高

5、斯白噪聲，對解調(diào)前后信號進行信噪比的對比分析，估計am調(diào)制解調(diào)系統(tǒng)的性能。調(diào)幅是中短波廣播中一種主要的調(diào)制方式,對其解調(diào)也有兩種方法。通過對同步檢波和包絡(luò)檢波兩種解調(diào)方法的分析,探討不同解調(diào)方法下系統(tǒng)的抗噪聲性能,從而從另一種角度指出了中短波廣播發(fā)射電臺必須通過提高發(fā)射功率和調(diào)幅度的方式來提升系統(tǒng)的抗噪聲性能。關(guān)鍵詞：am信號，調(diào)制，解調(diào)，信噪比，matlab abstractmodulation plays an important role in the communication system. through the modulation, no

6、t only can move frequency, the modulated signal spectrum to the desired position, which will convert the modulated signalsuitable for transmission of the modulated signal, but also has greatinfluence on the reliability of the system transmi

7、ssion efficiency andtransmission, modulation often determines the performance of acommunication system the. matlab software is widely used in digital signal analysis, system identification, analysis and modeling of time series, has a wide range

8、of application of neural network, dynamic simulation. the simulation of am modulation and demodulation systemwith matlab software, using 300hz sine wave and square wave, sine wave modulation of the 30khz, to observe the modulation

9、 signal,waveform and spectrum distribution of the modulated signal and demodulation signal, and the gauss white noise in the demodulation,compares the snr of demodulator and signal, performance estimation of am modulation and demo

10、dulation system. amplitude modulation is one of the main modulation in shortwave radio, there are two methods for the demodulation. through the analysis of synchronous detection and envelope detection of two kinds of demodulation

11、methods, to explore theanti noise performance of different demodulation system, which from another point of view in the short wave broadcast transmitting stationmust be to improve the anti noise performance of the system by i

12、ncreasing the transmission power and modulation mode.keyword: am signal, modulation, demodulation, noise ratio signal, matlab1 課程設(shè)計內(nèi)容與要求1.1 課程設(shè)計內(nèi)容調(diào)制信號：分別為300hz正弦信號和矩形信號；載波頻率：30khz；解調(diào)方式：分別為包絡(luò)解調(diào)和同步解調(diào)；1.2設(shè)計要求利線用matlab軟件進行振幅調(diào)制和解調(diào)程序設(shè)計，畫出以下三種情況下調(diào)制信號、已調(diào)信號、解調(diào)信號的波形、頻譜以及解調(diào)器輸入輸出信噪比的關(guān)系曲線；1）調(diào)制信號幅度

13、=0.8×載波幅度；2）調(diào)制信號幅度=載波幅度；3）調(diào)制信號幅度=1.5×載波幅度2 調(diào)制和解調(diào)原理2.1調(diào)制與解調(diào) 調(diào)制在通信系統(tǒng)中的作用至關(guān)重要。調(diào)制就是把信號轉(zhuǎn)換成適合在信道中傳輸?shù)男问降囊环N過程。廣義的調(diào)制分為基帶調(diào)制和帶通調(diào)制（也稱載波調(diào)制）。在無線通信中和其他大多數(shù)場合，調(diào)制一詞均指載波調(diào)制。 載波調(diào)制，就是用調(diào)制信號去控制載波的參數(shù)的過程，使載波的某一個或某幾個參數(shù)按照調(diào)制信號的規(guī)律而變化。調(diào)制信號是指來自信源的消息信號（基帶信號），這些信號可以是模擬的，也可以是數(shù)字的。未受到調(diào)制的周期性振蕩信號稱為載波，它可以是正弦波，也可以是非正弦波（如周期性脈沖序列）。

14、載波調(diào)制后稱為已調(diào)信號，它含有調(diào)制信號的全部特征。解調(diào)（也稱檢波）則是調(diào)制的逆過程，起作用是將已調(diào)信號中的調(diào)制信號恢復(fù)出來。對于振幅調(diào)制信號，解調(diào)(demodulation )就是從它的幅度變化上提取調(diào)制信號的過程?？衫贸朔e型同步檢波器實現(xiàn)振幅的解調(diào)，讓已調(diào)信號與本地恢復(fù)載波信號相乘并通過低通濾波可獲得解調(diào)信號。2.2 幅度調(diào)制的原理標準調(diào)幅就是常規(guī)雙邊帶調(diào)制，簡稱調(diào)幅（am）。假設(shè)調(diào)制信號m(t)的平均值為0，將其疊加一個直流偏量a0 后與載波相乘，即可形成調(diào)幅信號。其時域表示式為 式中：a0 為外加的直流分量；m(t)可以是確知信號，也可以是隨機信號。 當滿足條件： |m(t)|max

15、a0 （2.2.1）時，am波的包絡(luò)檢波與調(diào)制信號m(t)的形狀完全一樣，因此，用包絡(luò)檢波的方法很容易就恢復(fù)出原始調(diào)制信號；如果上述條件沒有滿足，就會出現(xiàn)“過調(diào)幅”現(xiàn)像，這時用包絡(luò)檢波將會產(chǎn)生失真。但是可以采用其他的解調(diào)方法，如同步檢波。am信號的頻譜由載頻分量、上邊帶、下邊帶三部分組成。上邊帶的頻譜結(jié)構(gòu)與原調(diào)制信號的頻譜結(jié)構(gòu)相同，下邊帶是上邊帶的鏡像。因此，am信號是帶有載波分量的雙邊帶信號，它的帶寬是基帶信號帶寬fh 的2倍，即 bam =2fh (2.2.1) 圖2.2.1 am調(diào)制模型3 同步解調(diào)與包絡(luò)解調(diào)解調(diào)是調(diào)制的逆過程，其作用是從接收的已調(diào)信號中恢復(fù)原基帶信號（即調(diào)制信號）。解調(diào)

16、的方法分為兩類：相干解調(diào)和非相干解調(diào)（包絡(luò)檢波）。3.1同步解調(diào)同步解調(diào)也叫相干解調(diào)。解調(diào)與調(diào)制的實質(zhì)一樣，均是頻譜搬移。調(diào)制是把基帶信號的譜搬到了載頻位置，這一過程可以通過一個相乘器與載波相乘來實現(xiàn)。解調(diào)則是調(diào)制的反過程，即把在載頻位置的已調(diào)信號的頻譜搬回到原始基帶位置，因此同樣可以用相乘器與載波相乘來實現(xiàn)。相干解調(diào)器的一般模型示于圖 中。lpf×sm(t)sp(t)sd(t)c(t)=coswct圖3.1.1同步解調(diào)器的一般模型相干解調(diào)時，為了無失真的恢復(fù)原基帶信號，接受端必須提供一個與接收的已調(diào)載波嚴格同步（同頻同相）的本地載波（稱為相干載波），它與接收的已調(diào)信號相乘后，經(jīng)低通

17、濾波器取出低頻分量，即可得到原始的基帶調(diào)制信號。將已調(diào)信號乘上一個與調(diào)制器同頻同相的載波，得 (3.1.1) 由上式可知，只要用一個低通濾波器，就可以將第1項與第2項分離，無失真的恢復(fù)出原始的調(diào)制信號 (3.1.2) 同步解調(diào)的關(guān)鍵是必須產(chǎn)生一個與調(diào)制器同頻同相位的載波。如果同頻同相位的條件得不到滿足，則會破壞原始信號的恢復(fù)。相干解調(diào)適用于所有線性調(diào)制信號的解調(diào)。只是信號的解調(diào)結(jié)果中含有直流成分，這時在解調(diào)后加一個簡單的隔直流電容即可。3.2 包絡(luò)解調(diào)信號在滿足|m(t)|max a0的條件下，其包絡(luò)與調(diào)制信號m(t)的形狀完全一樣。因此，am信號除了可以采用相干解調(diào)外，一般都采用簡單的包絡(luò)檢

18、波法來恢復(fù)信號。 包絡(luò)檢波器通常由半波或全波整流器和低通濾波器組成，它屬于非相干解調(diào)，因此不需要相干載波，廣播接收機中多采用此法。一個二極管峰值包絡(luò)檢波器如圖 所示，它由二極管vd和rc低通濾波器組成。 圖3.2.1包絡(luò)解調(diào)器設(shè)輸入信號是am信號 sam(t)=a0+m(t)coswct (3.2.1) 在大信號檢波時（一般大于0.5v），二極管處于受控的開關(guān)狀態(tài)。選擇rc滿足如下關(guān)系 fh<<1/rc<<fc (3.2.2) 式中：fh是調(diào)制信號的最高頻率；fc載波的頻率。在滿足上式的條件下，檢波器的輸出為 sd(t)=a0+m(t) (3.2.3)隔去直流后即可得到

19、原信號m(t)。 可見，包絡(luò)解調(diào)器就是直接從已調(diào)波的幅度中提取原調(diào)制信號。其解調(diào)輸出是相干解調(diào)輸出的2倍。因此，am信號幾乎無例外地采用包絡(luò)解調(diào)。4 噪聲與抗噪聲性能模型分析4.1 噪聲的分類噪聲的種類可廣義的分為人為噪聲、自然噪聲和內(nèi)部噪聲；也可以按噪聲對線性譜的影響是加性的還是乘性的區(qū)分，乘性噪聲又稱為相關(guān)噪聲；從信號分布來說，噪聲主要可以分為以下幾類4-5：(1) 單頻噪聲(2) 脈沖噪聲(3) 起伏噪聲(4) 平穩(wěn)非平穩(wěn)噪聲4.2 本文噪聲模型通常認為背景噪聲模型為具有雙邊譜密度的高斯隨機過程。以概率隨機取個數(shù)值中的某個值(代表時刻點)。若取，則噪聲模型可以簡化為譜密度為的平穩(wěn)加性高斯

20、白噪聲。信道的參數(shù)不可能是一直恒定的，它有可能會發(fā)生突變，這體現(xiàn)在時間下標上。但通常認為在考察時間段內(nèi)，這種突變發(fā)生的概率是很小的，即認為信道在考察時間段內(nèi)是平穩(wěn)的awgn信道3。因此本文認為信道中的噪聲都是平穩(wěn)加性高斯白噪聲。圖4.2.1 高斯白噪聲波形與功率譜 從圖4.2.1可以看出，高斯白噪聲在時域的顯著特征是波形比較雜亂，沒有任何的規(guī)律可言；在頻域的顯著特征是頻譜非常平坦，其功率在整個頻帶范圍內(nèi)均勻分布。4.3 抗噪聲性能的分析模型 各種線性已調(diào)信號在傳輸過程中不可避免地要受到噪聲的干擾，為了討論問題的簡單起見，我們這里只研究加性噪聲對信號的影響。因此，接收端收到的信號是發(fā)送信號與加性

21、噪聲之和8。 由于加性噪聲只對已調(diào)信號的接收產(chǎn)生影響，因而調(diào)制系統(tǒng)的抗噪聲性能主要用解調(diào)器的抗噪聲性能來衡量。為了對不同調(diào)制方式下各種解調(diào)器性能進行度量，通常采用信噪比增益g（又稱調(diào)制制度增益）來表示解調(diào)器的抗噪聲性能，即9 （4.3.1）有加性噪聲時解調(diào)器的數(shù)學模型如圖圖 4.3.1 am信號的解調(diào)原理圖圖中為已調(diào)信號，n(t)為加性高斯白噪聲。 和n(t)首先經(jīng)過一帶通濾波器，濾出有用信號，濾除帶外的噪聲。經(jīng)過帶通濾波器后到達解調(diào)器輸入端的信號為，噪聲為高斯窄帶噪聲，顯然解調(diào)器輸入端的噪聲帶寬與已調(diào)信號的帶寬是相同的。最后經(jīng)解調(diào)器解調(diào)輸出的有用信號為，噪聲為。4.4相干解調(diào)的抗噪聲性能各種

22、線性調(diào)制系統(tǒng)的相干解調(diào)模型如下圖所示。 圖 4.4.1 am信號的相干解調(diào)原理圖圖中可以是各種調(diào)幅信號，如am、dsb、ssb和vsb，帶通濾波器的帶寬等于已調(diào)信號帶寬10。下面討論各種線性調(diào)制系統(tǒng)的抗噪聲性能11。am信號的時域表達式為 (4.4.1) 通過分析可得am信號的平均功率為 (4.4.2)又已知輸入功率, 其中b表示已調(diào)信號的帶寬。由此可得am信號在解調(diào)器的輸入信噪比為 (4.4.3)am信號經(jīng)相干解調(diào)器的輸出信號為 (4.4.4)因此解調(diào)后輸出信號功率為 (4.4.5)在上圖中輸入噪聲通過帶通濾波器之后，變成窄帶噪聲，經(jīng)乘法器相乘后的輸出噪聲為 (4.4.6)經(jīng)lpf后， (4

23、.4.7)因此解調(diào)器的輸出噪聲功率為 (4.4.8)可得am信號經(jīng)過解調(diào)器后的輸出信噪比為 (4.4.9)由上面分析的解調(diào)器的輸入、輸出信噪比可得am信號的信噪比增益為 (4.4.10)5 am基于matlab的調(diào)制與解調(diào)5.1載波頻率為30kz的載波信號波形和頻譜的matlab程序及相應(yīng)波形圖如下：t=0:0.0000001:1; a0=10; %載波信號振幅fc=30000; %載波信號頻率uc=a0*cos(2*pi*fc*t); %載波信號figure(1);subplot(2,1,1);plot(t,uc);xlabel('t'),title('載頻信號波形&

24、#39;);axis(0,0.001,-15,15);subplot(2,1,2);y1=fft(uc); %對載波信號進行傅里葉變換plot(abs(y1);xlabel('w'),title('載波信號頻譜');axis(28500,31500,0,1000000); 圖5.1.15.2頻率為300hz的正弦波調(diào)制 設(shè)所疊加的直流分量幅值為10，以使調(diào)制信號幅度為8、10、15時，可分別觀察到欠調(diào)幅、滿調(diào)幅和過調(diào)幅現(xiàn)象。5.2.1調(diào)制信號幅度=0.8×載波幅度由于載波幅度為10，故調(diào)制信號的幅度為8，此時調(diào)制信號和已調(diào)信號的matlab程序及相應(yīng)波

25、形圖如下。t=0:0.00001:1;a0=10; %載波信號振幅a1=8; %調(diào)制信號振幅fc=30000; %載波信號頻率f=300; %調(diào)制信號頻率mes=a1*cos(2*pi*f*t);%調(diào)制信號uam=a0*(mes+10).*cos(2*pi*fc*t);%am已調(diào)信號subplot(221);plot(t,mes);xlabel('t'),title('調(diào)制信號'); axis(0,0.01,-20,+20) subplot(222);plot(t,uam);grid on;title('am已調(diào)信號波形');axis(0,0.0

26、1,-200,200) ; subplot(223);y1=fft(mes); % 對調(diào)制信號進行傅里葉變換plot(abs(y1);title('調(diào)制信號頻譜');axis(-5000,110000,0,500000); subplot(224);uam=fft(uam); % 對am已調(diào)信號進行傅里葉變換plot(abs(uam);title('am已調(diào)信號頻譜');axis(27000,77000,0,2000000) 圖5.2.15.2.2調(diào)制信號幅度=載波幅度時的調(diào)制信號由于載波幅度為10，故調(diào)制信號幅度為10, 此時有a1=10。程序略，相應(yīng)波形圖如

27、下。 圖5.2.25.2.3調(diào)制信號幅度=1.5×載波幅度時的調(diào)制信號由于載波幅度為10，故調(diào)制信號幅度為15, 此時有a1=15。程序略，相應(yīng)波形圖如下。 圖5.2.3 此時不滿足條件： |m(t)|max a0 ，由已調(diào)波形可以看出，包絡(luò)波形出現(xiàn)重疊，產(chǎn)生過調(diào)制現(xiàn)象。5.3頻率為300hz的正弦波解調(diào)與輸入輸出信噪比曲線5.3.1 調(diào)制信號幅度=0.8×載波信號幅度由于載波幅度為10，故調(diào)制信號的幅度為8，此時解調(diào)信號的matlab程序及相應(yīng)波形圖如下。fs=100000; %抽樣頻率 fc=30000; %載波頻率n=1000; %fft長度n=0:n-1;t=n/f

28、s; %截止時間和步長a1=8; %調(diào)制信號振幅a0=10; %直流分量x=a1*sin(2*pi*300*t); %基帶調(diào)制信號 y=ammod(x,fc,fs,0,10);%am已調(diào)信號yn=awgn(y,4); %加入高斯白噪聲yn1=awgn(y,10);yn2=awgn(y,15);yn3=awgn(y,20);yn4=awgn(y,25);y1=demod(y,fc,fs,'am'); %無噪聲已調(diào)信號解調(diào)yyn=demod(yn,30000,fs,'am'); %無噪聲已調(diào)信號解調(diào)yyn1=demod(yn1,30000,fs,'am

29、9;);yyn2=demod(yn2,30000,fs,'am');yyn3=demod(yn3,30000,fs,'am');yyn4=demod(yn4,30000,fs,'am');dy1=yn-y; %高斯白噪聲snr1=var(y)/var(dy1); %輸入信噪比dy2=yyn-y1; %解調(diào)后噪聲snr2=var(y1)/var(dy2); %輸出信噪比 dy11=yn1-y; snr11=var(y)/var(dy11); dy21=yyn1-y1; snr21=var(y1)/var(dy21); dy12=yn2-y; snr

30、12=var(y)/var(dy12); dy22=yyn2-y1; snr22=var(y1)/var(dy22); dy13=yn3-y; snr13=var(y)/var(dy13); dy23=yyn3-y1; snr23=var(y1)/var(dy23); dy14=yn4-y; snr14=var(y)/var(dy14); dy24=yyn4-y1; snr24=var(y1)/var(dy24); in=snr1,snr11,snr12,snr13,snr14; out=snr2,snr21,snr22,snr23,snr24; figure(1)plot(in,out,&#

31、39;*')hold onplot(in,out)xlabel('輸入信噪比')ylabel('輸入信噪比') figure(2)subplot(311)plot(t,yyn)xlabel('加噪聲解調(diào)信號波形')ff=fft(yyn,n);mag=(abs(ff);f=(0:length(ff)-1).*fs/length(ff);subplot(312)plot(f,mag)axis(0 2000 0 10000)xlabel('解調(diào)信號頻譜')subplot(313)plot(t,y1)xlabel('無噪聲

32、解調(diào)信號波形')figure(3) a=1; %調(diào)制度env=abs(hilbert(y);dem=2*(env-a0)/a;dem=fft(dem);subplot(211);plot(t,dem);title('包絡(luò)檢波已調(diào)信號');xlabel('時間')subplot(212);plot(f,abs(fftshift(dem);title('解調(diào)信號頻譜')xlabel('頻率') 圖5.3.1 圖5.3.25.3.2調(diào)制信號幅度=載波信號幅度由于載波幅度為10，故調(diào)制信號幅度為10, 此時有a1=10。程序略，相

33、應(yīng)波形圖如下。 圖5.3.3 圖5.3.45.3.3調(diào)制信號幅度=1.5×載波信號幅度由于載波幅度為10，故調(diào)制信號幅度為15, 此時有a1=15。程序略，相應(yīng)波形圖如下。圖5.3.5 圖5.3.6由圖5.3.5可以看出，由于此時產(chǎn)生過調(diào)幅現(xiàn)象，包絡(luò)檢波不能正確解調(diào)出調(diào)制信號的包絡(luò)。圖5.3.75.4調(diào)制信號為300hz矩形脈沖波5.4.1調(diào)制信號幅度=0.8×載波信號幅度由于載波幅度為10，故調(diào)制信號的幅度為8，此時調(diào)制信號和已調(diào)信號的matlab程序及相應(yīng)波形圖如下。fs=100000;n=10000;n=0:n-1;t=n/fs;a0=10;%載波信號振幅a1=8;

34、%調(diào)制信號振幅fc=30000; %載波信號頻率f=300; %調(diào)制信號幅度mes=a1*square(2*pi*f*t,50);%調(diào)制信號uam=a0*(mes+10).*cos(2*pi*fc*t);%am已調(diào)信號subplot(221);plot(t,mes);xlabel('t'),title('調(diào)制信號'); axis(0,0.01,-20,+20) subplot(222);plot(t,uam);grid on;title('am已調(diào)信號波形');axis(0,0.01,-300,300) ;subplot(223);y2=fft(

35、mes); % 對調(diào)制信號進行傅里葉變換plot(abs(y2);title('調(diào)制信號頻譜');axis( -10,500,0,60000);subplot(224);uam=fft(uam); % 對am已調(diào)信號進行傅里葉變換plot(abs(uam);title('am已調(diào)信號頻譜');axis(0,11000,0,400000) 圖5.4.15.4.2調(diào)制信號幅度=載波信號幅度由于載波幅度為10，故調(diào)制信號幅度為10, 此時有a1=10。程序略，相應(yīng)波形圖如下。 圖5.4.25.4.3調(diào)制信號幅度=1.5×載波信號幅度由于載波幅度為10，故調(diào)制

36、信號幅度為15, 此時有a1=15。程序略，相應(yīng)波形圖如下。 圖5.4.35.5 頻率為300hz的矩形波解調(diào) 及輸入輸出信噪比的關(guān)系曲線5.5.1 調(diào)制信號幅度=0.8×載波信號幅度由于載波幅度為10，故調(diào)制信號的幅度為8，此時解調(diào)信號的matlab程序及相應(yīng)波形圖如下。fs=100000; %抽樣頻率 fc=30000; %載波頻率n=1000; %fft長度f0=300;n=0:n-1;t=n/fs; %截止時間和步長a1=8; %調(diào)制信號振幅a0=10; x=a1*square(2*pi*f0*t,50); %基帶調(diào)制信號 y=ammod(x,fc,fs,0,10); %am

37、已調(diào)信號yn=awgn(y,4); %加入高斯白噪聲yn1=awgn(y,10);yn2=awgn(y,15);yn3=awgn(y,20);yn4=awgn(y,25);y1=demod(y,fc,fs,'am'); %無噪聲已調(diào)信號解調(diào)yyn=demod(yn,30000,fs,'am'); %無噪聲已調(diào)信號解調(diào)yyn1=demod(yn1,30000,fs,'am');yyn2=demod(yn2,30000,fs,'am');yyn3=demod(yn3,30000,fs,'am');yyn4=demod(

38、yn4,30000,fs,'am');dy1=yn-y; %高斯白噪聲snr1=var(y)/var(dy1); %輸入信噪比dy2=yyn-y1; %解調(diào)后噪聲snr2=var(y1)/var(dy2); %輸出信噪比 dy11=yn1-y; snr11=var(y)/var(dy11); dy21=yyn1-y1; snr21=var(y1)/var(dy21); dy12=yn2-y; snr12=var(y)/var(dy12); dy22=yyn2-y1; snr22=var(y1)/var(dy22); dy13=yn3-y; snr13=var(y)/var(dy

39、13); dy23=yyn3-y1; snr23=var(y1)/var(dy23); dy14=yn4-y; snr14=var(y)/var(dy14); dy24=yyn4-y1; snr24=var(y1)/var(dy24); in=snr1,snr11,snr12,snr13,snr14; out=snr2,snr21,snr22,snr23,snr24; figure(1)plot(in,out,'*')hold onplot(in,out)xlabel('輸入信噪比')ylabel('輸入信噪比') figure(2)subplot(311)plot(t,yyn)xlabel('加噪聲解調(diào)信號波形')ff=fft(yyn,n);mag=(abs(ff);f=(0:length(ff)-1).*fs/length(ff);subplo