完整版)Mathematica求解方程(组)、级数

1、方程（组）与级数的 Mathematica求解学习目标1. 能用Mathematica求各种方程（组）的数值解和近似解；2. 能对常见函数进行幕级数的展开。求解简单方程（组）数学里的方程是带有变量的等式。一般地说，一个或一组方程总是对于方程中出现的变量的可能取值范围增加了一些限制。所谓求解方程就是设法把方程对于变量取值的限制弄清 楚，最好的结果是用不含变量的表达式把变量的值表示出来。在这个系统里，方程也用含有变量的等式表示，要注意的是在这里等号用连续的两个等号（=）表示。方程的两端可以是任何数学表达式。用户可以自己操作 Mathematica系统去求解方程，例如使用移项一类的等价变换规则对 方

2、程加以变形、对方程的两端进行整理、把函数作用于方程的两端等等。系统也提供了一些用于求解方程的函数。1、求方程的代数解最基本的方程求解函数是Solve，它可以用于求解方程（主要是多项式方程）或方程组。Solve有两个参数，第一个参数是一个方程，或者是由若干个方程组的表（表示一个方程组）;第二个参数是要求解的变量或变量表。例如，下面的式子对于变量X求解方程432X4 X3 6x210:In 1:=SolvexA4-xA3-6xA2+仁=0,x输入了这个表达式，系统立刻就能计算出方程的四个根，求出的解都是精确解（代数根）。对于一般的多项式，这样得出的解常常是用根式描述的复数。方程的解被表示成一个表，

3、表中是几个子表，每一个子表的形式都是x-.，箭头后面是方程的一个解。Solve也可以求解多变量的方程或者方程组：In 2:=Solvex-2y=0,xA2-y=1,x,y这个表达式求解方程组x 2y 0 x2 y 1.有时求解方程会得到非常复杂的解。例如将上面的第一个方程稍加变形，所得到的解的表达式就会变得很长：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 类似，它的第二个参数用于说明希望消去的变量。另一 个函数 Reduce

5、用于化简复杂的方程或方程组，它试图用一组比较简单的逻辑关系来描述由 原来方程所描述的变量之间的关系。它的使用形式与 Solve, Eliminate 一样，这里不举例字 了。2、 求方程的数值解理论上已经证明，对于五次以上的多项式方程没有求代数解的一般方法， MATHEMATICA 也求不出那些不能分解因式的五次以上的多项式方程的解，例如：IN5 :=SOLVEXA5+5XA3-2= =0,X它返回一个带有函数 TORUOES 的表达式。可以把函数 N 作用到这个结果表达式上， 求出方程的数值解：IN6:=N%可以看到系统同时求出了方程的五个根的时候可以直接用函数N 和 SOLVE 结合完成工

6、作：IN7:=NSOLVEXA6+4XA2-31= =0,X在系统里直接提供了一个函数 NSOLVE 做这件事。对于更复杂的方程（或方程组） ，用 SOLVE 求不出根，使用函数 N 也解决不了问题。 对于这样的方程，用户可以使用 REDUCE,ELIMINATE 等函数去处理，设法把方程描述的 变量之间的关系搞清楚。如果需要的就是方程的根，那么只要用求数值根的函数 FINDROOT 。函数 FINDROOT 求数值根所采用的方法与人们一般用计算机求数值根的方法 一样。但是，由于 MATHEMATIC 有求导函数的能力，在这里计算有导函数的表达式的数 值根就非常简单。 不管表达式多么复杂， 系

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)。DSolveeqn1 ,

