信號(hào)與通信綜合設(shè)計(jì)項(xiàng)目報(bào)告-AM調(diào)制與解調(diào)系統(tǒng)的設(shè)計(jì)

信號(hào)與通信綜合設(shè)計(jì)項(xiàng)目報(bào)告信號(hào)與通信綜合設(shè)計(jì)項(xiàng)目報(bào)告--17-設(shè)計(jì)目的及內(nèi)容設(shè)計(jì)一個(gè)模擬系統(tǒng)，實(shí)現(xiàn)AM調(diào)制與解調(diào)，通過MATLAB仿真了解AM的工作原理，再根據(jù)給定的技術(shù)指標(biāo)通過程序設(shè)計(jì)實(shí)現(xiàn)系統(tǒng)仿真。用MATLAB產(chǎn)生一個(gè)頻率為20Hz的信號(hào)源，讓載波頻率為150Hz，A=2，仿真出：已調(diào)信號(hào)、正常調(diào)幅狀態(tài)、臨界調(diào)幅狀態(tài)、過調(diào)幅狀態(tài)及解調(diào)信號(hào)。1.1課程設(shè)計(jì)目的本次課程設(shè)計(jì)的內(nèi)容是利用MATLAB研究AM調(diào)制與解調(diào)模擬系統(tǒng)的理論設(shè)計(jì)和軟件仿真方法。通過完成本次課程設(shè)計(jì)，主要達(dá)到以下幾個(gè)目的：1.掌握模擬系統(tǒng)AM調(diào)制與解調(diào)的原理；2.掌握AM調(diào)制與解調(diào)模擬系統(tǒng)的理論設(shè)計(jì)方法；3.掌握應(yīng)用MATLAB分析系統(tǒng)時(shí)域、頻域特性的方法，進(jìn)一步鍛煉應(yīng)用MATLAB進(jìn)行編程仿真的能力。1.2課程設(shè)計(jì)內(nèi)容振幅調(diào)制，也可簡稱為調(diào)幅，AM（AmplitudeModulation）,通過改變輸出信號(hào)的振幅，來實(shí)現(xiàn)傳送信息的目的，一般在調(diào)制端輸出的高信號(hào)的幅度變化與原始信號(hào)成一定的函數(shù)關(guān)系，在解調(diào)端進(jìn)行解調(diào)并輸出原始低頻信號(hào)。并且利用MATLAB仿真出不同狀態(tài)下的波形與頻譜。分為以下三種狀態(tài)：正常調(diào)幅狀態(tài)下的波形與頻譜（Am<A0）；臨界調(diào)幅狀態(tài)下的波形與頻譜（Am=A0）；過調(diào)幅狀態(tài)下的波形與頻譜（Am>A0）。

