利用MATLAB进行s域分析

实验八利用MATLAB进行s域分析
一、实验目的
1、利用MATLAB实现部分分时展开式；
2、利用MATLAB计算H(S)的零极点与系统特性。

二、实验原理
MATLAB提供了函数residue可以得到复杂s域的部分分式展开式；MATLAB提供了函数roots可以计算H(S)的零极点。

三、程序示例见课本P255~258
四、实验内容
1、熟悉并验证示例程序
2、完成课本P265 M7-5；写实验报告（程序及波形）
num=[1 2];
den=[1 2 2 1];
[r,p]=residue(num,den);
sys=tf(num,den);
poles=roots(den)
figure(1);
pzmap(sys);
t=0:0.02:10;
h=impulse(num,den,t);
figure(2);
plot(t,h)
[H,w]=freqs(num,den);
figure(3);
plot(w,abs(H))
m=step(num,den,t);
figure(4);
plot(t,m)
运行结果为：
r =
1.0000
-0.5000 - 0.8660i
-0.5000 + 0.8660i
p =
-1.0000
-0.5000 + 0.8660i
-0.5000 - 0.8660i
poles =
-1.0000
-0.5000 + 0.8660i
-0.5000 - 0.8660i
-2-1.8-1.6-1.4-1.2-1
-0.8-0.6-0.4-0.20
P ole-Zero Map
Real Axis I m a g i n a r y A x i s
（a ）系统函数的零点分布图
012345678910
-0.20
0.2
0.4
0.6
0.8
1
1.2
(b)系统的冲激响应
012345678910
00.5
1
1.5
2
2.5
（c ）系统的频率响应
012345678910
00.5
1
1.5
2
2.5
（d ）系统的阶跃响应。