10、eqn2，, y 1x , y2x， , x求方程组的通解。DSolveequ1 ，y1x o= =y 10, , y 1x , y2x， , x 求方程组的特解。说明：应当特别注意，方程及各项参数的表述方式很严格，容易出现输入错误。微分方 程的表示法只有通过例题才能说清楚。例1解下列常微分方程(组)(1) y2yx 15(x 1护,(2) y1 y2(x x3)y(3)y z(4)y的通解及满足初始条件y (0) =0, z (0) =1的特解。z y解：In1 : =DSolvey x= =2yx/(x+1) + (x+1) a (5/2 ), X-Out1=2/7x02 - 31ln2=

11、DSolvey凶=(1+yxA2 )/(x+xA3xOut2=yx112c11xxr， yxyx) , yx,ln3=DSolvey x= =zx, z x= = -yxyx,zx , 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求解下列微分方程:(1) y3y 3y y (x 5)e x, (2) x2 (y )2xy。解：ln1:=DSolve y x +3yx +3y X + yx=(x - 5 ) Exp-x,yx,xOut1=yx 丄e xx225xx 5x2 e x - 23x_3ln2Out2=ln3Out3=ln4Out4=1 x 5x3 e2 3=Simplify%yx24e=

13、DSolvexA2 + yC1e xxC 2 ex2C3x(20x3x4xA2 = = 124C1 24xC2 24x2C3),yx , x2 ArcSi nx x、12 2 yx1x2如舸凶2 2=DSolveSqrty x = = x yxyxC1,C1,yx , xyx J x3C1说明：由以上可以看出对方程的类型并无限制,但是输出的答案未必符合习惯，例如第x一个方程的答案需要化简，有时即使化简后也未必与教材上的答案一致。例3 求微分方程xy + y - e = 0在初始条件y| x=1 = 2e下的特解。解：In1 : =DSolvex*y x+yx-EAx= =0, y1= =2E

14、, yx , xOut1= yxIn2的结果是用y表示解函数的名字，因此In3顺利画出解曲线如图13-43所示。2、常微分方程（组）的数值解函数NDSolve用于求给定初值条件或边界条件的常微分方程 （组）的近似解，其调用格 式如下：NDSolveeqns ，yi, y2, , x , xmin , xmax求常微分方程(组)的近似解。其中微分方程和初值条件的表示法如同DSolve，未知函数仍有带自变量和不带自变量两种形式，通常使用后一种更方便。初值点xo可以取在区间xmin , xmax上的任何一点处,得到插值函数InterpolatingFunctiondomain, table类型的近似

15、解，近似解的定义域domain 般为domain , table，也有可能缩小。例4 求常微分方程y = x2 + y2,满足初始条件y (0) = 0的数值解。解：In1 : =s1=NDSolvey x= =xA2+yxA2, y0= =0,y, x , -2 , 2Out1=yt InterpolatingFunction-2., 2. , In2: = y=y / . s11Out2=lnterpolatingFunction-2., 2., In3: =Plotyx , x , -2 , 2 , AspectRatio t Automatic ,PlotRa nge图13-43微分方

16、程的解曲线Out3= -Graphics-上例中包含许多值得学习的实用内容，其中第二项参数使用y而不是yx,比用yx好。如果求解区间改为x , -3 , 3,就会出现警告提示，实际得不到-3 , 3上的解。Out1表明返回的解放在一个表中，不便使用，实际的解就是插值函数：,2., In terpolati ngFunction -2.例5求常微分方程组:13x yxx3yx满足初始条件x（ 0）=0，y（ 0）=1的数值解。解：In1 : =s仁NDSolvex t= = yt -(xtA3/3 - xt ),t= = - xt,x0= =0 , y0= =1,x, y , t , -15 ,

17、15Out1=xt In terpolat ingFunction -15.,15., ,t in terpolat ingFun cti on -15.15. , In 2: = x=x / . s11, 1y=y / . s11, 2Out2=I nterpolatin gFu nctio n-15.Out3=I nterpolat in gFu nctio n-15.,15., ,15., In4=ParametricPlotxt,yt, t , -15, 15,AspectRatio图13-44 解的相轨线Out3= -Graphics-说明：上例是求一个著名方程组的近似解，其中In2

18、也可以改用一个赋值式x , y=x ,y / . Flattens1, 一次得到两个函数。通过求数值解容易得到它的相图，In4绘制了解的相轨线如图13-44所示，图中表明原点是奇点，极限环的形状也已经得到。为了应付复杂的情况，需要设置可选参数:Worki ngPrecisi on参见数值积分部分的介绍。AccuracyGoal计算结果的绝对误差。Precisi on Goal计算结果的相对误差。MaxSteps最大步数。Starti ngStepSize初始步长。以上可选参数的默认值都为Automatic ,其中 AccuracyGoal 和 Precisi on Goal 的默认值比 Wor

19、kingPrecision 小10,当解趋于 0时应将 AccuracyGoal取成Infinity 。对于常微 分方程，最大步长默认值为 1000。这个函数也可以解偏微分方程，最大步长默认值为200。例6解下列微分方程（组）：1（1）y |，满足初始条件y （0）=1的特解；4yx 3x 3y（2）y xz 26.5x y，满足初始条件x（0）=z（0）=0，y（0）=1的特解。z xy z解：In1 : =NDSolvey x= =l/4yx, y0= =1, y, x , 1,AccuracyGoal20, PrecisionGoal 20, WorkingPrecision 25Out

20、1=yt In terpolat in gFu nctio n0 , 1.000000000000000000000000000 , In 2: =y1 / . %Out2=0.968912424710644784118519+ 0.2474039592545229296234109ii x,y , z, t , 0 ,20, MaxStepst 3000Out3=xt In terpolati ngFu nctio n0.,20., ,in4t in terpolat ingFun cti on 0.t in terpolat ingFun cti on 0.=ParametricPlot3

21、DEvaluatext,20., ,20. , ,yt , zt / . %in 3:=NDSolvex t= = -3（xt -yt）,yt = = -xt zt+36.5xt -yt,zt = = xt yt- zt,x0 = = z0 = = 0,y0= =1t, 0, 20 , PlotPointst 1000图13-45 3维相轨线Out3= -Graphics3D-说明：以上范例中In1取高精度，而且是复系数方程。In2是求解在x=1时的近似1 1ii值，求精确解能得到准确值e4，读者可以求e4的近似值与 Out2的结果比较，验证近似解的精确度确实很高。In3在求解时增大步数，成功

22、地得到了由In4绘制的如图13-45所示的解的相轨线。In4所示的绘图语句与前面例子中的不同，现在只要会模仿使用它们 就行了，要想弄清原理请参阅相关Mathematica书籍。三、级数1. 求和与求积求有限或无穷和、积的函数是：i maxSumf , i , imin , imax求 f(i)，其中 imin 可以是-汽i i minimax可以是g (即+8),但是必须满足imin imax。基本输入模板中也有求和专用的符号， 使用模板输入更方便。Sumf , i , imin , imax , j , jmin , jmax，求多重和，也可以使用基本输入模板连续多次输入求和符号得到。i m

23、axProductf , i , imin , imax求 f (i)，基本输入模板中i i min也有求积符号。Productf , i , imin , imax , j , jmin , jmax，求多重积，也可以使用基本输入模板连续多次输入求积符号得到。NSum 和NProduct得到数值解。2. 将函数展开为幕级数Out1=3xx 一32x51517x7 62x910ox3152835ln2:=SeriesSi nx /x,x , 0, 9Out2=24“ x x16 1206x50408x10ox 362880将函数展开为幕级数的函数调用格式如下:Seriesf,x , X0, n将函数f (x)在X0处展成幂级数直到 n次项为止。Seriesf,x , X。，n, y ,y。，m将函数f(x, y)先对y后对x展开。例7展开下列函数为幂级数(1)y=tgx, (2)ysin xx(3) y =xyf (x), (4) y = e 。解：In1 : =SeriesTanx , x , 0, 9In 3:=Seriesfx , x , 1, 7Out3=f1 f 1(x