第9节多项式插值技巧.ppt

Hermite插值的计算 Hermite插值结果，保留了f’(x)的特征 Hermite插值与多项式插值比较 多项式插值结果，没有保留了f’(x)的特征 Hermite插值与多项式插值比较 clear all close all clc t = [1.3, 1.6, 1.9]; y1 = [0.6200, 0.4554, 0.2818]; dy = [-0.5220, -0.5699, -0.5812]; A = zeros(6); for iii = 1 : 3 for jjj = 1 : 6 A(iii, jjj) = t(iii) ^ (6 - jjj); end end for iii = 4 : 6 for jjj = 1 : 5 A(iii, jjj) = (6 - jjj) * t(iii - 3) ^ (5 - jjj); end end y = [y1, dy] p = A \ y.; x = 1.3 :0.001 : 1.9; yy = polyval(p,x); plot(x, yy) hold on plot(t, y1, o) dyy = diff(yy) * 1000; figure plot(x(1 : 600), dyy) hold on plot(t, dy, o) figure yy2 = polyval(polyfit(t, y1, 2), x) plot(x, yy2, r) hold on plot(t, y1, o) dyy2 = diff(yy2) * 1000; figure plot(x(1 : 600), dyy2) hold on plot(t, dy, o) Hermite插值的计算 sin函数的三阶多项式插值和Hermite插值 Hermite插值的计算 clear all close all clc t = 1 : 3; t = t / 3 * pi / 0.5; y1 = sin(t) dy = cos(t) A = zeros(6); for iii = 1 : 3 for jjj = 1 : 6 A(iii, jjj) = t(iii) ^ (6 - jjj); end end for iii = 4 : 6 for jjj = 1 : 5 A(iii, jjj) = (6 - jjj) * t(iii - 3) ^ (5 - jjj); end end y = [y1, dy] p = A \ y. x = 1 : 30; x = x / 30 * pi / 0.5; yy = polyval(p,x) plot(x, yy) hold on plot(t, y1, o) yy2 = polyval(polyfit(t, y1, 2), x) plot(x, yy2, r) plot(x, sin(x), g) 拉格朗日插值程序 %lagrangen.m function y=lagrangen(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i);s=0; for k=1:n L=1; for j=1:n if j~=k L=L*(z-x0(j))/(x0(k)-x0(j)); end end s=s+L*y0(k); end y(i)=s; end y; 拉格朗日插值程序 %Chazhibijiao.m x=-5:0.1:5;z=0*x;y=1./(1+x.^2); plot(x,z,k,x,y,r) axis([-5 5 -1.5 2]);pause,hold on for n=2:2:10 x0=linspace(-5,5,n+1); y0=1./(1+x0.^2); x=-5:0.1:5; y1=lagrangen(x0,y0,x); plot(x,y1), pause end y2=1./(1+x0.^2);y=interp1(x0,y2,x); plot (x,y,k),hold off gtext(n=2),gtext(n=4),gtext(n=6) gtext(n=8),gtext(n=10) gtext(f(x)=1/(1+x^2)) 第四章 函数的数值逼近 第一节 多项式插值 数值逼近（Approximation） Approxi