实验四+方程求实验.docVIP

实验 一. 实验目的 （1），学会解决某些实际问题，比较方法解题的不同处。 （2）熟悉Matlab编程环境，利用Matlab具体的。 二. 实验要求 用Matlab软件实现，并用实例在计算机上计算。 实验内容 1. 实验题目 （1）并给出了高精度的实根，试用Newton法和弦截法进行验证，要求精度，并绘制方程的图形。 （2）取，用迭代法求方程的根，然后用Aitken方法加速，要求精度为结果有4位有效数字。 （3）用快速弦截法求解方程，要求精度为，取作为开始值，并绘制的图形。 2. 设计思想 要求针对上述题目，详细分析每种算法的设计思想。 3. 对应程序 列出每种算法的程序。 4. 实验结果 列出相应的运行结果需要给出相应的图形。 实验体会 对实验过程进行分析总结，指出每种的设计要点及应注意的事项，以及自己通过实验所获得的对的理解。3. 对应程序（1）Newton法 function [x,k]=Newton(f,x0,emg) [f1,d1]=feval(f,x0); k=1; x(1)=x0; x(2)=x(1)-f1/d1; while abs(f1)emg u=1;k=k+1; [f1,d1]=feval(f,x(k)); x(k+1)=x(k)-u*f1/d1; [f2,d2]=feval(f,x(k+1)); end end 弦截法 function [x,k]=xianjie(f,x1,x2,emg) k=1; f1=feval(f,x1); f2=feval(f,x2); x(k)=x2-f2*(x2-x1)/(f2-f1); fk=feval(f,x(k)); k=k+1; x(k)=x(k-1)-(x(k-1)-x2)*f(k-1)/(fk-1-f2); while abs(x(k)-x(k-1))emg fk=feval(f,x(k)); x(k+1)=x(k)-(x(k)-x(1))*fk/(fk-x1); k=k+1; end （2） 迭代法 function [k,Deviation,xdDeviation,xk]=iteration(x0,k) x(1)=x0; for i=1:k Deviation=abs(x(i+1)-x(i)); xdDeviation=Deviation/(abs(x(i+1))+eps); i=i+1;xk=x(i); [(i-1) Deviation xdDeviation xk] end p=[(i-1) Deviation xdDeviation xk]; function y=f(x) y=exp(-x); Aitken加速 function [k,xk,yk,p]=Aitken(x0,emg,ddmax) x(1)=x0; for i=1:ddmax x1(i+1)=f(x(i)); x2(i+1)=f(x1(i+1)); x(i+1)=x2(i+1)-(x2(i+1)-x1(i+1))^2/(x2(i+1)-2*x1(i+1)+x(i)); Deviation=abs(x(i+1)-x(i)); xdDeviation=piancha/(abs(x(i+1))+eps); i=i+1;xk=x(i);yk=f(x(i)); if (Deviationemg)|(xdDeviationemg) k=i-1;xk=x(i);yk=f(x(i)); m=[0,1:i-1];p=[m,x1,x2,x]; return; end end if iddmax disp(μü′ú′?êy3?1y???¨μ?×?′ó?μ) k=i-1;xk=x(i);yk=f(x(i)); m=[0,1:i-1];p=[m,x1,x2,x]; return; end m=[0,1:i-1];p=[m,x1,x2,x]; function y=f(x) y=exp(-x); （3） 用快速弦截法 function [x,k]= fastxianjie(f,x1,x2,emg); k=1; y1=feval(f,x1); y2=feval(f,x2); x(k)=x2-y2*(x2-x1)/(y2-y1); y(k)=feval(f,x(k)); k=k+1; x(k)=x(k