预览加载中，请您耐心等待几秒...

昆工数值分析上机作业.doc

预览
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲，该文档总共38页，到这已经超出免费预览范围，如果喜欢就直接下载吧～

下载提示 文本预览

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

昆工数值分析上机作业（完整版）实用资料 (可以直接使用，可编辑完整版实用资料，欢迎下载） 国土资源工程学院测绘工程 课题一迭代格式的比较 要求: 编制一个程序进行运算，最后打印出每种迭代格式的敛散情况； =1\*GB3\*MERGEFORMAT①、建立迭代程序的M文件： function[k,pc,xk]=diedail(x0,k) x(1)=x0 fori=1:k x(i+1)=funl(x(i)); pc=abs(x(i+1)-x(i));%偏差 i=i+1; xk=x(i);%第k次迭代的结果 [(i-1)pcxk] end p=[(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.7000 x= 0.70006.3265 ans= 1.00005.62656.3265 x= 0.70006.32650.4992 ans= 2.00005.82740.4992 x= 0.70006.32650.499210.0231 ans= 3.00009.523910.0231 x= 0.70006.32650.499210.02310.3093 ans= 4.00009.71380.3093 x= 0.70006.32650.499210.02310.309320.1559 ans= 5.000019.846620.1559 x= 0.70006.32650.499210.02310.309320.15590.1513 ans= 6.000020.00460.1513 x= 0.70006.32650.499210.02310.309320.15590.151363.5112 ans= 7.000063.359963.5112 x= 0.70006.32650.499210.02310.309320.15590.151363.51120.0475 ans= 8.000063.46380.0475 x= 0.70006.32650.499210.02310.309320.15590.151363.51120.0475506.6983 ans= 9.0000506.6508506.6983 x= 0.70006.32650.499210.02310.309320.15590.151363.51120.0475506.69830.0059 ans= 10.0000506.69240.0059 k= 10 pc= 506.6924 xk= 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.6000 x= 0.6000-0.2613 ans= 1.00000.8613-0.2613 x= 0.6000-0.2613-0.3393 ans= 2.00000.0779-0.3393 x= 0.6000-0.2613-0.3393-0.3464 ans= 3.00000.0071-0.3464 x= 0.6000-0.2613-0.3393-0.3464-0.3472 ans= 4.00000.0008-0.3472 x= 0.6000-0.2613-0.3393-0.3464-0.3472-0.3473 ans= 5.00000.0001-0.3473 k= 5 pc= 9.9907e-005 xk= -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.6000 x= 0.6000-0.2613 ans= 1.00000.8613-0.2613 x= 0.6000-0.2613-0.3393 ans= 2.00000.0779-0.3393 x= 0.6000-0.2613-0.3393-0.3464 ans= 3.00000.0071-0.3464 x= 0.6000-0.2613-0.3393-0.3464-0.3472