作业帮 > 综合 > 作业

求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办?

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/11/06 01:52:39
求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办?
程序是
clear all;
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[T,Y]=ode45(@guanxi,[0 12],[0 1 1],options);
对应 guanxi.m 是
function dy=guanxi(t,y)
dy = zeros(2,1);
Dy(1)=Dy(2)*4*y(1)/3*((-2)*16*exp(-0.5*t)+85.5-y(2)) - Dy(2)*4*(-0.5)*16*exp(-0.5*t)/3*((-32)*exp(-0.5*t)+85.5-y(2));
自己认为可能的问题所在市我的微分式中有指数函数这个东西,所以不行~只是感觉~
初值不够呀,应该有四个,题目能不能贴出来?