拉普拉斯逆变换,必须有一个更好的方法?

38视图(30天)
我有一个传递函数 我应用阶跃输入 。我的终极目标是找到时间响应
用手解决 ,我知道输出
然后我知道部分分式扩张
从上面我可以作拉普拉斯逆变换和发现
我的目标是获得 尽可能少的按键。
我发现最快的方法是使用残留执行部分分式扩张():
num = 9;
分母项= [1 9 9 0];
(r、磷、钾)=残渣(num,分母项);
使结果:
r =
0.1708
-1.1708
1.0000
p =
-7.8541
-1.1459
0
k =
[]
我可以写:
信谊年代
F = 0.1708 / (-1.1459 + 7.8541) / (s + 1.1708) + 1 / s;
ilaplace (F)
使结果:
ans =
(427 * exp (10000 - (78541 * t) /)) / 2500 - (11459 * exp (2500 - (2927 * t) /)) / 10000 + 1
这是我想要的答案!
但是必须有一个更好的方法!有人可以请告知?

接受的答案

明星黾
明星黾 2019年11月6日
试试这个:
信谊s t
num = 9;
分母项= [1 9 9 0];
% (r、磷、钾)=残渣(num,分母项);
num = num;
洞穴= poly2sym(分母项,s);%创造象征性的多项式
F = num /洞穴* 1 / s%与阶跃输入传递函数
基维辛迪= partfrac (F)%部分分式扩张
基维辛迪f = ilaplace)%时域表达式
f =重写(f,“经验”)指数条款%重写
f = vpa (f, 4)%将分数转换成十进制等价物
flatex =乳胶(f)%乳胶表达式
生产:
2的评论
明星黾
明星黾 2019年11月7日
像往常一样,这是我的荣幸!
partfrac 函数中引入 R2015a ,似乎是自动调用在过去的两个或三个版本。它不是自动在那之前,所以我特别包括在这里。

登录置评。

更多的答案(0)

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!