信號與系統(tǒng)_課程設(shè)計(jì)基于Matlab的AM調(diào)制系統(tǒng)仿真

1、計(jì)算機(jī)科學(xué)系信號與系統(tǒng)課程設(shè)計(jì)目目 錄錄實(shí)驗(yàn)題目：基于Matlab的AM調(diào)制系統(tǒng)仿真成員：xx指導(dǎo)教師：xx2010-2011年度第二學(xué)期一、實(shí)驗(yàn)類型設(shè)計(jì)性實(shí)驗(yàn)設(shè)計(jì)性實(shí)驗(yàn)二、實(shí)驗(yàn)?zāi)康?.掌握振幅調(diào)制和解調(diào)原理。掌握振幅調(diào)制和解調(diào)原理。2.學(xué)會學(xué)會Matlab仿真軟件在振幅調(diào)制和解調(diào)中的應(yīng)用。仿真軟件在振幅調(diào)制和解調(diào)中的應(yīng)用。3.掌握參數(shù)設(shè)置方法和性能分析方法。掌握參數(shù)設(shè)置方法和性能分析方法。4.通過實(shí)驗(yàn)中波形的變換，學(xué)會分析實(shí)驗(yàn)現(xiàn)象。通過實(shí)驗(yàn)中波形的變換，學(xué)會分析實(shí)驗(yàn)現(xiàn)象。三、 實(shí)驗(yàn)內(nèi)容1.設(shè)計(jì)AM信號實(shí)現(xiàn)的Matlab程序，輸出調(diào)制信號、載波信號以及已調(diào)號波形以及頻譜圖，并改變參數(shù)觀察信號變

2、化情況，進(jìn)行實(shí)驗(yàn)分析。2.設(shè)計(jì)AM信號解調(diào)實(shí)現(xiàn)的Matlab程序，輸出并觀察解調(diào)信號波形，分析實(shí)驗(yàn)現(xiàn)象。四、 實(shí)驗(yàn)要求利用利用Matlab軟件進(jìn)行振幅調(diào)制和解調(diào)程序設(shè)計(jì)，輸出顯軟件進(jìn)行振幅調(diào)制和解調(diào)程序設(shè)計(jì)，輸出顯示調(diào)制信號、載波信號以及已調(diào)信號波形，并輸出顯示示調(diào)制信號、載波信號以及已調(diào)信號波形，并輸出顯示三種信號頻譜圖。對產(chǎn)生波形進(jìn)行分析，并通過參數(shù)的三種信號頻譜圖。對產(chǎn)生波形進(jìn)行分析，并通過參數(shù)的改變，觀察波形變化，分析實(shí)驗(yàn)現(xiàn)象。改變，觀察波形變化，分析實(shí)驗(yàn)現(xiàn)象。五、振幅調(diào)制原理五、振幅調(diào)制原理5.1振幅調(diào)制產(chǎn)生原理振幅調(diào)制產(chǎn)生原理 所謂調(diào)制，就是在傳送信號的一方將所要傳送的信號附加在高

3、頻振蕩所謂調(diào)制，就是在傳送信號的一方將所要傳送的信號附加在高頻振蕩上，再由天線發(fā)射出去。這里高頻振蕩波就是攜帶信號的運(yùn)載工具，也叫上，再由天線發(fā)射出去。這里高頻振蕩波就是攜帶信號的運(yùn)載工具，也叫載波。振幅調(diào)制，就是由調(diào)制信號去控制高頻載波的振幅，直至隨調(diào)制信載波。振幅調(diào)制，就是由調(diào)制信號去控制高頻載波的振幅，直至隨調(diào)制信號做線性變化。在線性調(diào)制系列中，最先應(yīng)用的一種幅度調(diào)制是全調(diào)幅或號做線性變化。在線性調(diào)制系列中，最先應(yīng)用的一種幅度調(diào)制是全調(diào)幅或常規(guī)調(diào)幅，簡稱為調(diào)幅（常規(guī)調(diào)幅，簡稱為調(diào)幅（AM）。在頻域中已調(diào)波頻譜是基帶調(diào)制信號頻）。在頻域中已調(diào)波頻譜是基帶調(diào)制信號頻譜的線性位移；在時(shí)域中，已

