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

조회 수: 7 (최근 30일)
Preeti Prajapati 2015년 7월 13일
댓글: Star Strider 2015년 7월 19일
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이전 댓글 4개 표시이전 댓글 4개 숨기기
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 2015년 7월 19일
Yes.
I also posted the inverse Laplace transform of the partial fraction decomposition.

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

### 카테고리

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