资源描述:
《昆工数值分析上机作业》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、国土资源工程学院测绘工程课题一迭代格式的比较要求:1、编制一个程序进行运算,最后打印出每种迭代格式的敛散情况;①.建立迭代程序的M文件:function[k,pc,xk]=diedail(xO,k)x(1)=x0fori=1:kx(i+1)=funl(x(i));pc=abs(x(i+l)-x(i));%偏差i=i+l;xk=x(i);%第)<次迭代的结果[(i-1)pcxk]endp=[(i-1)pcxk];%输出迭代次数、偏差、第k次迭代的结果②、对于不同的迭代式子建立不同的funl.m文件(1)X=^^~建立迭代式子(1)的M文件:functiony1=funi(x)
2、yl=(3*x+l)/xA2;在matlab命令窗口中运行以下命令»[k,pc,xk]二dieda订(0.7,10)%初值是0.7,迭代10次X二0.7000X=0.70006.3265ans1.00005.62656.3265X=0.70006.32650.4992ans2.00005.82740.4992X二0.70006.32650.499210.0231ans3.00009.523910.02310.70006.32650.499210.02310.3093ans—4.00009.71380.3093x—0.70006.32650.499210.02310.3093
3、20.1559ans二5.000019.846620.1559X0.70006.32650.499210.02310.309320.15590.1513ans=6.000020.00460.1513x=0.70006.32650.499210.02310.309320.15590.151363.5112ans=7.000063.359963.5112x—0.70006.32650.499210.02310.309320.15590.151363.51120.0475ans=8.000063.46380.0475x—0.70006.32650.499210.02310.309
4、320.15590.151363.51120.0475506.6983ans二9.0000506.6508506.6983x—0.70006.32650.499210.02310.309320.15590.151363.51120.0475506.69830.0059ans=10.0000506.69240.0059k二10pc=506.6924xk=0.0059由以上结果可知迭代式(1)是发散的。建立迭代式子(2)的M文件:functionyl=fun1(x)y1=(xA3-l)/3;在matlab命令窗口中运行以下命令»[k,pc,xk]二diedail(O.6,5)%
5、初值0.6,迭代5次0.60000.6000-0.26131.00000.8613-0.26130.6000-0.2613-0.33932.00000.0779-0.33930.6000-0.2613-0.3393-0.34643.00000.0071-0.34640.6000-0.2613-0.3393-0.34644.00000.0008-0.34720.6000-0.2613-0.3393-0.34645.00000.0001-0.3473ansXansxansxansxansX=-0.3472-0.3472-0-3473k=5pc二9.9907e-005xk=-0.
6、3473由以上结果可知迭代式(3)X=(3x+1)%(2)是收敛的。建立迭代式子(3)的M文件:functiony1=funi(x)yl=(3*x+l)A(1/3);在matlab命令窗口中运彳丁以下命令:»[k,pc,xk]=diedail(0.9,12)%初值0.9,迭代12次结果为0.60000.6000—0.26131.00000.8613-0.2613X=0.6000-0.2613-0.3393ans二2.00000.0779-0.3393x二0.6000-0.2613-0.3393-0.3464ans—3.00000.0071-0.3464x=0.6000-0.
7、2613-0.3393-0.3464ans—4.00000.0008-0.3472x=0.6000-0.2613-0.3393-0.3464ans二5.00000.0001-0.3473ans=-0.3472-0.3472k二-0.3473pc二9.9907e-005xk=-0.3473»[k,pc,xk]二dieda订(0.9,12)%初值0.9,迭代12次0.90000.90001.54671.00000.64671.54670.90001.54671.78002.00000.23341.78000.90001