实验三卷积积分与卷积和运算.docVIP

下载本文档

45
0
约3.03千字
约 9页
2017-03-24 发布于重庆
举报
版权申诉

实验三卷积积分与卷积和运算.doc

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

实验三卷积积分与卷积和运算

a=1000; t1=-5:1/a:5; f1=stepfun(t1,0); f2=stepfun(t1,-1/a)-stepfun(t1,1/a); subplot(231); plot(t1,f1); axis([-5,5,0,1.2]); xlabel(时间(t));ylabel(幅值f1(t));title(f1); subplot(232); plot(t1,f2); ylabel(幅值f2(t)); title(f2); y=conv(f1,f2); r=2*length(t1)-1; t=-10:1/a:10; subplot(233); plot(t,y); axis([-5,5,0,1.2]); title(f1和f2的卷积); ylabel(f(t)); f11=conv(f1,f1);f22=conv(f2,f2); subplot(234); plot(t,f11); title(f1和f1的卷积); ylabel(f11(t)); axis([-5,5,0,5000]); subplot(235); plot(t,f22); title(f2和f2的卷积); ylabel(f22(t)); s=0.01; k1=0:s:2; k2=k1; f1=3*k1; f2=3*k2; f=conv(f1,f2); f=f*s; k0=k1(1)+k2(1); k3=length(f1)+length(f2)-2; k=k0:s:k3*s; subplot(3,1,1); plot(k1,f1); title(f1(t)); subplot(3,1,2); plot(k2,f2); title(f2(t)); subplot(3,1,3); plot(k,f); title(f(t)) clear all; T=0.1; t=0:T:10; f=sin(t); h=0.5*(exp(-t)+exp(-3*t)); Lf=length(f); Lh=length(h); for k=1:Lf+Lh-1 y(k)=0; for i=max(1,k-(Lh-1)):min(k,Lf) y(k)=y(k)+f(i)*h(k-i+1) ; end yzsappr(k)=T*y(k) ; end subplot(3,1,1); plot(t,f); title(f(t)); subplot(3,1,2); plot(t,h); title(h(t)); subplot(3,1,3); plot(t,yzsappr(1:length(t))); title(卷积近似计算结果); xlabel(时间); a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel(n);ylabel(幅值); k1=-5:15; f1=[zeros(1,5),ones(1,16)]; subplot(3,1,1); stem(k1,f1); title(f1(k)); k2=k1; f2=[zeros(1,5),ones(1,3),zeros(1,13)]; subplot(3,1,2); stem(k2,f2); title(f2(k)); k3=k1(1)+k2(1):k1(end)+k2(end); f3=conv(f1,f2); subplot(3,1,3); stem(k3,f3); title(f3(k)); a=1000; t1=-5:1/a:5; f1=stepfun(t1,-1/a)-stepfun(t1,1/a); f2=stepfun(t1,0); subplot(331); plot(t1,f1); axis([-5,5,0,1.2]); xlabel(时间(t));ylabel(幅值f1(t));title(单位冲激f1); subplot(332); plot(t1,f2); axis([-5,5,0,1.2]); xlabel(时间(t)); ylabel(幅值f2(t)); title(单位阶跃f2); t=0:0.0001:8; t2=0;t3=4; u1=stepfun(t,t2); u2=stepfun(t,t3); subplot(333); f3=u1-u2;plot(t,u); xlabel(时间(t));ylabel(幅值f3(t));title(u(t)-u(t-4)); y=conv(f1,f2); r=2*length(t1)-1; t=-10: