MATLAB实例:Hermite插值多项式
作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/
问题描述:清华大学出版社李庆扬《数值分析》第五版教材P49习题14:
求次数小于等于3的多项式P(x),使其满足条件 P(0)=0, P’(0)=1, P(1)=1, P’(1)=2
1. MATLAB程序
Hermite_kailugaji.m
function f = Hermite_kailugaji(t,y,y_1,t0)
%t y为坐标向量 y_1为一阶导的值 t0为插值点的t坐标|| f0为t0对应的值
syms x;
f = 0.0;
if(length(t) == length(y))
if(length(y) == length(y_1))
n = length(t);
else
disp('y和y的导数的维数不相等!');
return;
end
else
disp('x和y的维数不相等!');
return;
end
for i=1:n
h = 1.0;
a = 0.0;
for j=1:n
if( j ~= i)
h = h*(x-t(j))^2/((t(i)-t(j))^2);
a = a + 1/(t(i)-t(j));
end
end
f = f + h*((t(i)-x)*(2*a*y(i)-y_1(i))+y(i));
if(i==n)
if(nargin == 4)
f = subs(f,'x',t0);
else
f = vpa(f,6);
end
end
end
demo.m
clear
clc
% 求次数小于等于3的多项式P(x),使其满足条件
% P(0)=0, P’(0)=1, P(1)=1, P’(1)=2
x=[0 1];
y=[0 1];
y_1=[1 2];
%x, y为坐标向量, y_1为一阶导的值, x0为插值点的x坐标|| f0为x0对应的值
f=Hermite_kailugaji(x,y,y_1);
f=collect(f); % Matlab 合并同类项,将多项式化为一般式
fprintf('Hermite插值多项式为: y=%s
', f);
f=Hermite_kailugaji(x,y,y_1,2);
fprintf('在x=2处的Hermite插值为: %f
', f);
2. 结果
Hermite插值多项式为: y=x^3 - 1.0*x^2 + 1.0*x
在x=2处的Hermite插值为: 6.000000
提示:MATLAB多项式整理小技巧:collect()是化为一般式,factor()是因式分解,合并同类项。