二分法matlab实验报告总结,matlab计算方法迭代法牛顿法二分法实验报告.doc

matlab计算方法迭代法牛顿法二分法实验报告.doc

姓名实验报告成绩评语指导教师(签名)年月日说明指导教师评分后，实验报告交院(系)办公室保存。实验一方程求根一、实验目的用各种方法求任意实函数方程在自变量区间A，B上，或某0XF一点附近的实根。并比较方法的优劣。二、实验原理1、二分法对方程在A，B内求根。将所给区间二分，在分点0XF判断是否；若是，则有根。否则，继续判断是否2ABXF2ABX,若是，则令,否则令。否则令。重复此过程直FXBX至求出方程在A,B中的近似根为止。0XF(2)、迭代法将方程等价变换为Ψ()形式，并建立相应的迭代公式XFXΨ()。1KX(3)、牛顿法若已知方程的一个近似根，则函数在点附近可用一阶泰勒多项0X0X式来近似，因此方程可近似表示为 001XFXFPF0XF设,则。取作为原方程新的近似根 0F 0F0X 0FX，然后将作为代入上式。迭代公式为。1X1X01K0 KXF三、实验设备MATLAB70软件四、结果预测(1)009033(2)009052(3)0,090521X5X2X五、实验内容(1)、在区间0,1上用二分法求方程的近似根，要求误差021XE不超过。3105(2)、取初值，用迭代公式，求方程的0X1KX0 KXF021XE近似根。要求误差不超过。35(3)、取初值，用牛顿迭代法求方程的近似根。要求0X02XE误差不超过。315六、实验步骤与实验程序(1)二分法第一步在MATLAB70软件，建立一个实现二分法的MATLAB函数文件AGUI_BISECTM如下FUNCTIONXAGUI_BISECTFNAME,A,B,EFNAME为函数名，A,B为区间端点，E为精度FAFFNAME,A把A端点代入函数，求的FAFBFFNAME,B把B端点代入函数，求的FBIFFAFB0ERROR 两端函数值为同号 END如果FAFB0，则输出两端函数值为同号K0XAB/2WHILEBA2E循环条件的限制FXFFNAME,X把X代入代入函数，求的FXIFFAFXFUNINLINE EXPX10X2 XAGUI_BISECTFUN,0,1,05103第三步得到计算结果，且计算结果为KX00500000000000001025000000000000201250000000000030062500000000004009375000000000500781250000000060085937500000007008984375000000800917968750000090090820312500001000903320312500011009033203125000(2)迭代法第一步第一步在MATLAB70软件，建立一个实现迭代法的MATLAB函数文件AGUI_MAINM如下FUNCTIONXAGUI_MAINFNAME,X0,EFNAME为函数名DFNAME的函数FNAME的导数，X0为迭代初值E为精度，N为最大迭代次数默认为100N100XX0把X0赋给X，再算X2E赋给X0X0X2EK0WHILEABSX0XEXX0把X0赋给X，再算X2E赋给X0X0X2EK0WHILEABSX0XEKFUNINLINE EXPX10X2 DFUNINLINE EXPX10 XAGUI_NEWTONFUN,DFUN,0,05103第三步得出结果，且结果为KX100909090909090920090525108583393009052510858339以下是结果的屏幕截图七、实验结果(1)009033(2)009052(3)0,090521X5X2X八、实验分析与结论由上面的对二分法、迭代法、牛顿法三种方法的三次实验结果，我们可以得出这样的结论二分法要循环K11次，迭代法要迭代K5次，牛顿法要迭代K2次才能达到精度为的要求，而且方程3105的精确解经计算，为00905250,计算量从大到小依次是二021XE分法,迭代法,牛顿法。由此可知，牛顿法和迭代法的精确度要优越于二分法。而这三种方法中，牛顿法不仅计算量少，而且精确度高。从而可知牛顿迭代法收敛速度明显加快。可是迭代法是局部收敛的，其收敛性与初值X0有关。二分法收敛虽然是速度最慢，但也有自己的优势，可常用于求精度不高的近似根。迭代法是逐次逼近的方法,原理简单,但存在收敛性和收敛速度的问题。对与不同的题目，可以从三种方法的优缺点考虑用哪一种方法比较好。