matlab数值计算法程序,数值计算方法matlab程序.doc

数值计算方法matlab程序.doc

数值计算方法matlab程序二分法function[x0,k]=bisect1(fun1,a,b,ep)ifnargin0 x0=[fa,fb];k=0;return;endk=1;whileabs(b-a)/2>epx=(a+b)/2;fx=f(fun1,x);iffx*fa>fun1=inline( x^3-x-1 );>>[x0,k]=bisect1(fun1,1.3,1.4,1e-4)x0=1.3247k=7>>简单迭代法function[x0,k]=iterate1(fun1,x0,ep,N)ifnarginep>>[x0,k]=iterate1(fun1,1.5)x0=1.3247k=7>>fun1=inline( x^3-1 );>>[x0,k]=iterate1(fun1,1.5)x0=Infk=9>>Steffesen加速迭代(简单迭代法的加速)function[x0,k]=steffesen1(fun1,x0,ep,N)ifnarginep>>[x0,k]=steffesen1(fun1,1.5)x0=1.3247k=3>>fun1=inline( x^3-1 );>>[x0,k]=steffesen1(fun1,1.5)x0=1.3247k=6Newton迭代function[x0,k]=Newton7(fname,dfname,x0,ep,N)ifnarginep>>dfname=inline( 1+sin(x) );>>[x0,k]=Newton7(fname,dfname,pi/4,1e-8)x0=0.7391k=4数值计算方法matlab程序非线性方程求根的Matlab函数调用举例：1.求多项式的根：求f(x)=x^3-x-1=0的根：>>roots([10-1-1])ans=1.3247-0.6624+0.5623i-0.6624-0.5623i2.求一般函数的根>>fun=inline( x*sin(x^2-x-1) , x )fun=Inlinefunction:fun(x)=x*sin(x^2-x-1)>>fplot(fun,[-20.1]);gridon>>x=fzero(fun,[-2,-1])x=-1.5956>>x=fzero(fun,[-1-0.1])x=-0.6180[x,f,h]=fsolve(fun,-1.6)x=-1.5956f=数值计算方法matlab程序1.4909e-009h=1(h>0表示收敛，h>A=[234;352;4330];>>b=[6,5,32] b=6532>>[A,x]=gauss3(A,b)A=2.00003.00004.00006.000000.5000-4.0000-4.000000-2.0000-4.0000 x=-1382列选主元的高斯消元法：数值计算方法matlab程序function[A,x]=gauss5(A,b)%本算法用列选主元的高斯消元法求解线性方程组n=length(b);A=[A,b];fork=1:n-1%选主元[ap,p]=max(abs(A(k:n,k)));p=p+k-1;ifp>kt=A(k,:);A(k,:)=A(p,:);A(p,:)=t;end%消元A((k+1):n,(k+1):(n+1))=A((k+1):n,(k+1):(n+1))-A((k+1):n,k)/A(k,k)*A(k,(k+1):(n+1));A((k+1):n,k)=zeros(n-k,1);end%回代x=zeros(n,1);x(n)=A(n,n+1)/A(n,n);fork=n-1:-1:1x(k)=(A(k,n+1)-A(k,(k+1):n)*x((sk+1:n)))/A(k,k);end>>A=[234;352;4330];b=[6,5,32] ;>>[A,x]=gauss5(A,b)A=4.00003.000030.000032.000002.7500-20.5000-19.0000000.18180.3636x=-1382三角分解法：Doolittle分解function[L,U]=doolittle1(A)n=length(A);U=zeros(n);L=eye(n);U(1,:)=A(1,:);L(2:n,1)=A(2:n,1)/U(1,1);fork=2:nU(k,k:n)=A(k,k:n)-L(k,1:k-1)*U(1:k-1,k:n);数值计算方法matlab程序L(k+1:n,k)=A(k+1:n,k)-L(k+1:n,1:k-1)*U(1:k-1,n)/U(k,k);Endy=zeros(n,1);x=y;y(1)=b(1);fori=2:ny(i)=b(i)-L(i,1:i-1)*y(1:i-1);endx(n)=y(n)/U(n,n);fori=n-1:-1:1x(i)=(y(i)-U(i,i+1:n)*x(i+1:n))/U(i,i);end>>A=[123;252;315];b=[141820] ;>>[L,U,x]=doolittle1(A,b)L=1002103-81U=12301-400-36x=2.83331.33332.8333平方根法：function[L,x]=choesky3(A,b)n=length(A);L=zeros(n);L(:,1)=A(:,1)/sqrt(A(1,1));fork=2:nL(k,k)=A(k,k)-L(k,1:k-1)*L(k,1:k-1) ;L(k,k)=sqrt(L(k,k));fori=k+1:nL(i,k)=(A(i,k)-L(i,1:k-1)*L(k,1:k-1) )/L(k,k);endendy=zeros(n,1);x=y;y(1)=b(1)/L(1,1);数值计算方法matlab程序fori=2:ny(i)=(b(i)-L(i,1:i-1)*y(1:i-1))/L(i,i);endx(n)=y(n)/L(n,n);fori=n-1:-1:1x(i)=(y(i)-L(i+1:n,i) *x(i+1:n))/L(i,i);end>>A=[4-11;-14.252.75;12.753.5]A=4.0000-1.00001.0000-1.00004.25002.75001.00002.75003.5000>>b=[467.25] b=4.00006.00007.2500[L,x]=choesky3(A,b)L=2.000000-0.50002.000000.50001.50001.0000 x=111>>迭代法求方程组的解Jacobi迭代法：function[x,k]=jacobi2(a,b,x0,ep,N)%本算法用Jacobi迭代求解ax=b,用分量形式n=length(b);k=0;ifnarginep>>[