Mathematica求解方程组、级数

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

2、方程的两端进行整理、把函数作用于方程的两端等等.系统也提供了一些用于求解方程的函数.1、求方程的代数解最根本的方程求解函数是Solve,它可以用于求解方程主要是多项式方程或方程组.Solve有两个参数,第一个参数是一个方程,或者是由假设干个方程组的表表示一个方程组;第二个参数是要求解的变量或变量表.例如,下面的式子对于变量X求解方程x4 x3 6x2 10:In1:=SolvexA4-xA3-6xA2+1=0,x输入了这个表达式,系统马上就能计算出方程的四个根,求出的解都是精确解代数根.对于一般的多项式,这样得出的解常常是用根式描述的复数.方程的解被表示成一个表,表中是几个子表,每一个子表的形

3、式都是 x->.,箭头后面是方程的一个解.Soke也可以求解多变量的方程或者方程组：In2:=Solvex-2y=0,xA2-y=1,x,y这个表达式求解方程组x 2y 027x y 1.有时求解方程会得到非常复杂的解.例如将上面的第一个方程稍加变形,所得到的解的表达式就会变得很长：In3:=SolvexA4-xA3-6xA2=2=0,x这个表达式求出的解的表达式非常长,以至一个计算机屏幕显示不下.使用MS-DOS系统上的Mathematica的读者可以用键盘上的PgUP键和PgDn键把计算机屏幕上已经卷出的表达式翻回来阅读,附录B里提供了使用这类计算机的有关操作的更详细的说明.对于使用

4、图形界面提供的功能去翻阅前面的结果.在被求解的方程里还可以有其他符号参数,可以要求系统对于这一个或者那一个变量求解方程.对于Mathematica系统来说,方程中的符号变量无论使用什么变量名都是一样 的.对于处理复杂的方程,MATHEMETICA 系统还提供了例外两个有用的函数.函数Eliminate用于从方程组消去一个或几个变量,例如下面的表达式消去方程组里的变量Y:IN4:=EliminateXA2-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结合完成工作：IN7:=NSOLVEXA6+4XA2-31= =0,X在系统里

6、直接提供了一个函数NSOLVE做这件事.对于更复杂的方程或方程组,用SOLVE求不出根,使用函数 N也解决不了问题. 对于这样的方程,用户可以使用REDUCE,ELIMINATE 等函数去处理,设法把方程描述的变量之间的关系搞清楚.如果需要的就是方程的根,那么只要用求数值根的函数 FINDROOT.函数FINDROOT求数值根所采用的方法与人们一般用计算机求数值根的方法 一样.但是,由于 MATHEMATIC 有求导函数的水平,在这里计算有导函数的表达式的数 值根就非常简单.不管表达式多么复杂,系统都能自动的求出它的导函数.求数值根使用的也是牛顿法,用户必须给 FINDROOT提供一个初始值.

7、下面一个简单的例子：IN8:=FINDROOTSINXEXP2X-COSX= =0,X,0.5对于求不出导函数的表达式,例如用户自己定义的一个复杂计算函数,使用FindRoot提供函数值取不同符号正负号的两个点用表的形式放在上面初始值0.5的位置,形式是：IN9:= FindRoot FUN1X= =0,X,0,1这里假使FUN1是用户定义的一个函数.使用计算机求数值根的第一个问题是确定初始点,假设初始值选取得不好将给求根带来困难.再一个麻烦是用户要自己求出函数的导函数.在MATHEMATICA系统里处理的对象是表达式,一个表达式可以效劳于不同的用途,可以作为求值的对象,作为画图的对象,也可以

8、作为演算的对象.当需要求一个表达式的数值根的时候,表达式的这样的多种功 能,或者说MATHEMATICA系统对于表达式的多方面的操作水平就表现出很大的优越性. 一个代数表达式,无论多么复杂,MATHEMA TICA系统都可以直接求出它的导函数,可以作出它的图形.从图形上我们很容易熟悉这个函数表达式在某一个区间的大致性质,包括它的根的出现和分布情况.对表达式的这些熟悉为人们确定如何取初始值、如何求根提供了很有价值的线索.这样,某些比拟难以处理的问题可能就容易解决了.二、求解常微分方程组1、常微分方程组的精确解Mathematica能求常微分方程组的准确解,能求解的类型大致覆盖了人工求解的范围,功

