使用Solve解决线性系统问题

Question

我有一个由变量表示的4个线性方程组，我从解决以前的系统得到的，但Solve函数不返回输出，尽管它看起来是一个非常简单的系统来解决：

Solve[{
 -d5c2 dn5t1 - d5c3 dn5t1 - a3 n3t1 - 
  (d4c1 n4t1 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2)
    ==  dn3t1,
 -a3 n3t2 - (d5c2 dn5t1 n5t2)/n5t1 - 
  (d5c3 dn5t1 n5t2)/n5t1 -
  (d4c1 n4t2 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2) 
    == dn3t2,
 -a3 n3t3 - 
  (d4c1 n4t3 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2) - 
  (d5c2 dn5t1 n5t3)/n5t1 - (d5c3 dn5t1 n5t3)/n5t1 
    ==  dn3t3,
 -a3 n3t4 - 
  (d4c1 n4t4 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2) -
  (d5c2 dn5t1 n5t4)/n5t1 - (d5c3 dn5t1 n5t4)/n5t1
    == dn3t4
       }, {a3, d5c2, d5c3, d4c1}]

返回空白输出：

{}

我是这门语言的新手; Solve可以处理的非数值表达式的大小是否存在某种限制？

Answer 1

您的系统没有解决方案。展示如何将其置于规范的线性代数形式可能是有用的：

sys={
 -d5c2 dn5t1 - d5c3 dn5t1 - a3 n3t1 - 
  (d4c1 n4t1 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2)
    ==  dn3t1,
 -a3 n3t2 - (d5c2 dn5t1 n5t2)/n5t1 - 
  (d5c3 dn5t1 n5t2)/n5t1 -
  (d4c1 n4t2 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2) 
    == dn3t2,
 -a3 n3t3 - 
  (d4c1 n4t3 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2) - 
  (d5c2 dn5t1 n5t3)/n5t1 - (d5c3 dn5t1 n5t3)/n5t1 
    ==  dn3t3,
 -a3 n3t4 - 
  (d4c1 n4t4 (dn4t2 n5t1 - dn4t1 n5t2))/(n4t2 n5t1 - n4t1 n5t2) -
  (d5c2 dn5t1 n5t4)/n5t1 - (d5c3 dn5t1 n5t4)/n5t1
    == dn3t4
       }

lhs = sys[[All, 1]];
rhs = sys[[All, 2]];

(m = Transpose[Coefficient[lhs, #] & /@ {a3, d5c2, d5c3, d4c1}]) // MatrixForm

此时您可以尝试LinearSolve[m,rhs]，但在这种情况下，它会报告

  

遇到没有解决方案的线性方程

我们认为这是因为行列式为零。

  Det[m] 

  

0

从根本上说，你的未知数d5c2和d5c3在每个方程中具有相同的系数，因此你实际上有四个方程式，只有三个未知数。