插值与数值积分与matlab.pptVIP

实验3-插值与数值积分 1. 插值方法 插值问题 已知n+1个节点(xj,yj) (j=0,1,…,n), 其中xj互不相同, 不妨设 a=x0x1…xn=b, 求任一插值点 处的插值 . 求解思路 构造一个相对简单的函数y=f(x), 称为插值函数, 使f通过全部节点f(xj)=yj,再用f(x)计算插值 . 实验3-插值与数值积分 1. 插值方法 常用的插值方法：Lagrange多项式插值; 分段线性插值; 三次样条插值. ① Lagrange多项式插值 构造插值函数f(x)=Ln(x)= , 其中 , i, j=0,1,…,n 实验3-插值与数值积分 1. 插值方法 ① Lagrange多项式插值 例: 用g(x)=sin(pi*x/2)产生等距节点, 做插值多项式Ln(x)与g(x)比较. function y=lagr(x0,y0,x) %函数M文件lagr.m n=length(x0);m=length(x); for i=1:m %得到X每个点相应的L函数值 z=x(i);s=0; for k=1:n %Lag函数 p=1; for j=1:n %Lag基函数 if j~=k, p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end 实验3-插值与数值积分 1. 插值方法 ① Lagrange多项式插值 m=4:2:10; %应用程序 for k=1:4 n=m(k);x0=linspace(-5,5,n); %产生节点x坐标 y0=sin(pi*x0/2); %产生对应节点y坐标 x=-5:0.1:5;y=sin(pi*x/2); %计算比照函数 y1=lagr(x0,y0,x); %调用函数计算插值 subplot(2,2,k),plot(x,y,x,y1,r--), hold on end gtext(L_4(x)), gtext(L_6(x)), gtext(L_8(x)), gtext(L_1_0(x)) 实验3-插值与数值积分 1. 插值方法 ① Lagrange多项式插值 用Lagrange多项式近似g(x), 随着节点个数的增加, Ln(x)的次数变大, 多数情况下误差会变小. 但n增加时, Ln(x)的起伏增加, 有时会出现很大的振荡, Runge现象的例子: g(x)=1/(1+x^2). 用前面的程序计算: 实验3-插值与数值积分 1. 插值方法 ① Lagrange多项式插值 m=4:2:10; %Runge 的例子: g(x)=1/(1+x^2) for k=1:4 n=m(k);x0=linspace(-5,5,n); %产生节点x坐标 y0=1./(1+x0.^2); %产生对应节点y坐标 x=-5:0.1:5;y=1./(1+x.^2); %计算比照函数 y1=lagr(x0,y0,x); %调用函数计算插值 subplot(2,2,k),plot(x,y,x,y1,r--), hold on end gtext(L_4(x)), gtext(L_6(x)), gtext(L_8(x)), gtext(L_1_0(x)) 实验3-插值与数值积分 1. 插值方法 ② 分段线性插值 构造插值函数f(x)=In(x)= , 其中 即将每两个相邻的节点用直线连起来, 如此形成的一条折线就是分段线性插值函数. In(x)有良好的收敛性, 实验3-插值与数值积分 1. 插值方法 ② 分段线性插值 Runge的例子用分段线性插值 function