要解一个挺大的方程组,代码太长我就不贴完了,关键就是这么几行:
eq1 = Dt[Subscript[q, m] == \[Rho] * A * V]
(*这里有 eqxxx 的定义,方程组*)
SolveResult =
Solve[eqs, {Dt[Ma], Dt[V], Dt[p], Dt[\[Rho]], Dt[T], Dt[Subscript[p, t]], Dt[s], Dt[F]}]
其中一个解出来的结果初步处理之后是这样的:
-(((2 + (-1 + k) Ma^2) (Subscript[q, m] -
A k Ma^2 V (-1 + y) \[Rho]))/(2 A (-1 + Ma^2) V \[Rho]))
然而我想进一步化简它,此过程中需要传递给Simply一个Assuming:
Subscript[q, m] == \[Rho] * A * V
但是一旦执行了这一句,上面最初定义的eq1,再次计算的时候就直接变成了一个`eq1 = True`,再也没法拿来当方程用了
这该怎么办呢?
eq1 = Dt[Subscript[q, m] == \[Rho] * A * V]
(*这里有 eqxxx 的定义,方程组*)
SolveResult =
Solve[eqs, {Dt[Ma], Dt[V], Dt[p], Dt[\[Rho]], Dt[T], Dt[Subscript[p, t]], Dt[s], Dt[F]}]
其中一个解出来的结果初步处理之后是这样的:
-(((2 + (-1 + k) Ma^2) (Subscript[q, m] -
A k Ma^2 V (-1 + y) \[Rho]))/(2 A (-1 + Ma^2) V \[Rho]))
然而我想进一步化简它,此过程中需要传递给Simply一个Assuming:
Subscript[q, m] == \[Rho] * A * V
但是一旦执行了这一句,上面最初定义的eq1,再次计算的时候就直接变成了一个`eq1 = True`,再也没法拿来当方程用了
这该怎么办呢?