課程設(shè)計(jì)基本原理調(diào)制就是對(duì)信號(hào)源的信息進(jìn)行處理，使其變?yōu)檫m合于信道傳輸?shù)男问降倪^程，一般來說，信號(hào)源的信息（也稱為信源）含有直流分量和頻率較低的頻率分量，稱為基帶信號(hào)?；鶐盘?hào)往往不能作為傳輸信號(hào)，因此必須把基帶信號(hào)轉(zhuǎn)變?yōu)橐粋€(gè)相對(duì)基帶頻率而言頻率非常高的信號(hào)以適合于信道傳輸。這個(gè)信號(hào)叫做已調(diào)信號(hào)，而基帶信號(hào)叫做調(diào)制信號(hào)，調(diào)制是通過改變高頻載波即消息的載體信號(hào)的幅度、相位或者頻率，使其隨著基帶信號(hào)幅度的變化而變化來實(shí)現(xiàn)的。而解調(diào)則是將基帶信號(hào)從載波中提取出來以便預(yù)定的接受者（也成為信宿）處理和理解的過程。2.1AM調(diào)制原理AM調(diào)制采用集電極調(diào)幅，所謂的集電極調(diào)幅就是用調(diào)制信號(hào)來改變高頻功率放大器的集電極直流電壓，以實(shí)現(xiàn)調(diào)幅。其中的載波可以是正弦波或脈沖序列，以正弦型信號(hào)作載波的調(diào)制叫做連續(xù)波調(diào)制，調(diào)制后的載波就載有調(diào)制信號(hào)所包含的信息，稱為已調(diào)波。幅度調(diào)制使用調(diào)制信號(hào)去控制高頻載波的振幅，使其按照調(diào)制信號(hào)的規(guī)律而變化的過程。幅度調(diào)制的一般模型如REF_Ref137022466\h圖2.1所示。圖STYLEREF1\s2.SEQ圖\*ARABIC\s11幅度調(diào)制的一般模型圖STYLEREF1\s2.SEQ圖\*ARABIC\s12AM信號(hào)的波形與頻譜AM信號(hào)的典型波形和頻譜如REF_Ref137022473\h圖圖STYLEREF1\s2.SEQ圖\*ARABIC\s12AM信號(hào)的波形與頻譜2.2基本原理（1）AM調(diào)制的基本原理：圖STYLEREF1\s2.SEQ圖\*ARABIC\s13AM調(diào)制器的模型AM調(diào)制是指對(duì)信號(hào)進(jìn)行幅度調(diào)制，而幅度調(diào)制是正弦載波的幅度隨調(diào)制信號(hào)作線性變化的過程，如REF_Ref137022482\h圖2.3所示為AM調(diào)制器的模型，將調(diào)制信號(hào)加上一個(gè)直流分量，保證信號(hào)的最小值大于零，然后再和載波相乘，得到已調(diào)信號(hào)。在波形上，幅度已調(diào)信號(hào)的幅度隨基帶信號(hào)的規(guī)律而成正比的變化；在頻譜結(jié)構(gòu)上，它的頻譜完全是基帶信號(hào)頻譜在頻域內(nèi)的簡單搬移（精確到常數(shù)因子）。幅度調(diào)制信號(hào)可以表示成m(t)。對(duì)于AM來說m(t)是帶有直流分量的基帶信號(hào)，可以表示成m圖STYLEREF1\s2.SEQ圖\*ARABIC\s13AM調(diào)制器的模型（2）AM解調(diào)的基本原理：AM調(diào)幅解調(diào)的時(shí)候可以采用相干解調(diào)，解調(diào)就是調(diào)制的逆過程，實(shí)現(xiàn)頻譜搬移，將已調(diào)信號(hào)乘以載波后通過低通濾波器并在幅度上做一定調(diào)整即可恢復(fù)出原來的調(diào)制信號(hào)。另外AM信號(hào)在滿足m1(t)圖STYLEREF1\s2.SEQ圖\*ARABIC\s14相干解調(diào)器的一般模型相干解調(diào)也叫同步檢波，它由乘法器和低通濾波器組成，如REF_Ref137022489\h圖STYLEREF1\s2.SEQ圖\*ARABIC\s14相干解調(diào)器的一般模型同時(shí)相干解調(diào)適用于所有線性調(diào)制信號(hào)的解調(diào)。實(shí)現(xiàn)相干解調(diào)的關(guān)鍵是在接收端恢復(fù)出一個(gè)調(diào)制載波嚴(yán)格同步的相干載波，相干解調(diào)是指利用乘法器輸入一路與載頻相干（同頻同相）的參考信號(hào)與載頻相乘，恢復(fù)出載波性能的好壞，直接關(guān)系到接收機(jī)解調(diào)性能的優(yōu)勢(shì)。如REF_Ref136970057\h圖2.5所示為解調(diào)電路的組成框圖。圖STYLEREF1\s2.SEQ圖\*ARABIC\s15解調(diào)電路的組成框圖在現(xiàn)代科技中，乘法器是非常重要的一種電子元件，它可以實(shí)現(xiàn)數(shù)字信號(hào)的乘法運(yùn)算，廣泛應(yīng)用于各種計(jì)算機(jī)、通信、控制等領(lǐng)域。乘法器是一種數(shù)字電路，它的主要作用是實(shí)現(xiàn)數(shù)字信號(hào)的乘法運(yùn)算。在乘法器中，輸入信號(hào)經(jīng)過一系列的邏輯門電路處理后，輸出結(jié)果為兩個(gè)輸入信號(hào)的乘積。乘法器的基本原理可以用以下公式表示：A×B=C其中，A和B為輸入信號(hào)，C為輸出結(jié)果。乘法器的工作原理可以分為兩種類型：串行和并行。串行乘法器是將兩個(gè)輸入信號(hào)分別進(jìn)行位移和加法運(yùn)算，最終得到輸出結(jié)果。而并行乘法器則是將兩個(gè)輸入信號(hào)分別進(jìn)行分解和加法運(yùn)算，最終得到輸出結(jié)果。兩種類型的乘法器均采用邏輯門電路實(shí)現(xiàn)，具體實(shí)現(xiàn)方式有多種。根據(jù)乘法器的不同實(shí)現(xiàn)方式，可以將其分為以下幾種類型：a.串行乘法器串行乘法器是一種最簡單的乘法器，它采用逐位相乘的方式實(shí)現(xiàn)。串行乘法器的輸入信號(hào)經(jīng)過位移和加法運(yùn)算后，得到輸出結(jié)果。串行乘法器的優(yōu)點(diǎn)是結(jié)構(gòu)簡單，適合于低速應(yīng)用。但是，由于其逐位相乘的方式，其速度較慢，不適合于高速應(yīng)用。b.并行乘法器并行乘法器是一種較為復(fù)雜的乘法器，它采用分解和加法運(yùn)算的方式實(shí)現(xiàn)。并行乘法器的輸入信號(hào)分別進(jìn)行分解，然后進(jìn)行加法運(yùn)算，得到輸出結(jié)果。并行乘法器的優(yōu)點(diǎn)是速度快，適合于高速應(yīng)用。但是，由于其結(jié)構(gòu)較為復(fù)雜，相對(duì)于串行乘法器來說，成本較高。c.Booth乘法器Booth乘法器是一種改進(jìn)的串行乘法器，它采用位移和加減運(yùn)算的方式實(shí)現(xiàn)。Booth乘法器的輸入信號(hào)經(jīng)過位移和加減運(yùn)算后，得到輸出結(jié)果。Booth乘法器的優(yōu)點(diǎn)是速度快，適合于高速應(yīng)用。但是，由于其結(jié)構(gòu)較為復(fù)雜，相對(duì)于串行乘法器來說，成本較高。d.Wallace樹乘法器Wallace樹乘法器是一種改進(jìn)的并行乘法器，它采用分解和加法運(yùn)算的方式實(shí)現(xiàn)。Wallace樹乘法器的輸入信號(hào)分別進(jìn)行分解，然后進(jìn)行加法運(yùn)算，得到輸出結(jié)果。Wallace樹乘法器的優(yōu)點(diǎn)是速度快，適合于高速應(yīng)用。但是，由于其結(jié)構(gòu)較為復(fù)雜，相對(duì)于串行乘法器來說，成本較高?？傊?，乘法器是一種非常重要的電子元件，廣泛應(yīng)用于各種計(jì)算機(jī)、通信、控制等領(lǐng)域。隨著科技的不斷發(fā)展，乘法器的結(jié)構(gòu)和性能也在不斷改進(jìn)和提高，為各種應(yīng)用提供了更好的支持和保障。在本次課設(shè)中，應(yīng)用了串行乘法器。