4、調(diào)波包絡(luò)與調(diào)制信號波形呈線性關(guān)系。譜的線性位移；在時(shí)域中，已調(diào)波包絡(luò)與調(diào)制信號波形呈線性關(guān)系。設(shè)正弦載波為設(shè)正弦載波為：式中，式中，A為載波幅度；為載波幅度； 為載波角頻率；為載波角頻率； 為載波初始相位為載波初始相位(假設(shè)假設(shè) =0).調(diào)制信號（基帶信號）為調(diào)制信號（基帶信號）為 。根據(jù)調(diào)制的定義，振幅調(diào)制信號（已調(diào)。根據(jù)調(diào)制的定義，振幅調(diào)制信號（已調(diào)信號）一般可以表示為信號）一般可以表示為 ： 設(shè)調(diào)制信號設(shè)調(diào)制信號 的頻譜為的頻譜為 ，則已調(diào)信號，則已調(diào)信號 的頻譜的頻譜 ： 五、振幅調(diào)制原理五、振幅調(diào)制原理5.2調(diào)幅電路方案分析調(diào)幅電路方案分析標(biāo)準(zhǔn)調(diào)幅波（標(biāo)準(zhǔn)調(diào)幅波（AM）產(chǎn)生原理）產(chǎn)生

5、原理 調(diào)制信號是只來來自信源的調(diào)制信號（基帶信號），這些信號可以是模擬調(diào)制信號是只來來自信源的調(diào)制信號（基帶信號），這些信號可以是模擬的，亦可以是數(shù)字的。為首調(diào)制的高頻振蕩信號可稱為載波，它可以是正的，亦可以是數(shù)字的。為首調(diào)制的高頻振蕩信號可稱為載波，它可以是正弦波，亦可以是非正弦波弦波，亦可以是非正弦波(如周期性脈沖序列如周期性脈沖序列)。載波由高頻信號源直接產(chǎn)。載波由高頻信號源直接產(chǎn)生即可，然后經(jīng)過高頻功率放大器進(jìn)行放大，作為調(diào)幅波的載波，調(diào)制信生即可，然后經(jīng)過高頻功率放大器進(jìn)行放大，作為調(diào)幅波的載波，調(diào)制信號由低頻信號源直接產(chǎn)生，二者經(jīng)過乘法器后即可產(chǎn)生雙邊帶的調(diào)幅波。號由低頻信號源直接

6、產(chǎn)生，二者經(jīng)過乘法器后即可產(chǎn)生雙邊帶的調(diào)幅波。設(shè)載波信號的表達(dá)式為設(shè)載波信號的表達(dá)式為 ，調(diào)制信號的表達(dá)式為，調(diào)制信號的表達(dá)式為 ，則調(diào)幅信號的表達(dá)式為則調(diào)幅信號的表達(dá)式為 圖5.1 標(biāo)準(zhǔn)調(diào)幅波示意圖五、振幅調(diào)制原理五、振幅調(diào)制原理5.3信號解調(diào)思路信號解調(diào)思路 從高頻已調(diào)信號中恢復(fù)出調(diào)制信號的過程稱為解調(diào)從高頻已調(diào)信號中恢復(fù)出調(diào)制信號的過程稱為解調(diào)(demodulation )，又稱為檢波又稱為檢波(detection )。對于振幅調(diào)制信號，解調(diào)。對于振幅調(diào)制信號，解調(diào)(demodulation )就是就是從它的幅度變化上提取調(diào)制信號的過程。解調(diào)從它的幅度變化上提取調(diào)制信號的過程。解調(diào)(de

7、modulation )是調(diào)制的是調(diào)制的逆過程。逆過程。 可利用乘積型同步檢波器實(shí)現(xiàn)振幅的解調(diào)，讓已調(diào)信號與本地恢復(fù)載可利用乘積型同步檢波器實(shí)現(xiàn)振幅的解調(diào)，讓已調(diào)信號與本地恢復(fù)載波信號相乘并通過低通濾波可獲得解調(diào)信號。波信號相乘并通過低通濾波可獲得解調(diào)信號。六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.1載波信號與調(diào)制信號分析載波信號與調(diào)制信號分析% =載波信號載波信號=t=-1:0.00001:1; A0=10; %載波信號振幅載波信號振幅f=6000; %載波信號頻率載波信號頻率w0=f*pi;Uc=A0*cos(w0*t); %載波信號載波信號figure(1); subpl

