How graph second order differential equations with Matlab?

170 ビュー (過去 30 日間)
jake thompson
jake thompson 2018 年 1 月 18 日
回答済み: Phloykan Tepwarin 2022 年 1 月 18 日
Looking to get some help on how to use matlab to solve the following equation problem 6, all help is appreciated! Thank You! I haven't used matlab in 2 years very rusty, image is reference to original problem statement.
4y"-20y'+25y=0
%%PROBLEM 6
clear; syms y(x)
ode = 4*diff(y,x,2) - 20*diff(y,x) + 25*y == 0;
%ySol(x) = dsolve(ode)
%ezplot('ode')
ySol = dsolve(ode, y(0) == 0)
figure
subplot(2,1,1)
ezplot(real(ySol(1)))
subplot(2,1,2)
ezplot(imag(ySol(1)))
ERROR
Warning: Contour not rendered for constant ZData
> In contour (line 52)
In ezplot>ezimplicit (line 315)
In ezplot (line 153)
In sym/ezplot (line 61)
  1 件のコメント
Rena Berman
Rena Berman 2019 年 12 月 12 日
(Answers Dev) Restored edit

サインインしてコメントする。

採用された回答

Star Strider
Star Strider 2018 年 1 月 18 日
The problem is that you need to define an initial condition for the first derivative.
Try this:
syms y(x)
Dy = diff(y);
D2y = diff(y,2);
ode = 4*D2y - 20*Dy + 25*y == 0;
ySol = dsolve(ode, y(0) == 0, Dy(0) == 1) % Define The Initial Condition For ‘Dy(0)’ To Be ‘Some Value’
figure
ezplot(ySol)
  3 件のコメント
Star Strider
Star Strider 2018 年 1 月 18 日
As always, my pleasure!
Hamid Ghassri
Hamid Ghassri 2019 年 7 月 8 日
編集済み: Hamid Ghassri 2019 年 7 月 8 日
how can you set a range for x

サインインしてコメントする。

その他の回答 (2 件)

Phloykan Tepwarin
Phloykan Tepwarin 2022 年 1 月 18 日
yms s t Y u ; u = 10*exp(-t); ic = [0 0]; U = laplace(u, t, s); Y1 = s*Y; Y2 = s*Y1; S = Y2+9*Y-U; Y(s) = solve(S, Y); y = ilaplace(Y(s)); B = simplify(y); pretty(B); h = ezplot(y, [0 10]); set(h,'linewidth',2.5); title('function y(t)'); xlabel('t'); ylabel('y(t)');

Phloykan Tepwarin
Phloykan Tepwarin 2022 年 1 月 18 日
syms s t Y u ;
u = 10*exp(-t);
ic = [0 0];
U = laplace(u, t, s);
Y1 = s*Y; Y2 = s*Y1;
S = Y2+9*Y-U;
Y(s) = solve(S, Y);
y = ilaplace(Y(s));
B = simplify(y);
pretty(B);
h = ezplot(y, [0 10]);
set(h,'linewidth',2.5);
title('function y(t)');
xlabel('t');
ylabel('y(t)');

カテゴリ

Help Center および File ExchangeCalculus についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by