向量行列式矩阵与线性方程组实验-Mathematica.docVIP

向量、行列式、矩阵与线性方程组实验 [学习目标] 会用Mathematica进行向量的计算； 能用Mathematica进行行列式的计算； 会利用Mathematica进行矩阵的运算与初等变换； 能利用Mathematica解线性方程组。 线性代数的数值计算程序并不稀奇，早有大量的算法和软件。然而这里是进行准确的符号运算，学习了本节以后，就可以摆脱冗繁的矩阵运算了。本节介绍用Mathematica实现线性代数运算的各种专用函数，它们基本上满足了线性代数计算的需求。读者将会看到，以下的一些计算功能是十分出色的。但从我国的教材来看，还有个别计算功能没有涉及，留有继续开发的余地。 矩阵的输入与输出 在Mathematica中向量和矩阵就是一个表。 {a1，a2，…，an} 表示一个向量。 {{a11，a12，…，a1n}，{a21，a22，…，a2n}，…，{am1，am2，…，amn}} 表示一个m行n列的矩阵，其中每一个子表表示矩阵的一行。 直接输入矩阵 直接输入矩阵的方法有3种，如下所述。 按表的形式输入矩阵 既然矩阵和向量都是表，表的一般操作对于矩阵和向量仍然适用。但是，按表的格式键入矩阵和向量，会让人很不习惯。因此，Mathematica也提供了矩阵和向量的常规形式的输入、输出方法。 由模板输入矩阵 基本输入模板中有输入2阶方阵的模板，单击该模板输入一个空白的2阶方阵。按“Ctrl +”使矩阵增加一列，按“Ctrl + Enter”使矩阵增加一行。如果矩阵不大，此法较方便。 由菜单输入矩阵 如果输入行、列数较多的矩阵，可以打开主菜单的Input项，其中Create Table/Matrix/Palette可用于建立一个矩阵，单击该项出现一个的对话框。选择Make：Matrix，再输入行数和列数，单击OK按钮，于是一个空白矩阵被输入到工作区窗口。 空白矩阵的每个小方块代表一个元素的位置，光标所在的小方块与众不同，可以用Tab键将光标从一个方块跳到下一个方块，也可以用鼠标选中一个方块。 以矩阵形式输出矩阵 不管用何种方法输入矩阵，矩阵总是按表的形式输出。这既违背常规，又难于阅读。因此，Mathematica提供了以矩阵形式输出矩阵的函数： MatrixForm[list] 将表list按矩阵的形式输出。 观察下面矩阵的输出。 In[1]：= a={{1，2，3}，{4，5，6}} Out[1]= {{1，2，3}，{4，5，6}} In[2]：=MatrixForm[a] Out[2] / / MatrixForm= In[3]：= a={{1，2，3}，{4，5，6}} / / MatrixForm Out[3] / / MatrixForm= In[4]：= Out[4]={{1，2}，{3，4}，{5，6}} In[5]：= % // MatrixForm Out[5] / / MatrixForm= 说明：由上例可以看出，不管输入的形式是否为矩阵，必须使用MatrixForm才能使输出为矩阵形式。这不合习惯且费事，解决的方法是自制一个模板：// MatrixForm，以便快速输入。 使用函数MatrixForm又会出现另一个问题，可以通过以下例子来说明。 观察下面矩阵的输出。 In[1]：= a= Out[1]={{1，2}，{3，4}} In[2]：= b= // MatrixForm Out[2] // MatrixForm= In[3]：=Inverse[a] / / MatrixForm Out[3] / / MatrixForm= In[4]：=Inverse[b] // MatrixForm Out[4] / / MatrixForm= Inverse 说明：以上In[3]和In[4]是求逆矩阵，Mathematica求出a的逆矩阵，对b却失败！变量a形式上是表，但能被Mathematica作为矩阵处理。而变量b虽然表示常规形式的矩阵，但不能对b进行各种矩阵计算，务必注意。 技巧：使用括号能够改变表达式的含义，解决上述问题。 观察下面矩阵的输出 In[1]：=// MatrixForm