8、ot(2,1,1);plot(t,Uc);title(載頻信號波形載頻信號波形);axis(0,0.01,-15,15);subplot(2,1,2);Y1=fft(Uc); %對載波信號進(jìn)行傅里葉變換對載波信號進(jìn)行傅里葉變換 plot(abs(Y1);title(載波信號頻譜載波信號頻譜);axis(5800,6200,0,1000000); 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.1載波信號與調(diào)制信號分析載波信號與調(diào)制信號分析 (載波信號載波信號) 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.1載波信號與調(diào)制信號分析載波信號與調(diào)制信號分析 % =調(diào)制信號調(diào)制

9、信號=t=-1:0.00001:1; A1=5; %調(diào)制信號振幅調(diào)制信號振幅f=6000; %載波信號頻率載波信號頻率w0=f*pi;mes=A1*cos(0.001*w0*t); %調(diào)制信號調(diào)制信號subplot(2,1,1);plot(t,mes);xlabel(t),title(調(diào)制信號調(diào)制信號); subplot(2,1,2);Y2=fft(mes); % 對調(diào)制信號進(jìn)行傅里葉變換對調(diào)制信號進(jìn)行傅里葉變換 plot(abs(Y2);title(調(diào)制信號頻譜調(diào)制信號頻譜);axis(198000,202000,0,1000000); 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)

10、6.1載波信號與調(diào)制信號分析載波信號與調(diào)制信號分析 (調(diào)制信號調(diào)制信號) 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.1載波信號與調(diào)制信號分析載波信號與調(diào)制信號分析 % =AM已調(diào)信號已調(diào)信號=t=-1:0.00001:1; A0=10; %載波信號振幅載波信號振幅A1=5; %調(diào)制信號振幅調(diào)制信號振幅A2=3; %已調(diào)信號振幅已調(diào)信號振幅f=3000; %載波信號頻率載波信號頻率w0=2*f*pi; m=0.15; %調(diào)制度調(diào)制度mes=A1*cos(0.001*w0*t); %消調(diào)制信號消調(diào)制信號Uam=A2*(1+m*mes).*cos(w0).*t); %AM 已調(diào)信號已

11、調(diào)信號 subplot(2,1,1); plot(t,Uam);grid on; title(AM調(diào)制信號波形調(diào)制信號波形);subplot(2,1,2); Y3=fft(Uam); % 對對AM已調(diào)信號進(jìn)行傅里葉變換已調(diào)信號進(jìn)行傅里葉變換 plot(abs(Y3),grid; title(AM調(diào)制信號頻譜調(diào)制信號頻譜);axis(5950,6050,0,500000); 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.1載波信號與調(diào)制信號分析載波信號與調(diào)制信號分析 (AM已調(diào)信號已調(diào)信號) 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.2 設(shè)計(jì)設(shè)計(jì)FIR數(shù)字低通濾波器數(shù)

12、字低通濾波器 FIR濾波器比鞥采用間接法，常用的方法有窗函數(shù)法、頻率采樣法和濾波器比鞥采用間接法，常用的方法有窗函數(shù)法、頻率采樣法和切比雪夫等波紋逼近法。對于線性相位濾波器，經(jīng)常采用切比雪夫等波紋逼近法。對于線性相位濾波器，經(jīng)常采用FIR濾波器。濾波器。 對于數(shù)字高通、帶通濾波器的設(shè)計(jì)，通用方法為雙線性變換法?？梢越柚鷮τ跀?shù)字高通、帶通濾波器的設(shè)計(jì)，通用方法為雙線性變換法。可以借助于模擬濾波器的頻率轉(zhuǎn)換設(shè)計(jì)一個(gè)所需類型的過渡模擬濾波器，再經(jīng)過雙于模擬濾波器的頻率轉(zhuǎn)換設(shè)計(jì)一個(gè)所需類型的過渡模擬濾波器，再經(jīng)過雙線性變換將其轉(zhuǎn)換策劃那個(gè)所需的數(shù)字濾波器。具體設(shè)計(jì)步驟如下：線性變換將其轉(zhuǎn)換策劃那個(gè)所需