課程設(shè)計(jì)與仿真3.1課程設(shè)計(jì)（1）AM的調(diào)制設(shè)載波ct=Ac*cos?(20πt)，Ac=2，信源mt=Am*cos?(πt)sS式中，ωc為載波角頻率，H若濾波器Hω為全通網(wǎng)絡(luò)，即ht=δ(t)，并假設(shè)調(diào)制信號(hào)mt的平均值為0.將sS式中，m(t)可以是確知信號(hào)，也可以是隨機(jī)信號(hào)（此時(shí)，已調(diào)信號(hào)的頻域表示必須用功率譜描述），設(shè)載波信號(hào)為uct=Ucmcosωct=UcmcosU（2）AM解調(diào)設(shè)輸入解調(diào)器的信號(hào)為：S與相干載波相乘后，得：S==經(jīng)隔直電容和低通濾波器后，輸出信號(hào)為：m信號(hào)功率為：S相干解調(diào)的關(guān)鍵是必須產(chǎn)生一個(gè)與調(diào)制器同頻同相位的載波，如果同頻同相位的條件得不到滿足，則會(huì)破壞原始信號(hào)的恢復(fù)。3.2仿真結(jié)果分析圖STYLEREF1\s3.SEQ圖\*ARABIC\s11調(diào)制信號(hào)的波形與頻譜如REF_Ref136970105\h圖3.1所示是調(diào)制信號(hào)的波形與頻譜，采樣頻率為50Hz。幅度調(diào)制的特點(diǎn)是載波的頻率始終保持不變，它的振幅卻是變化的。對(duì)于已調(diào)信號(hào)則是由原信號(hào)與載波信號(hào)相乘的結(jié)果。由于正弦函數(shù)幅度均勻，故與載波相乘后，幅度分布也比較均勻。此時(shí)再對(duì)已調(diào)信號(hào)求取頻譜，結(jié)果如圖所示，證明了調(diào)制過程中有頻譜的搬移。圖STYLEREF1\s3.SEQ圖\*ARABIC\s12正常調(diào)幅狀態(tài)下的波形與頻譜圖STYLEREF1\s3.SEQ圖\*ARABIC\s13臨界調(diào)幅狀態(tài)下的波形與頻譜圖STYLEREF1\s3.SEQ圖\*ARABIC\s14過調(diào)幅狀態(tài)下的波形與頻譜圖STYLEREF1\s3.SEQ圖\*ARABIC\s15相干解調(diào)后的波形與頻譜如REF_Ref136970230\h圖3.2所示是正常調(diào)幅狀態(tài)下的波形與頻譜，載波ct=Ac*cos?(20πt)，設(shè)定Ac=2，Am=2，自定義直流分量A0，令A(yù)0=2.1，則A0>Am，利用程序畫出正常調(diào)幅狀態(tài)下的波形；如REF_Ref136970238\h圖3.3所示是臨界調(diào)幅狀態(tài)下的波形與頻譜，令A(yù)0=2，可以看出臨界調(diào)幅狀態(tài)下產(chǎn)生了失真；如REF_Ref136970244\h圖3.4所示是過調(diào)幅狀態(tài)下的波形與頻譜，令A(yù)0=0.5，波形圖中產(chǎn)生了頻譜混疊形成了失真。如REF_Ref136970253\h圖3.5所示是相干解調(diào)后的波形與頻譜，即將已調(diào)信號(hào)乘以同頻同相的本地載波通過低通濾波器作相干解調(diào)，即可恢復(fù)出原始信號(hào)，在圖中可以看出解調(diào)后的信號(hào)與原始信號(hào)相同。

