[已答复] 非线性微分方程的数值解

[复制链接]
2345@2345 发表于 2022-5-12 12:07:24
求助大神非线性微分方程组的数值解

5QOIYMSYX(A2@BM8[YXD2TD.png





function [ f ] = egl_21( t,i)
f=[2*i(2)*i(1)-1*i(1);-2*i(2)*i(1)];
[T,I]=ode23('egl_21',[0,62],[0.00000032;0.99999968]);
figure(1)
plot(T,I(:,1),'g-',T,I(:,2),'r--')


end
初值条件是t=0时,i(1)=0.00000032,i(2)=0.99999968,这里我把i(1)作为i的解,i(2)作为s的解
用这个程序不知道哪里出错了,求大神指点

1 条回复


谢中华 发表于 2022-5-12 16:24:20
  1. fun = @(t,s)[2*s(2)*s(1)-s(1);-2*s(2)*s(1)];
  2. [t,s] = ode45(fun,[0,62],[0.00000032;0.99999968]);
  3. plot(t,s(:,1),'g-',t,s(:,2),'r--')
复制代码
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表