13、的數(shù)字濾波器。具體設(shè)計(jì)步驟如下：（1）確定所需類型數(shù)字濾波器的技術(shù)指標(biāo)。）確定所需類型數(shù)字濾波器的技術(shù)指標(biāo)。（2）將所需類型數(shù)字濾波器的邊界頻率轉(zhuǎn)換成相應(yīng)的模擬濾波器的邊界）將所需類型數(shù)字濾波器的邊界頻率轉(zhuǎn)換成相應(yīng)的模擬濾波器的邊界頻率，轉(zhuǎn)換公式為頻率，轉(zhuǎn)換公式為=2/T tan(0.5)(3)將相應(yīng)類型的模擬濾波器技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器技術(shù)指標(biāo)。將相應(yīng)類型的模擬濾波器技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器技術(shù)指標(biāo)。（4）設(shè)計(jì)模擬低通濾波器。）設(shè)計(jì)模擬低通濾波器。（5）通過頻率變換將模擬低通轉(zhuǎn)換成相應(yīng)類型的過渡模擬濾波器。）通過頻率變換將模擬低通轉(zhuǎn)換成相應(yīng)類型的過渡模擬濾波器。（6）采用雙線性變換

14、法將相應(yīng)類型的過渡模擬濾波器轉(zhuǎn)換成所需類型的）采用雙線性變換法將相應(yīng)類型的過渡模擬濾波器轉(zhuǎn)換成所需類型的數(shù)字濾波器。數(shù)字濾波器。 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.2 設(shè)計(jì)設(shè)計(jì)FIR數(shù)字低通濾波器數(shù)字低通濾波器 我們知道，脈沖響應(yīng)不變法的主要缺點(diǎn)是會產(chǎn)生頻譜混疊現(xiàn)象，使數(shù)我們知道，脈沖響應(yīng)不變法的主要缺點(diǎn)是會產(chǎn)生頻譜混疊現(xiàn)象，使數(shù)字濾波器的頻響偏離模擬濾波器的頻響特性。為了克服之一缺點(diǎn)，可以采字濾波器的頻響偏離模擬濾波器的頻響特性。為了克服之一缺點(diǎn)，可以采用雙線性變換法。用雙線性變換法。下面我們介紹用窗函數(shù)法設(shè)計(jì)下面我們介紹用窗函數(shù)法設(shè)計(jì)FIR濾波器的步驟。如下：濾波器

15、的步驟。如下：（1）根據(jù)對阻帶衰減及過渡帶的指標(biāo)要求，選擇串窗數(shù)類型（矩形窗、）根據(jù)對阻帶衰減及過渡帶的指標(biāo)要求，選擇串窗數(shù)類型（矩形窗、三角窗、漢寧窗、哈明窗、凱塞窗等），并估計(jì)窗口長度三角窗、漢寧窗、哈明窗、凱塞窗等），并估計(jì)窗口長度N。先按照阻帶。先按照阻帶衰減選擇窗函數(shù)類型。原則是在保證阻帶衰減滿足要求的情況下，盡量選衰減選擇窗函數(shù)類型。原則是在保證阻帶衰減滿足要求的情況下，盡量選擇主瓣的窗函數(shù)。擇主瓣的窗函數(shù)。（2）構(gòu)造希望逼近的頻率響應(yīng)函數(shù)。）構(gòu)造希望逼近的頻率響應(yīng)函數(shù)。（3）計(jì)算）計(jì)算h(n).。（4）加窗得到設(shè)計(jì)結(jié)果。）加窗得到設(shè)計(jì)結(jié)果。 六六.AM基于基于matlab的調(diào)制與

