文档介绍:实验目的:
1 学会运用MATLAB求拉普拉斯变换与逆变换,拉普拉斯变换法求解微分方程;
2 学会运用MATLAB求解系统函数的零极点,分析其分布和时域特性关系;
3 学会运用MATLAB分析系统函数的极点分布与系统稳定性的关系。
实验内容:(请将实验的题目内容、实验及过程代码、实验结果(必要时可以进行拷屏)、实验体会等填写到此处。页面空间不够,可另附页或另附文件。
。
ilaplace函数求的逆拉普拉斯变换。
。
,
起始条件为求系统的零状态响应,零输入响应和完全响应。
。
,同时分析系统函数极点对时域波形的影响。
程序代码
1. f=sym('exp(-t)*sin(a*t)');
L=laplace(f)
2. F=sym('s^2/(s^2+1)');
ft=ilaplace(F)
3. format rat;
B=[1,2];
A=[1,4,3,0];
[r,p]=residue(B,A)
4. syms t s
Yzis=(3*s+13)/(s*2+3*s+2);
yzi=ilaplace(Yzis)
xt=4*exp(-2*t)*Heaviside(t);
Xs=laplace(xt);
Yzss=Xs/(s*2+3*s+2);
yzs=ilaplace(Yzss)
yt=simplify(yzi+yzs)
5. b=[1,-2];
a=[1,4,5];
sys=tf(b,a);
pzmap(sys)
axis([-3,3,-2,2]);
6. b1=[1];
a1=[1 0];
sys1=tf(b1,a1);
subplot(121)
pzmap(sys1)
axis([-2,2,-2,2])
subplot(122)
impulse(b1,a1)
figure
b2=[1];
a2=[1 -1];
sys2=tf(b2,a2);
subplot(121)
pzmap(sys2)
axis([-2,2,-2,2])
subplot(122)
impulse(b2,a2)
figure
b3=[1];
a3=[1 1];
sys3=tf(b3,a3);
subplot(121)
pzmap(sys3)
axis([-2,2,-2,2])
subplot(122)
impulse(b3,a3)
figure
b4=[1];
a4=[1 -2 50];
sys4=tf(b4,a4);
subplot(121)
pzmap(sys4)
axis([-2,2,-8,8])
subplot(122)
impulse(b4,a4)
figure
b5=