迭代格式的比较.docVIP

课题一 迭代格式的比较 一、问题提出 设方程f(x)=x - 3x–1=0 有三个实根x*1=1.8793, x*2=0.34727,x*3 =-1.53209 现采用下面六种不同计算格式，求f(x)=0的根 x*1 或x*2 1. 2. 3. 4. 5. 6. 二、要求 1、编制一个程序进行运算，最后打印出每种迭代格式的敛散情况； 2、用事后误差估计来控制迭代次数，并且打印出迭代的次数； 3、初始值的选取对迭代收敛有何影响； 4、分析迭代收敛和发散的原因。 三、目的和意义 1、通过实验进一步了解方程求根的算法； 2、认识选择计算格式的重要性； 3、掌握迭代算法和精度控制； 4、明确迭代收敛性与初值选取的关系。1．六种不同计算格式运行程序及运算结果 1.1 运行程序 （1） function iteration_1 x0=input(输入初始值x0 = ); k=input(输入迭代次数k = ); x(1)=x0 ; i=1; while i = k x(i+1)=fun(x(i)) i=i+1; xk=x(i); end function y =fun(x) y=(3*x+1)/x^2; （2） function iteration_2 x0=input(输入初始值x0 = ); k=input(输入迭代次数k = ); x(1)=x0 ; i=1; while i = k x(i+1)=fun(x(i)) i=i+1; xk=x(i); end function y =fun(x) y=(x^3-1)/3; （3） function iteration_3 x0=input(输入初始值x0 = ); k=input(输入迭代次数k = ); x(1)=x0 ; i=1; while i = k x(i+1)=fun(x(i)) i=i+1; xk=x(i); end function y =fun(x) y=(3*x+1)^(1/3); （4） function iteration_4 x0=input(输入初始值x0 = ); k=input(输入迭代次数k = ); x(1)=x0 ; i=1; while i = k x(i+1)=fun(x(i)) i=i+1; xk=x(i); end function y =fun(x) y=1/(x^2-3); （5） function iteration_5 x0=input(输入初始值x0 = ); k=input(输入迭代次数k = ); x(1)=x0 ; i=1; while i = k x(i+1)=fun(x(i)) i=i+1; xk=x(i); end function y =fun(x) y=(3+1/x)^(1/2); （6） function iteration_6 x0=input(输入初始值x0 = ); k=input(输入迭代次数k = ); x(1)=x0 ; i=1; while i = k x(i+1)=fun(x(i)) i=i+1; xk=x(i); end function y =fun(x) y=x-(1/3)*((x^3-3*x-1)/(x^2-1)); 1.2运行结果 在matlab的command window中调用写好的m文件，运行方式如下 iteration 输入初始值x0 = 输入迭代次数k = （1） iteration_1 输入初始值x0 = 1 输入迭代次数k = 5 计算格式（1）最终运行结果如下： x = 1.0000 4.0000 0.8125 5.2071 0.6130 7.5549 由此知，计算格式（1）发散。 （2） iteration_2 输入初始值x0 = 1 输入迭代次数k = 5 计算格式（2）最终运行结果如下： x = 1.0000 0 -0.3333 -0.3457 -0.3471 -0.3473 由此知，计算格式（2）收敛。 （3） iteration_3 输入初始值x0 = 1 输入迭代次数k = 8 计算格式（2）最终运行结果如下： x = 1.0000 1.5874 1.7928 1.8545 1.8723 1.8774 1.8788 1.8792 1.8793 由此知，计算格式（3）收敛。 （4） iteration_4 输入初始值x0 = 1 输入迭代次数k = 5 计算格式（2）最终运行结果如下： x = 1.0000 -0.5000 -0.3636 -0.348