8112 Lagrange插值算法及应用 已知x2,y点,可求出c向量上各点处的 插值为: 插值算法为 P(a)=>yi j=1,j≠讠 MATLAB函数调用格式 y=lagrange(o, yo, c) 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 8.1.1.2 Lagrange插值算法及应用 已知 点，可求出 向量上各点处的 插值为： 插值算法为 MATLAB函数调用格式
Lagrange插值算法的 MATLAB实现 function y=lagrange(xo, yO, x) ii=1: length(x0); y=zeros(size(x)) for i=ii ij=find (ii"=i); y1=1 for j=1: length(ij) y1=y1.*(x-x0(ij(j); end y=y+y1*yo(i)/prod(xo(i)-xo(ij)) end 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 Lagrange插值算法的MATLAB实现
4例8.3 对f(x)=1/(1+25x2),-1≤x≤1,进彳 Lagrange插值 >>x0=-1+2*[0:10]/10;y0=1./(1+25*x0.^2); x=-1: 01: 1; y=lagrange(xo, yo, x) ya=1./(1+25*x.^2);plot(x,ya,x,y,’:’) 调用 interp1()函数: s >>y1=interp1(x, yo, x,,cubic') y2=interp1(xO, yO, x,",); plot(x, ya, x,y1,':',x, y2 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 对 ，进行 Lagrange插值 调用interp1()函数： 例 8.3
812已知样本点的定积分计算 编写 MATLAB函数 function y=quadspln(xo, yo, a, b) f=@(x, xO, yO)interp1(xo, yo, x,'spline) y=quadI(f, a,b, 1e-8,[, xO, yO) 函数调用格式: 1=quadspln(ao, y0, a, b) 2/20/2021星期六, 高等应用数学问题的 MATLAB求解 2008-9-6,13:1107 Slide 1(of 11) 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 8.1.2 已知样本点的定积分计算 编写MATLAB函数： 函数调用格式：
4例84 利用样条插值算法求解 T sin dr →用30个采样点求解: LAE >>x0=0: pi/30: pi; yO=sin(x0); I=quadspln(xO, yo, 0, pi) 2/20/2021星期六, 高等应用数学问题的 MATLAB求解 2008-9-6,13:1107 Slide 1(of 11) 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 例 8.4 利用样条插值算法求解 用30个采样点求解：