查看: 97|回复: 7|关注: 0
打印 上一主题 下一主题

[已解决] 常微分方程组的数值解

[复制链接]

新手

7 麦片

财富积分


050


2

主题

7

帖子

0

最佳答案
求帮助啊 !!!

新手

5 麦片

财富积分


050


0

主题

2

帖子

0

最佳答案
2#
发表于 4?天前 | 只看该作者
函数结尾有end吗?

新手

7 麦片

财富积分


050


2

主题

7

帖子

0

最佳答案
3#
?楼主| 发表于 4?天前 | 只看该作者
您好,感谢回帖,有没有end不影响,问题在于那个 H 该怎么输入 不是没有定义H, 就是输入参数不足 望回复 万分感谢

新手

18 麦片

财富积分


050


9

主题

240

帖子

1

最佳答案
4#
发表于 4?天前 | 只看该作者 |此回复为最佳答案
function dy =syg(t,y)
H=sqrt(0.3909-0.1303*y(1)-2*y(2)-y(3));
dy(1)=1.79*106^(-4)*(1-y(1))*H;
dy(2)=1.79*106^(-4)*(1-y(1))-1.855*10^(-4)*y(2)*H;
dy(3)=1.855*10^(-4)*y(2)*H-6.49*10^(-4)*(y(3)*H-0.1303*(3*y(1)-2*y(2)-y(3)/5*H));
dy=[dy(1);dy(2);dy(3)]
end

新手

7 麦片

财富积分


050


2

主题

7

帖子

0

最佳答案
5#
?楼主| 发表于 4?天前 | 只看该作者
function dy =syg(t,y)
H=sqrt(0.3909-0.1303*(3*y(1)-2*y(2)-y(3)));
dy(1)=1.79*10^(-4)*(1-y(1))*H;
dy(2)=(1.79*10^(-4)*(1-y(1))-1.855*10^(-4)*y(2))*H;
dy(3)=1.855*10^(-4)*y(2)*H-6.49*10^(-4)*(y(3)*H-0.1303*(y(1)-y(2)-y(3))*(3*y(1)-2*y(2)-y(3))/(5*H));
dy=[dy(1);dy(2);dy(3)];
end

这个确实可以运行了??不过我执行后y(3)列结果为零
[T,Y]=ode45(@syg,[0:0.001:0.02],[0,0,0])
真看不出哪错了??有眼尖的大神吗

新手

7 麦片

财富积分


050


2

主题

7

帖子

0

最佳答案
6#
?楼主| 发表于 4?天前 | 只看该作者
tyut2002489 发表于 2019-9-20 18:03
function dy =syg(t,y)
H=sqrt(0.3909-0.1303*y(1)-2*y(2)-y(3));
dy(1)=1.79*106^(-4)*(1-y(1))*H;

function dy =syg(t,y)
H=sqrt(0.3909-0.1303*(3*y(1)-2*y(2)-y(3)));
dy(1)=1.79*10^(-4)*(1-y(1))*H;
dy(2)=(1.79*10^(-4)*(1-y(1))-1.855*10^(-4)*y(2))*H;
dy(3)=1.855*10^(-4)*y(2)*H-6.49*10^(-4)*(y(3)*H-0.1303*(y(1)-y(2)-y(3))*(3*y(1)-2*y(2)-y(3))/(5*H));
dy=[dy(1);dy(2);dy(3)];
end
[T,Y]=ode45(@syg,[0:0.001:0.02],[0,0,0]
大神啊??可以运行了??不过结果不对??我更迷茫了??能再指导一下吗

论坛优秀回答者

权威

3783 麦片

财富积分



3

主题

3962

帖子

845

最佳答案
  • 关注者: 178
7#
发表于 昨天?08:59 | 只看该作者
syg123456 发表于 2019-9-20 19:32
function dy =syg(t,y)
H=sqrt(0.3909-0.1303*(3*y(1)-2*y(2)-y(3)));
dy(1)=1.79*10^(-4)*(1-y(1))*H;

数值太小了吧?
看微分方程、时间、初值。
106^(-4)约为7.9209e-09

新手

7 麦片

财富积分


050


2

主题

7

帖子

0

最佳答案
8#
?楼主| 发表于 昨天?09:03 | 只看该作者
对 确实是数值太小了??书上的数据错了??感谢回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表