求助，用 MATLAB 的 ode45 求解微分方程组
答案:1 悬赏:10 手机版
解决时间 2021-03-27 13:46
- 提问者网友：最美的风景
- 2021-03-26 19:29
求助，用 MATLAB 的 ode45 求解微分方程组
最佳答案
- 五星知识达人网友：等灯
- 2021-03-26 20:16
例如：求解下列微分方程组
求解步骤：
①自定义函数 rigid。m
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
②在执行窗口下，执行下列命令
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[t,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
③求解结果
求解步骤：
①自定义函数 rigid。m
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
②在执行窗口下，执行下列命令
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[t,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
③求解结果
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
大家都在看
推荐资讯