数值分析-MALAB相关算法.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语言程序 #include 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 x 1 ; 2、幂法求特征值 function [m x biaozhi] mifa A,jingdu,cishu %幂法求矩阵最大特征值，其中 %m为绝对值最大的特征值，x为对应最大特征值的特征向量 %biaozhi表明迭代是否成功 if nargin 3 cishu 100; end if nargin 2 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 Dy x+y; 5、最小二乘 M文件： function c zxrc x,y,m %x是数据点横坐标，y数据点纵坐标 %m要构造的多项式的系数，c是多项式由高到低的系数所组成的向量 n length x ; b zeros 1:m+1 ; f zeros n,m+1 ; for k 1:m+1 f :,k x.^ k-1 ; end a f*f; b f*y; c a\b; c flipud c ; 6、矩阵相关的算法 （1）.求矩阵的