《mathematica在高等代数与微积分中的应用》由会员分享，可在线阅读，更多相关《mathematica在高等代数与微积分中的应用（27页珍藏版）》请在金锄头文库上搜索。
1、MATHEMATICA在高等代数与微积分中的应用1 高等代数运算1.1 矩阵的输入、表输入：例：输入矩阵 命令：A=1,2,3,4,5,6,7,8,9不过，我们看到输出的结果不是矩阵形式，如果希望得到矩阵形式，可再使用函数MatrixForm，如： 或者：、二阶方阵可直接用模板输入单击输入面板上的“ ”，再输入矩阵的元素即可，例如，求矩阵的逆：求矩阵逆的函数是：Inverse，或：两种不同的使用方式或：、菜单来输入操作：“输入”“创建表单矩阵面板T” 对话框选择“矩阵” 输入行数和列数 空白矩阵计算结果如下图示：例：、增加行与列按Ctrl+ Shift +“，”； 增加行，Ctrl+“”增加列
2、。、输入任意矩阵例：输入任意矩阵，可用命令：Arraya,2,2 / MatrixFormArraya，m，n 创建m行、n列的矩阵，元素为ai，j 、创建一个n阶单位矩阵：IdentityMatrixn、创建一个对角线上为表list的元素的方阵：DiagonalMatrix list 例： a1=1,2,3,4,5DiagonalMatrixa1 / MatrixForm1.2 MATHEMATICA的矩阵运算命令(1) a=a1,a2,an功能：定义一个一维向量（），这里是数或字母(2) a=Tablefj，j,n 例： (3) a=a11,a12,a1n,a21,a22,a2n,am1,
3、am2,amn功能: 定义一个矩阵: 例：(4) a=Tablefi,j，i，m，j，n功能: 定义一个分量可以用fi,j计算的矩阵，其中f是关于i和j的函数，给出矩阵在第i行第j列的元素值例：(5) MatrixForma 功能：把a按通常的矩阵或向量形式输出，其中a是矩阵或向量(6) DiagonalMatrixlist功能：使用列表中list的元素生成一个对角矩阵.例：(7) IdentityMatrixn 功能：生成n阶单位阵(8) A+B 功能：求A与B的和, 这里A与B都是矩阵或都是向量(9) A-B 功能：求A与B的差这里A与B都是矩阵或都是向量(10) k*A 功能：求常数k与
4、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).TransposeA 功能：求矩阵A的转置矩阵.(15). InverseA 功能：求矩阵A的逆矩阵(16). MatrixPowerA，n 功能:计算方阵A的n次幂(17). DetA功能:求方阵A的行列式(18) ai, j 功能：取矩阵a的位于第i
5、行，第j列的元素.(19). ai 功能：取矩阵a的第i行的所有元素或取向量a的第i个分量.(20) Transposeaj 功能：取矩阵a的第j列的所有元素.1.3 多项式运算命令 PolynomialGCDf,g 功能：求多项式f、g的最大公因式。例： f=4 x4-2 x3-16 x2+5 x+9;g=2 x3-x2-5 x+4;PolynomialGCDf,g PolynomialQuotientf,g,x 功能: 求g除f的商，x为变量。 PolynomialRemainderf,g,x 功能：求g除f的余式，x为变量。 Lengthq1功能：求表q1中元素的个数。 Expandu功
6、能：Expandf 把分式u的分子展开,分母不变且被看成单项。例： f=4 x4-2 x3-16 x2+5 x+9;g=2 x3-x2-5 x+4;Expandf/g Collectexpr，x，y 将expr表示成x的多项式，再把多项式的每一项系数表示成y的多项式。 PolynomialLCMp1，p2， 求多项式p1，p2，的最小公倍式。 PowerExpandexpr 将(xy)n分解成 xnyn 的形式。 用mathematica进行分式运算Denominatorf 提取分式f的分母Numeratorf 提取分式f的分子ExpandDenominatorf 展开分式f的分母Expand
7、Numeratorf 展开分式f的分子Expandf 把分式f的分子展开,分母不变且被看成单项ExpandAllf 把分式f的分母和分子全部展开ExpandAllf, x 只展开分式f中与x匹配的项Togetherf 把分式f的各项通分后再合并成一项Apartf 把分式f拆分成多个分式的和的形式Apartf, x 对指定的变量x（x以外的变量作为常数）,把分式f拆分成多个分式的和的形式Cancelf 把分式f的分子和分母约分Factorf 把分式f的分母和分子因式分解Factorf 把f因式分解。1.3 应用举例：例A 求多项式的最大公因式及相应的、。在高等代数中有以下结论： 最大公因式设F是
8、一个数域，Fx是F上一元多项式环。 定义1 令f(x)和g(x)是Fx的两个多项式.若是Fx的一个多项式h(x)同时整除 f(x)和g(x),那么h(x)叫做f(x)与g(x)的一个公因式.定义2 设d(x)是多项式f(x)与g(x)的一个公因式。若是d(x)能被f(x)与g(x)的每一个公因式整除，那么d(x)叫做f(x)与g(x)的一个最大公因式. 最大公因式的求法 辗转相除法设f(x)和g(x)是Fx的两个多项式且g(x)0，用g(x)除 f(x)，得到商式q(x)及余式r(x)，如果r(x) 0，那么再以r(x)除g(x), 得商式q(x)及余式r(x )，如果r(x ) 0，再以r(
9、x ) 除r(x)，如此继续下去，因为余式的次数每次降低，所以此过程必在有限次后得到这样一个余式r(x)：它整除前一个余式r(x). 这样我们就得到一串等式： f(x)=g(x) q(x)+ r(x) g(x)= r(x) q(x)+ r(x ) r(x)= r(x )q(x)+r(x) （1）r(x)=r(x) q(x)+ r(x)r(x)= r(x) q(x)+ r(x)r(x)= r(x)q(x)则r(x)就是f(x)与g(x)的一个最大公因式. 由(1)的倒数第二个等式得： 令： ， （2） 则： （3）由(1)的倒数第三个等式得： （4）把（4）代入（3），并令 ， （5）即得： （
10、6）一直下去，最后可得到与，使得： （7） 算法描述：根据上述推导及结论，可以得到以下算法： 输入f(x)与g(x)； 辗转相除： 定义数组q1：存储每次带余除法所得的商式。 定义数组r1：存储每次带余除法所得的余式。 While(1) f(x)=g(x) q(x)+ r(x); if(r(x)=0) Break; 添加q(x)到q1，添加r(x)到r1; f(x)=g(x); g(x)=r(x); r1中的最后一个元素就是所求最大公因式. 求与 k= q1中元素的个数；u(x)=1;v(x)=-q1k;for(i=1,ik,i+) w(x)=u(x); u(x)=v(x); v(x)=w(x)-v(x)*q1k-i; 输出结果； mathematic程序f=4 x4-2 x3-16 x2+5 x+9;g=2 x3-x2-5 x+4;f1=f;g1=g;d=PolynomialGCDf,g; (*求多项式f，g的最大公因式*)q1=;WhileTrue, q=PolynomialQuotientf,g,x; (*求g除f的商,x为变量*) r=PolynomialRemainderf,g,x ;(*求g除f的余式,x为变量*) Ifr=0,Break; AppendToq