16、解調(diào)的調(diào)制與解調(diào)6.2 設(shè)計(jì)設(shè)計(jì)FIR數(shù)字低通濾波器數(shù)字低通濾波器 接下來，我們根據(jù)語音信號的特點(diǎn)給出有關(guān)濾波器的技術(shù)指標(biāo)：接下來，我們根據(jù)語音信號的特點(diǎn)給出有關(guān)濾波器的技術(shù)指標(biāo)：低通濾波器的性能指標(biāo)：低通濾波器的性能指標(biāo)：通帶邊界頻率通帶邊界頻率fp=300Hz，阻帶截止頻率，阻帶截止頻率fc=320Hz，阻帶最小衰減，阻帶最小衰減As=100db , 通帶最大衰減通帶最大衰減Ap=1dB在在Matlab中，可以利用函數(shù)中，可以利用函數(shù)fir1設(shè)計(jì)設(shè)計(jì)FIR濾波器，利用函數(shù)濾波器，利用函數(shù)butter,cheby1和和ellip設(shè)計(jì)設(shè)計(jì)IIR濾波器，濾波器，利用利用Matlab中的函數(shù)中的函

17、數(shù)freqz畫出各步步器的頻率響應(yīng)。畫出各步步器的頻率響應(yīng)。hn=fir1(M，wc，window)，可以指定窗，可以指定窗函數(shù)向量函數(shù)向量window。如果缺省。如果缺省window參數(shù)，則參數(shù)，則fir1默認(rèn)為哈明窗。其中可選的窗函數(shù)有默認(rèn)為哈明窗。其中可選的窗函數(shù)有Rectangular Barlrtt Hamming Hann Blackman窗，其相應(yīng)的都有實(shí)現(xiàn)函數(shù)。窗，其相應(yīng)的都有實(shí)現(xiàn)函數(shù)。MATLAB信號處理工具箱函數(shù)信號處理工具箱函數(shù)buttp buttor butter是巴特沃斯濾波器設(shè)計(jì)函數(shù)，其有是巴特沃斯濾波器設(shè)計(jì)函數(shù)，其有5種調(diào)用格種調(diào)用格式，本課程設(shè)計(jì)中用到的是式，本

18、課程設(shè)計(jì)中用到的是N,wc=butter(N,wc,Rp,As,s),該格式用于計(jì)算巴特沃斯模擬濾波該格式用于計(jì)算巴特沃斯模擬濾波器的階數(shù)器的階數(shù)N和和3dB截止頻率截止頻率wc。MATLAB信號處理工具箱函數(shù)信號處理工具箱函數(shù)cheblap,cheblord和和cheeby1是切比雪夫是切比雪夫I型濾波器設(shè)計(jì)函數(shù)。我型濾波器設(shè)計(jì)函數(shù)。我們用到的是們用到的是cheeby1函數(shù)，其調(diào)用格式如下：函數(shù)，其調(diào)用格式如下：B,A=cheby1(N,Rp,wpo,ftypr)B,A=cheby1(N,Rp,wpo,ftypr,s) 函數(shù)函數(shù)butter,cheby1和和ellip設(shè)計(jì)設(shè)計(jì)IIR濾波器時(shí)都

19、是默認(rèn)的雙線性變換法，所以在設(shè)計(jì)濾波器時(shí)只濾波器時(shí)都是默認(rèn)的雙線性變換法，所以在設(shè)計(jì)濾波器時(shí)只需要代入相應(yīng)的實(shí)現(xiàn)函數(shù)即可。下面我們將給出需要代入相應(yīng)的實(shí)現(xiàn)函數(shù)即可。下面我們將給出FIR數(shù)字濾波器的主要程序。數(shù)字濾波器的主要程序。 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.2 設(shè)計(jì)設(shè)計(jì)FIR數(shù)字低通濾波器數(shù)字低通濾波器%=FIR低通濾波器低通濾波器=Ft=2000; %采樣頻率采樣頻率fpts=100 120; %通帶邊界頻率通帶邊界頻率fp=100Hz，阻帶截止頻率，阻帶截止頻率fs=120Hzmag=1 0; dev=0.01 0.05; %通帶波動(dòng)通帶波動(dòng)1%，阻帶波動(dòng)，阻

