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-
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/27 16:06:19
![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-](/uploads/image/z/5170020-60-0.jpg?t=matlab%E6%B1%82%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E6%95%B0%E5%80%BC%E8%A7%A3ode45%E5%87%BA%E9%94%99%2C%E6%80%8E%E4%B9%88%E5%8A%9E%3F%E5%BB%BA%E7%AB%8Bequation.m%E6%96%87%E4%BB%B6%2C%E9%87%8C%E9%9D%A2%E5%86%99%E5%85%A5%EF%BC%9Afunction+dy+%3Dequation%28t%2Cy%29dy%3Dzeros%283%2C1%29%3Bdy%281%29%3D10%5E9%2A0.6%2A%28pi%2F2%2Batan%28%280.5%2Ay%281%29%2B0.5%2Ay%282%29%2B0.5%2Ay%283%29%2B1%29%5E0.3%2A%281.5%2Ay%281%29%2B1.5%2Ay%282%29%2B1.5%2Ay%283%29%2B1%29%5E0.3%2A%28y%281%29%5E0.6%2Ay%282%29%5E0.4%29%5E0.3-)
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-
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]);
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-
这里面的初始值可真够大的,估计许用超出内存了