编程论坛 → 数据结构与算法

关于ode45求解二阶二元微分方程

nicksasaki 发布于 2016-12-12 18:12， 1578 次点击

求大神帮忙看看我的编程啊，始终不对

只有本站会员才能查看附件，请登录

2 回复

azzbcc2016-12-12 19:17

所以说，代码呢？

nicksasaki2016-12-12 20:22

回复 2楼 azzbcc

clc;clear;
t0=0;tf=5;
[t,y]=ode45('f1',[t0,tf],[0,0,0,0]);
figure;
plot(t,y(1),'r',t,y(3),'b');
hold on

%%%%%%%%
function dy=f1(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=-y(2)-10*y(1)-0.25*y(4)-2.5*y(3)
dy(4)=-554/256*y(4)-5540/256*y(3)-100/256*y(2)-1000/256*y(1)

求指导