线性代数/主题：计算机代数系统/解决方案

解决方案

本主题的答案使用 Maple 作为计算机代数系统。特别是，所有这些都在运行于 MS-DOS NT 版本 4.0 下的 Maple V 上进行了测试。（在所有这些中，加载线性代数包的预备命令以及 Maple 对 Enter 键的响应已被省略。）其他系统也有类似的命令。

其他答案将为 Wolfram Mathematica 13.0 添加。

{{TextBox|1=

问题 1

使用计算机解决本章开头提出的两个问题。

1. 这是静力学问题。

   ${\begin{array}{rl}40h+15c&=100\\25c&=50+50h\end{array}}$

2. 这是化学问题。

   ${\begin{array}{rl}7h&=7j\\8h+1i&=5j+2k\\1i&=3j\\3i&=6j+1k\end{array}}$

答案

1. 命令
   > A:=array( [[40,15],
   [-50,25]] );
   > u:=array([100,50]);
   > linsolve(A,u);
   
   得出答案 $[1,4]$. Mathematica 答案

   eqns = {40 h + 15 c == 100, 25 c == 50 + 50 h};
   {b, m} = CoefficientArrays[eqns, {c, h}]
   LinearSolve[m, -b]
   {c, h} /. Solve[eqns, {c, h}]
   

   返回 {c,h} 的向量 {1,4}。
2. 这里有一个自由变量
   > A:=array( [[7,0,-7,0],
   [8,1,-5,2],
   [0,1,-3,0],
   [0,3,-6,-1]] );
   > u:=array([0,0,0,0]);
   > linsolve(A,u);
   
   提示回复 $[\_t_{1},3\_t_{1},\_t_{1},3\_t_{1}]$. Mathematica 代码

   eqns = {7 h == 7 j, 8 h + i == 5 j + 2 k, 1 i == 3 j, 
      3 i == 6 j + 1 k};
   {b, m} = CoefficientArrays[eqns, {h, i, j, k}]
   LinearSolve[m, -b]
   {h, i, j, k} /. Solve[eqns, {h, i, j, k}]
   

   返回 LinearSolve 的 {0,0,0,0} 以及最后一行的 $\{\{h,3h,3h,3h\}\}$。

}}

问题 2

使用计算机解决第一小节中的这些系统，或得出结论“多个解决方案”或“无解决方案”。

1. ${\begin{array}{*{2}{rc}r}2x&+&2y&=&5\\x&-&4y&=&0\end{array}}$
2. ${\begin{array}{*{2}{rc}r}-x&+&y&=&1\\x&+&y&=&2\end{array}}$
3. ${\begin{array}{*{3}{rc}r}x&-&3y&+&z&=&1\\x&+&y&+&2z&=&14\end{array}}$
4. ${\begin{array}{*{2}{rc}r}-x&-&y&=&1\\-3x&-&3y&=&2\end{array}}$
5. ${\begin{array}{*{3}{rc}r}&&4y&+&z&=&20\\2x&-&2y&+&z&=&0\\x&&&+&z&=&5\\x&+&y&-&z&=&10\end{array}}$
6. ${\begin{array}{*{4}{rc}r}2x&&&+&z&+&w&=&5\\&&y&&&-&w&=&-1\\3x&&&-&z&-&w&=&0\\4x&+&y&+&2z&+&w&=&9\end{array}}$

Maple 答案

这些很容易输入。例如，第一个

> A:=array( [[2,2],
[1,-4]] );
> u:=array([5,0]);
> linsolve(A,u);

得到预期答案 $[2,1/2]$。其他输入方式类似。

1. 答案是 $x=2$ 和 $y=1/2$。
2. 答案是 $x=1/2$ 和 $y=3/2$。
3. 该系统有无穷多个解。在第一小节中，以 $z$ 为参数，我们得到 $x=(43-7z)/4$ 和 $y=(13-z)/4$。Maple 返回 $[-12+7\_t_{1},\_t_{1},13-4\_t_{1}]$，出于某种原因，它更喜欢 $y$ 作为参数。
4. 该系统没有解。当给出数组 $A$ 和向量 $u$ 并要求 Maplelinsolve(A,u)时，它根本没有返回结果，也就是说，它没有给出任何解。
5. 解是 $(x,y,z)=(5,5,0)$。
6. 有很多解。Maple 给出 $[1,-1+\_t_{1},3-\_t_{1},\_t_{1}]$。

