目錄1課程設(shè)計(jì)的目的、意義……………………（3）2課題任務(wù)……………………（3）3設(shè)計(jì)思路與方案……………（4）4設(shè)計(jì)內(nèi)容、步驟及要求……………………（4）5設(shè)計(jì)步驟及結(jié)果分析………………………（4）5.1必做部分………………（6）5.1.1Matlab程序及運(yùn)行結(jié)果…………（6）1.普通AM調(diào)制與解調(diào)……………………（6）2.抑制雙邊帶調(diào)制與解調(diào)………………（10)3.單邊帶調(diào)制與解調(diào)……………………(14)5.1.2Simulink仿真及運(yùn)行結(jié)果…………（16）1.普通AM調(diào)制與解調(diào)……………………(16)1.1單音普通調(diào)制解調(diào)…………………(16)1.2復(fù)音普通調(diào)制解調(diào)…………………(18)2.抑制雙邊帶調(diào)制解調(diào)…………………（20)2.1單音雙邊帶調(diào)制解調(diào)………………（20）2.2復(fù)音抑制雙邊帶調(diào)制解調(diào)…………（21）3.單邊帶調(diào)制解調(diào)………………………（22）3.1單音單邊帶調(diào)制解調(diào)………………（22）3.2復(fù)音單邊帶調(diào)制解調(diào)………………（24)5.2拓展部分………………（26）5.2.1單音普通AM調(diào)制解調(diào)……………（26）5.2.2單音抑制雙邊帶調(diào)制解調(diào)…………（27）5.2.3單音單邊帶調(diào)制解調(diào)………………（27）6總結(jié)……………（29）7參考文獻(xiàn)………………………（30）8意見(jiàn)、建議……………………（31）摘要：本課程設(shè)計(jì)主要利用MATLAB集成環(huán)境下的Simulink仿真平臺(tái)及Labview虛擬儀器仿真研究AM調(diào)制與解調(diào)模擬系統(tǒng)的理論設(shè)計(jì)和軟件仿真方法。從而實(shí)現(xiàn)單音調(diào)制的普通調(diào)幅方式（AM）、抑制載波的雙邊帶調(diào)制（DSB-SC）和單邊帶調(diào)制（SSB）的系統(tǒng)設(shè)計(jì)及仿真，并顯示仿真結(jié)果，根據(jù)仿真顯示結(jié)果分析所設(shè)計(jì)的系統(tǒng)性能。在課程設(shè)計(jì)中，幅度調(diào)制是用調(diào)制信號(hào)去控制高頻載波的振幅，使其按調(diào)制信號(hào)的規(guī)律變化，其它參數(shù)不變。同時(shí)也是使高頻載波的振幅載有傳輸信息的調(diào)制方式。關(guān)鍵詞：Simulink,GUI友好界面，調(diào)制與解調(diào)，Labview1、本課題的目的與意義1.1目的：本課程設(shè)計(jì)課題主要研究AM調(diào)制與解調(diào)模擬系統(tǒng)的理論設(shè)計(jì)和軟件仿真方法。通過(guò)完成本課題的設(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)行編程仿真的能力；4．熟悉基于Simulink的動(dòng)態(tài)建模和仿真的步驟和過(guò)程；5．了解基于LabVIEW虛擬儀器的特點(diǎn)和使用方法，熟悉采用LabVIEW進(jìn)行仿真的方法。1.2意義：通過(guò)本次課程設(shè)計(jì)使我們了解了幅度調(diào)制與解調(diào)的基本原理。在進(jìn)行了專業(yè)基礎(chǔ)知識(shí)課程教學(xué)的基礎(chǔ)上，設(shè)計(jì)分析一些簡(jiǎn)單的仿真系統(tǒng)，有助于加深對(duì)所學(xué)知識(shí)的鞏固和理解。課題任務(wù)設(shè)計(jì)AM調(diào)制與解調(diào)模擬系統(tǒng)，仿真實(shí)現(xiàn)相關(guān)功能。包括:可實(shí)現(xiàn)單音調(diào)制的普通調(diào)幅方式（AM）、抑制載波的雙邊帶調(diào)制（DSB-SC）和單邊帶調(diào)制（SSB）的系統(tǒng)設(shè)計(jì)及仿真，要求給出系統(tǒng)的設(shè)計(jì)框圖、源程序代碼及仿真結(jié)果，并要求給出程序的具體解釋說(shuō)明，記錄系統(tǒng)的各個(gè)輸出點(diǎn)的波形和頻譜圖。具體內(nèi)容為：1）設(shè)計(jì)實(shí)現(xiàn)AM(包括普通AM、DSB-SC和SSB)調(diào)制與解調(diào)的模擬系統(tǒng)，給出系統(tǒng)的原理框圖,對(duì)系統(tǒng)的主要參數(shù)進(jìn)行設(shè)計(jì)說(shuō)明。2）采用Matlab語(yǔ)言設(shè)計(jì)相關(guān)程序,實(shí)現(xiàn)1）中所設(shè)計(jì)模擬系統(tǒng)的功能，要求采用兩種方式進(jìn)行仿真，即直接采用Matlab語(yǔ)言編程的靜態(tài)仿真方式、采用Simulink進(jìn)行動(dòng)態(tài)建模和仿真的方式。要求采用兩種以上調(diào)制信號(hào)源（如正弦波、三角波和方波）進(jìn)行仿真，并記錄系統(tǒng)的各個(gè)輸出點(diǎn)的波形和頻譜圖。3)設(shè)計(jì)圖形用戶界面。采用Matlab語(yǔ)言，利用GUI設(shè)計(jì)友好的圖形用戶界面，完成AM調(diào)制與解調(diào)的功能。4)采用LabVIEW進(jìn)行仿真設(shè)計(jì),實(shí)現(xiàn)系統(tǒng)的功能,要求給出系統(tǒng)的前面板和框圖，采用兩種以上調(diào)制信號(hào)源（如正弦波、三角波和方波）進(jìn)行仿真，并記錄仿真結(jié)果。5)要求分析上述三種實(shí)現(xiàn)方式（直接采用Matlab語(yǔ)言編程的靜態(tài)仿真方式、采用Simulink進(jìn)行動(dòng)態(tài)建模和仿真的方式和采用LabVIEW進(jìn)行仿真設(shè)計(jì)）進(jìn)行對(duì)比分析，并與理論設(shè)計(jì)結(jié)果進(jìn)行比較分析。6)對(duì)系統(tǒng)功能進(jìn)行綜合測(cè)試，整理數(shù)據(jù)，撰寫設(shè)計(jì)報(bào)告。3、設(shè)計(jì)方案及論證：模擬調(diào)制方式是載頻信號(hào)的幅度、頻率或相位隨著欲傳輸?shù)哪M輸入基帶信號(hào)的變化而相應(yīng)發(fā)生變化的調(diào)制方式，包括：幅度調(diào)制（AM）、頻率調(diào)制（FM）、相位調(diào)制（PM）三種。幅度調(diào)制是用調(diào)制信號(hào)去控制高頻載波的振幅，使其按調(diào)制信號(hào)的規(guī)律變化，其它參數(shù)不變。是使高頻載波的振幅載有傳輸信息的調(diào)制方式。振幅調(diào)制分為三種方式：普通調(diào)幅方式（AM）、抑制載波的雙邊帶調(diào)制（DSB-SC）和單邊帶調(diào)制（SSB）。所得的已調(diào)信號(hào)分別稱為調(diào)幅波信號(hào)、雙邊帶信號(hào)和單邊帶信號(hào)。先在Matlab環(huán)境下寫出三種振幅調(diào)制方式的程序，運(yùn)行得出結(jié)果；再在Simulink環(huán)境下做出仿真圖，再次得到結(jié)果；最后在Labview集成環(huán)境下進(jìn)行仿真設(shè)計(jì)，賦值運(yùn)行，得出結(jié)果。對(duì)所得的三次結(jié)果進(jìn)行對(duì)比分析。設(shè)計(jì)內(nèi)容、步驟及要求4.1必選部分(1)設(shè)計(jì)實(shí)現(xiàn)AM調(diào)制與解調(diào)的模擬系統(tǒng)，給出系統(tǒng)的原理框圖,對(duì)系統(tǒng)的主要參數(shù)進(jìn)行設(shè)計(jì)說(shuō)明,具體參數(shù)包括：載波頻率、調(diào)制信號(hào)頻率、載波大小、調(diào)制信號(hào)大小、濾波器參數(shù)等。并對(duì)所設(shè)計(jì)的系統(tǒng)進(jìn)行理論分析計(jì)算。(2)根據(jù)所設(shè)計(jì)的AM調(diào)制與解調(diào)的模擬系統(tǒng)，進(jìn)行基于Matlab語(yǔ)言的靜態(tài)仿真設(shè)計(jì)。分別實(shí)現(xiàn)單音調(diào)制的普通調(diào)幅方式（AM）、抑制載波的雙邊帶調(diào)制（DSB-SC）和單邊帶調(diào)制（SSB）的系統(tǒng)仿真設(shè)計(jì)，要求給出系統(tǒng)的Matlab編程仿真程序及結(jié)果，并要求寫出程序的具體解釋說(shuō)明，記錄系統(tǒng)的各個(gè)輸出點(diǎn)的波形和頻譜圖。要求調(diào)制信號(hào)分別采用不同類型的信號(hào)進(jìn)行仿真，至少給出兩種以上調(diào)制信號(hào)源，如：?jiǎn)我粜盘?hào)、合成復(fù)雜音信號(hào)、直接錄制的模擬語(yǔ)音信號(hào)。(3)根據(jù)所設(shè)計(jì)的AM調(diào)制與解調(diào)的模擬系統(tǒng)，采用Simulink進(jìn)行動(dòng)態(tài)建模仿真設(shè)計(jì)。分別實(shí)現(xiàn)普通調(diào)幅方式（AM）、抑制載波的雙邊帶調(diào)制（DSB-SC）和單邊帶調(diào)制（SSB）的系統(tǒng)動(dòng)態(tài)仿真設(shè)計(jì)，要求包括調(diào)制和解調(diào)的部分，并給出采用Simulink進(jìn)行動(dòng)態(tài)建模仿真的系統(tǒng)方框圖，同時(shí)記錄系統(tǒng)的各個(gè)輸出點(diǎn)的波形和頻譜圖。要求采用兩種以上調(diào)制信號(hào)源進(jìn)行仿真，具體參數(shù)自定。載波信號(hào)頻率根據(jù)設(shè)計(jì)情況設(shè)定。(4)根據(jù)所設(shè)計(jì)的AM調(diào)制與解調(diào)的模擬系統(tǒng)，采用Matlab語(yǔ)言，利用GUI設(shè)計(jì)友好的圖形用戶界面。完成4.1必選部分（2）所要求的功能。(5)根據(jù)仿真結(jié)果,對(duì)系統(tǒng)的時(shí)域、頻域特性進(jìn)行分析，并與理論設(shè)計(jì)結(jié)果進(jìn)行比較分析。4.2選作部分(1)根據(jù)所設(shè)計(jì)的AM調(diào)制與解調(diào)的模擬系統(tǒng)，說(shuō)明具體的參數(shù)，進(jìn)行基于LabVIEW環(huán)境的仿真，分別實(shí)現(xiàn)普通調(diào)幅方式（AM）、抑制載波的雙邊帶調(diào)制（DSB-SC）和單邊帶調(diào)制（SSB）的系統(tǒng)仿真設(shè)計(jì)，要求包括調(diào)制和解調(diào)的部分，給出系統(tǒng)的前面板和框圖，并記錄仿真結(jié)果。(2)要求調(diào)制信號(hào)采用不同類型的信號(hào)源，進(jìn)行進(jìn)一步的仿真，給出系統(tǒng)的前面板和框圖，并記錄仿真結(jié)果，觀察分析頻譜的變化情況。5、設(shè)計(jì)步驟及結(jié)果分析5.1必做部分：（1）有關(guān)參數(shù)說(shuō)明：Wp(或wp)為通帶截止頻率，Ws(或ws)為阻帶截止頻率，Rp為通帶衰減，As為阻帶衰減（2）有關(guān)函數(shù)說(shuō)明：%butterworth低通濾波器原型設(shè)計(jì)函數(shù)，要求Ws>Wp>0，As>Rp>0function[b,a]=afd_butt(Wp,Ws,Rp,As)N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws)));%上條語(yǔ)句為求濾波器階數(shù)：N為整數(shù);%ceil朝正無(wú)窮大方向取整;fprintf('\nButterworthFilterOrder=%2.0f\n',N)OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N)))%求對(duì)應(yīng)于N的3db截止頻率;[b,a]=u_buttap(N,OmegaC);%傅里葉變換函數(shù)function[Xk]=dft(xn,N)n=[0:1:N-1];k=[0:1:N-1];WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^(nk);Xk=xn*WNnk;5.1.1Matlab程序代碼及程序運(yùn)行結(jié)果（一）普通AM調(diào)制與解調(diào)：%單音普通調(diào)幅波的調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ma),要求fs>2fc;%x為調(diào)制信號(hào),t為調(diào)制信號(hào)自變量,2t0為采樣區(qū)間，fs為采樣頻率;%fc為載波頻率，Vm0為未經(jīng)調(diào)制的輸出載波電壓振幅，ma調(diào)幅度;t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;x=4*cos(150*pi*n);%調(diào)制信號(hào)y2=Vm0*cos(2*pi*fc*n);%載波信號(hào)figure(1)subplot(2,1,1);plot(n,y2);axis([-0.01,0.01,-5,5]);title('載波信號(hào)');N=length(x);Y2=fft(y2);subplot(2,1,2);plot(n,Y2);title('載波信號(hào)頻譜');%畫出頻譜波形y=Vm0*(1+ma*x/Vm0).*cos(2*pi*fc*n);figure(2)subplot(2,1,1);plot(n,x);title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');X=fft(x);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(3)subplot(2,1,1);plot(w,abs(X))axis([0,pi/4,0,2000]);title('調(diào)制信號(hào)頻譜');subplot(2,1,2);plot(w,abs(Y));axis([pi/6,pi/4,0,1200]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;Rp=1;As=15;ws=60/N*pi;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(4)subplot(2,1,1);plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2);plot(w,abs(Y))axis([0,pi/6,0,1000]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171%復(fù)雜音普通調(diào)幅波調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ma)，要求fs>2fc%x為調(diào)制信號(hào)，t為調(diào)制信號(hào)自變量,2t0為采樣區(qū)間，fs為采樣頻率，%fc為載波頻率，Vm0為未經(jīng)調(diào)制的輸出載波電壓振幅，ma為調(diào)幅度;t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;x=4*cos(150*pi*n)-3*cos(200*pi*n);%調(diào)制信號(hào)n=-t0/2:1/fs:t0/2;y=Vm0*(1+ma*x/Vm0).*cos(2*pi*fc*n);figure(1)subplot(2,1,1)plot(n,x)title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');N=length(x);X=fft(x);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(2)subplot(2,1,1)plot(w,abs(X))axis([0,pi/4,0,2000]);title('調(diào)制信號(hào)頻譜');%畫出頻譜波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,1200]);title('已調(diào)波信號(hào)頻譜');y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(3)subplot(2,1,1)plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,1000]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171%三角波普通調(diào)幅波調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ma)，要求fs>2fc%x為調(diào)制信號(hào)，t為調(diào)制信號(hào)自變量,2t0為采樣區(qū)間，fs為采樣頻率，%fc為載波頻率，Vm0為未經(jīng)調(diào)制的輸出載波電壓振幅，ma為調(diào)幅度;t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;x=sawtooth(150*pi*n,0.5)%調(diào)制信號(hào)y2=Vm0*cos(2*pi*fc*n);%載波信號(hào)figure(1)subplot(2,1,1);plot(n,y2);axis([-0.01,0.01,-5,5]);title('載波信號(hào)');N=length(x);Y2=fft(y2);subplot(2,1,2);plot(n,Y2);title('載波信號(hào)頻譜');%畫出頻譜波形y=Vm0*(1+ma*x/Vm0).*cos(2*pi*fc*n);figure(2)subplot(2,1,1);plot(n,x);title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');X=fft(x);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(3)subplot(2,1,1);plot(w,abs(X))axis([0,pi/4,0,2000]);title('調(diào)制信號(hào)頻譜');subplot(2,1,2);plot(w,abs(Y));axis([pi/6,pi/4,0,1200]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;Rp=1;As=15;ws=60/N*pi;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(4)subplot(2,1,1);plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2);plot(w,abs(Y))axis([0,pi/6,0,1000]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171(二)抑制載波的雙邊帶調(diào)制%單音抑制載波雙邊帶調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ma)，要求fs>2fc,%x為調(diào)制信號(hào)，2t0為采樣區(qū)間，fs為采樣頻率，%fc為載波頻率，Vm0為未經(jīng)調(diào)制的輸出載波電壓振幅，ma為調(diào)幅度;t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;x=4*cos(150*pi*n);%調(diào)制信號(hào)y=Vm0*x.*cos(2*pi*fc*n);%載波信號(hào)figure(1)subplot(2,1,1)plot(n,x)title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');N=length(x);X=fft(x);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(2)subplot(2,1,1)plot(w,abs(X))axis([0,pi/4,0,2000]);title('調(diào)制信號(hào)頻譜');%畫出頻譜波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2200]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(2000*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(3)subplot(2,1,1)plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,5000]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171%復(fù)音抑制載波雙邊帶調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ma)要求fs>2fc%x調(diào)制信號(hào)，t調(diào)制信號(hào)自變量,t0采樣區(qū)間，fs采樣頻率，%fc載波頻率，Vm0輸出載波電壓振幅，ma調(diào)幅度t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;x=4*cos(150*pi*n)-3*cos(200*pi*n);%調(diào)制信號(hào)y=Vm0*x.*cos(2*pi*fc*n);%載波信號(hào)figure(1)subplot(2,1,1)plot(n,x)title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');N=length(x);X=fft(x);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(2)subplot(2,1,1)plot(w,abs(X))axis([0,pi/4,0,3000]);title('調(diào)制信號(hào)頻譜');%畫出頻譜波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2200]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(2000*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(3)subplot(2,1,1)plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,5500]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171%三角波抑制載波雙邊帶調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ma)要求fs>2fc%x調(diào)制信號(hào)，t調(diào)制信號(hào)自變量,t0采樣區(qū)間，fs采樣頻率，%fc載波頻率，Vm0輸出載波電壓振幅，ma調(diào)幅度t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;x=sawtooth(150*pi*n,0.5)%調(diào)制信號(hào)y=Vm0*x.*cos(2*pi*fc*n);%載波信號(hào)figure(1)subplot(2,1,1)plot(n,x)title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');N=length(x);X=fft(x);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(2)subplot(2,1,1)plot(w,abs(X))axis([0,pi/4,0,2000]);title('調(diào)制信號(hào)頻譜');%畫出頻譜波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2200]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(2000*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(3)subplot(2,1,1)plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,5000]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171（三）單邊帶調(diào)制與解調(diào)%單音單邊帶調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ma),要求fs>2fc%x為調(diào)制信號(hào)，2t0為采樣區(qū)間，fs為采樣頻率，%fc為載波頻率，Vm0為未經(jīng)調(diào)制的輸出載波電壓振幅，ma為調(diào)幅度;t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;N=length(n);x1=4*cos(150*pi*n);%調(diào)制信號(hào)x2=hilbert(x1,N);y=(Vm0*x1.*cos(2*pi*fc*n)-Vm0*x2.*sin(2*pi*fc*n))/2;figure(1)subplot(2,1,1)plot(n,x1)title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');X=fft(x1);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(2)subplot(2,1,1)plot(w,abs(X))axis([0,pi/4,0,3000]);title('調(diào)制信號(hào)頻譜');%畫出頻譜波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2500]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(3)subplot(2,1,1)plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,2500]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171%復(fù)雜音單邊帶調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ka),要求fs>2fc%x為調(diào)制信號(hào)，2t0為采樣區(qū)間，fs為采樣頻率，%fc為載波頻率，Vm0為未經(jīng)調(diào)制的輸出載波電壓振幅，ma為調(diào)幅度;t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;N=length(n);x1=4*cos(150*pi*n)-3*cos(200*pi*n);%調(diào)制信號(hào)x2=hilbert(x1,N);y=(Vm0*x1.*cos(2*pi*fc*n)-Vm0*x2.*sin(2*pi*fc*n))/2;figure(1)subplot(2,1,1)plot(n,x1)title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');X=fft(x1);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(2)subplot(2,1,1)plot(w,abs(X))axis([0,pi/4,0,3000]);title('調(diào)制信號(hào)頻譜');%畫出頻譜波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2500]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(3)subplot(2,1,1)plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,2500]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171%三角波單邊帶調(diào)制y=amod(x,t,fs,t0,fc,Vm0,ka),要求fs>2fc%x為調(diào)制信號(hào)，2t0為采樣區(qū)間，fs為采樣頻率，%fc為載波頻率，Vm0為未經(jīng)調(diào)制的輸出載波電壓振幅，ma為調(diào)幅度;t0=0.1;fs=12000;fc=1000;Vm0=2.5;ma=0.25;n=-t0/2:1/fs:t0/2;N=length(n);x1=sawtooth(150*pi*n,0.5)%調(diào)制信號(hào)x2=hilbert(x1,N);y=(Vm0*x1.*cos(2*pi*fc*n)-Vm0*x2.*sin(2*pi*fc*n))/2;figure(1)subplot(2,1,1)plot(n,x1)title('調(diào)制信號(hào)');subplot(2,1,2)plot(n,y)title('已調(diào)波信號(hào)');X=fft(x1);Y=fft(y);w=0:2*pi/(N-1):2*pi;figure(2)subplot(2,1,1)plot(w,abs(X))axis([0,pi/4,0,3000]);title('調(diào)制信號(hào)頻譜');%畫出頻譜波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2500]);title('已調(diào)波信號(hào)頻譜');%畫出頻譜波形y1=y-2*cos(100*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)頻譜處wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T);y=filter(b,a,y2);figure(3)subplot(2,1,1)plot(n,y)title('解調(diào)波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,2500]);title('解調(diào)信號(hào)頻譜');%畫出頻譜波形程序運(yùn)行結(jié)果：ButterworthFilterOrder=6OmegaC=0.1171基于C8051F單片機(jī)直流電動(dòng)機(jī)反饋控制系統(tǒng)的設(shè)計(jì)與研究基于單片機(jī)的嵌入式Web服務(wù)器的研究MOTOROLA單片機(jī)MC68HC（8）05PV8/A內(nèi)嵌EEPROM的工藝和制程方法及對(duì)良率的影響研究基于模糊控制的電阻釬焊單片機(jī)溫度控制系統(tǒng)的研制基于MCS-51系列單片機(jī)的通用控制模塊的研究基于單片機(jī)實(shí)現(xiàn)的供暖系統(tǒng)最佳啟停自校正（STR）調(diào)節(jié)器單片機(jī)控制的二級(jí)倒立擺系統(tǒng)的研究基于增強(qiáng)型51系列單片機(jī)的TCP/IP協(xié)議棧的實(shí)現(xiàn)基于單片機(jī)的蓄電池自動(dòng)監(jiān)測(cè)系統(tǒng)基于32位嵌入式單片機(jī)系統(tǒng)的圖像采集與處理技術(shù)的研究基于單片機(jī)的作物營(yíng)養(yǎng)診斷專家系統(tǒng)的研究基于單片機(jī)的交流伺服電機(jī)運(yùn)動(dòng)控制系統(tǒng)研究與開(kāi)發(fā)基于單片機(jī)的泵管內(nèi)壁硬度測(cè)試儀的研制基于單片機(jī)的自動(dòng)找平控制系統(tǒng)研究基于C8051F040單片機(jī)的嵌入式系統(tǒng)開(kāi)發(fā)基于單片機(jī)的液壓動(dòng)力系統(tǒng)狀態(tài)監(jiān)測(cè)儀開(kāi)發(fā)模糊Smith智能控制方法的研究及其單片機(jī)實(shí)現(xiàn)一種基于單片機(jī)的軸快流CO〈,2〉激光器的手持控制面板的研制基于雙單片機(jī)沖床數(shù)控系統(tǒng)的研究基于CYGNAL單片機(jī)的在線間歇式濁度儀的研制基于單片機(jī)的噴油泵試驗(yàn)臺(tái)控制器的研制基于單片機(jī)的軟起動(dòng)器的研究和設(shè)計(jì)基于單片機(jī)控制的高速快走絲電火花線切割機(jī)床短循環(huán)走絲方式研究基于單片機(jī)的機(jī)電產(chǎn)品控制系統(tǒng)開(kāi)發(fā)基于PIC單片機(jī)的智能手機(jī)充電器基于單片機(jī)的實(shí)時(shí)內(nèi)核設(shè)計(jì)及其應(yīng)用研究基于單片機(jī)的遠(yuǎn)程抄表系統(tǒng)的設(shè)計(jì)與研究基于單片機(jī)的煙氣二氧化硫濃度檢測(cè)儀的研制基于微型光譜儀的單片機(jī)系統(tǒng)單片機(jī)系統(tǒng)軟件構(gòu)件開(kāi)發(fā)的技術(shù)研究基于單片機(jī)的液體點(diǎn)滴速度自動(dòng)檢測(cè)儀的研制基于單片機(jī)系統(tǒng)的多功能溫度測(cè)量?jī)x的研制基于PIC單片機(jī)的電能采集終端的設(shè)計(jì)和應(yīng)用基于單片機(jī)的光纖光柵解調(diào)儀的研制氣壓式線性摩擦焊機(jī)單片機(jī)控制系統(tǒng)的研制基于單片機(jī)的數(shù)字磁通門傳感器基于單片機(jī)的旋轉(zhuǎn)變壓器-數(shù)字轉(zhuǎn)換器的研究基于單片機(jī)的光纖Bragg光柵解調(diào)系統(tǒng)的研究單片機(jī)控制的便攜式多功能乳腺治療儀的研制基于C8051F020單片機(jī)的多生理信號(hào)檢測(cè)儀基于單片機(jī)的電機(jī)運(yùn)動(dòng)控制系統(tǒng)設(shè)計(jì)Pico專用單片機(jī)核的可測(cè)性設(shè)計(jì)研究基于MCS-51單片機(jī)的熱量計(jì)基于雙單片機(jī)的智能遙測(cè)微型氣象站MCS-51單片機(jī)構(gòu)建機(jī)器人的實(shí)踐研究基于單片機(jī)的輪軌力檢測(cè)基于單片機(jī)的GPS定位儀的研究與實(shí)現(xiàn)基于單片機(jī)的電液伺服控制系統(tǒng)用于單片機(jī)系統(tǒng)的MMC卡文件系統(tǒng)研制基于單片機(jī)的時(shí)控和計(jì)數(shù)系統(tǒng)性能優(yōu)化的研究基于單片機(jī)和CPLD的粗光柵位移測(cè)量系統(tǒng)研究單片機(jī)控制的后備式方波UPS提升高職學(xué)生單片機(jī)應(yīng)用能力的探究基于單片機(jī)控制的自動(dòng)低頻減載裝置研究基于單片機(jī)控制的水下焊接電源的研究基于單片機(jī)的多通道數(shù)據(jù)采集系統(tǒng)基于uPSD3234單片機(jī)的氚表面污染測(cè)量?jī)x的研制基于單片機(jī)的紅外測(cè)油儀的研究96系列單片機(jī)仿真器研究與設(shè)計(jì)基于單片機(jī)的單晶金剛石刀具刃磨設(shè)備的數(shù)控改造基于單片機(jī)的溫度智能控制系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)基于MSP430單片機(jī)的電梯門機(jī)控制器的研制基于單片機(jī)的氣體測(cè)漏儀的研究基于三菱M16C/6N系列單片機(jī)的CAN/USB協(xié)議轉(zhuǎn)換器基于單片機(jī)和DSP的變壓器油色譜在線監(jiān)測(cè)技術(shù)研究基于單片機(jī)的膛壁溫度報(bào)警系統(tǒng)設(shè)計(jì)基于AVR單片機(jī)的低壓無(wú)功補(bǔ)償控制器的設(shè)計(jì)基于單片機(jī)船舶電力推進(jìn)電機(jī)監(jiān)測(cè)系統(tǒng)基于單片機(jī)網(wǎng)絡(luò)的振動(dòng)信號(hào)的采集系統(tǒng)基于單片機(jī)的大容量數(shù)據(jù)存儲(chǔ)技術(shù)的應(yīng)用研究基于單片機(jī)的疊圖機(jī)研究與教學(xué)方法實(shí)踐基于單片機(jī)嵌入式Web服務(wù)器技術(shù)的研究及實(shí)現(xiàn)基于AT89S52單片機(jī)的通用數(shù)據(jù)采集系統(tǒng)基于單片機(jī)的多道脈沖幅度分析儀研究機(jī)器人旋轉(zhuǎn)電弧傳感角焊縫跟蹤單片機(jī)控制系統(tǒng)基于單片機(jī)的控制系統(tǒng)在PLC虛擬教學(xué)實(shí)驗(yàn)中的應(yīng)用研究基于單片機(jī)系統(tǒng)的網(wǎng)絡(luò)通信研究與應(yīng)用基于PIC16F877單片機(jī)的莫爾斯碼自動(dòng)譯碼系統(tǒng)設(shè)計(jì)與研究基于單片機(jī)的模糊控制器在工業(yè)電阻爐上的應(yīng)用研究基于雙單片機(jī)沖床數(shù)控系統(tǒng)的研究與開(kāi)發(fā)基于Cygnal單片機(jī)的μC/OS-Ⅱ的研究基于單片機(jī)的一體化智能差示掃描量熱儀系統(tǒng)研究基于TCP/IP協(xié)議的單片機(jī)與Internet互聯(lián)的研究與實(shí)現(xiàn)變頻調(diào)速液壓電梯單片機(jī)控制器的研究基于單片機(jī)γ-免疫計(jì)數(shù)器自動(dòng)換樣功能的研究與實(shí)現(xiàn)基于單片機(jī)的倒立擺控制系統(tǒng)設(shè)計(jì)與實(shí)