信息与通信数字信号处理课设报告_第1页
信息与通信数字信号处理课设报告_第2页
信息与通信数字信号处理课设报告_第3页
已阅读5页，还剩26页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

%普通调幅波调制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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论