发布时间 : 星期一 文章（完整版）Mathematica求解方程（组）、级数更新完毕开始阅读2dc88859ce7931b765ce0508763231126edb77a8
方程（组）与级数的Mathematica求解
[学习目标]
1. 能用Mathematica求各种方程（组）的数值解和近似解； 2. 能对常见函数进行幂级数的展开。
一、 求解简单方程（组）
数学里的方程是带有变量的等式。一般地说，一个或一组方程总是对于方程中出现的变量的可能取值范围增加了一些限制。所谓求解方程就是设法把方程对于变量取值的限制弄清楚，最好的结果是用不含变量的表达式把变量的值表示出来。在这个系统里，方程也用含有变量的等式表示，要注意的是在这里等号用连续的两个等号(==)表示。方程的两端可以是任何数学表达式。
用户可以自己操作Mathematica系统去求解方程，例如使用移项一类的等价变换规则对方程加以变形、对方程的两端进行整理、把函数作用于方程的两端等等。系统也提供了一些用于求解方程的函数。
1、 求方程的代数解
最基本的方程求解函数是Solve，它可以用于 求解方程(主要是多项式方程)或方程组。Solve有两个参数，第一个参数是一个方程，或者是由若干个方程组的表(表示一个方程组);第二个参数是要求解的变量或变量表。例如，下面的式子对于变量X求解方程
x4?x3?6x2?1?0：
In[1]:=Solve[x^4-x^3-6x^2+1==0,x] 输入了这个表达式，系统立刻就能计算出方程的四个根，求出的解都是精确解(代数根)。对于一般的多项式，这样得出的解常常是用根式描述的复数。方程的解被表示成一个表，表中是几个子表，每一个子表的形式都是{x->...}，箭头后面是方程的一个解。Solve也可以求解多变量的方程或者方程组:
In[2]:=Solve[{x-2y==0,x^2-y==1},{x,y}]
这个表达式求解方程组: ??x?2y?02x??y?1.
有时求解方程会得到非常复杂的解。例如将上面的第一个方程稍加变形，所得到的解的
表达式就会变得很长:
In[3]:=Solve[x^4-x^3-6x^2=2==0,x]
这个表达式求出的解的表达式非常长，以至一个计算机屏幕显示不下。使用MS-DOS系统上的Mathematica的读者可以用键盘上的PgUP键和PgDn键把计算机屏幕上已经卷出的表达式翻回来阅读，附录B里提供了使用这类计算机的有关操作的更详细的说明。对于使用图形界面提供的功能去翻阅前面的结果。
在被求解的方程里还可以有其他符号参数，可以要求系统对于这一个或者那一个变量求解方程。对于Mathematica系统来说，方程中的符号变量（无论使用什么变量名）都是一样的。
对于处理复杂的方程，MATHEMETICA系统还提供了例外两个有用的函数。函数
Eliminate用于从方程组消去一个或几个变量，例如下面的表达式消去方程组里的变量Y:
IN[4]:=Eliminate[{X^2-2Y= =1,X+2Y= =4},Y]
Eliminate 的使用形式与Solve类似，它的第二个参数用于说明希望消去的变量。另一个函数Reduce用于化简复杂的方程或方程组，它试图用一组比较简单的逻辑关系来描述由原来方程所描述的变量之间的关系。它的使用形式与Solve, Eliminate一样，这里不举例字了。
2、 求方程的数值解
理论上已经证明，对于五次以上的多项式方程没有求代数解的一般方法，MATHEMATICA也求不出那些不能分解因式的五次以上的多项式方程的解，例如：
IN[5]:=SOLVE[X^5+5X^3-2= =0,X]
它返回一个带有函数TORUOES的表达式。可以把函数N作用到这个结果表达式上，求出方程的数值解：
IN[6]:=N[%]
可以看到系统同时求出了方程的五个根的时候可以直接用函数N和SOLVE结合完成工作： IN[7]:=N[SOLVE[X^6+4X^2-31= =0,X]]
在系统里直接提供了一个函数NSOLVE做这件事。 对于更复杂的方程（或方程组），用SOLVE求不出根，使用函数N也解决不了问题。对于这样的方程，用户可以使用REDUCE,ELIMINATE等函数去处理，设法把方程描述的变量之间的关系搞清楚。如果需要的就是方程的根，那么只要用求数值根的函数FINDROOT。函数FINDROOT求数值根所采用的方法与人们一般用计算机求数值根的方法一样。但是，由于MATHEMATIC有求导函数的能力，在这里计算有导函数的表达式的数值根就非常简单。不管表达式多么复杂，系统都能自动的求出它的导函数。求数值根使用的也是牛顿法，用户必须给FINDROOT提供一个初始值。下面一个简单的例子：
IN[8]:=FINDROOT[SIN[X]EXP[2X]-COS[X]= =0,{X,0.5}]
对于求不出导函数的表达式，例如用户自己定义的一个复杂计算函数，使用FindRoot提供函数值取不同符号（正负号）的两个点（用表的形式放在上面初始值0.5的位置），形式是：
IN[9]:= FindRoot [FUN1[X]= =0,{X,{0,1}}] 这里假使FUN1是用户定义的一个函数。
使用计算机求数值根的第一个问题是确定初始点，若初始值选取得不好将给求根带来困难。再一个麻烦是用户要自己求出函数的导函数。在MATHEMATICA系统里处理的对象是表达式，一个表达式可以服务于不同的用途，可以作为求值的对象，作为画图的对象，也可以作为演算的对象。当需要求一个表达式的数值根的时候，表达式的这样的多种功
能，或者说MATHEMATICA系统对于表达式的多方面的操作能力就表现出很大的优越性。一个代数表达式，无论多么复杂，MATHEMATICA系统都可以直接求出它的导函数，可以作出它的图形。从图形上我们很容易认识这个函数表达式在某一个区间的大致性质，包括它的根的出现和分布情况。对表达式的这些认识为人们确定如何取初始值、如何求根提供了很有价值的线索。这样，某些比较难以处理的问题可能就容易解决了。
二、 求解常微分方程（组） 1、常微分方程（组）的精确解
Mathematica能求常微分方程（组）的准确解，能求解的类型大致覆盖了人工求解的范
围，功能很强。但不如人灵活（例如在隐函数和隐方程的处理方面），输出的结果与教材上的答案可能在形式上不同。另外，Mathematica求数值解也很方便，且有利于作出解的图形。 求准确解的函数调用格式如下：
DSolve[eqn，y[x]，x] 求方程eqn的通解y（x），其中自变量是x。 DSolve[{eqn，y[x0]= =y0}，y[x]，x] 求满足初始条件y（x0）= y0的特解y（x）。 DSolve[{eqn1，eqn2，…}，{y1[x]，y2[x]，…}，x] 求方程组的通解。 DSolve[{equ1，…，y1[x0]= =y10，…}，{y1[x]，y2[x]，…}，x] 求方程组的特解。 说明：应当特别注意，方程及各项参数的表述方式很严格，容易出现输入错误。微分方程的表示法只有通过例题才能说清楚。
例1
解下列常微分方程（组）：
5?y??z1?y22y2?(x?1)， （1）y??（2）y??， （3） ， ?3(x?x)yx?1?z???y?y??z （4）?的通解及满足初始条件y（0）=0，z（0）=1的特解。
?z??y? 解：In[1]：=DSolve[y′[x]= =2y[x]/（x+1）+（x+1）^（5/2）， y[x]，x] Out[1]=??y[x]?????2??(1?x)7/2?(1?x)2c[1]?? 3?? In[2]：=DSolve[y′[x]= =（1+y[x]^2）/((x+x^3）y[x])，y[x]，x]
?1? Out[2]={{y[x]??11?c[1]?1??c[1]22xx}， {y[x]?}} 111?21?2xx In[3]：=DSolve[{y′[x]= =z[x]，z′[x]= = -y[x]}， {y[x]，z[x]}，x]
Out[3]={{y[x]→C[1]Cos[x]+ C[2]Sin[x]， z[x]→C[2]Cos[x]- C[1]Sin[x]}}
In[4]：=DSolve[{y′[x]= =z[x]，z′[x]= = -y[x]，y[0]= =0，z[0]= =1}， {y[x]，z[x]}，x]
Out[4]={{y[x]→Sin[x]，z[x]→Cos[x]}}
提示：认真观察上例，可以从中学习输入格式，未知函数总带有自变量，等号用连续键入两个等号表示，这两点由于不习惯会出错！导数符号用键盘上的撇号，连续两撇表示二阶导数，这与习惯相同。自变量、未知量、初始值的表示法与普通变量相同。
说明：输出结果总是尽量用显式解表出，有时反而会使表达式变得复杂，这与教科书的习惯不同。当求显式解遇到问题时，会给出提示。通解中的任意常数用C[1]，C[2]，…表示。
例2
求解下列微分方程：
?x （1）y????3y???3y??y?(x?5)e，（2）x?(y?)?1，（3）
22y??xy。
解：In[1]：=DSolve[y???[x]+3y″[x] +3y′[x] + y[x] = =（x - 5）Exp[-x]， y[x]，x]
1?x2?x2??x?5x2x3? Out[1]={{y[x]?ex??5x???ex???? ????22?3???25x3x4??x?x?x2 1e?x?}} ???????eC[1]?exC[2]?exC[3]2?34? In[2]：=Simplify[%] Out[2]={{y[x]?1?xe(?20x3?x4?24C[1]?24xC[2]?24x2C[3])}} 24 In[3]：=DSolve[x^2 + y′[x]^2 = = 1，y[x]，x]
1ArcSin[x]x1?x2??C[1]}， 221ArcSin[x]2 {y[x]?x1?x??C[1]}}
22 Out[3]={{y[x]?? In[4]：=DSolve[Sqrt[y′[x]] = = x y[x]，y[x]，x] Out[4]={{y[x]??3}}
x3?C[1] 说明：由以上可以看出对方程的类型并无限制，但是输出的答案未必符合习惯，例如第一个方程的答案需要化简，有时即使化简后也未必与教材上的答案一致。
例3
求微分方程xy+ y - e= 0在初始条件y|x=1 = 2e下的特解。
′′
x
解：In[1]：=DSolve[x*y[x]+y[x]-E^x= =0，y[1]= =2E，y[x]，x]
e?ex
Out[1]= {{y[x]→}}
x