Mathematica 答案

RowReduceAugmentedMatrix[matrix_] := MatrixForm[RowReduce[matrix]]

将此函数应用于 MatrixForm 中的 RowReduce 对 ${\begin{pmatrix}-1&1&1\\1&1&2\end{pmatrix}}$ 进行行化简，得到 ${\begin{pmatrix}1&0&1/2\\0&1&3/2\end{pmatrix}}$。

问题 3

使用计算机来解决第二小节中的这些方程组。

1. ${\begin{array}{*{2}{rc}r}3x&+&6y&=&18\\x&+&2y&=&6\end{array}}$
2. ${\begin{array}{*{2}{rc}r}x&+&y&=&1\\x&-&y&=&-1\end{array}}$
3. ${\begin{array}{*{3}{rc}r}x_{1}&&&+&x_{3}&=&4\\x_{1}&-&x_{2}&+&2x_{3}&=&5\\4x_{1}&-&x_{2}&+&5x_{3}&=&17\end{array}}$
4. ${\begin{array}{*{3}{rc}r}2a&+&b&-&c&=&2\\2a&&&+&c&=&3\\a&-&b&&&=&0\end{array}}$
5. ${\begin{array}{*{4}{rc}r}x&+&2y&-&z&&&=&3\\2x&+&y&&&+&w&=&4\\x&-&y&+&z&+&w&=&1\end{array}}$
6. ${\begin{array}{*{4}{rc}r}x&&&+&z&+&w&=&4\\2x&+&y&&&-&w&=&2\\3x&+&y&+&z&&&=&7\end{array}}$

答案

与上一个问题一样，输入这些方程组很容易。

1. 这个方程组有无穷多个解。在第二小节中，我们给出了解集为

   $\{{\begin{pmatrix}6\\0\end{pmatrix}}+{\begin{pmatrix}-2\\1\end{pmatrix}}y\,{\big |}\,y\in \mathbb {R} \}$

   Maple 的响应结果是 $[6-2\_t_{1},\_t_{1}]$。
2. 解集只有一个元素

   $\{{\begin{pmatrix}0\\1\end{pmatrix}}\}$

   Maple 毫不费力地找到了它 $[0,1]$.
3. 这个系统的解集是无限的

   $\{{\begin{pmatrix}4\\-1\\0\end{pmatrix}}+{\begin{pmatrix}-1\\1\\1\end{pmatrix}}x_{3}\,{\big |}\,x_{3}\in \mathbb {R} \}$

   Maple 给出了 $[\_t_{1},-\_t_{1}+3,-\_t_{1}+4]$.
4. 存在唯一的解

   $\{{\begin{pmatrix}1\\1\\1\end{pmatrix}}\}$

   Maple 给出了 $[1,1,1]$.
5. 该系统有无数个解；在第二小节中，我们用两个参数描述了解集

   $\{{\begin{pmatrix}5/3\\2/3\\0\\0\end{pmatrix}}+{\begin{pmatrix}-1/3\\2/3\\1\\0\end{pmatrix}}z+{\begin{pmatrix}-2/3\\1/3\\0\\1\end{pmatrix}}w\,{\big |}\,z,w\in \mathbb {R} \}$

   Maple 也给出了 $[3-2\_t_{1}+\_t_{2},\_t_{1},\_t_{2},-2+3\_t_{1}-2\_t_{2}]$.
6. 解集为空，Maple 对该命令的回复没有返回解。linsolve(A,u)命令。

问题 4

计算机对一般 $2\!\times \!2$ 系统的解给出什么？

${\begin{array}{*{2}{rc}r}ax&+&cy&=&p\\bx&+&dy&=&q\end{array}}$

答案

响应此提示

> A:=array( [[a,c],
[b,d]] );
> u:=array([p,q]);
> linsolve(A,u);

Maple 可能思考了二十秒，然后给出了以下回复。

${\bigl [}-{\frac {-d\,p+q\,c}{-b\,c+a\,d}},{\frac {-b\,p+a\,q}{-b\,c+a\,d}}{\bigr ]}$

Mathematica 解决方案

RowReduce[( {
    {a, c, p},
    {b, d, q}
   } )] // MatrixForm

返回 ${\begin{pmatrix}1&0&{\frac {dp-cq}{-bc+ad}}\\0&1&{\frac {bp-aq}{bc-ad}}\end{pmatrix}}$，大约在 2 毫秒内。