《Mathematica求解方程(组)、级数》由会员分享，可在线阅读，更多相关《Mathematica求解方程(组)、级数（16页珍藏版）》请在金锄头文库上搜索。
1、方程（组）与级数的 Mathematica求解学习目标1. 能用Mathematica求各种方程（组）的数值解和近似解；2. 能对常见函数进行幕级数的展开。求解简单方程（组）数学里的方程是带有变量的等式。一般地说，一个或一组方程总是对于方程中出现的变量的可能取值范围增加了一些限制。所谓求解方程就是设法把方程对于变量取值的限制弄清 楚，最好的结果是用不含变量的表达式把变量的值表示出来。在这个系统里，方程也用含有变量的等式表示，要注意的是在这里等号用连续的两个等号（=）表示。方程的两端可以是任何数学表达式。用户可以自己操作 Mathematica系统去求解方程，例如使用移项一类的等价变换规则对 方
2、程加以变形、对方程的两端进行整理、把函数作用于方程的两端等等。系统也提供了一些用于求解方程的函数。1、求方程的代数解最基本的方程求解函数是Solve，它可以用于求解方程（主要是多项式方程）或方程组。Solve有两个参数，第一个参数是一个方程，或者是由若干个方程组的表（表示一个方程组）;第二个参数是要求解的变量或变量表。例如，下面的式子对于变量X求解方程432x-x3 _6x21=0 :In 1:=SolvexA4-xA3-6xA2+仁=0,x输入了这个表达式，系统立刻就能计算出方程的四个根，求出的解都是精确解（代数根）。对于一般的多项式，这样得出的解常常是用根式描述的复数。方程的解被表示成一个
3、表，表中是几个子表，每一个子表的形式都是 x-.，箭头后面是方程的一个解。Solve也可以求解多变量的方程或者方程组：In 2:=Solvex-2y=0,xA2-y=i,x,y这个表达式求解方程组：x _ 2y = 0 x2 -y 有时求解方程会得到非常复杂的解。例如将上面的第一个方程稍加变形，所得到的解的表达式就会变得很长：In 3:=SolvexA4-xA3-6xA2=2=0,x这个表达式求出的解的表达式非常长，以至一个计算机屏幕显示不下。使用MS-DOS系统上的Mathematica的读者可以用键盘上的 PgUP键和PgDn键把计算机屏幕上已经卷出 的表达式翻回来阅读，附录B里提供了使用
4、这类计算机的有关操作的更详细的说明。对于使用图形界面提供的功能去翻阅前面的结果。在被求解的方程里还可以有其他符号参数，可以要求系统对于这一个或者那一个变量求解方程。对于Mathematica系统来说，方程中的符号变量（无论使用什么变量名）都是一样 的。对于处理复杂的方程，MATHEMETICA 系统还提供了例外两个有用的函数。函数Eliminate 用于从方程组消去一个或几个变量，例如下面的表达式消去方程组里的变量 Y:IN4:=Elimi nateXA2-2Y= =1,X+2Y= =4,YEliminate 的使用形式与 Solve 类似，它的第二个参数用于说明希望消去的变量。另一 个函数
5、Reduce 用于化简复杂的方程或方程组，它试图用一组比较简单的逻辑关系来描述由 原来方程所描述的变量之间的关系。它的使用形式与 Solve, Eliminate 一样，这里不举例字 了。2、 求方程的数值解理论上已经证明，对于五次以上的多项式方程没有求代数解的一般方法， MATHEMATICA 也求不出那些不能分解因式的五次以上的多项式方程的解，例如：IN5 :=SOLVEXA5+5XA3-2= =0,X它返回一个带有函数 TORUOES 的表达式。可以把函数 N 作用到这个结果表达式上， 求出方程的数值解：IN6:=N%可以看到系统同时求出了方程的五个根的时候可以直接用函数N 和 SOLV
6、E 结合完成工作：IN7:=NSOLVEXA6+4XA2-31= =0,X在系统里直接提供了一个函数 NSOLVE 做这件事。 对于更复杂的方程（或方程组） ，用 SOLVE 求不出根，使用函数 N 也解决不了问题。 对于这样的方程，用户可以使用 REDUCE,ELIMINATE 等函数去处理，设法把方程描述的 变量之间的关系搞清楚。如果需要的就是方程的根，那么只要用求数值根的函数 FINDROOT 。函数 FINDROOT 求数值根所采用的方法与人们一般用计算机求数值根的方法 一样。但是，由于 MA THEMA TIC 有求导函数的能力，在这里计算有导函数的表达式的数 值根就非常简单。 不管
7、表达式多么复杂， 系统都能自动的求出它的导函数。 求数值根使用的 也是牛顿法，用户必须给 FINDROOT 提供一个初始值。下面一个简单的例子：IN8:=FINDROOTSINXEXP2X-COSX= =0,X,0.5对于求不出导函数的表达式，例如用户自己定义的一个复杂计算函数，使用FindRoot提供函数值取不同符号（正负号）的两个点（用表的形式放在上面初始值0.5 的位置），形式是：IN9:= FindRoot FUN1X= =0,X,0,1这里假使 FUN1 是用户定义的一个函数。 使用计算机求数值根的第一个问题是确定初始点， 若初始值选取得不好将给求根带来困 难。再一个麻烦是用户要自己
8、求出函数的导函数。 在 MATHEMATICA 系统里处理的对象是 表达式，一个表达式可以服务于不同的用途， 可以作为求值的对象，作为画图的对象， 也可 以作为演算的对象。当需要求一个表达式的数值根的时候，表达式的这样的多种功 能，或者说 MATHEMATICA 系统对于表达式的多方面的操作能力就表现出很大的优越性。 一个代数表达式，无论多么复杂， MATHEMA TICA 系统都可以直接求出它的导函数，可以 作出它的图形。 从图形上我们很容易认识这个函数表达式在某一个区间的大致性质， 包括它 的根的出现和分布情况。 对表达式的这些认识为人们确定如何取初始值、 如何求根提供了很 有价值的线索。
9、这样，某些比较难以处理的问题可能就容易解决了。二、 求解常微分方程（组）1、常微分方程（组）的精确解Mathematica 能求常微分方程（组）的准确解，能求解的类型大致覆盖了人工求解的范围，功能很强。但不如人灵活(例如在隐函数和隐方程的处理方面)，输出的结果与教材上的答案可能在形式上不同。 另外，Mathematica求数值解也很方便， 且有利于作出解的图形。 求准确解的函数调用格式如下：DSolveeqn, yx ，x求方程eqn的通解y (x)，其中自变量是 x。DSolveeqn ，yx 0= =yo，yx ，x求满足初始条件 y(xo) = yo 的特解 y (x)。DSolveeq
10、n1 , eqn2，, y 1x , y2x， , x求方程组的通解。DSolveequ1 ，y1x 0= =y 10, , y 1x , y2x， , x求方程组的特解。说明：应当特别注意，方程及各项参数的表述方式很严格，容易出现输入错误。微分方 程的表示法只有通过例题才能说清楚。例1解下列常微分方程(组)(1)y = 2y. (x 1)2, (2) y =y；x+1(x + x )y(3)y = zz = _y(4)y = zy的通解及满足初始条件y (0) =0 , zZ = y(0) =1的特解。解: ln1=DSolvey x= =2yx/(x+1) + (x+1)A (5/2 ),
11、yx,xOut1=r 2WTj+x)7/2+(1+x5ln2=DSolvey 凶=(1+yxF2 )/(x+xA3 )Out2=yx 匚T yx ln3=DSolvey x= =zxz，X= = -yxyx,zx , xyx) , yx , x- 1 - 12c1 1 2 xOut3=yxt C1Cosx+ C2Si nx,zxT C2Cosx- C1Si nxln4=DSolvey x= =zx, z x= = -yx,y0= =0 , z0= =1,yx,zx , xOut4=yxt Sinx , zx t Cosx提示：认真观察上例，可以从中学习输入格式，未知函数总带有自变量，等号用连续
12、键入两个等号表示，这两点由于不习惯会出错！导数符号用键盘上的撇号，连续两撇表示二阶 导数，这与习惯相同。自变量、未知量、初始值的表示法与普通变量相同。说明：输出结果总是尽量用显式解表出， 有时反而会使表达式变得复杂， 这与教科书的 习惯不同。当求显式解遇到问题时，会给出提示。通解中的任意常数用C1 , C2，表示。例2求解下列微分方程：门)y 3y 3y y =(x -5)e(2)x2 (y)2 =1, (3) . y 二 xy。(x - 5 ) Exp-x,解：In1 : =DSolve y x +3y x+3y x + yx=yxOut1=,x2、-5x + x2丿5x22yx； 1 ex
13、22竺 xeC1 exC2 ex2C3.3 4xIn 2: =Simplify%1Out2= yxe(-20x3x4 24C1 24xC2 24x2C3) 24In3: =DSolvexA2 + y xA2 = = 1, yx , xxxOut3=In 4yx1 - x221y x x2Out4=ArcSi nx2C1,ArcSi nxC12=DSolveSqrty x = = x yx , yx , xyx -33x -C1xx说明：由以上可以看出对方程的类型并无限制,但是输出的答案未必符合习惯，例如第xx一个方程的答案需要化简，有时即使化简后也未必与教材上的答案一致。例3 求微分方程xy + y - e x = 0在初始条件y| x=1 = 2e下的特解。解：In1 : =DSolvex*y x+yx-EAx= =0, y1= =2E , yx , xOut1= yxx例5求常微分方程组:2、常微分方程（组）的数值解函数NDSolve用于求给定初值条件或边界条件的常微分方程（组）的近似解，其调用格式如下：NDSolveeqns ，yi, y2, , x , xmin , xmax 求常微分方程（组）的近似解。其中微分方程
