不同方法求方程的根matlab 代码.docVIP

1、%使用二分法找到方程 600 x^4 -550 x^3 +200 x^2 -20 x -1 =0 在区间[0.1,1]上的根， %误差限为 e=10^-4 disp(二分法) a=0.2;b=0.26; tol=0.0001; n0=10; fa=600*(a.^4)-550*(a.^3)+200*(a.^2)-20*a-1; for i=1:n0 p=(a+b)/2; fp=600*(p.^4)-550*(p.^3)+200*(p.^2)-20*p-1; if fp==0||(abs((b-a)/2)tol) disp(用二分法求得方程的根p=) disp(p) disp(二分迭代次数为：) disp(i) break; end if fa*fp0 a=p; else b=p; end end if i==n0~(fp==0||(abs((b-a)/2)tol)) disp(n0) disp(次二分迭代后没有求出方程的根) end 2、%使用牛顿法找到方程 600 x^4 -550 x^3 +200 x^2 -20 x -1 =0 在区间[0.1,1]上的根， %误差限为 e=10^-4 disp(牛顿法) p0=0.3; for i=1:n0 p=p0-(600*(p0.^4)-550*(p0.^3)+200*(p0.^2)-20*p0-1)./(2400*(p0.^3)-1650*p0.^2+400*p0-20); if(abs(p-p0)tol) disp(用牛顿法求得方程的根p=) disp(p) disp(牛顿迭代次数为：) disp(i) break; end p0=p; end if i==n0~(abs(p-p0)tol) disp(n0) disp(次牛顿迭代后没有求出方程的根) end 3、%使用割线法找到方程 600 x^4 -550 x^3 +200 x^2 -20 x -1 =0 在区间[0.1,1]上的根， %误差限为 e=10^-4 disp(割线法) 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; for i=2:n0 p=p1-q1*(p1-p0)/(q1-q0); if abs(p-p1)tol disp(用割线法求得方程的根p=) disp(p) disp(割线法迭代次数为：) disp(i) break; end p0=p1; q0=q1; pp=p1; p1=p; q1=600*(p.^4)-550*(p.^3)+200*(p.^2)-20*p-1; end if i==n0~(abs(p-pp)tol) disp(n0) disp(次割线法迭代后没有求出方程的根) end 4、%使用试位法找到方程 600 x^4 -550 x^3 +200 x^2 -20 x -1 =0 在区间[0.1,1]上的根， %误差限为 e=10^-4 disp(试位法) 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; for i=2:n0 p=p1-q1*(p1-p0)/(q1-q0); if abs(p-p1)tol disp(用试位法求得方程的根p=) disp(p) disp(试位法迭代次数为：) disp(i) break; end q=600*(p.^4)-550*(p.^3)+200*(p.^2)-20*p-1; if q*q10 p0=p1; q0=q1; end pp=p1; p1=p; q1=q; end if i==n0~(abs(p-pp)tol) disp(n0) disp(次试位法迭代后没有求出方程的根) end 5、%使用muller方法找到方程 600 x^4 -550