mathematica求解方程(组)、级数

《mathematica求解方程(组)、级数》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

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

2、用移项一类的等价变换规则对方程加以变形、对方程的两端进行整理、把函数作用于方程的两端等等。系统也提供了一些用于求解方程的函数。1、求方程的代数解最基本的方程求解函数是Solve，它可以用于求解方程(主要是多项式方程)或方程组。Solve有两个参数，第一个参数是一个方程，或者是由若干个方程组的表(表示一个方程组);第二个参数是要求解的变量或变量表。例如，下面的式子对于变量X求解方程：In[1]:=Solve[x^4-x^3-6x^2+1==0,x]输入了这个表达式，系统立刻就能计算出方程的四个根，求出的解都是精确解(代数根)。对于一般的多项式，这样得出的解常常是用根

3、式描述的复数。方程的解被表示成一个表，表中是几个子表，每一个子表的形式都是{x->...}，箭头后面是方程的一个解。Solve也可以求解多变量的方程或者方程组:In[2]:=Solve[{x-2y==0,x^2-y==1},{x,y}]这个表达式求解方程组:有时求解方程会得到非常复杂的解。例如将上面的第一个方程稍加变形，所得到的解的表达式就会变得很长:In[3]:=Solve[x^4-x^3-6x^2=2==0,x]这个表达式求出的解的表达式非常长，以至一个计算机屏幕显示不下。使用MS-DOS系统上的Mathematica的读者可以用键盘上的PgUP键和PgDn键

4、把计算机屏幕上已经卷出的表达式翻回来阅读，附录B里提供了使用这类计算机的有关操作的更详细的说明。对于使用图形界面提供的功能去翻阅前面的结果。在被求解的方程里还可以有其他符号参数，可以要求系统对于这一个或者那一个变量求解方程。对于Mathematica系统来说，方程中的符号变量（无论使用什么变量名）都是一样的。对于处理复杂的方程，MATHEMETICA系统还提供了例外两个有用的函数。函数Eliminate用于从方程组消去一个或几个变量，例如下面的表达式消去方程组里的变量Y:IN[4]:=Eliminate[{X^2-2Y==1,X+2Y==4},Y]Eliminat

5、e的使用形式与Solve类似，它的第二个参数用于说明希望消去的变量。另一个函数Reduce用于化简复杂的方程或方程组，它试图用一组比较简单的逻辑关系来描述由原来方程所描述的变量之间的关系。它的使用形式与Solve,Eliminate一样，这里不举例字了。2、求方程的数值解理论上已经证明，对于五次以上的多项式方程没有求代数解的一般方法，MATHEMATICA也求不出那些不能分解因式的五次以上的多项式方程的解，例如：IN[5]:=SOLVE[X^5+5X^3-2==0,X]它返回一个带有函数TORUOES的表达式。可以把函数N作用到这个结果表达式上，求出方程的数值解：

6、IN[6]:=N[%]可以看到系统同时求出了方程的五个根的时候可以直接用函数N和SOLVE结合完成工作：IN[7]:=N[SOLVE[X^6+4X^2-31==0,X]]在系统里直接提供了一个函数NSOLVE做这件事。对于更复杂的方程（或方程组），用SOLVE求不出根，使用函数N也解决不了问题。对于这样的方程，用户可以使用REDUCE,ELIMINATE等函数去处理，设法把方程描述的变量之间的关系搞清楚。如果需要的就是方程的根，那么只要用求数值根的函数FINDROOT。函数FINDROOT求数值根所采用的方法与人们一般用计算机求数值根的方法一样。但是，由于MATH

7、EMATIC有求导函数的能力，在这里计算有导函数的表达式的数值根就非常简单。不管表达式多么复杂，系统都能自动的求出它的导函数。求数值根使用的也是牛顿法，用户必须给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

8、是用户定义