数值分析第二章上机报告.docxVIP

数值分析上机报告姓名：学号第二章上机题1.程序：（1）定义函数f(x)function Fx=fx(x)Fx=(x^3)/3-x;（2）定义函数df(x)function Fx=dfx(x)Fx=x^2-1;（3）Newton法程序clearef=10^-6; %???¨èYDí?ó2?10^-6k=0;x0=input(Please input initial value Xo:);disp(k Xk);fprintf(0 %f\n,x0); flag=1;while flag==1 k=10^3 x1=x0-fx(x0)/dfx(x0);if abs(x1-x0)ef flag=0;end k=k+1; x0=x1;fprintf(%d %f\n,k,x0); end（4）寻找最大delta值程序clearflag=1;k=1;x0=0;while flag==1 delta=k*10^-6; x0=delta; k=k+1; m=0; flag1=1;while flag1==1 m=10^3 x1=x0-fx(x0)/dfx(x0);if abs(x1-x0)10^-6 flag1=0;end m=m+1; x0=x1;endif flag1==1||abs(x0)=10^-6 flag=0; endendfprintf(The maximun delta is %f\n,delta); 2.运行结果：（1）运行search.m文件结果为：The maximum delta is 0.774597即得最大的δ为0.774597，Newton迭代序列收敛于根=0的最大区间为（-0.774597，0.774597）；（2）运行Newton.m文件在区间上各输入若干个数，计算结果如下：区间上取-1000，-100，-50，-30，-10，-8，-7，-5，-3，-1.5结果显示，以上初值迭代序列均收敛于-1.732051，即根在区间即区间（-1，-0.774597）上取-0.774598，-0.8，-0.85，-0.9，-0.99，计算结果如下：计算结果显示，迭代序列局部收敛于-1.732051，即根，局部收敛于1.730251，即根。在区间即区间（-0.774597，0.774597）上，由search.m的运行过程表明，在整个区间上均收敛于0，即根。在区间即区间（0.774597，1）上取0.774598，0.8，0.85，0.9，0.99，计算结果如下：计算结果显示，迭代序列局部收敛于-1.732051，即根，局部收敛于1.730251，即根。区间上取100，60，20，10，7，6，4，3，1.5，计算结果如下:]结果显示，以上初值迭代序列均收敛于1.732051，即根。综上所述：(-∞,-1)区间收敛于-1.73205,(-1,δ)区间局部收敛于1.73205,局部收敛于-1.73205,(-δ,δ)区间收敛于0,(δ,1)区间类似于(-1,δ)区间，(1,∞)收敛于1.73205。通过本上机题，明白了对于多根方程，Newton法求方程根时，迭代序列收敛于某一个根有一定的区间限制，在一个区间上，可能会局部收敛于不同的根。