数值分析9-牛顿法和线切割法编程.docxVIP

作业十：1、在区间[0,1.5]上分别用二分法、牛顿法和线切割法编程求下面的函数的零点，精度要求10-10。F(x)= cos2(2x)-x2 2、求下面的函数在区间[0,13]上的所有零点，精度要求10-10。提示：先扫描得到解所在区间，再用迭代法求解。F(x)=2-x-cosx1、二分法：程序clear;clc;a=1;b=1.5;N=0;x=1f =@(x)(cos(2*x))^2-x^2;c=(a+b)/2;q=f(c)*f(b)while (abs(b-a)1e-10)||abs(f(a))1e-10;if f(c)*f(b)0 a=c;else b=c;end c=(a+b)/2; x=c; N=N+1;endfprintf(\n x = %.5f, f(x) = %.5f,12?-?·%d′? \n, x, f(x),N);程序结果牛顿法：程序clear;clc;f =@(x)(cos(2*x))^2-x^2;a=0;b=1.5;symsxg=diff((cos(2*x))^2-x^2);x1=b;N=0; x0=(a+b)/2;while abs(f(x0))1e-10; G=subs(g,x,x1);x0=x1-f(x1)/G; if abs(x0-x1)10^(-10); x1=x0;if x1ax1=a+(b+x0);endelsebreakend N=N+1;endfprintf(\n x = %.5f, f(x) = %.5f,共循环%d次\n, x0, f(x0),N);程序结果线切割法：程序clear;clc;a=0;b=1.5;N=0;f =@(x)(cos(2*x))^2-x^2;x0=a;x1=(a+b)/2;x=0:0.01:1.5;figureplot(x,(cos(2*x)).^2-x.^2);while abs(f(x1))1e-10;x2=x1-(f(x1)*(x1-x0))/(f(x1)-f(x0));if abs(x1-x2)1e-10x1=x2;elsebreakend N=N+1;endfprintf(\n x = %.5f, f(x) = %.5f,共循环%d次 \n, x2, f(x2),N);程序结果2、扫描法：clear;clc;a=0;b=13;N=0;f =@(x)2^(-x)-cos(x);for i=a:0.001:b-0.001%扫描确定解区间if f(i)==0fprintf(\n x = %.5f, f(x) = %.5f\n, i, f(i)); N=N+1;elseif f(i)*f(i+0.001)0 e=i; g=i+0.001; c=(e+g)/2;while (abs(e-g)1e-10)||abs(f(c))1e-10;%二分法求解if f(c)*f(e)0 g=c;else e=c;end c=(e+g)/2; x=c;endfprintf(\n x = %.5f, f(x) = %.5f \n, x, f(x)); N=N+1; elsecontinueendendend;fprintf(\n 共有%d个解 \n, N);程序结果