计算方法简明教程王能超 引论分解38页PPT
- 格式:ppt
- 大小:2.57 MB
- 文档页数:38
王能超计算方法第一章插值方法 1.1 Lagrange插值 1.2 逐步插值1.3 分段三次Hermite插值 1.4 分段三次样条插值第二章数值积分 2.1 Simpson公式 2.2 变步长梯形法 2.3 Romberg加速算法 2.4 三点Gauss公式第三章常微分方程德差分方法 3.1 改良的Euler方法 3.2 四阶Runge-Kutta方法 3.3 二阶Adams预报校正系统3.4 改良的四阶Adams预报校正系统第四章方程求根 4.1 二分法 4.2 开方法4.3 Newton下山法 4.4 快速弦截法第五章线性方程组的迭代法 5.1 Jacobi迭代5.2 Gauss-Seidel迭代 5.3 超松弛迭代 5.4 对称超松弛迭代第六章线性方程组的直接法 6.1 追赶法6.2 Cholesky方法 6.3 矩阵分解方法6.4 Gauss列主元消去法1第一章插值方法 1.1 Lagrange插值计算Lagrange插值多项式在x=x0处的值. MATLAB文件:〔文件名:Lagrange_eval.m〕 function [y0,N]= Lagrange_eval(X,Y,x0) %X,Y是插值点坐标 %x0是插值点%y0是Lagrange插值多项式在x0处的值 %N是Lagrange插值函数的权系数m=length(X); N=zeros(m,1); y0=0; for i=1:m N(i)=1; for j=1:m if j~=i; N(i)=N(i)*(x0-X(j))/(X(i)-X(j)); end endy0=y0+Y(i)*N(i); end用法》X=[…];Y=[…];》x0= ;》[y0,N]= Lagrange_eval(X,Y,x0)1.2 逐步插值计算逐步插值多项式在x=x0处的值. MATLAB文件:〔文件名:Neville_eval.m〕 function y0=Neville_eval(X,Y,x0)%X,Y是插值点坐标 %x0是插值点%y0是Neville逐步插值多项式在x0处的值 m=length(X); P=zeros(m,1); P1=zeros(m,1); P=Y; for i=1:m P1=P; k=1;for j=i+1:m k=k+1;2P(j)=P1(j-1)+(P1(j)-P1(j-1))*(x0-X(k-1))/(X(j)-X(k-1)); endif abs(P(m)-P(m-1))=X(i)&& x0=X(1,j)& x0。