9、能很强.但不如人灵活(例如在隐函数和隐方程的处理方面),输出的结果与教材上的答案可能在形式上不同.另外,Mathematica求数值解也很方便,且有利于作出解的图形.求准确解的函数调用格式如下：DSolveeqn, yx, x求方程eqn的通解y(x),其中自变量是x.DSolveeqn , yxo= =yo, yx , x 求满足初始条件 y(xo)= yo的特解 y(x).DSolveeqn1 , eqn2, y 1x , y2x, x求方程组的通解.DSolveequ1 ,yx o= =y 10, y 1x , y2x, x 求方程组的特解.说明：应当特别注意,方程及各项参数的表述方式很

10、严格,容易出现输入错误.微分方程的表示法只有通过例题才能说清楚.例1解以下常微分方程(组)：,八 2V1 vy z y(x 1)2, y -,(3)x 1(x x )vz yy z ,(4)7的通解及满足初始条件 y (0) =0, z (0) =1的特解.z y解：In1 : =DSolvey ' x= =2yx/(x+1) + (x+1) a (5/2),.X,Out1=2132/7)1d2)In2=DSolveyx= =11+yxA2 )/(x+xA3yx) , yx,xOut2=yx112 C11 yx112 c1-x1 ；In3=DSolvey ' x= =zx, z

11、 ' x= = -yxyx,zx , xOut3=yxC1Cosx+ C2Sinx ,zx- C2Cosx- C1SinxIn4: =DSolvey ' x= =zx, z' x= = -yx, y0= =0 , z0= =1,yx, zx , xOut4=yxSinx , zx - Cosx提示：认真观察上例,可以从中学习输入格式,未知函数总带有自变量,等号用连续键入两个等号表示,这两点由于不习惯会出错！导数符号用键盘上的撇号,连续两撇表示二阶导数,这与习惯相同.自变量、未知量、初始值的表示法与普通变量相同.说明：输出结果总是尽量用显式解表出,有时反而会使表达式变得复

12、杂,这与教科书的习惯不同.当求显式解遇到问题时,会给出提示.通解中的任意常数用C1 , C2,表示.例2求解以下微分方程：x22(1) y 3y 3y y (x 5)e , (2) x (y)1 , (3) Jyxy.(x - 5 ) Exp-x,解：In1 : =DSolve y x+3y x +3y ' x + yx=yx, xOut1=In2Out2=In3Out3=In41x2yx 2e x341 x 5x x-e 2342 x 5x2x5x2e xC1 e xxC2 e xx2C3=Simplify%1yx金 ex( 20x3x424C1 24xC2=DSolvexA2 +

13、y ' xA2 = = 1, yx , x1:2 ArcSinxyx 厂1 x2C1,12 ArcSinxyx -x 1 xC1=DSolveSqrty ' x = = x yx , yx , xOut4= yxx3 C1»24x2C3) 说明：由以上可以看出对方程的类型并无限制,但是输出的答案未必符合习惯,例如第一个方程的答案需要化简,有时即使化简后也未必与教材上的答案一致.例3求微分方程xy + y - e x = 0在初始条件y| x=1 = 2e下的特解.一 一 '一解：In1 : =DSolvex*y x+yx-EAx= =0, y1= =2E ,

14、yx , xOut1= yx2、常微分方程组的数值解函数NDSolve用于求给定初值条件或边界条件的常微分方程组的近似解,其调用格式如下：NDSolveeqns , y 1, y2, x , xmin , xmax求常微分方程组的近似解.其中微分方程和初值条件的表示法如同DSolve,未知函数仍有带自变量和不带自变量两种形式,通常使用后一种更方便.初值点xo可以取在区间xmin , xmax上的任何一点处,得到插值函数 InterpolatingFunctiondomain , table类型的近似解,近似解的定义域 domain 一般为domain , table,也有可能缩小.例4求常微分

15、方程y' = x2+ y2,满足初始条件y 0 = 0的数值解.解：In1 : =s1=NDSolvey ' x= =xA2+yxA2, y0= =0,y, x , -2 , 2Out1=y f InterpolatingFunction-2. , 2. , < >In2: = y=y / . s11Out2=InterpolatingFunction-2., 2., < >f -1.5 , 1.5In3: =Plotyx , x , -2 , 2 , AspectRatio fAutomatic ,PlotRange图13-43微分方程的解曲线Out3

16、= -Graphics-上例中包含许多值得学习的实用内容,其中第二项参数使用y而不是yx,比用yx好.如果求解区间改为x, -3, 3,就会出现警告提示,实际得不到-3 , 3上的解.Out1,2., < >说明返回的解放在一个表中,不便使用,实际的解就是插值函数：InterpolatingFunction-2.In2的结果是用y表示解函数的名字,因此In3顺利画出解曲线如图13-43所示.例5求常微分方程组:1 3x y - x x3y x满足初始条件x (0) =0, y (0) =1的数值解.解：In1 : =s1=NDSolvex ' t= = yt -(xtA3/

