Second order diffrential equation solve using ode 45

4 次查看（过去 30 天）
Abhik Saha
Abhik Saha 2022-8-22
评论： Abhik Saha 2022-8-22
I want to solve a second order differential equation of the form
d^2y/dt^2+dy/dt=sin(100t)
I want the value of dy/dt and y.
I am copied my code below where easily I can get the dy/dt as a function of t but I am not able find y as a function of t. I am new to Matlab. Please help me regarding this issue. I copied the code below
---------------------------------------------
main program
---------------------------------------------
t=linspace(0,500,1000); %time span%
y0=[0];
%%solving using ode45
[tsol, ysol]=ode45(@(t,y0) firstodefun2(t,y0), t, y0);
---------------------------------------------------
function definition
----------------------------------------------------
function dy=firstodefun2(t,y0)
dy=zeros(1,1);
dy(1)=sin(100*t)-y0(1);
end

请先登录，再回答此问题。

采纳的回答

Sam Chak
Sam Chak 2022-8-22
Ran in:
Hi @Abhik Saha
The code is fixed below. This is a 2nd-order system, and so there are two state variables and two initial values.
t = [0, 10];
y0 = [0 0];
% solving using ode45
[tsol, ysol] = ode45(@(t, y) firstodefun2(t, y), t, y0);
plot(tsol, ysol(:, 1)), grid on, xlabel('t, [sec]'), ylabel('y(t)')
function dy = firstodefun2(t, y)
dy = zeros(2, 1);
dy(1) = y(2);
dy(2) = sin(100*t) - y(2);
end

更多回答（1 个）

Torsten
Torsten 2022-8-22
Set x(1) = y and x(2) = dy/dt.
Then you get a system of differential equations
dx(1)/dt = x(2)
dx(2)/dt = -x(2) + sin(100*t)
Of course, you also will have to supply two initial conditions (one for y, the second for dy/dt at t=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!

Translated by