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

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积分器吗？