17、3 - xt),y' t= = - xt , x0= =0 , y0= =1,x,y , t , -15 ,15Out1=xfInterpolatingFunction-15.fInterpolatingFunction-15.,15., < >,15. , < >In2: = x=x / . s11, 1y=y / . s11, 2Out2=InterpolatingFunction-15.Out3=InterpolatingFunction-15.,15., < >,15., < >In4=ParametricPlotxt,yt ,

18、t , -15, 15,AspectRatiofAutomatic图13-44解的相轨线Out3= -Graphics-说明：上例是求一个著名方程组的近似解,其中In2也可以改用一个赋值式x,y=x ,y / . Flattens1 , 一次得到两个函数.通过求数值解容易得到它的相图,In4绘制了解的相轨线如图13-44所示,图中说明原点是奇点,极限环的形状也已经得到.为了应付复杂的情况,需要设置可选参数:WorkingPrecisionAccuracyGoalPrecisionGoalMaxStepsStartingStepSize以上可选参数的默认值都为值比 WorkingPrecisio

19、n参见数值积分局部的介绍.计算结果的绝对误差.计算结果的相对误差.最大步数.初始步长.分方程,最大步长默认值为例6Automatic ,其中 AccuracyGoal 和 PrecisionGoal小10,当解趋于 0时应将 AccuracyGoal取成Infinity .1000.这个函数也可以解偏微分方程,最大步长默认值为解以下微分方程组：的默认对于常微200.,八1.(1) y ,满足初始条件y (0) =1的特解；4yx 3x 3y(2) y xz 26.5x y ,满足初始条件 x (0) =z (0) =0, y (0) =1z xy z解：In1 : =NDSolvey '

20、; x= =I/4yx, y0= =1 , y, x , 1,AccuracyGoal-20, PrecisionGoal -20, WorkingPrecisionOut1=y f InterpolatingFunction0, 1.000000000000000000000000000 , < >In2: =y1 / . %Out2=0.968912424710644784118519+ 0.2474039592545229296234109In3: =NDSolvex ' t= = -3(xt -yt),y' t = = -xt zt+36.5xt -yt,z

21、' t = = xt yt- zt,的特解. 25x0 = = z0 = = 0,y0= =1x,y, z , t , 0, 20,MaxSteps-3000Out3=xfInterpolatingFunction0.fInterpolatingFunction0.fInterpolatingFunction0.,20., < >,20., < >,20. , < >In4=ParametricPlot3DEvaluatext,yt , zt / . %t, 0, 20 , PlotPoints - 1000图13-45 3维相轨线Out3= -Gr

22、aphics3D-说明：以上范例中In1取高精度,而且是复系数方程.In2是求解在x=1时的近似1 1ii值,求精确解能得到准确值e4,读者可以求e4的近似值与 Out2的结果比拟,验证近似解的精确度确实很高.In3在求解时增大步数,成功地得到了由 In4绘制的如图13-45所示的解的相轨线.In4所示的绘图语句与前面例子中的不同,现在只要会模仿使用它们就行了,要想弄清原理请参阅相关Mathematica书籍.三、级数1. 求和与求积求有限或无穷和、积的函数是：i maxSumf , i , imin , imax求 f ,其中 imin 可以是-8,i i minimax可以是00 (即+&

23、#176;°),但是必须满足imin w imax.根本输入模板中也有求和专用的符号,使用模板输入更方便.Sumf , i , imin , imax , j , jmin , jmax, 求多重和,也可以使用根本输入模板连续屡次输入求和符号得到.i maxProductf , i , imin , imax求 f (i),根本输入模板中i i min也有求积符号.Productf,i , imin , imax , j , jmin , jmax,求多重积,也可以使用根本输入模板连续屡次输入求积符号得到.NSum 和NProduct得到数值解.2. 将函数展开为哥级数将函数展开为哥

24、级数的函数调用格式如下:Seriesf , x , X0, n将函数f (x)在X0处展成哥级数直到n次项为止.Seriesf , x , x0, n, y , y0, m 将函数 f (x, y)先对 y 后对 x 展开.例7展开以下函数为哥级数：sinx,、.,、,、xy(1) y=tgx, (2)y,(3)y =f(x),(4)y = e.x解：In1:=SeriesTanx, x, 0,9Out1=35x 2xx - 3157917x7 62x9ioox3152835In2=SeriesSinx /x,x , 0, 9Out2=246( x x x1 612050408x36288010oxIn3Out3=Seriesfx , x , 1, 7f1 f 1(x 1) 1f 1(x 1)2 1f (3)1(x 1)3 26f (4)1(x 1