總結(jié)本次課設(shè)研究了AM信號(hào)的調(diào)制與解調(diào)技術(shù)，包括其原理、實(shí)現(xiàn)方法與仿真結(jié)果測(cè)試。在實(shí)驗(yàn)中，通過改變A0的值，可以得到正常調(diào)幅、臨界調(diào)幅和過調(diào)幅三種不同狀態(tài)下的波形。接著進(jìn)行了解調(diào)實(shí)驗(yàn)，通過解調(diào)器將AM調(diào)制信號(hào)進(jìn)行解調(diào)，恢復(fù)出原始的低頻信號(hào)，同時(shí)，當(dāng)解調(diào)參數(shù)設(shè)置合理時(shí)，解調(diào)后的信號(hào)波形和頻譜與原始低頻信號(hào)基本一致。并且，我們深入了解了AM信號(hào)調(diào)制與解調(diào)的原理和方法，同時(shí)也了解到AM調(diào)制與解調(diào)的一些特點(diǎn)和應(yīng)用場(chǎng)景。例如：AM調(diào)制信號(hào)具有較好的傳輸性能和抗干擾能力，適用于廣播、電視等通信領(lǐng)域；而AM解調(diào)技術(shù)也廣泛應(yīng)用于語音信號(hào)處理、音頻信號(hào)處理等領(lǐng)域。通過這次課設(shè)，讓我感受到了MATLAB的強(qiáng)大，也鍛煉了實(shí)際運(yùn)用的能力，也感謝老師的支持與幫助，但由于經(jīng)驗(yàn)的匱乏，難免會(huì)有許多考慮不周全的地方，之后會(huì)繼續(xù)改善。

