数值分析matlab代码.doc

《数值分析matlab代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1、%用牛顿法求f(x)=x-sinx的零点，e=10^(-6)disp('牛顿法');i=1;n0=180;p0=pi/3;tol=10^(-6);fori=1:n0p=p0-(p0-sin(p0))/(1-cos(p0));ifabs(p-p0)<=10^(-6)disp('用牛顿法求得方程的根为')disp(p);disp('迭代次数为：')disp(i)break;endp0=p;endifi==n0&&~(abs(p-p0)<=10^(-6))disp(n0)disp('次牛顿迭代后无法求出方程的解')en

2、d2、disp('Steffensen加速');p0=pi/3;fori=1:n0p1=0.5*p0+0.5*cos(p0);p2=0.5*p1+0.5*cos(p1);p=p0-((p1-p0).^2)./(p2-2.*p1+p0);ifabs(p-p0)<=10^(-6)disp('用Steffensen加速求得方程的根为')disp(p);disp('迭代次数为：')disp(i)break;endp0=p;endifi==n0&&~(abs(p-p0)<=10^(-6))disp(n0)disp('次Stef

3、fensen加速后无法求出方程的解')end1、%使用二分法找到方程600x^4-550x^3+200x^2-20x-1=0在区间[0.1,1]上的根，%误差限为e=10^-4disp('二分法')a=0.2;b=0.26;tol=0.0001;n0=10;fa=600*(a.^4)-550*(a.^3)+200*(a.^2)-20*a-1;fori=1:n0p=(a+b)/2;fp=600*(p.^4)-550*(p.^3)+200*(p.^2)-20*p-1;iffp==0

4、

5、(abs((b-a)/2)

6、disp('用二分法求得方程的根p=')disp(p)disp('二分迭代次数为：')disp(i)break;endiffa*fp>0a=p;elseb=p;endendifi==n0&&~(fp==0

7、

8、(abs((b-a)/2)

9、0.^4)-550*(p0.^3)+200*(p0.^2)-20*p0-1)./(2400*(p0.^3)-1650*p0.^2+400*p0-20);if(abs(p-p0)

10、[0.1,1]上的根，%误差限为e=10^-4disp('割线法')p0=0.2;p1=0.25;q0=600*(p0.^4)-550*(p0.^3)+200*(p0.^2)-20*p0-1;q1=600*(p1.^4)-550*(p1.^3)+200*(p1.^2)-20*p1-1;fori=2:n0p=p1-q1*(p1-p0)/(q1-q0);ifabs(p-p1)

11、1;pp=p1;p1=p;q1=600*(p.^4)-550*(p.^3)+200*(p.^2)-20*p-1;endifi==n0&&~(abs(p-pp)

12、00*(p1.^4)-550*(p1.^3)+200*(p1.^2)-20*p1-1;fori=2:n0p=p1-q1*(p1-p0)/(q1-q0);ifabs(p-p1)