Unable to understand the output of the inverse laplace.
3 次查看(过去 30 天)
显示 更早的评论
I am attaching the program file. The output of the inverse Laplace is as follows:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/168577/image.jpeg)
I do want to know what does the term r4 mean?
0 个评论
采纳的回答
Walter Roberson
2018-1-8
r4 is a variable that was introduced to make the solution more compact.
Notice under the first summation it says "r4 in #1", and then it gives a definition for #1 in terms of RootOf(). What that is telling you is that the RootOf() has multiple roots (four of them in this case) and that a list of those roots is to be calculated and that the summation should then proceed with r4 being substituted with each member of that list of roots.
The RootOf() is being applied to a quartic equation. There are closed form solutions for that quartic, but they are typically quite long to write out, and if you were to write it all out then the result would be very difficult to make sense of as a human.
So in this case, the inverse laplace has a definite value, and it would be possible to list it all out, but your mind would probably go on strike if you tried to understand it in the expanded form.
As I look at those coefficients in the image, it looks to me as if you attempted to do an symbolic inverse laplace on a system that involved floating point values. Symbolic inverse laplace attempts to find a exact closed form solution if possible, but finding exact solutions does not make sense when the original coefficients are uncertain.
For example, something "simple" like 1.1E-6 might look like an exact coefficient, but internally it is 1.10000000000000005610204141526242693771564518101513385772705078125e-06 and from a scientific point of view, writing 1.1E-6 represents "some value, X, such that 1.05E-6 <= X < 1.15E-6"
If possible, you should convert all of your coefficients to symbolic as they are introduced, so that you can reduce round-off error. For example, instead of 1.1E-6 you might use sym(11)/sym(10^7)
0 个评论
更多回答(1 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!