二阶微分方程组的MATLAB解法

二阶微分方程组的MATLAB解法
请告诉一种用ode45来解二阶微分方程组的方法,最好把算法写出来
0714angel 1年前 已收到1个回答 举报

kk艳阳天 幼苗

共回答了12个问题采纳率:91.7% 举报

一个二阶微分方程:
y''+y'+y=sin(t)
初始条件为y(0)=5,y'(0)=6.
过程:
先降阶为一阶微分方程组
y'=z
z'=-z-y+sin(t)
编制如下函数m文件
function dy=weifen(t,x)
dy=zeros(2,1);
%y=x(1)
%z=x(2)
dy(1)=x(2);
dy(2)=sin(t)-x(2)-x(1);
然后用ode45解方程
[t,y]=ode45(@weifen,[0 20],[5 6])
plot(t,y)就可以画出y和y'的图像

1年前

1
可能相似的问题
Copyright © 2024 YULUCN.COM - 雨露学习互助 - 17 q. 2.017 s. - webmaster@yulucn.com