數(shù)值分析中求解非線性方程的MATLAB求解程序(6種)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)數(shù)值分析中求解非線性方程的MATLAB求解程序（6種）1.求解不動點function k,p,err,P=fixpt(g,p0,tol,max1)%求解方程 x=g(x) 的近似值，初始值為 p0 %迭代式為 Pn+1=g(Pn)%迭代條件為：在迭代范圍內(nèi)滿足|k|1(根及附近且包含初值)k為斜率P(1)=p0;for k=2:max1 P(k)=feval(g,P(k-1); err=abs(P(k)-P(k-1); relerr=err/(abs(P(k)+eps)

2、; p=P(k); if (errtol)|(relerr0 break;endmax1=1+round(log(b-a)-log(delta)/log(2);for k=1:max1 c=(a+b)/2; yc=feval(f,c); if yc=0 a=c; b=c; elseif yb*yc0 b=c; yb=yc; else a=c; ya=yc; end if b-a0 disp(Note:f(a)*f(b)0);endfor k=1:max1 dx=yb*(b-a)/(yb-ya); c=b-dx; ac=c-a; yc=feval(f,c); if yc=0 break; els

3、eif yb*yc0 b=c; yb=yc; else a=c; ya=yc; end dx=min(abs(dx),ac); if abs(dx)delta|abs(yc)epsilon break; endendc;err=abs(b-a)/2;yc=feval(f,c);4.求解非線性方程根的近似位置function R=approot(X,epsilon)%求解根近似位置%為了粗估算方程f(x)=0在區(qū)間a,b的根的位置，%使用等間隔采樣點(xk,f(xk)和如下的評定準(zhǔn)則：%f(xk-1)與f(xk)符號相反，%或者|f(xk)|足夠小且曲線y=f(x)的斜率在%(xk,f(xk)附

4、近改變符號。Y=f(X);yrange=max(Y)-min(Y);epsilon2=yrange*epsilon;n=length(X);m=0;X(n+1)=X(n);Y(n+1)=Y(n);for k=2:n if Y(k-1)*Y(k)0 m=m+1; R(m)=(X(k-1)+X(k)/2; end s=(Y(k)-Y(k-1)*(Y(k+1)-Y(k); if (abs(Y(k)epsilon2)&(s=0) m=m+1; R(m)=X(k); endend5.牛頓拉夫森迭代function p0,err,k,y=newton(f,df,p0,delta,epsilon,max1)

5、%牛頓拉夫森迭代求解非線性方程%df需要自己f手動求得%p0為迭代初始值，max1為最大迭代次數(shù)for k=1:max1 p1=p0-feval(f,p0)/feval(df,p0); err=abs(p1-p0); relerr=2*err/(abs(p1)+delta); p0=p1; y=feval(f,p0); if(errdelta)|(relerrdelta)|(abs(y)epsilon) break endend6.割線法function p1,err,k,y=secant(f,p0,p1,delta,epsilon,max1)%割線法迭代求解非線性方程for k=1:max1 p2=p1-feval(f,p1)*(p1-p0)/(feval(f,p1)-feval(f,p0); err=abs(p2-p1); relerr=2*err/(abs(p2)+delta); p0