关于ode45求解二阶二元微分方程 - 数据结构与算法

编程论坛 → 开发语言 → 『数据结构与算法』 → 关于ode45求解二阶二元微分方程

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

只看楼主

nicksasaki

等　级：新手上路
帖　子：2
专家分：0
注　册：2016-12-12
结帖率：0

楼主

已结贴√ 问题点数：20 回复次数：2

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

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

2016-12-12 18:12

azzbcc

来　自：江西财经大学
等　级：贵宾
威　望：81
帖　子：3293
专家分：12919
注　册：2012-11-4

第 2 楼

得分:20

所以说，代码呢？

[fly]存在即是合理[/fly]

2016-12-12 19:17

nicksasaki

等　级：新手上路
帖　子：2
专家分：0
注　册：2016-12-12

第 3 楼

得分:0

回复 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)

求指导

2016-12-12 20:22

参与讨论请移步原网站贴子：https://bbs.bccn.net/thread-472325-1-1.html