Question

我试图使用Mathematica解决微分方程系统并发生错误：

DSolve :: deqx：提供的方程式不是给定函数的微分方程。

这是我的代码。

f = 100000;
s = 1;
DSolve[
 {
  X'[t] == -s*X[t]*w[t] + s*Z[t]*y[t] - 2 s*X[t]^2 - f*Iz[t]*X[t] + 
    2 f*Ix[t] + 3 f*Ix[t]*Z[t],
  Y'[t] == -s*Y[t]*x[t] + s*W[t]*z[t] - f*Iw[t]*Y[t] - 2 s*Y[t]^2 + 
    2 f*Iy[t] + 3 f*Iy[t]*W[t],
  Z'[t] == -s*Z[t]*y[t] + s*X[t]*w[t] - 2 s*Z[t]^2 + 2 f*Iz[t] + 
    3 f*Iz[t]*X[t] - f*Ix[t]*Z[t],
  W'[t] == -s*W[t]*z[t] - 2 s*W[t]^2 + 2 f*Iw[t] + 3 f*Iw[t]*Y[t] + 
    s*Y[t]*x[t] - f*Iy[t]*W[t],
  x'[t] == -s*X[t]*x[t] - s*Y[t]*x[t] + s*emp1[t],
  y'[t] == -s*Y[t]*y[t] - s*Z[t]*y[t] + s*emp2[t],
  z'[t] == -s*Z[t]*z[t] - s*W[t]*z[t] + s*emp3[t],
  w'[t] == -s*W[t]*w[t] - s*X[t]*w[t] + s*emp4[t],
  Ix'[t] == s*X[t]^2 - f*Ix[t] - f*Ix[t]*Z[t],
  Iy'[t] == s*Y[t]^2 - f*Iy[t] - f*Iy[t]*W[t],
  Iz'[t] == s*Z[t]^2 - f*Iz[t] - f*Iz[t]*X[t],
  Iw'[t] == s*W[t]^2 - f*Iw[t] - f*Iw[t]*Y[t],
  E0'[t] == -f*E0[t]*x[t] + f*E1[t]*y[t] - f*E0[t]*z[t] + f*E1[t]*w[t],
  E1'[t] == -E0'[t],

  emp1'[t] == 0,
  emp2'[t] == 0,
  emp3'[t] == 0,
  emp4'[t] == 0,
  emp1[0] == 100,
  emp2[0] == 100,
  emp3[0] == 100,
  emp4[0] == 100,
  X[0] == 1,
  Y[0] == 1,
  Z[0] == 0,
  W[0] == 0,
  x[0] == 0,
  y[0] == 0,
  z[0] == 0,
  w[0] == 0,
  Ix[0] == 0,
  Iy[0] == 0,
  Iz[0] == 0,
  Iw[0] == 0,
  E0[0] == 1,
  E1[0] == 0

  },
 {X[t], Y[t], Z[t], W[t], x[t], y[t], z[t], w[t], Ix[t], Iy[t], Iz[t],
   Iw[t], E0[t], E1[t]},
 t]

我在使用DSolve之前尝试清除所有变量，但它根本不起作用。

有人可以帮助我吗？谢谢！

Answer 1

需要将

emp1 [t]到emp4 [t]添加到倒数第二行的花括号中的函数列表中。

DSolve :: deqx：提供的方程式不是给定函数的微分方程

1 个答案: