>> impulse(sys)
0.6 0.5 0.4 Am
plit0.3 ude
0.2 0.1 0 -0.1 0
Impulse Response
1 2 3 4 5 6 7 8 9 10
Time (seconds)
>> a=[1 sqrt(2) 1]; >> b=[1 0 0]; >> sys=tf(b,a); >> impulse(sys)
0.2 0 -0.2 -0.4 Am
plit-0.6 ude
-0.8 -1 -1.2 -1.4 -1.6 0
1
2
3
4
5
6
7
8
9
Time (seconds)
>> a=[1 1 1]; >> b=[1 0]; >> sys=tf(b,a); >> impulse(sys)
1 0.8 0.6 Am
plit0.4 ude
0.2 0 -0.2 -0.4 0
2 4 6
Time (seconds)
8 10 12
>> a=[1 1 1]; >> b=[1 0 1]; >> sys=tf(b,a); >> impulse(sys) 0.4 0.2 0 Am
plit-0.2 ud
e -0.4 -0.6 -0.8 -1 0
2
4
6
Time (seconds)
8
10
12
2、已知某系统可以由微分方程描述
y(n)+y(n-1)+y(n-2)=x(n)
①请利用MATLAB绘出该系统冲激响应和阶跃响应的时域波形； ②根据冲激响应的时域波形分析系统的稳定性； ③如果系统的输入为，求系统的零状态响应;
>> a=[1 1 6]; >> b=[1]; >> sys=tf(b,a); >> impulse(sys) >> step(sys)
0.5 0.4 0.3 Am
plit0.2 ude
0.1 0 -0.1 -0.2 0
2 4 6 8 10 12 14
Time (seconds)
0.35 0.3 0.25 Am
plit0.2 ude
0.15 0.1 0.05 0 0
Step Response
2 4 6 8 10 12 14
Time (seconds)
是稳定系统
>> t=0:0.01:10; >> x=exp(-t); >>lsim(sys,x,t)
1.2 1 0.8 Am
plit0.6 ude
0.4 0.2 0 -0.2 0
Linear Simulation Results
1 2 3 4 5 6 7 8 9 10
Time (seconds)
3、已知描述离散系统的差分方程如下，试采用MATLAB绘出各系统的单位抽样响应，并根据单位抽样响应的时域波形分析系统的稳定性。
① y(n)+3y(n-1)+2y(n-2)=x(n)
>> a=[1 3 2]; >> b=[1]; >> impz(b,a)
不是稳定系统
5、采用MATLAB计算如下两个序列的卷积，并绘出图形
x1(n)?[1,2,1,1]?x2(n)?1,(?2?n?2)
>> x1=[1 2 1 1]; >> x2=[1 1 1 1 1]; >> x=conv(x1,x2) x =
1 3 4 5 5 4 2 1