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

5 visualizaciones (últimos 30 días)
Preeti Prajapati
Preeti Prajapati el 13 de Jul. de 2015
Comentada: Star Strider el 19 de Jul. de 2015
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 comentarios
Mostrar 4 comentarios más antiguos
Preeti Prajapati
Preeti Prajapati el 19 de Jul. de 2015
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 el 19 de Jul. de 2015
Yes.
I also posted the inverse Laplace transform of the partial fraction decomposition.

Iniciar sesión para comentar.

Iniciar sesión para responder a esta pregunta.

Respuestas (0)

Categorías

Más información sobre Calculus en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by