昆工数值分析上机作业（完整版）实用资料(可以直接使用，可编辑完整版实用资料，欢迎下载）国土资源工程学院测绘工程课题一迭代格式的比较要求:编制一个程序进行运算，最后打印出每种迭代格式的敛散情况；=1\*GB3\*MERGEFORMAT①、建立迭代程序的M文件：function[k,pc,xk]=diedail(x0,k)x(1)=x0fori=1:kx(i+1)=funl(x(i));pc=abs(x(i+1)-x(i));%偏差i=i+1;xk=x(i);%第k次迭代的结果[(i-1)pcxk]endp=[(i-1)pcxk];%输出迭代次数、偏差、第k次迭代的结果=2\*GB3\*MERGEFORMAT②、对于不同的迭代式子建立不同的funl.m文件建立迭代式子（1）的M文件：functiony1=funl(x)y1=(3*x+1)/x^2;在matlab命令窗口中运行以下命令>>[k,pc,xk]=diedail(0.7,10)%初值是0.7，迭代10次x=0.7000x=0.70006.3265ans=1.00005.62656.3265x=0.70006.32650.4992ans=2.00005.82740.4992x=0.70006.32650.499210.0231ans=3.00009.523910.0231x=0.70006.32650.499210.02310.3093ans=4.00009.71380.3093x=0.70006.32650.499210.02310.309320.1559ans=5.000019.846620.1559x=0.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.309320.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）建立迭代式子（2）的M文件：functiony1=funl(x)y1=(x^3-1)/3;在matlab命令窗口中运行以下命令>>[k,pc,xk]=diedail(0.6,5)%初值0.6，迭代5次x=0.6000x=0.6000-0.2613ans=1.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.2613-0.3393-0.3464-0.3472ans=4.00000.0008-0.3472x=0.6000-0.2613-0.3393-0.3464-0.3472-0.3473ans=5.00000.0001-0.3473k=5pc=9.9907e-005xk=-0.3473由以上结果可知迭代式（2）是收敛的。（3）建立迭代式子（3）的M文件：functiony1=funl(x)y1=(3*x+1)^(1/3);在matlab命令窗口中运行以下命令：>>[k,pc,xk]=diedail(0.9,12)%初值0.9，迭代12次结果为x=0.6000x=0.6000-0.2613ans=1.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.2613-0.3393-0.3464-0.3472