- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
. WORD格式整理. .
. .专业知识分享. .
【研讨题目2】 信号与系统时域分析专题研讨
【目的】
1.研究用离散方法近似计算连续信号的卷积积分;
2.通过分析近似计算卷积积分过程中出现的问题,锻炼学生分析问题和解决问题的能力;
【知识点】
信号时域分析,卷积积分,卷积和
【研讨题目】连续信号卷积积分的数值近似计算
两个连续信号的卷积积分定义为
为了能用数值方法进行计算,需对连续信号进行抽样。记x[k]=x(k?), h[k]=h(k),为进行数值计算所选定的抽样间隔,可以证明连续信号卷积积分可近似的表示为
(1)
由式(1)可知,可以利用Matlab提供的conv函数近似计算连续信号的卷积积分。
一、(*)理论分析
为了对近似计算的结果进行分析,用解析的方法计算下列卷积积分,推出卷积积分的解析表达式;
(1) 时限信号卷积积分
x1(t)=u(t)?u(t?1),y1(t)=x1(t)?x1(t);
卷积结果为:y1(t)= x1(t)?x1(t)=r(t)-2*r(t-1)+r(t-2)
(2) 分段常数信号卷积积分
x2(t)= x1(t)+2 x1(t?1)+ x1(t?2),h2(t)= x1(t)? x1(t?1), y2(t)=x2(t)?h2(t);
卷积结果为:y2(t)= x2(t)?h2(t)
=y1(t)+y1(t-1)-y1(t-2)-y1(t-3)
=r(t)-r(t-1)-2*r(t-2)+2*r(t-3)+r(t-4)-r(t-5)
(3) 非时限信号卷积积分
x3(t)=u(t),h3(t)=e?tu(t), y3(t)=x3(t)?h3(t)
卷积结果为:y3= x3(t)?h3(t) =[1-exp(-t)]*u(t)
二、(*)时限信号卷积积分的近似计算
取不同的△值,用Matlab函数conv近似计算卷积积分y1(t)并画出其波形,讨论?的取值对计算结果的影响。
上图中,绿线为间隔0.01的结果,蓝线是间隔0.1结果,红线为实际结果,
由此可见:时间间隔越小,与实际结果越接近。
附程序代码:
t1=[0:0.01:5];
t2=[0:0.1:5];
t=[0:0.1:5];
x1=1.*(t1=0)-1.*(t1=1);
x2=1.*(t2=0)-1.*(t2=1);
y1=convn(x1,x1);
y2=convn(x2,x2);
y=t.*[t=0]-2*(t-1).*[t=1]+(t-2).*[t=2]
N1=length(y1); %length函数取y1的长度%
N2=length(y2);
plot(t,y,r);
hold on;
plot(0:0.01:(N1-1).*0.01,y1*0.01,g);
plot(0:0.1:(N2-1).*0.1,y2*0.1,b);
axis([0 5 0 1])
三、(**)分段常数信号卷积积分的Matlab计算
(1)若x2[k]={1,2 ,1,0; k=0,1,2}, h2[k]= {1,1 ; k=0,1},计算离散卷积y2[k]=x2[k]?h2[k];
y2[k]=x2[k]?h2[k]结果如下:
附程序代码:x2=[1,2,1,0];h2=[1,-1];
y2=conv(x2,h2);
N=length(y2);
stem(0:N-1,y2);
axis([0 8 -1 1])
(2)比较y2(t)和y2[k],你发现了什么?
y2(t)的图像如下:
附程序代码:t=[0:0.1:5];
y=t.*[t=0]-(t-1).*[t=1]-2*(t-2).*[t=2]+2*(t-3).*[t=3]+(t-4).*[t=4]-(t-5).*[t=5]
plot(t,y);
hold on;
axis([0 8 -1 1])
y2(t)和y2[k]图像比较:
附程序代码:x2=[1,2,1,0];h2=[1,-1];
y2=conv(x2,h2);
t=[0:0.1:5];
y=t.*[t=0]-(t-1).*[t=1]-2*(t-2).*[t=2]+2*(t-3).*[t=3]+(t-4).*[t=4]-(t-5).*[t=5]
N=length(y2);
stem(0:N-1,y2);
hold on;
axis([0 8 -1 1])
plot(t,y);
hold on;
axis([0 8 -
文档评论(0)