第三章Maple代數.docVIP

第三章 线性代数 教学目的：。 教学：。 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 多项式相除的余式(rem)/商式(quo) 计算p1除以p2的余式, 将商式赋值给q的命令格式为：rem(p1, p2, x, q); 计算p1除以p2的商式, 将余式赋值给r的命令格式为：quo(p1, p2, x, r); 余式和商式满足：p1=p2*q+r, 其中degree(r, x)degree(p2, x) rem(x^5+x^3+x, x+1, x,