Wolfram mathematica 用Solve求解线性系统问题
我有一个由4个线性方程组组成的系统,其变量是我从求解以前的系统中获得的,但尽管求解函数看起来是一个非常简单的系统,但它不会返回输出:Wolfram mathematica 用Solve求解线性系统问题,wolfram-mathematica,Wolfram Mathematica,我有一个由4个线性方程组组成的系统,其变量是我从求解以前的系统中获得的,但尽管求解函数看起来是一个非常简单的系统,但它不会返回输出: 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)
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可以处理的非数值表达式的大小是否有某种限制?您的系统没有解决方案。演示如何将其转换为标准线性代数形式可能会很有用:
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
在每个方程中具有相同的系数,因此您实际上有四个方程,只有三个未知数。空白输出可能意味着“无解”,即您的一些方程相互矛盾。尝试求解[{a==2,a==3},{a}]
,得到相同的结果。