第9节多项式插值答辩.pptVIP

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,