Why am i getting an unknown variable while taking the inverse Laplace of a 4th order system?

조회 수: 7 (최근 30일)
I am calculating the inverse Laplace of a 4th order transfer function with the following code:
syms s t
B=((0.8236*s^2+0.4877*s+1.315)/(s^4+1.385*s^3+1.978*s^2+1.317*s));
y=ilaplace(B,s,t)
Here is the ilaplace i get:
y =
1315/1317 - (3682969*sum((r3*exp(r3*t))/(2*(1500*r3^2 + 1385*r3 + 989)), r3 in RootOf(s3^3 + (277*s3^2)/200 + (989*s3)/500 + 1317/1000, s3)))/6585 - (1315000*sum((r3^2*exp(r3*t))/(2*(1500*r3^2 + 1385*r3 + 989)), r3 in RootOf(s3^3 + (277*s3^2)/200 + (989*s3)/500 + 1317/1000, s3)))/1317 - (19587691*sum(exp(r3*t)/(2*(1500*r3^2 + 1385*r3 + 989)), r3 in RootOf(s3^3 + (277*s3^2)/200 + (989*s3)/500 + 1317/1000, s3)))/13170
After that i am calculating integral squared error with the following code:
h=(1/s);
Y=ilaplace(h,s,t);
Ise=int((Y-y)^2,t,0,100)
and getting :
Ise =
((5260*exp(100*r3))/1734489 + (1729225*exp(200*r3))/3468978 - 193305/385442)/r3 + 400/1734489
Now, the problem is that i am getting this unknown variable 'r3' due to which i am unable to calculate a numerical value of integral squared error when i am using :ISE=double(Ise)
Can anyone please let me know how to get a numerical value of integral squared error. Thank you in advance.
  댓글 수: 6
Preeti Prajapati
Preeti Prajapati 2015년 7월 19일
Sir, i tried the partfrac command, but it showed an error 'Undefined function 'partfrac' for input arguments of type 'sym''. I am working in 2013a. Does it only work in 2015a version?
Star Strider
Star Strider 2015년 7월 19일
Yes.
I also posted the inverse Laplace transform of the partial fraction decomposition.

댓글을 달려면 로그인하십시오.

답변 (0개)

카테고리

Help CenterFile Exchange에서 Calculus에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by