云南大学信息学院信号系统实验二.docx
《云南大学信息学院信号系统实验二.docx》由会员分享，可在线阅读，更多相关《云南大学信息学院信号系统实验二.docx（9页珍藏版）》请在冰豆网上搜索。
云南大学信息学院信号系统实验二
《信号与系统》
实验报告
姓名:
学号:
20081专业:
信息学院08级通信工程
实验目的
1、利用MATLAB求连续系统的冲激响应与阶跃响应、求离散系统的单位响应与阶跃响应;
2、利用MATLAB求离散时间卷积和及连续时间卷积的数值近似;
3、比较手工计算结果与用MATLAB计算的结果异同,进一步加深对信号与系统理论知识的深入理解。
实验内容
8.1（2、8.3（3、8.5、8.6（3、8.7（a、8.8（2、8.9（2（5
实验设备软件硬件:
电脑
软件:
MATLAB
实验原理
1、对应知识点:
（1LTI连续时间系统的响应（2
（32、在编写
实验记录8.1
MATLABb=[1,3];
lsim（sys,f,t
012345678910
00.51
1.52
2.5
3LinearSimulationResults
Time（sec
Amplitude
8.3
（3y''（t+4y'（t+5y（t=f'（tMATLAB源程序如下:
a=[1,4,5];b=[1,0];
y1=step（b,a,0:
1:
10y2=impulse（b,a,0:
1:
10
subplot（2,1,1step（b,a,10subplot（2,1,2
impulse（b,a,10
8.6
（3y（n+y（n-1+（1/4y（n-2=x（n
a=[1,0,1/4];b=[1];n=1:
20;
y=impz（b,a,20subplot（2,1,1
stem（n,y,'filled'z=dstep（b,a,20subplot（2,1,2
stem（n,z,'filled'
00.5Time（sec0
1
234
5678910
00.5Time（sec
012345678910
LinearSimulationResults
Time（sec
Amplitude
y=10101010101010101010z=00112233445566778899
有单位序列响应的时域波形可见,该系统不是稳定系统。
8.7（a
x1=[0,1,2,1,0,0];x2=[0,1,1,1,1,0,0];n1=-2:
3;n2=-1:
5;
x=conv（x1,x2;ns=n1（1+n2（1;
leg=length（x1+length（x2-2;
n=ns:
（ns+leg;
stem（n,x,'filled'
8.9（2
建立M文件
functionf=gggfconv（f1,f3,t1,t2
f=conv（f1,f3;
f=f*0.01;
ts=t1（1+t2（1;
leg=length（t1+length（t2-2;
t=ts:
0.01:
（ts+leg*0.01;
subplot（2,2,1
plot（t1,f1
axis（[min（t1,max（t1,min（f1-min（f1*0.2,max（f1+max（f1*0.2]title（'f1（t'
xlabel（'t'
subplot（2,2,2
plot（t2,f2
axis（[min（t2,max（t2,min（f2-min（f2*0.2,max（f2+max（f2*0.2]title（'f2（t'
xlabel（'t'
subplot（2,2,3
plot（t,f;
axis（[min（t,max（t,min（f-min（f*0.2,max（f+max（f*0.2]
p=get（gca,'position';
p（3=2.4*p（3;
set（gca,'position',p;
title（'f（t=f1（t*f2（t'
xlabel（'t'
t1=0:
0.01:
2;
t2=1:
0.01:
3;
t3=-2:
0.01:
2;
f1=Heaviside（t1-Heaviside（t1-2;
f2=Heaviside（t2-1-Heaviside（t2-3;
f3=（0.5.*t3+1.*（t3<0
+（-0.5.*t3+1.*（t3>=0;gggfconv（f1,f2,t1,t2
t1=0:
0.01:
2;
t2=1:
0.01:
3;
t3=-2:
0.01:
2;
f1=Heaviside（t1-Heaviside（t1-2;
f2=Heaviside（t2-1-Heaviside（t2-3;
f3=（0.5.*t3+1.*（t3<0+（-0.5.*t3+1.*（t3>=0;gggfconv（f1,f3,t1,t3
（5
建立M文件
functionf=gggfconv（f3,f4,t3,t4
f=conv（f3,f4;
f=f*0.01;
ts=t3（1+t2（4;
leg=length（t3+length（t4-2;
t=ts:
0.01:
（ts+leg*0.01;
subplot（2,2,1
plot（t3,f3
axis（[min（t3,max（t3,min（f3-min（f3*0.2,max（f3+max（f3*0.2]title（'f3（t'
xlabel（'t'
subplot（2,2,2
plot（t4,f4
axis（[min（t4,max（t4,min（f4-min（f4*0.2,max（f4+max（f4*0.2]title（'f4（t'
xlabel（'t'
subplot（2,2,3
plot（t,f;
axis（[min（t,max（t,min（f-min（f*0.2,max（f+max（f*0.2]
p=get（gca,'position';
p（3=2.4*p（3;
set（gca,'position',p;
title（'f（t=f3（t*f4（t'
xlabel（'t'
t1=0:
0.01:
2;
t2=1:
0.01:
3;
t4=-2:
0.01:
2;
f1=Heaviside（t1-Heaviside（t1-2;
f2=Heaviside（t2-1-Heaviside（t2-3;
f4=（0.5.*t4+1.*（t4<0+（-0.5.*t4+1.*（t4>=0;gggfconv（f3,f4,t3,t4
f1（t11f2（t0.50.5000.51t1.52011.52t2.53f（t=f1（t*f2（t21.510.5011.522.53t3.544.55t1=0:
0.01:
2;t2=1:
0.01:
3;t4=-2:
0.01:
2;f1=Heaviside（t1-Heaviside（t1-2;f2=Heaviside（t2-1-Heaviside（t2-3;f4=（0.5.*t4+1.*（t4<0+（-0.5.*t4+1.*（t4>=0;gggfconv（f1,f4,t1,t4实验分析及实验小实验分析及实验小结分析及实验通过本次试验，我对MATLAB有了更深入的了解，试验中也遇到了些问题，有的程序编了又编才运行出来，实验过后了解了连续时间系统的时域分析。
升级会员 