MATHEMATICA在高等代数和微积分中应用.doc

MATHEMATICA在高等代数与微积分中的应用 1 高等代数运算 1.1 矩阵的输入 ①、表输入： 例：输入矩阵 命令：A={{1,2,3},{4,5,6},{7,8,9}} 不过，我们看到输出的结果不是矩阵形式，如果希望得到矩阵形式，可再使用函数MatrixForm，如： 或者： ②、二阶方阵可直接用模板输入——单击输入面板上的“ ”，再输入矩阵的元素即可，例如，求矩阵的逆： 求矩阵逆的函数是：Inverse， 或： 两种不同的使用方式 两种不同的使用方式 或： ③、菜单来输入．操作：“输入”→“创建表单／矩阵／面板［T］…” ? 对话框→选择“矩阵”→ 输入行数和列数→ ? 空白矩阵．计算结果如下图示： 例： ④、增加行与列 按Ctrl+ Shift +“，”； 增加行，Ctrl+“?”增加列。 ⑤、输入任意矩阵 例：输入任意矩阵，可用命令：Array[a,{2,2}] // MatrixForm Array[a，{m，n}] 创建m行、n列的矩阵，元素为a[i，j]． Array[a，{m，n}] 创建m行、n列的矩阵，元素为a[i，j]． ⑦、创建一个对角线上为表list的元素的方阵： DiagonalMatrix[ list ] 例： a1={1,2,3,4,5} DiagonalMatrix[a1 ] // MatrixForm 1.2 MATHEMATICA的矩阵运算命令 (1) a={a1,a2,…,an} 功能：定义一个一维向量（），这里是数或字母． (2) a=Table[f[j]，{j,n}] 例： (3) a={{a11,a12,…,a1n},{a21,a22,…,a2n},…,{am1,am2,…,amn}} 功能: 定义一个矩阵: 例： (4) a=Table[f[i,j]，{i，m}，{j，n}] 功能: 定义一个分量可以用f[i,j]计算的矩阵，其中f是关于i和j的函数，给出矩阵在第i行第j列的元素值． 例： (5) MatrixForm[a] ???? 功能：把a按通常的矩阵或向量形式输出，其中a是矩阵或向量． (6) DiagonalMatrix[list] 功能：使用列表中list的元素生成一个对角矩阵. 例： (7) IdentityMatrix[n] ?? 功能：生成n阶单位阵 (8) A+B ??? 功能：求A与B的和, 这里A与B都是矩阵或都是向量． (9) A-B ?? 功能：求A与B的差．这里A与B都是矩阵或都是向量． (10) ?k*A ??? 功能：求常数k与A的数乘，这里A是矩阵或向量． (11) A.B ??? ?功能：求矩阵A与矩阵B的乘积，注意A与B之间的乘号“.”必须使用数字键盘上的小数点． (12)? a.b ???? 功能：求向量a与向量b的内积，注意a与b之间的乘号“.”必须使用数字键盘上的小数点． (13) ?A.b或b.A ???? 功能：求矩阵A与向量b的乘积，注意A与b之间的乘号“.”必须使用数字键盘上的小数点． (14).?Transpose［A］ ????? 功能：求矩阵A的转置矩阵. (15). Inverse［A］ ??? ??功能：求矩阵A的逆矩阵 (16). MatrixPower[A，n] ???? 功能:计算方阵A的n次幂． (17). Det[A] 功能:求方阵A的行列式 (18) ?a［［i, j］］ ?? ? 功能：取矩阵a的位于第i行，第j列的元素. (19). a［［i］］ ??? ?功能：取矩阵a的第i行的所有元素或取向量a的第i个分量. (20) Transpose［a］［［j］］ ??? ?功能：取矩阵a的第j列的所有元素. 1.3 多项式运算命令 ① PolynomialGCD[f,g] ??? 功能：求多项式f、g的最大公因式。 例： f=4 x^4-2 x^3-16 x^2+5 x+9; g=2 x^3-x^2-5 x+4; PolynomialGCD[f,g] ② PolynomialQuo