怎样用matlab解非线性微分方程请问用matlab怎么解dx/dt=2y^2-2dy/dt=2yz-x-1dz/zt=2z^2-2y-4x(3)=1,y(3)=0,z(3)=2 只知道一个边值求在[0,3]上的解,有表达式更好,数值解也行,或者能得出[0,3]上的曲线也行,
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/30 03:03:48
![怎样用matlab解非线性微分方程请问用matlab怎么解dx/dt=2y^2-2dy/dt=2yz-x-1dz/zt=2z^2-2y-4x(3)=1,y(3)=0,z(3)=2 只知道一个边值求在[0,3]上的解,有表达式更好,数值解也行,或者能得出[0,3]上的曲线也行,](/uploads/image/z/3687858-18-8.jpg?t=%E6%80%8E%E6%A0%B7%E7%94%A8matlab%E8%A7%A3%E9%9D%9E%E7%BA%BF%E6%80%A7%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E8%AF%B7%E9%97%AE%E7%94%A8matlab%E6%80%8E%E4%B9%88%E8%A7%A3dx%2Fdt%3D2y%5E2-2dy%2Fdt%3D2yz-x-1dz%2Fzt%3D2z%5E2-2y-4x%283%29%3D1%2Cy%283%29%3D0%2Cz%283%29%3D2+%E5%8F%AA%E7%9F%A5%E9%81%93%E4%B8%80%E4%B8%AA%E8%BE%B9%E5%80%BC%E6%B1%82%E5%9C%A8%5B0%2C3%5D%E4%B8%8A%E7%9A%84%E8%A7%A3%2C%E6%9C%89%E8%A1%A8%E8%BE%BE%E5%BC%8F%E6%9B%B4%E5%A5%BD%2C%E6%95%B0%E5%80%BC%E8%A7%A3%E4%B9%9F%E8%A1%8C%2C%E6%88%96%E8%80%85%E8%83%BD%E5%BE%97%E5%87%BA%5B0%2C3%5D%E4%B8%8A%E7%9A%84%E6%9B%B2%E7%BA%BF%E4%B9%9F%E8%A1%8C%2C)
怎样用matlab解非线性微分方程请问用matlab怎么解dx/dt=2y^2-2dy/dt=2yz-x-1dz/zt=2z^2-2y-4x(3)=1,y(3)=0,z(3)=2 只知道一个边值求在[0,3]上的解,有表达式更好,数值解也行,或者能得出[0,3]上的曲线也行,
怎样用matlab解非线性微分方程
请问用matlab怎么解
dx/dt=2y^2-2
dy/dt=2yz-x-1
dz/zt=2z^2-2y-4
x(3)=1,y(3)=0,z(3)=2 只知道一个边值
求在[0,3]上的解,有表达式更好,数值解也行,或者能得出[0,3]上的曲线也行,
怎样用matlab解非线性微分方程请问用matlab怎么解dx/dt=2y^2-2dy/dt=2yz-x-1dz/zt=2z^2-2y-4x(3)=1,y(3)=0,z(3)=2 只知道一个边值求在[0,3]上的解,有表达式更好,数值解也行,或者能得出[0,3]上的曲线也行,
%by dynamic
%see also http://www.matlabsky.com
%contact me matlabsky@gmail.com
%2009.2
%
%关于【微分法的Matlab求解详解】参见http://www.matlabsky.com/thread-526-1-2.html其中对微分方程的Matlab讲解分析的很透彻
%
%选取状态变量
%x1=x,x2=y,x3=z
%各状态变量的一阶导数
%x1'=x'=2*x(2)^2-2
%x2'=y'=2*x(2)*x(3)-x(1)-1
%x3'=z'=2*x(3)^2-2*x(2)-4
%
%微分方程表达式
odefun=@(t,x)[2*x(2).^2-2
2*x(2).*x(3)-x(1)-1
2*x(3).^2-2*x(2)-4];
%初值,就是各个状态变量的初值
x0=[1 0 2];
%微分变量区间
tspan=[3 0]; %注意这里必须写成[3 0]写成[0 3]是错误的,一楼就是错误的
%调用ode45函数
[t,x]=ode45(odefun,tspan,x0);
%绘制[x,y,z]曲线
plot3(x(:,1),x(:,2),x(:,3))
x