matlab 微分方程组解不出来
来源:学生作业帮 编辑:搜搜考试网作业帮 分类:数学作业 时间:2024/07/12 12:42:13
matlab 微分方程组解不出来
m文件
function xdot=peng(t,x)a=87183;b=0.81;c=9.237;xdot=[x(2)+0.1*x(1)*sin(x(2));-(a/b)*x(4)*x(5)-(a/b)*c*x(5)^2*sign(x(1));x(4)-10*x(3)*x(5);-a*x(4)*x(5)-a*c*x(5)^2*sign(x(1))-x(3)*x(5)^2;x(1)^2/x(5)^2+x(3)^2/x(5)^2+(x(4)/x(5)^2+c*sign(x(1)/x(5)))^2];
命令窗口
[t,x]=ode45('peng',[0,10],[1,0,1.5,-9.24,1]);
m文件
function xdot=peng(t,x)a=87183;b=0.81;c=9.237;xdot=[x(2)+0.1*x(1)*sin(x(2));-(a/b)*x(4)*x(5)-(a/b)*c*x(5)^2*sign(x(1));x(4)-10*x(3)*x(5);-a*x(4)*x(5)-a*c*x(5)^2*sign(x(1))-x(3)*x(5)^2;x(1)^2/x(5)^2+x(3)^2/x(5)^2+(x(4)/x(5)^2+c*sign(x(1)/x(5)))^2];
命令窗口
[t,x]=ode45('peng',[0,10],[1,0,1.5,-9.24,1]);
![matlab 微分方程组解不出来](/uploads/image/z/18451411-43-1.jpg?t=matlab+%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E8%A7%A3%E4%B8%8D%E5%87%BA%E6%9D%A5)
xdot前加一句xdot=zeros(5,1)
sign符号函数不一定解得出来
再问: ����ʹ������ô�찡
sign符号函数不一定解得出来
再问: ����ʹ������ô�찡