(完整版)Mathematica数值分析和数值计算

1、 第五章数值分析和数值计算1如何求插值多项式给定n个点（令,片）,（i=l,2,.,n）,构造一个次数不超过n-1的多项式函数f（x）,使得f（xi）=yi,则称f（x）为拉格朗日插值多项式。11可以证明该多项式函数由公式（x一x）（x一x）（x一x）（x一x）（x一x）（x一x）y=y23n+y13n+1（x-x）（x-x）（x-x）1（x-x）（x-x）（x-x）12131n21232n（x一x）（x一x）（x一x）+y12n-1n（x一x）（x一x）（x一x）n1n2nn-1唯一给定。Mathematica提供了根据插值点数据计算拉格朗日插值多项式的函数InterpolatingPoly

2、nomial,下面是其调用格式：InterpolatingPolynomialdata,var作出以data为插值点数据，以var为变量名的插值多项式。例：iri；=data=lf16f2r46f4.60f（5f0;f=InterpolatingPolynomia1daxOuti=16+30+-3-C-4+x）3（-2|硝：=Expand黑140233k吨F下在多数情况下，我们构造插值函数的目的在于计算函数f(x)的值，而并不在意插值多项式的具体表示形式。对于拉格朗日插值多项式，当n较大时，得到的高次插值多项式由于截断误差和舍入误差的影响,往往误差较大。此时在实际应用中,一般采用分段插值。Ma

3、thematica提供了分段插值函数Interpolation，其使用格式为：Interpolationdata,InterpolationOrder-n这里InterpolationOrder-n指定插值多项式的次数，默认值为3。此外数据data中还可以包括插值点处的导数，格式为：xl,yl,dyl,x2,y2,dy2,.例：已知f(0)=0,f(l)=2f(0)=l,F(l)=l,求3次插值多项式f(x)，并计算f(0.72)和画出函数f(x)在0,1区间上的图形。nil=data=OrOf1lr2f1;f=InterpolationdaInterpolationOrder-3uti=工i

4、rtq匸po丄mlziELqEu.口crtio门Cl#1fln2:=f0.72utp=1.5287ln5;=PlotfXf(Xf0f1ut5=-Graphics-注：虽然Interpolation函数生成一个没有显式表达式的函数，但我们可以计算该函数在某点的值，也可以画出它在插值区间内的图形。2如何进行曲线拟合曲线拟合是根据已知的某函数在某些点处的函数值，按照最小二乘法的原理，得到该函数的近似函数的过程。一般来说，拟合的过程是先选定一组函数，通过它们的线性组合来近似原函数。在Mathematica中，使用函数Fit来进行拟合，其一般使用格式是：Fit拟合数据，用于拟合的函数列表，变量面是其常用

5、的几种形式：拟合函数形式意义Fitdata,l,x,x作线性拟合a+bxFitdata,l,x,x2,x构造二次拟合a+bx+cx2Fitdata,Tablexi,i,O,n,xn次多项式拟合ExpFitLogdata,l,x,x拟合曲线为ea+bx你也可以根据你的需要自己选择合适的拟合函数。例：ini：=data=TablePrimexf20Fitdataflfxfdcuti=-767363十277368in2：=FitdataflrxA2fxut2=-1.92363+2.2055x+0.0746753xSolve和FindRoot函数都能求方程的根，它们之间的区别在于，Solve主要用于求

6、多项式方程或方程组的所有根，而FindRoot则给出任意方程或方程组的一个根(注意不是全部),FindRoot的求根方法是牛顿切线法(或弦截法)，其使用格式为：函数意义FindRoot方程，x,x0从x=x0开始，计算方程的一个数值解FindRoot方程，x,x0,xmin,xmax求根区间范围(xmin,xmax)FindRoot方程组,x,x0,y,y0,.计算方程组的一个数值解3如何求解超越方程注：初值xO在一些情况下对寻找根是有影响的，如果找不到，请改变初值再试。初值的选取可以结合所求函数的图形来分析。例：ini：=FindRootk+Sinx=0.7fx1uti=kt0.836278

7、in2：=FindRootExpx*y=xA2rxA2十2=2#匚lfy.1autp=xt1.34403fyt0.435964如何求函数的极小值函数的极小值用函数FindMinimum计算，和FindRoot样，该函数也只能给出初值附近的一个极小值，你可以根据函数的图形来选择初值。下面是其使用格式：函数意义FindMinimumf,x,x0以x-x0为起始点计算f的一个局部极小值FindMinimumf,x,xs,x0,xlxs为起始点，在区间(x0,xl)上求f的极小值FindMinimumf,x,x0,y,y0,.计算多元函数的极值例：In8:=PlotSillX4-X/2fxf-10f1

8、0K-Hi-5a-4outpe=-Graphics-in9：=FindMinimumSinx4-x/2rxf51.22B37,4.1BB79)5如何求常微分方程的数值解求常微分方程或方程组的数值解的函数是NDSolve，由于是求数值解，所以初始条件也必须一起给出，和方程列在一起。函数使用的一般形式为：NDSolve方程1,方程2,.,yx,x,xmin,xmax对常微分方程，求函数y关于自变量x在区间xmin,xmax内的数值解；NDSolve方程1,方程2,.,ylx,y2x,.,x,xmin,xmax对常微分方程组，求函数yl,y2,关于自变量x在区间xmin,xmax内的数值解;例：(1

9、)求解y=x2+y2,y(0)=0,并计算y(1)和画出y(x)的图形(2)求解方程组x(t)=-y(t)-x2(t),y(t)=2x(t)-y(t),x(0)=y(0)=1解：ini：=NDSolveyTx=+ycy0=0fyxfg21Outi=y-IntsrpolatingFun-Cti2.f|丐：=yy=y3cf.outp=Interpolat.ingFunction-2.2.fxlnp:=yy/X-1out3=U350248in4：=Plotyyfxf-J2io厂i；12/-丄巾-out4j=-Graphics-in8：=NDSolvek1t=-yt凸yTt=2xt-ytfx0=y0=

10、1ftbyt,t,0,3ut6=x七-TnterpolatingFunction03丁tfyt-InterpolatingFunction0.,3.t-)mic：=fi=xt/.鋼I1fy=yt/.61outi&=InterpolatingFunction0.f3.ftfInterpolatingFunctionf3.ftIn18:=fx/t-1outi8=-0.1810441ln19:=Plotfyftf0f3uti9=-Graphics6线性规划函数ConstrainedMax和ConstrainedMin可求解任意线性规划在限定区域上的目标函数的最大值或最小值问题。其使用格式为：ConstrainedMin目标函数，限定区域,变量列表ConstrainedMax目标函数，限定区域,变量列表例：求解线性规划问题maxf=3x+2y-1s.t.0=x=10=y=2解：ini：=Cons七iiiedMQiK3x十2一1丁x1,y2注：函数中的变量限于取非负值，对不等式也是约定可以取等号；关于这方面的内容可以参考有关线性规划方面的书籍。习题1.已知函数f(x)过点(-1,3),(0,0.5),(0.5,0),(1,1),试求拉格郎日插值多项式并画出其图形;2给出离散数据如下：X1923303540Y19.028.547.068.290.0试