參考文獻(xiàn)張輝,曹麗娜.現(xiàn)代通信原理與技術(shù)[M].西安電子科技大學(xué)出版社,2018樊昌信.通信原理[M].(第1版至第5版).北京；北京國防工業(yè)出版社,2005胡廣書.現(xiàn)代信號(hào)處理[M].北京；北京清華大學(xué)出版社,2004胡宴如,耿蘇燕.高頻電子線路[M].北京；高等教育出版社,2015羅明.數(shù)字通信信號(hào)的自動(dòng)識(shí)別與參數(shù)估計(jì)研究[D].西安；西安電子科技大學(xué),2005陳后金.信號(hào)與系統(tǒng)[M].北京；高等教育出版社,2006張智星.《MATLAB程序設(shè)計(jì)與應(yīng)用》.北京；清華大學(xué)出版社.2002吳艷君.AM系統(tǒng)的建模與仿真[J].電子技術(shù).上海；2012.

附錄仿真程序代碼clc;clear;closeall;%AM調(diào)制過程figure(1)%打開圖形窗口并命名a0=3;%增加直流成分f0=20;%信號(hào)頻率fc=50;%載波頻率fs=1000;%采樣頻率snr=100;%信噪比figure(1);t=(-20:0.001:20);%時(shí)間段f=(0:40000)*fs/40001-fs/2;%頻率段W0=2*pi*f0;%信號(hào)角頻率Wc=2*pi*fc;%載波角頻率am1=cos(2*pi*f0*t);%信號(hào)波形（信號(hào)幅度直接取1）AM1=fft(am1);%信號(hào)傅里葉變換am=a0+am1;%調(diào)制后總的AM幅度（直流加信號(hào)）t1=cos(2*pi*fc*t);%載波波形T1=fft(t1);%載波傅里葉變換s_am=am.*t1;%調(diào)制信號(hào)S_AM=fft(s_am);%調(diào)制后波形傅里葉變換subplot(3,2,1);plot(t(19801:20200),am1(19801:20200));title('信號(hào)時(shí)域圖');subplot(3,2,2);plot(f,fftshift(abs(AM1)));title('信號(hào)頻譜圖');subplot(3,2,3);plot(t(19801:20200),t1(19801:20200));title('載波時(shí)域圖');subplot(3,2,4);plot(f,fftshift(abs(T1)));title('載波頻譜圖');subplot(3,2,5);plot(t(19801:20200),s_am(19801:20200));title('已調(diào)信號(hào)時(shí)域圖');subplot(3,2,6);plot(f,fftshift(abs(S_AM)));title('已調(diào)信號(hào)頻譜圖');%產(chǎn)生噪聲figure(2);y=awgn(s_am,snr);%添加加性高斯白噪聲subplot(2,2,1);plot(t(19801:20200),y(19801:20200));title('添加噪聲后信號(hào)波形圖');Y=fft(y);%添加噪聲后信號(hào)傅里葉變換subplot(2,2,2);plot(f,fftshift(abs(Y)));title('添加噪聲后信號(hào)頻譜圖');%設(shè)計(jì)濾波器%Wp=2*pi*100;%帯通頻率（載波頻率）%Ws=2*pi*300;%帶阻頻率%Rp=3;%帯通頻率衰減%As=15;%帶阻頻率衰減%[N,wc]=buttord(Wp,Ws,Rp,As,'s');%計(jì)算巴特沃斯低通濾波器階數(shù)N，3db截止頻率Wc%[B,A]=butter(N,wc,'s');%計(jì)算巴特沃斯模擬低通濾波器系統(tǒng)函數(shù)系數(shù)%tam=filtfilt(B,A,y);%信號(hào)經(jīng)過濾波器%TAM=fft(tam);%濾波后信號(hào)傅里葉變換%設(shè)計(jì)巴特沃斯低通濾波器，接收端輸入前濾波a=[35,65];%通帶范圍b=[30,70];%阻帶范圍Wp=a/(fs/2);%帶通頻率Ws=b/(fs/2);%帶阻頻率Rp=5;%帯通頻率衰減Rs=15;%帶阻頻率衰減[N,w3]=buttord(Wp,Ws,Rp,Rs);%計(jì)算巴特沃斯低通濾波器階數(shù)N，3db截止頻率w3[B,A]=butter(N,w3,'s');%計(jì)算巴特沃斯模擬低通濾波器系統(tǒng)函數(shù)系數(shù),s表示模擬濾波器tam=filtfilt(B,A,y);%信號(hào)經(jīng)過濾波器subplot(2,2,3);plot(t(19801:20200),tam(19801:20200));title('帶通濾波后信號(hào)波形圖');TAM=fft(tam);%濾波后信號(hào)傅里葉變換subplot(2,2,4);plot(f,fftshift(abs(TAM)));title('帶通濾波后信號(hào)頻譜圖');%相干解調(diào)figure(3);jtam=tam.*t1;%相干解調(diào)（調(diào)制已接收信號(hào)tam乘以載波信號(hào)t1）%設(shè)計(jì)巴特沃斯低通濾波器，對(duì)相干解調(diào)后的信號(hào)濾波Wt=15/(fs/2);%帶通頻率Wz=40/(fs/2);%帶通頻率Rt=3;%帯通頻率衰減Rz=20;%帶阻頻率衰減[N1,W3db]=buttord(Wt,Wz,Rt,Rz);%計(jì)算巴特沃斯低通濾波器階數(shù)N，3db截止頻率W3db[B,A]=butter(N1,W3db,'low');%計(jì)算巴特沃斯模擬低通濾波器系統(tǒng)函數(shù)系數(shù)outam=filtfilt(B,A,jtam);%相干相乘后信號(hào)jtam經(jīng)過巴特沃斯低通濾波器，輸出為outam。filtfilt為0相移濾波器subplot(2,1,1);plot(t(19801:20200),outam(19801:20200));title('解調(diào)信號(hào)時(shí)域圖');M0OUTAM=fft(outam);%解調(diào)輸出結(jié)果的傅里葉變換subplot(2,1,2);plot(f,fftshift(abs(M0OUTAM)));title('解調(diào)信號(hào)頻譜圖');%設(shè)置信源m(t)=2*cos(pi*t)，載波c(t)=Ac*cos(20*pi*t),做已調(diào)信號(hào)的波形和頻譜%設(shè)定Ac=2，自定義直流分量A0，分別做出正常調(diào)幅(Am<A0)、臨界調(diào)幅(Am=A0)、過調(diào)幅(Am>A0)%Am=2，AM(t)=(A0+mt)*ct，Am為調(diào)制系數(shù)%正常調(diào)幅——A0=3clearall;closeall;Fs=100;%采樣頻率，要大于時(shí)域信號(hào)中最大頻率(是多少？)的兩倍Ts=1/Fs;%采樣周期，或者說是時(shí)域信號(hào)的時(shí)間步長，因?yàn)槟闾幚淼男盘?hào)不會(huì)是連續(xù)的L=5000;%信號(hào)長度，或者說你選取的信號(hào)離散點(diǎn)的個(gè)數(shù)，就是總共L個(gè)時(shí)間點(diǎn)由此知，頻率分辨率為（T*L）hzT=L*Ts;df1=Fs/L;t=(0:L-1)*Ts;%時(shí)間向量，可以看成時(shí)域信號(hào)的橫坐標(biāo)，類似于1dt,2dt,3dt，直到（L-1）t%這意味著，你的Fs越大，或者T越小，得到的結(jié)果越準(zhǔn)確。這一段是自己假設(shè)一個(gè)輸入時(shí)域信號(hào)f0=1;%頻率值%直流分量A0=3;Ac=2;%基帶信號(hào)——調(diào)制信號(hào)m_t=2*cos(pi*f0*t);%載波c_t=Ac*cos(20*pi*t);%調(diào)制y=(A0+m_t).*c_tS=y;%原函數(shù)Y=fft(S);%調(diào)用fft(x)，不用fft(x,n)，據(jù)說前者更準(zhǔn)確P2=abs(Y/L);%每個(gè)量除以數(shù)列長度L%**********(為什么要除以L？)********P1=P2(1:L/2+1);%取交流部分，這個(gè)L/2+1是對(duì)的，已經(jīng)驗(yàn)證過了P1(2:end-1)=2*P1(2:end-1);%交流部分模值乘以2,(雙邊轉(zhuǎn)單邊)f=Fs*(0:(L/2))/L;%計(jì)算頻域中傅里葉信號(hào)實(shí)際的頻率的橫坐標(biāo)，類似于1df,2df,3df,直到Fs/2figure(2)%建立幕布1subplot(2,1,1);%分成三行一列的區(qū)域，占第一個(gè)區(qū)域plot(t,S)%原函數(shù)作圖，分別為x，y軸xlabel('t(s)')%x軸標(biāo)簽ylabel('AM_1');axis([020-1212]);title('正常調(diào)幅(Am<A0)');subplot(2,1,2);%分成三行一列的區(qū)域，占第二個(gè)區(qū)域plot(f,P1)%傅里葉作圖這個(gè)振幅才是真實(shí)的振幅，對(duì)Y已經(jīng)乘以2/L了，所以振幅是真實(shí)的title('正常調(diào)幅頻譜')%標(biāo)題,單邊振幅圖xlabel('f(Hz)')ylabel('|P1(f)|')axis([01508]);%坐標(biāo)軸的范圍，[xmin,xmax,ymin,ymax]%臨界調(diào)幅——A0=2Fs=100;%采樣頻率，要大于時(shí)域信號(hào)中最大頻率(是多少？)的兩倍Ts=1/Fs;%采樣周期，或者說是時(shí)域信號(hào)的時(shí)間步長，因?yàn)槟闾幚淼男盘?hào)不會(huì)是連續(xù)的L=5000;%信號(hào)長度，或者說你選取的信號(hào)離散點(diǎn)的個(gè)數(shù)，就是總共L個(gè)時(shí)間點(diǎn)由此知，頻率分辨率為（T*L）hzT=L*Ts;df1=Fs/L;t=(0:L-1)*Ts;%時(shí)間向量，可以看成時(shí)域信號(hào)的橫坐標(biāo)，類似于1dt,2dt,3dt，直到（L-1）t%這意味著，你的Fs越大，或者T越小，得到的結(jié)果越準(zhǔn)確。這一段是自己假設(shè)一個(gè)輸入時(shí)域信號(hào)f0=1;%頻率值A(chǔ)0=2;Ac=2;m_t=2*cos(pi*f0*t);c_t=Ac*cos(20*pi*t);y=(A0+m_t).*c_tS=y;%原函數(shù)Y=fft(S);%調(diào)用fft(x)，不用fft(x,n)，據(jù)說前者更準(zhǔn)確P2=abs(Y/L);%每個(gè)量除以數(shù)列長度L%**********(為什么要除以L？)********P1=P2(1:L/2+1);%取交流部分，這個(gè)L/2+1是對(duì)的，已經(jīng)驗(yàn)證過了P1(2:end-1)=2*P1(2:end-1);%交流部分模值乘以2,(雙邊轉(zhuǎn)單邊)f=Fs*(0:(L/2))/L;%計(jì)算頻域中傅里葉信號(hào)實(shí)際的頻率的橫坐標(biāo)，類似于1df,2df,3df,直到Fs/2figure(3)%建立幕布1subplot(2,1,1);%分成三行一列的區(qū)域，占第一個(gè)區(qū)域plot(t,S)%原函數(shù)作圖，分別為x，y軸xlabel('t(s)')%x軸標(biāo)簽ylabel('AM_2');axis([020-1010]);title('臨界調(diào)幅(Am=A0)');subplot(2,1,2);%分成三行一列的區(qū)域，占第二個(gè)區(qū)域plot(f,P1)%傅里葉作圖這個(gè)振幅才是真實(shí)的振幅，對(duì)Y已經(jīng)乘以2/L了，所以振幅是真實(shí)的title('臨界調(diào)幅頻譜')%標(biāo)題,單邊振幅圖xlabel('f(Hz)')ylabel('|P1(f)|')axis([01505]);%坐標(biāo)軸的范圍，[xmin,xmax,ymin,ymax]%過調(diào)幅——A0=0.5Fs=100;%采樣頻率，要大于時(shí)域信號(hào)中最大頻率(是多少？)的兩倍Ts=1/Fs;