matlab求微分方程数值解ode45出错,怎么办?

来源：学生作业帮编辑：作业帮分类：综合作业时间：2024/07/01 14:18:27

matlab求微分方程数值解ode45出错,怎么办?
建立equation.m文件,里面写入：
function dy =equation(t,y)
dy=zeros(3,1);
dy(1)=10^9*0.6*(pi/2+atan((0.5*y(1)+0.5*y(2)+0.5*y(3)+1)^0.3*(1.5*y(1)+1.5*y(2)+1.5*y(3)+1)^0.3*(y(1)^0.6*y(2)^0.4)^0.3-3))-y(1)*0.5*(3-(pi/2+atan((0.5*y(1)+0.5*y(2)+0.5*y(3)+1)^0.3*(1.5*y(1)+1.5*y(2)+1.5*y(3)+1)^0.3*(y(1)^0.6*y(2)^0.4)^0.3)));
dy(2)=0.25*(y(1)^0.6*y(2)^0.4)-10^(-7)*y(2)*y(3);
dy(3)=0.1*y(3)*(y(2)*1000/20-y(3));
end
然后在命令窗口输入：
[T,Y]=ode45('equation',[0 10],[1.636*10^7 13597 10^6]);
plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+');
运行一段时间后,出现如下错误：
Error using ==> horzcatOut of memory.Type HELP MEMORY for your options.
Error in ==> ode45 at 487 yout = [yout,zeros(neq,chunk,dataType)];
Error in ==> system_dynamic at 1[T,Y]=ode45('equation',[0 10],[1.636*10^7 13597 10^6]);

这里面的初始值可真够大的,估计许用超出内存了

求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办? MATLAB 常微分方程数值解出错 matlab用ode45 解微分方程 matlab数值解法求解二阶微分方程 ODE45函数 matlab ode45 求解二阶常微分方程关于matlab ode45解微分方程解不出来 matlab ode45解线性微分方程时的参数 matlab中的ode45怎样设置定步长解微分方程? MATLAB 求解微分方程数值解常微分方程数值解Matlab 求大神帮忙用matlab求解微分方程的数值解 matlab ode45 解方程