[信息与通信]数字信号处理课设报告.docVIP

下载本文档

5
0
约5.34千字
约 31页
2018-02-17 发布于浙江
举报
版权申诉

[信息与通信]数字信号处理课设报告.doc

1、本文档共31页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

[信息与通信]数字信号处理课设报告

%普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs2fc % 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(解调信号频谱); %普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs2fc % 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(解调信号频谱); %抑制载波双边带调制y=amod(x,t,fs,t0,fc,Vm0,ka)要求fs2fc % x调制信号， t0采样区间，fs采样频率， % fc载波频率，Vm0输出载波电压振幅，ka调幅度 t0=0.1; fs=10000; fc=1000; Vm0=2; ka=0.25; x=5*cos(2