第三章 Maple代数6.docVIP

PAGE 71 -第三章 线性代数教学目的：学习并掌握计算机代数系统Maple下进行多项式的各种运算，包括合并、分解因式、排序等，学习并掌握行列式、矩阵的各种符号计算，包括特征矩阵、特征值、特征向量等，学习并掌握线性方程组求解的符号计算方法，了解线性空间的有关内容。教学目标：掌握多项式、行列式、矩阵的各种运算，求解线性方程组。重点内容：多项式、行列式、矩阵的各种运算，线性方程组求解。难点内容：多项式的符号计算，矩阵的运算。1多项式多项式是代数学中最基本的对象之一, 它不但与高次方程的讨论有关, 而且是进一步学习代数以及其它数学分支的基础. 1.1 多项式生成及类型测试在Maple中, 多项式由名称、整数和其他Maple值, 通过+、-、*和^等组合而成. 例如： p1:=5*x^5+3*x^3+x+168;这是一个整系数单变量多项式. 多元多项式和定义在其他数域上的多项式可以类似构造： p2:=3*x*y^2*z^3+2*sqrt(-1)*x^2*y*z+2002;由此可以看出, Maple中多项式的生成与“赋值”命令相似. 另外, 还可以通过函数randpoly生成随机多项式, 生成一个关于vars的随机多项式的格式如下：randpoly(vars, opts);其中, vars表示变量或者变量列表或集合, opts为可选项方程或者指定属性的名称. 如： randpoly(x); ＃随机生成关于x的5次(默认)多项式 randpoly([x, y], terms=8); ＃随机生成关于[x, y]二元8项多项式 randpoly([x, sin(x), cos(x)]);而要随机生成关于[x, y, z]的密集的、均匀的、度为2的多项式的命令为： randpoly([x,y,z],dense,homogeneous,degree=2);用type命令可以测试多项式的类型： type(p1, polynom(integer, x)); #测试p1是否是一个关于x的整系数多项式 type(p2, polynom(complex, {x, y, z})); #测试p2是否是一个关于{x, y, z}的复系数多项式1.2 提取多项式系数coeff函数用来提取一元多项式的系数, 而多元多项式所有系数的提取用命令coeffs, 指定系数的提取用命令coftayl. (1) 提取多项式p中x^n的系数使用命令：coeff(p, x^n);或coeff(p, x, n); (2) 提取多项式p中变量x的所有系数并将相应的x幂存于变量t中：coeffs(p, x, ’t’); (3) 返回expr在x=a处的Taylor展式中(x-a)^k的系数: coeftayl(expr, x=a, k); p:=2*x^2+3*y^3*x-5*x+68; coeff(p, x); coeff(x^4-5*x^2-sin(a)*(x+1)^2, x^2); s:=3*x^2*y^2+5*x*y; coeffs(s); coeffs(s, x, t); t; coeftayl(exp(x), x=0, 10); p:=3*(x+1)^3+sin(Pi/3)*x^2*y+x*y^3+x-6; coeftayl(p, x=-1, 1); coeftayl(p, [x, y]=[0, 0], [1, 0]);返回默认为降序排列的多元多项式的首项和末项系数分别使用命令lcoeff、tcoeff： lcoeff(p, x); tcoeff(p, x);1.3 多项式的约数和根1.3.1多项式的最大公约因式(gcd)/最小公倍因式(lcm)求多项式的最大公约因式/最小公倍因式的命令与求两个整数最大公约数/最小公倍数命令一样, 都是gcd/lcm. 命令格式分别为：gcd(p1, p2, t, s); lcm(p1, p2, t, s);其中, 第3个参数t赋值为余因子p1/gcd(p1, p2), 第4个参数s赋值为余因子p2/gcd(p1, p2). p1:=x^4+x^3+2*x^2+x+1; p2:=x^2+x+1; gcd(p1, p2, t, s); t, s; lcm(p1, p2);1.3.2多项式的平方根(psqrt)和第n次方根(proot)求多项式p的平方根, 若不是完全平方, 则返回_NOSQRT：psqrt(p); 求多项式p的n次方根, 若不是完全n次方, 则返回_NOROOT：proot(p, n); p:=x^4+4*x^3+6*x^2+4*x+1; psqrt(p); proot(p, 4); proot(p, 8);1.3.3 多项