AM调制与解调课程设计

目录1课程设计的目的、意义……………………（3）2课题任务……………………（3）3设计思路与方案……………（4）4设计内容、步骤及要求……………………（4）5设计步骤及结果分析………………………（4）5.1必做部分………………（6）5.1.1Matlab程序及运行结果…………（6）1.普通AM调制与解调……………………（6）2.抑制双边带调制与解调………………（10)3.单边带调制与解调……………………(14)5.1.2Simulink仿真及运行结果…………（16）1.普通AM调制与解调……………………(16)1.1单音普通调制解调…………………(16)1.2复音普通调制解调…………………(18)2.抑制双边带调制解调…………………（20)2.1单音双边带调制解调………………（20）2.2复音抑制双边带调制解调…………（21）3.单边带调制解调………………………（22）3.1单音单边带调制解调………………（22）3.2复音单边带调制解调………………（24)5.2拓展部分………………（26）5.2.1单音普通AM调制解调……………（26）5.2.2单音抑制双边带调制解调…………（27）5.2.3单音单边带调制解调………………（27）6总结……………（29）7参考文献………………………（30）8意见、建议……………………（31）摘要：本课程设计主要利用MATLAB集成环境下的Simulink仿真平台及Labview虚拟仪器仿真研究AM调制与解调模拟系统的理论设计和软件仿真方法。从而实现单音调制的普通调幅方式（AM）、抑制载波的双边带调制（DSB-SC）和单边带调制（SSB）的系统设计及仿真，并显示仿真结果，根据仿真显示结果分析所设计的系统性能。在课程设计中，幅度调制是用调制信号去控制高频载波的振幅，使其按调制信号的规律变化，其它参数不变。同时也是使高频载波的振幅载有传输信息的调制方式。关键词：Simulink,GUI友好界面，调制与解调，Labview1、本课题的目的与意义1.1目的：本课程设计课题主要研究AM调制与解调模拟系统的理论设计和软件仿真方法。通过完成本课题的设计，拟主要达到以下几个目的：1．掌握模拟系统AM调制与解调的原理。2．掌握AM调制与解调模拟系统的理论设计方法；3．掌握应用MATLAB分析系统时域、频域特性的方法，进一步锻炼应用Matlab进行编程仿真的能力；4．熟悉基于Simulink的动态建模和仿真的步骤和过程；5．了解基于LabVIEW虚拟仪器的特点和使用方法，熟悉采用LabVIEW进行仿真的方法。1.2意义：通过本次课程设计使我们了解了幅度调制与解调的基本原理。在进行了专业基础知识课程教学的基础上，设计分析一些简单的仿真系统，有助于加深对所学知识的巩固和理解。课题任务设计AM调制与解调模拟系统，仿真实现相关功能。包括:可实现单音调制的普通调幅方式（AM）、抑制载波的双边带调制（DSB-SC）和单边带调制（SSB）的系统设计及仿真，要求给出系统的设计框图、源程序代码及仿真结果，并要求给出程序的具体解释说明，记录系统的各个输出点的波形和频谱图。具体内容为：1）设计实现AM(包括普通AM、DSB-SC和SSB)调制与解调的模拟系统，给出系统的原理框图,对系统的主要参数进行设计说明。2）采用Matlab语言设计相关程序,实现1）中所设计模拟系统的功能，要求采用两种方式进行仿真，即直接采用Matlab语言编程的静态仿真方式、采用Simulink进行动态建模和仿真的方式。要求采用两种以上调制信号源（如正弦波、三角波和方波）进行仿真，并记录系统的各个输出点的波形和频谱图。3)设计图形用户界面。采用Matlab语言，利用GUI设计友好的图形用户界面，完成AM调制与解调的功能。4)采用LabVIEW进行仿真设计,实现系统的功能,要求给出系统的前面板和框图，采用两种以上调制信号源（如正弦波、三角波和方波）进行仿真，并记录仿真结果。5)要求分析上述三种实现方式（直接采用Matlab语言编程的静态仿真方式、采用Simulink进行动态建模和仿真的方式和采用LabVIEW进行仿真设计）进行对比分析，并与理论设计结果进行比较分析。6)对系统功能进行综合测试，整理数据，撰写设计报告。3、设计方案及论证：模拟调制方式是载频信号的幅度、频率或相位随着欲传输的模拟输入基带信号的变化而相应发生变化的调制方式，包括：幅度调制（AM）、频率调制（FM）、相位调制（PM）三种。幅度调制是用调制信号去控制高频载波的振幅，使其按调制信号的规律变化，其它参数不变。是使高频载波的振幅载有传输信息的调制方式。振幅调制分为三种方式：普通调幅方式（AM）、抑制载波的双边带调制（DSB-SC）和单边带调制（SSB）。所得的已调信号分别称为调幅波信号、双边带信号和单边带信号。先在Matlab环境下写出三种振幅调制方式的程序，运行得出结果；再在Simulink环境下做出仿真图，再次得到结果；最后在Labview集成环境下进行仿真设计，赋值运行，得出结果。对所得的三次结果进行对比分析。设计内容、步骤及要求4.1必选部分(1)设计实现AM调制与解调的模拟系统，给出系统的原理框图,对系统的主要参数进行设计说明,具体参数包括：载波频率、调制信号频率、载波大小、调制信号大小、滤波器参数等。并对所设计的系统进行理论分析计算。(2)根据所设计的AM调制与解调的模拟系统，进行基于Matlab语言的静态仿真设计。分别实现单音调制的普通调幅方式（AM）、抑制载波的双边带调制（DSB-SC）和单边带调制（SSB）的系统仿真设计，要求给出系统的Matlab编程仿真程序及结果，并要求写出程序的具体解释说明，记录系统的各个输出点的波形和频谱图。要求调制信号分别采用不同类型的信号进行仿真，至少给出两种以上调制信号源，如：单音信号、合成复杂音信号、直接录制的模拟语音信号。(3)根据所设计的AM调制与解调的模拟系统，采用Simulink进行动态建模仿真设计。分别实现普通调幅方式（AM）、抑制载波的双边带调制（DSB-SC）和单边带调制（SSB）的系统动态仿真设计，要求包括调制和解调的部分，并给出采用Simulink进行动态建模仿真的系统方框图，同时记录系统的各个输出点的波形和频谱图。要求采用两种以上调制信号源进行仿真，具体参数自定。载波信号频率根据设计情况设定。(4)根据所设计的AM调制与解调的模拟系统，采用Matlab语言，利用GUI设计友好的图形用户界面。完成4.1必选部分（2）所要求的功能。(5)根据仿真结果,对系统的时域、频域特性进行分析，并与理论设计结果进行比较分析。4.2选作部分(1)根据所设计的AM调制与解调的模拟系统，说明具体的参数，进行基于LabVIEW环境的仿真，分别实现普通调幅方式（AM）、抑制载波的双边带调制（DSB-SC）和单边带调制（SSB）的系统仿真设计，要求包括调制和解调的部分，给出系统的前面板和框图，并记录仿真结果。(2)要求调制信号采用不同类型的信号源，进行进一步的仿真，给出系统的前面板和框图，并记录仿真结果，观察分析频谱的变化情况。5、设计步骤及结果分析5.1必做部分：（1）有关参数说明：Wp(或wp)为通带截止频率，Ws(或ws)为阻带截止频率，Rp为通带衰减，As为阻带衰减（2）有关函数说明：%butterworth低通滤波器原型设计函数，要求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)));%上条语句为求滤波器阶数：N为整数;%ceil朝正无穷大方向取整;fprintf('\nButterworthFilterOrder=%2.0f\n',N)OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N)))%求对应于N的3db截止频率;[b,a]=u_buttap(N,OmegaC);%傅里叶变换函数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程序代码及程序运行结果（一）普通AM调制与解调：%单音普通调幅波的调制y=amod(x,t,fs,t0,fc,Vm0,ma),要求fs>2fc;%x为调制信号,t为调制信号自变量,2t0为采样区间，fs为采样频率;%fc为载波频率，Vm0为未经调制的输出载波电压振幅，ma调幅度;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);%调制信号y2=Vm0*cos(2*pi*fc*n);%载波信号figure(1)subplot(2,1,1);plot(n,y2);axis([-0.01,0.01,-5,5]);title('载波信号');N=length(x);Y2=fft(y2);subplot(2,1,2);plot(n,Y2);title('载波信号频谱');%画出频谱波形y=Vm0*(1+ma*x/Vm0).*cos(2*pi*fc*n);figure(2)subplot(2,1,1);plot(n,x);title('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');subplot(2,1,2);plot(w,abs(Y));axis([pi/6,pi/4,0,1200]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;Rp=1;As=15;ws=60/N*pi;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2);plot(w,abs(Y))axis([0,pi/6,0,1000]);title('解调信号频谱');%画出频谱波形程序结果：ButterworthFilterOrder=6OmegaC=0.1171%复杂音普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ma)，要求fs>2fc%x为调制信号，t为调制信号自变量,2t0为采样区间，fs为采样频率，%fc为载波频率，Vm0为未经调制的输出载波电压振幅，ma为调幅度;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);%调制信号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('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');%画出频谱波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,1200]);title('已调波信号频谱');y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,1000]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171%三角波普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ma)，要求fs>2fc%x为调制信号，t为调制信号自变量,2t0为采样区间，fs为采样频率，%fc为载波频率，Vm0为未经调制的输出载波电压振幅，ma为调幅度;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)%调制信号y2=Vm0*cos(2*pi*fc*n);%载波信号figure(1)subplot(2,1,1);plot(n,y2);axis([-0.01,0.01,-5,5]);title('载波信号');N=length(x);Y2=fft(y2);subplot(2,1,2);plot(n,Y2);title('载波信号频谱');%画出频谱波形y=Vm0*(1+ma*x/Vm0).*cos(2*pi*fc*n);figure(2)subplot(2,1,1);plot(n,x);title('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');subplot(2,1,2);plot(w,abs(Y));axis([pi/6,pi/4,0,1200]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;Rp=1;As=15;ws=60/N*pi;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2);plot(w,abs(Y))axis([0,pi/6,0,1000]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171(二)抑制载波的双边带调制%单音抑制载波双边带调制y=amod(x,t,fs,t0,fc,Vm0,ma)，要求fs>2fc,%x为调制信号，2t0为采样区间，fs为采样频率，%fc为载波频率，Vm0为未经调制的输出载波电压振幅，ma为调幅度;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);%调制信号y=Vm0*x.*cos(2*pi*fc*n);%载波信号figure(1)subplot(2,1,1)plot(n,x)title('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');%画出频谱波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2200]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(2000*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,5000]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171%复音抑制载波双边带调制y=amod(x,t,fs,t0,fc,Vm0,ma)要求fs>2fc%x调制信号，t调制信号自变量,t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ma调幅度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);%调制信号y=Vm0*x.*cos(2*pi*fc*n);%载波信号figure(1)subplot(2,1,1)plot(n,x)title('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');%画出频谱波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2200]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(2000*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,5500]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171%三角波抑制载波双边带调制y=amod(x,t,fs,t0,fc,Vm0,ma)要求fs>2fc%x调制信号，t调制信号自变量,t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ma调幅度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)%调制信号y=Vm0*x.*cos(2*pi*fc*n);%载波信号figure(1)subplot(2,1,1)plot(n,x)title('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');%画出频谱波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2200]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(2000*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,5000]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171（三）单边带调制与解调%单音单边带调制y=amod(x,t,fs,t0,fc,Vm0,ma),要求fs>2fc%x为调制信号，2t0为采样区间，fs为采样频率，%fc为载波频率，Vm0为未经调制的输出载波电压振幅，ma为调幅度;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);%调制信号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('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');%画出频谱波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2500]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,2500]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171%复杂音单边带调制y=amod(x,t,fs,t0,fc,Vm0,ka),要求fs>2fc%x为调制信号，2t0为采样区间，fs为采样频率，%fc为载波频率，Vm0为未经调制的输出载波电压振幅，ma为调幅度;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);%调制信号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('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');%画出频谱波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2500]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(800*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,2500]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171%三角波单边带调制y=amod(x,t,fs,t0,fc,Vm0,ka),要求fs>2fc%x为调制信号，2t0为采样区间，fs为采样频率，%fc为载波频率，Vm0为未经调制的输出载波电压振幅，ma为调幅度;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)%调制信号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('调制信号');subplot(2,1,2)plot(n,y)title('已调波信号');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('调制信号频谱');%画出频谱波形subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,2500]);title('已调波信号频谱');%画出频谱波形y1=y-2*cos(100*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号频谱处wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计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('解调波');Y=fft(y);subplot(2,1,2)plot(w,abs(Y))axis([0,pi/6,0,2500]);title('解调信号频谱');%画出频谱波形程序运行结果：ButterworthFilterOrder=6OmegaC=0.1171基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协议的单片机与Internet互联的研究与实现变频调速液压电梯单片机控制器的研究基于单片机γ-免疫计数器自动换样功能的研究与实现基于单片机的倒立摆控制系统设计与实