基于matlab的数值分析( 非线性方程求根)上机实验报告

《基于matlab的数值分析( 非线性方程求根)上机实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数值分析试验报告非线性方程求根二分法解方程1.题目：用二分法求方程的所有根2.方法：二分法3.程序functionx=erfenfa(a,b)if(a*a*a-3*a-1)*(b*b*b-3*b-1)>0disp('区间错误，请重新调试！');breakelsex=(a+b)/2.0;k=0;while(k<=100&abs(x*x*x-3*x-1)>0.0001&b-a>0.0001)if(x*x*x-3*x-1)*(a*a*a-3*a-1)<0b=x;elsea=x;endx=(a+b)/2.0;k=k+1;k,xif

2、k==100display('计算方法错误，请重新调试！');endendenddisp(['k=',num2str(k)])4.结果>>formatcompact>>erfenfa(1,2)k=1x=1.75000000000000k=2x=1.87500000000000k=3x=1.93750000000000k=4x=1.90625000000000k=5x=1.89062500000000k=6x=1.88281250000000k=7x=1.87890625000000k=8x=1.88085937500000

3、k=9x=1.87988281250000k=10x=1.87939453125000k=10ans=1.87939453125000>>erfenfa(-1,1)k=1x=-0.50000000000000k=2x=-0.25000000000000k=3x=-0.37500000000000k=4x=-0.31250000000000k=5x=-0.34375000000000k=6x=-0.35937500000000k=7x=-0.35156250000000k=8x=-0.34765625000000k=9x=-

4、0.34570312500000k=10x=-0.34667968750000k=11x=-0.34716796875000k=12x=-0.34741210937500k=13x=-0.34729003906250k=13ans=-0.34729003906250>>erfenfa(-2,-1)k=1x=-1.75000000000000k=2x=-1.62500000000000k=3x=-1.56250000000000k=4x=-1.53125000000000k=5x=-1.54687500000000k=6x=

5、-1.53906250000000k=7x=-1.53515625000000k=8x=-1.53320312500000k=9x=-1.53222656250000k=10x=-1.53173828125000k=11x=-1.53198242187500k=12x=-1.53210449218750k=12ans=-1.532104492187505，拓展：在程序编写完成调试过程中发现如果给出错误区间程序仍能给出一个错误答案，于是找同学讨教，讨论得出了正确的方法不动点解方程1.题目：解方程的根1.方法：不动点迭代法2.

6、程序functionx=budongdian(x0)p=x0;fork=1:10x0=p;p=20/(x0*x0+4*x0+10);ifabs(p-x0)<0.0001break;endkpifk==10disp('计算错误')endend3.结果>>formatcompact>>budongdian(1)k=1p=1.3333k=2p=1.1688k=3p=1.2468k=4p=1.2091k=5p=1.2271k=6p=1.2184k=7p=1.2226k=8p=1.2206k=9p=1.2216k=10p=1.221

7、1k=11p=1.2213k=12p=1.2212p=1.2213Newton迭代法求方程的根1.题目：求方程的根2.方法：newton迭代法3.程序：functionx=newton(a)p=a;fork=1:100a=p;p=a-(3*a*a*a+5*a*a+10)/(9*a*a+10*a);ifabs(p-a)<0.00001break;endifk==100display('计算错误');endkpendp1.结果>>formatlong>>newton(-2)k=1p=-2.37500000000000k=2p=

8、-2.30147484094853k=3p=-2.29793295142995p=-2.29792492582976>>对于不同的函数只需将文中x*x*x-3*x-1你所需求解的函数就ok了原文已完。下文为附加文档，如不需要，下载后可以编辑删除，谢谢！村支部书记在乡镇会议上的汇报发言原文已完。下文为附加