作业帮 > 综合 > 作业

怎么用matlab编辑一个微分方程组,比如:dx=y-0.5x dy=z-x dz=-x+sgn(x);

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/23 08:48:53
怎么用matlab编辑一个微分方程组,比如:dx=y-0.5x dy=z-x dz=-x+sgn(x);
其中sgn(x)是符号函数
用ode45编辑的程序是怎么样的?谢谢哪个大侠给我说一下
function DX=YD(t,X)
x=X(1);
y=X(2);
z=X(3);
dx=y-0.5*x;
dy=z-x;
dz=-x+sgn(x);
DX=[dx,dy,dz];
end
再问: 嗯,这个我之前已经弄出来了,我现在又遇到一个问题了,比如编写了个M文件: function dxdt=fun(t,x) x1=x(1); x2=x(2); x3=x(3); f(x1)=tanh(20*x1+40)+tanh(20*x1)+tanh(20*x1-40); dxdt=[ 10*x2-10*f(x1); x1-x2+x3; -15*x2; ]; end 然后再命令窗口中输入了: >> t_end=1000; >> x0=[0;1;0]; 还需要对f(x1)这个函数赋初值,但是我弄的初值都是错误的,你能给我说怎么算这个的初值吗?
再答: 初值的含义是微分方程组的初值,只要初值确定了,解曲线就清楚了。f(x) 是定义的一个函数还是中间变量?如果是中间变量则不需要对f(x1)赋值。有点乱,我不知你在算什么?是rkf积分器吗?