20、帶波動(dòng)5%n21,wn21,beta,ftype=kaiserord(fpts,mag,dev,Ft);%kaiserord估計(jì)采估計(jì)采用凱塞窗設(shè)計(jì)的用凱塞窗設(shè)計(jì)的FIR濾濾 波器的參數(shù)波器的參數(shù)b21=fir1(n21,wn21,Kaiser(n21+1,beta); %由由fir1設(shè)計(jì)濾波器設(shè)計(jì)濾波器h,w=freqz(b21,1); %得到頻率響應(yīng)得到頻率響應(yīng)plot(w/pi,abs(h);grid ontitle(FIR低通濾波器低通濾波器); 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.2 設(shè)計(jì)設(shè)計(jì)FIR數(shù)字低通濾波器數(shù)字低通濾波器 (FIR低通濾波器低通濾波器) 六

21、六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.3 AM解調(diào)解調(diào)%=AM信號解調(diào)信號解調(diào)=t=-1:0.00001:1; A0=10; %載波信號振幅載波信號振幅A1=5; %調(diào)制信號振幅調(diào)制信號振幅A2=3; %已調(diào)信號振幅已調(diào)信號振幅f=3000; %載波信號頻率載波信號頻率w0=2*f*pi; m=0.15; %調(diào)制度調(diào)制度mes=A1*cos(0.001*w0*t); %調(diào)制信號調(diào)制信號Uam=A2*(1+m*mes).*cos(w0).*t); %AM 已調(diào)信號已調(diào)信號Dam=Uam.*cos(w0*t); %對對AM調(diào)制信號進(jìn)行解調(diào)調(diào)制信號進(jìn)行解調(diào)subplot(2,1,1

22、); plot(t,Dam); grid on;title(濾波前濾波前AM解調(diào)信號波形解調(diào)信號波形); subplot(2,1,2); Y5=fft(Dam); plot(abs(Y5),grid; title(濾波前濾波前AM解調(diào)信號頻譜解調(diào)信號頻譜);axis(187960,188040,0,200000); 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.3 AM解調(diào)解調(diào) (AM解調(diào)信號解調(diào)信號) 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.3 AM解調(diào)解調(diào)%=AM解調(diào)信號解調(diào)信號FIR濾波濾波=t=-1:0.00001:1;A0=10; %載波信號振幅載波信號

23、振幅A1=5; %調(diào)制信號振幅調(diào)制信號振幅A2=3; %已調(diào)信號振幅已調(diào)信號振幅 *f=6000; %載波信號頻率載波信號頻率w0=f*pi; m=0.15; %調(diào)制度調(diào)制度 *Uc=A0.*cos(w0*t); %載波信號載波信號subplot(5,2,1); plot(t,Uc);title(載波信號載波信號); axis(0,0.01,-15,15);T1=fft(Uc); %傅里葉變換傅里葉變換subplot(5,2,2); plot(abs(T1);title(載波信號頻譜載波信號頻譜); axis(5800,6200,0,1000000); 六六.AM基于基于matlab的調(diào)制與解

24、調(diào)的調(diào)制與解調(diào)6.3 AM解調(diào)解調(diào)mes=A1*cos(0.001*w0*t); %調(diào)制信號調(diào)制信號subplot(5,2,3); plot(t,mes);title(調(diào)制信號調(diào)制信號);T2=fft(mes); subplot(5,2,4);plot(abs(T2); title(調(diào)制信號頻譜調(diào)制信號頻譜);axis(198000,202000,0,2000000);Uam=A2*(1+m*mes).*cos(w0).*t); %AM 已調(diào)信號已調(diào)信號 *subplot(5,2,5); plot(t,Uam);title(已調(diào)信號已調(diào)信號);T3=fft(Uam);subplot(5,2,6

25、);plot(abs(T3);title(已調(diào)信號頻譜已調(diào)信號頻譜);axis(5950,6050,0,500000); 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.3 AM解調(diào)解調(diào)Dam=Uam.*cos(w0*t); %對對AM已調(diào)信號進(jìn)行解調(diào)已調(diào)信號進(jìn)行解調(diào)subplot(5,2,7); plot(t,Dam);title(濾波前的濾波前的AM解調(diào)信號波形解調(diào)信號波形);T4=fft(Dam); %求求AM信號的頻譜信號的頻譜subplot(5,2,8); plot(abs(T4);title(濾波前的濾波前的AM解調(diào)信號頻譜解調(diào)信號頻譜);axis(187960,1880

26、40,0,200000);z21=fftfilt(b21,Dam); %FIR低通濾波低通濾波subplot(5,2,9); plot(t,z21,r);title(濾波后的濾波后的AM解調(diào)信號波形解調(diào)信號波形);T5=fft(z21); %求求AM信號的頻譜信號的頻譜subplot(5,2,10); plot(abs(T5),r);title(濾波后的濾波后的AM解調(diào)信號頻譜解調(diào)信號頻譜);axis(198000,202000,0,200000); (AM解調(diào)信號解調(diào)信號FIR低通濾波低通濾波)(AM解調(diào)信號解調(diào)信號FIR低通濾波低通濾波)六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解

27、調(diào)6.3 AM解調(diào)解調(diào)%=AM解調(diào)信號解調(diào)信號FIR濾波濾波=加噪加噪=t=-1:0.00001:1;A0=10; A1=5; A2=3; f=6000; w0=f*pi; m=0.15; Uc=A0.*cos(w0*t); subplot(5,2,1); plot(t,Uc); title(載波信號載波信號); axis(0,0.01,-15,15); T1=fft(Uc); %傅里葉變換傅里葉變換subplot(5,2,2); plot(abs(T1);title(載波信號頻譜載波信號頻譜); axis(5800,6200,0,1000000);mes=A1*cos(0.001*w0*t)

28、; subplot(5,2,3); plot(t,mes);title(調(diào)制信號調(diào)制信號); T2=fft(mes); subplot(5,2,4); plot(abs(T2);title(調(diào)制信號頻譜調(diào)制信號頻譜); axis(198000,202000,0,2000000);Uam1=A2*(1+m*mes).*cos(w0).*t); %AM 已調(diào)信號已調(diào)信號subplot(5,2,5); plot(t,Uam1);title(已調(diào)信號已調(diào)信號); T3=fft(Uam1); subplot(5,2,6); plot(abs(T3);title(已調(diào)信號頻譜已調(diào)信號頻譜); axis(5

29、950,6050,0,500000); 六六.AM基于基于matlab的調(diào)制與解調(diào)的調(diào)制與解調(diào)6.3 AM解調(diào)解調(diào) (AM解調(diào)信號解調(diào)信號FIR濾波濾波=加噪加噪= 續(xù)續(xù))sn1=20; %信噪比信噪比db1=A12/(2*(10(sn1/10); %計(jì)算對應(yīng)噪聲方差計(jì)算對應(yīng)噪聲方差n1=sqrt(db1)*randn(size(t); %生成高斯白噪聲生成高斯白噪聲Uam=n1+Uam1; Dam=Uam.*cos(w0*t); %對對AM已調(diào)信號進(jìn)行解調(diào)已調(diào)信號進(jìn)行解調(diào)subplot(5,2,7); plot(t,Dam);title(濾波前的濾波前的AM解調(diào)信號波形解調(diào)信號波形); T4=fft(Dam); %求求AM信號的頻譜信號的頻譜subplot(5,2,8); plot(abs(T4);title(濾波前的濾波前的AM解調(diào)信號頻譜解調(diào)信號頻譜); axis(187960,188040,0,200000);z21=fftfilt(b21,Dam); %FIR低通濾波低通濾波subplot(5,2,9); plot(t,z2