symbolic Integration issue while solving differential equation

조회 수: 12 (최근 30일)
Santosh Bhagat
Santosh Bhagat 2018년 7월 4일
댓글: Torsten 2018년 7월 4일
I am trying to solve differential equation, but i instead of getting explicit solution, i am getting 'int ' expression as shown below. kindly please suggest how i can solve this. Thanks in advance.
yr1=-0.493*exp(-0.1099*t)+0.493*exp(0.0049*t)*cos(1.0005*t); Diffrential equation-> y1(t) = dsolve(10*diff(y1, t, 3)+diff(y1, t, 2)+10*diff(y1, t)+1.1*y1== 0.55*yr1*yr1, y1(0) == 0, Dy1(0) == 0, D2y1(0) == 0);
expression-> int(-(72185553*exp(-(1099*x)/5000)*exp((x*(60*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(1/3) + 900*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3) - 299))/(1800*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(1/3)))*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3)*(exp((287*x)/1250)/2 + (cos((2001*x)/1000)*exp((287*x)/1250))/2 - 2*cos((2001*x)/2000)*exp((287*x)/2500) + 1)*(2700*cos((3^(1/2)*x*(900*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3) + 299))/(1800*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(1/3)))*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3) - 897*cos((3^(1/2)*x*(900*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3) + 299))/(1800*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(1/3))) + 299*3^(1/2)*sin((3^(1/2)*x*(900*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3) + 299))/(1800*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(1/3))) + 900*3^(1/2)*sin((3^(1/2)*x*(900*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3) + 299))/(1800*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(1/3)))*((205957^(1/2)*5400000^(1/2))/5400000 - 259/6750)^(2/3)))/(800*(1390209750*3^(1/2) - 259*3^(1/2)*205957^(1/2)*5400000^(1/2))), x, 0, t, 'IgnoreSpecialCases', true, 'IgnoreAnalyticConstraints', true)

답변 (1개)

Torsten
Torsten 2018년 7월 4일
편집: Torsten 2018년 7월 4일
Seems the solution to your differential equation is given by this integral, but MATLAB is not able to find an analytical expression for this integral. Use a numerical ODE integrator instead (e.g. ODE45).
Best wishes
Torsten
  댓글 수: 2
Santosh Bhagat
Santosh Bhagat 2018년 7월 4일
when i put the solution in 'vpa' then vpa changes 'int' in 'numeric::int' with limits. so is there any chance to find the solution through numeric::int, and second how to use this numeric::int in symbolic matlab not in mupad.
Thanks
Torsten
Torsten 2018년 7월 4일
so is there any chance to find the solution through numeric::int
If you mean "analytical expression for the solution" by "solution" - then no. If you mean "numerical solution", then yes because vpa gives the numerical solution.
how to use this numeric::int in symbolic matlab not in mupad
You can't use this solution for symbolic computations.
Best wishes
Torsten.

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

카테고리

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

제품

Community Treasure Hunt

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

Start Hunting!

Translated by