版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
%普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs>2fc%x调制信号，t调制信号自变量,t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ka调幅度t0=0.1;fs=10000;fc=1000;Vm0=2;ka=0.25;n=-t0/2:1/fs:t0/2;w=0:1/(N-1):1;x=5*cos(200*pi*n);%调制信号y=Vm0*(1+ka*x/Vm0).*cos(2*pi*fc*n);h=Vm0*cos(2*pi*fc*n);figure(1)subplot(2,2,1)plot(n,x)title('调制信号');subplot(2,2,2)plot(n,y)title('已调波信号');subplot(2,2,3)plot(n,h)title('载波信号');axis([-0.01,0.01,-5,5]);H=fft(h);subplot(2,2,4)plot(w,abs(H))title('载波频谱');axis([0,1,0,500]);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,1,0,3000]);title('subplot(2,1,2)plot(w,abs(Y))axis([0,1,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,1,0,1000]);title(' 解 调 信 号 频 谱');调制信号50-5-0.05420-2
-0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05已调波信号-4-0.05
-0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05调制信号频谱30002000100000 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5已调波信号频谱100050000.55 0.6 0.65 0.7 0.75解调波43210-1-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05解调信号频谱100050000 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5%普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs>2fc%x调制信号，t调制信号自变量,t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ka调幅度t0=0.1;fs=10000;fc=1000;Vm0=2;ka=0.25;x=3*cos(200*pi*n)-2*cos(250*pi*n);n=-t0/2:1/fs:t0/2;y=Vm0*(1+ka*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,3000]);title('调制信号频谱');subplot(2,1,2)plot(w,abs(Y))axis([pi/6,pi/4,0,1200]);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,1000]);title('解调信号频谱');调制信号50-5-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05已调波信号420-2-4-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05调制信号频谱30002000100000 0.1 0.2 0.3 0.4 0.5 0.6 0.7已调波信号频谱100050000.55 0.6 0.65 0.7 0.75解调波210-1-2-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05解调信号频谱100050000 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5%抑制载波双边带调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs>2fc%x调制信号，t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ka调幅度t0=0.1;fs=10000;fc=1000;Vm0=2;ka=0.25;x=5*cos(200*pi*n);%n=-t0/2:1/fs:t0/2;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,2000]);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,3000]);title('解调信号频谱');调制信号50-5-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05已调波信号1050-5-10-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05调制信号频谱30002000100000 0.1 0.2 0.3 0.4 0.5 0.6 0.720001500100050001050-5-10
已调波信号频谱0.55 0.6 0.65 0.7 0.75解调波-15-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05解调信号频谱30002000100000 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5%抑制载波双边带调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs>2fc%x调制信号，t调制信号自变量,t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ka调幅度t0=0.1;fs=10000;fc=1000;Vm0=2;ka=0.25;x=3*cos(200*pi*n)-2*cos(250*pi*n);%调制信号n=-t0/2:1/fs:t0/2;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,1200]);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,2000]);title('解调信号频谱');调制信号50-5-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05已调波信号1050-5-10-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05调制信号频谱30002000100000 0.1 0.2 0.3 0.4 0.5 0.6 0.7已调波信号频谱100050000.55 0.6 0.65 0.7 0.75解调波1050-5-10-15-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05解调信号频谱20001500100050000 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5%单边带调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs>2fc%x调制信号，t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ka调幅度t0=0.1;fs=10000;fc=1000;Vm0=2;ka=0.25;n=-t0/2:1/fs:t0/2;N=length(n);x=5*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,3000]);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,3000]);title(‘解调波频谱’)调制信号50-5-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05已调波信号1050-5-10-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05调制信号频谱30002000100000 0.1 0.2 0.3 0.4 0.5 0.6 0.7300020001000050-5
已调波信号频谱0.55 0.6 0.65 0.7 0.75解调波-10-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05解调信号频谱30002000100000 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5%单边带调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs>2fc%x调制信号，t0采样区间，fs采样频率，%fc载波频率，Vm0输出载波电压振幅，ka调幅度t0=0.1;fs=10000;fc=1000;Vm0=2;ka=0.25;n=-t0/2:1/fs:t0/2;N=length(n);x1=3*cos(200*pi*n)-2*cos(250*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
温馨提示
- 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 深度解析(2026)《GBT 34263-2017工业用羟丙基甲基纤维素》
- 深度解析(2026)《GBT 34229-2017煤层气井试井方法 钻杆地层测试法》
- 2026年华南师范大学环境学院非事业编制人员招聘备考题库及参考答案详解一套
- 安钢总医院2026年度招聘备考题库参考答案详解
- 2025年佛山市南海区桂城街道映月中学教师招聘备考题库及1套完整答案详解
- 汕头大学2025年公开招聘非事业编工作人员备考题库及一套答案详解
- 福州市交通建设集团有限公司2025年度公开招聘备考题库有完整答案详解
- 北京医院（国家老年医学中心）2026年招聘104人备考题库（含答案详解）
- 平果市2025年度乡镇卫生院公开招聘乡村医生备考题库含答案详解
- 大理农林职业技术学院2026年公开招聘非编工作人员备考题库带答案详解
- 环泊酚的临床应用及注意事项
- 浅谈施工组织设计的优化
- 工程伦理与管理知到智慧树章节测试课后答案2024年秋山东大学
- 老房子维修安全协议书（2篇）
- 教育信息化促进乡村教育振兴
- DB45-T 2854-2024 侗族服饰标准规范
- CTD申报资料：创新药IND模块一-行政文件和药品信息
- 人教版四年级数学上册数与代数专项复习卷(含答案)
- 改善重症心肌梗死患者预后仍任重道远
- 呼吸道疾病健康知识讲座
- 离婚登记申请受理回执单模板
评论
0/150
提交评论