搜索
数值分析第九次实验报告姓名：董安葳学号：5123119题目：斯特芬森迭代法和牛顿法实验方法：本人根据课本中的解非线性方程的斯特芬森迭代法和牛顿法，先是深刻理解算法，在草稿纸上写出伪代码，然后上机写出程序代码，并测试，最后得到斯特芬森迭代法和牛顿法的程序，并用写的程序完成了课后题目。实验过程：sitefensen.mfunctionx=sitefensen(f,x0,min)%其中f为迭代函数，x0为迭代初值，min为设置精度X=zeros(100);%设置迭代100次Y=zeros(100);Z=zeros(100);f=inline(f);X(1)=x0;t=0;forii=1:100%设置迭代100次Y(ii)=f(X(ii));Z(ii)=f(Y(ii));X(ii+1)=X(ii)-((Y(ii)-X(ii))^2)/(Z(ii)-2*Y(ii)+X(ii));ifabs(X(ii+1)-X(ii))min%一旦达到精度就不用再迭代了t=1;x=X(ii+1);break;endendift~=1x=X(100);%迭代到了100次endniudunfa.mfunctionx=niudunfa(f,x0,min)%其中f为原方程，x0为迭代初值，min为设置精度X=zeros(100);%设置迭代100次X(1)=x0;df=diff(f);t=0;forii=1:100%设置迭代100次X(ii+1)=X(ii)-subs(f,X(ii))/subs(df,X(ii));ifabs(X(ii+1)-X(ii))min%一旦达到精度就不用再迭代了t=1;x=X(ii+1);break;endendift~=1%迭代到了100次x=X(100);end实验截图：课后5题课后7题