数值分析-ATLAB相关算法.docVIP

数值分析-MATLAB算法 刘 亚 1、四阶龙格库塔法： function yout=xin(bianliang) %定义输入输出 clear all x0=0;xn=1;y0=1;h=0.1; %设置初始值、区间和步长 [y,x]=lgkt4j(x0,xn,y0,h); %四阶龙格库塔法 n=length(x); fprintf( i x(i) y(i)\n); %输出格式 for i=1:n fprintf(%2d %3.3f %4.4f\n,i,x(i),y(i)); end function [y,x]=lgkt4j(x0,xn,y0,h) x=x0:h:xn; %设置区间 n=length(x); y1=x; y1(1)=y0; for i=1:n K1=f(x(i),y1(i)); K2=f(x(i)+h/2,y1(i)+h/2*K1); K3=f(x(i)+h/2,y1(i)+h/2*K2); K4=f(x(i)+h,y1(i)+h*K3); y1(i+1)=y1(i)+h/6*(K1+2*K2+2*K3+K4); end y=y1; function Dy=f(x,y) Dy=y-2*x/y; C语言程序 #includemath.h main() { float x=0,y0=1,h=0.2,y1,k1,k2,k3,k4; k1=y0-2*x/y0; k2=y0+h/2*k1-(2*x+h)/(y0+h/2*k1); k3=y0+h/2*k2-(2*x+h)/(y0+h/2*k2); k4=y0+h*k3-(2*x+2*h)/(y0+h*k3); y1=y0+h/6*(k1+2*k2+2*k3+k4); do { printf(%5.4f\n,y1); x=x+h; y0=y1; k1=y0-2*x/y0; k2=y0+h/2*k1-(2*x+h)/(y0+h/2*k1); k3=y0+h/2*k2-(2*x+h)/(y0+h/2*k2); k4=y0+h*k3-(2*x+2*h)/(y0+h*k3); y1=y0+h/6*(k1+2*k2+2*k3+k4); } while(x1); } 2、幂法求特征值 function [m x biaozhi]=mifa(A,jingdu,cishu) %幂法求矩阵最大特征值，其中 %m为绝对值最大的特征值，x为对应最大特征值的特征向量 %biaozhi表明迭代是否成功 if nargin3 cishu=100; end if nargin2 jingdu=1e-5; end n=length(A); x=ones(n,1); biaozhi=迭代失败！; k=0;m1=0; while k=cishu v=A*x; [vmax,k]=max(abs(v)); m=v(k); x=v/m; if abs(m-m1)jingdu biaozhi=迭代成功！;break; end m1=m; k=k+1; end 3、拉格朗日插值 function [c,l]=lglr(x,y) %x为n个节点的横坐标组成的向量，y为纵坐标组成的向量 %c为插值函数的系数组成的向量 %输出为差值多项式的系数 w=length(x); n=w-1; l=zeros(w,w); for k=1:n+1 v=1; for j=1:n+1 if k~=j v=conv(v,poly(x(j)))/(x(k)-x(j)); end end l(k,:)=v; end c=y*l; 举例 4、改进欧拉法 function yout=gaijinoula(f,x0,y0,xn,n) %定义输入输出 x=zeros(1,n+1);y=zeros(1,n+1); x(1)=x0;y(1)=y0; h=(xn-x0)/n; for i=1:n x(i+1)=x(i)+h; z0=y(i)+h*feval(f,x(i),y(i)); y(i+1)=y(i)+(feval(f,x(i),y(i))+feval(f,x(i+1),z0))*h/2; end shuchu=[x,y] fprintf( x(i) y(i)) function Dy=f(x,y) D