Maple教程 - 第4章 - 方程求解.pdfVIP

- 100 - 第四章 方程求解 1 代数方程(组)求解 1.1 常用求解工具—solve 求解代数方程或代数方程组, 使用Maple中的 solve函数. 求解关于 x的方程 eqn=0 的命令格式为： solve(eqn, x); 求解关于变量组 vars的方程组 eqns的命令为： solve(eqns, vars); eqn:=(x^2+x+2)*(x-1); := eqn ( ) + + x2 x 2 ( ) ? x 1 solve(eqn,x); , ,1 ? + 12 1 2 I 7 ? ? 1 2 1 2 I 7 当然, solve也可以求解含有未知参数的方程： eqn:=2*x^2-5*a*x=1; := eqn = ? 2 x2 5 a x 1 solve(eqn,x); , + 54 a 1 4 + 25 a 2 8 ? 54 a 1 4 + 25 a 2 8 solve 函数的第一个参数是有待求解的方程或方程的集合, 当然也可以是单个表达 式或者表达式的集合, 如下例： solve(a+ln(x-3)-ln(x),x); 3 e a ? + 1 e a 对于第二个参数, Maple的标准形式是未知变量或者变量集合, 当其被省略时, 函数 indets自动获取未知变量. 但当方程中含有参数时, 则会出现一些意想不到的情况： solve(a+ln(x-3)-ln(x)); - 101 - { }, = x x = a ? + ( )ln ? x 3 ( )ln x 很多情况下, 我们知道一类方程或方程组有解, 但却没有解决这类方程的一般解法, 或者说没有解析解. 比如, 一般的五次或五次以上的多项式, 其解不能写成解析表达式. Maple 具备用所有一般算法尝试所遇到的问题, 在找不到解的时候, Maple 会用 RootOf 给出形式解. x^7-2*x^6-4*x^5-x^3+x^2+6*x+4; ? ? ? + + + x7 2 x6 4 x5 x3 x2 6 x 4 solve(%); + 1 5 ? 1 5 ( )RootOf , ? ? _Z5 _Z 1 = index 1 ( )RootOf , ? ? _Z5 _Z 1 = index 2 ( )RootOf , ? ? _Z5 _Z 1 = index 3, , , , , ( )RootOf , ? ? _Z5 _Z 1 = index 4 ( )RootOf , ? ? _Z5 _Z 1 = index 5, solve(cos(x)=x,x); ( )RootOf ? _Z ( )cos _Z 对于方程组解的个数可用 nops命令获得, 如： eqns:={seq(x[i]^2=x[i],i=1..7)}; := eqns { }, , , , , , = x1 2 x1 = x2 2 x2 = x3 2 x3 = x4 2 x4 = x5 2 x5 = x6 2 x6 = x7 2 x7 nops({solve(eqns)}); 128 但是, 有时候, Maple甚至对一些“显而易见”的结果置之不理, 如： solve(sin(x)=3*x/Pi,x); ( )RootOf ? 3 _Z ( )sin _Z π 此方程的解为 0 , 6 π ± , 但 Maple 却对这个超越方程无能为力, 即便使用 allvalues 求解也只有下述结果： allvalues(%); ( )RootOf , ? 3 _Z ( )sin _Z π 0. 另外一个问题是, Maple 在求解方程之前,会对所有的方程或表达式进行化简, 而不 管表达式的类型, 由此而产生一些低级的错误： (x-1)^2/(x^2-1); ( ) ? x 1 2 ? x2 1 solve(%); 1 - 102 - 但是, 大量实验表明, solve 的确是一个实用的方程求解工具, 但是也不可盲目相信 它给出的一切结果, 特别是对于非线性方程而言, 对于给出的结果需要加以验证. 下面通过几个例子说明在Maple中非线性方程组的求解问题. 例：求解方程组： ? ? ? =? =+ yx yx 9 25 2 22 eqns:={x^2+y^2=25,y=x^2-5}; := eqns { }, = y ? x2 5 = + x2 y2 25 vars:={x,y}; :=