How to solve 2nd ODE equation in numerical and analytical method at same plot graph?
2 ビュー (過去 30 日間)
古いコメントを表示
Hello,
I am was looking for help to solve this equation
d^x/dt^2 + dx/dt + x =1
in Matlab, can I get help to solve this equation of plot the numerical solutions and analytical solutions in the same graph and compare them.
with ode45
can anyone help me with this code?
0 件のコメント
採用された回答
Star Strider
2022 年 9 月 28 日
Interesting problem!
I was curious enough to do the entire simulation to see what the results are —
syms x(t) x0 Dx0 T Y
Dx = diff(x);
D2x = diff(Dx);
DEqn = D2x + Dx + x == 1;
xs = dsolve(DEqn, x(0)==x0, Dx(0)==Dx0)
xs = subs(xs, {x0,Dx0},{1,1}) % Numerical Initial Conditions
[VF,Sbs] = odeToVectorField(DEqn)
DEfcn = matlabFunction(VF, 'Vars',{T,Y})
tspan = [0 10];
[t,x] = ode89(DEfcn, tspan, [1 1]);
figure
fplot(xs, tspan, 'DisplayName','Analytic')
hold on
plot(t, x(:,1), 'DisplayName','Numeric')
hold off
grid
xlabel('Time')
ylabel('x_1(t)')
legend('Location','best')
They match almost exactly.
.
4 件のコメント
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Ordinary Differential Equations についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!