艾尔米特插值的MATLAB实现,埃尔米特(Hermite)插值

Hermite插值满足在节点上等于给定函数值，而且在节点上的导数值也等于给定的导数值。对于高阶导数的情况，Hermite插值多项式比较复杂，在实际情况中，常常遇到的是函数值与一阶导数给定的情况。在此情况下，n个节点x1,x2,…,xn的Hermite插值多项式的表达式如下：

其中：

这样就很容易写出代码了，关键就是记住公式。

matlab中的实现：

Hermite.m

function f = Hermite(x,y,y_1,x0)

%求已知数据点的向后差分牛顿插值多项式

%已知数据点的x 坐标向量:x

%已知数据点的y 坐标向量:y

%已知数据点的导数向量:y_1

%求得的Hermite插值多项式或x0处的插值:f

syms t;

f = 0.0;

if(length(x) == length(y))

if(length(y) == length(y_1))

n = length(x);

else

disp('y和y的导数的维数不相等！');

return;

end

else

disp('x和y的维数不相等！');

return;

end

for i=1:n

h = 1.0;

a = 0.0;

%%计算hi和ai

for j=1:n

if( j ~= i)

h = h*(t-x(j))^2/((x(i)-x(j))^2);

a = a + 1/(x(i)-x(j));

end

f = f + h*((x(i)-t)*(2*a*y(i)-y_1(i))+y(i));

if(i==n)

if(nargin == 4)

f = subs(f,'t',x0);

else

f = vpa(f,6);

end

HermiteInsert.m

x=1:0.2:1.8;

y=[1 1.0954 1.1832 1.2649 1.3416];

y_1=[0.5 0.4564 0.4226 0.3953 0.3727];

f=Hermite(x,y,y_1)

f=Hermite(x,y,y_1,1.44)

x1=0:2*pi;

y1=sin(x1);

y1_1=cos(x1);

xx=0:0.5:2*pi;

yy=Hermite(x1,y1,y1_1,xx);

plot(x1,y1,'o',xx,yy,'r')

本文同步分享在博客“风翼冰舟”(CSDN)。

如有侵权，请联系 support@oschina.cn 删除。

本文参与“OSC源创计划”，欢迎正在阅读的你也加入，一起分享。

艾尔米特插值的MATLAB实现,埃尔米特(Hermite)插值相关推荐

c语言埃尔米特插值思路,【数学建模算法】(26)插值和拟合：埃尔米特(Hermite)插值和样条插值...
1.埃尔米特(Hermite)插值 1.1.Hermite插值多项式如果对插值函数,不仅要求它在节点处与函数同值,而且要求它与函数有相同的一阶.二阶甚至更高阶的导数值,这就是 Hermite 插值问 ...
插值与拟合 (一) ：拉格朗日多项式插值、Newton插值、分段线性插值、Hermite插值、样条插值、 B 样条函数插值、二维插值
插值:求过已知有限个数据点的近似函数. 拟合:已知有限个数据点,求近似函数,不要求过已知数据点,只要求在某种意义下它在这些点上的总偏差最小. 插值和拟合都是要根据一组数据构造一个函数作为近似,由于近似 ...
埃尔米特(Hermite)插值
Hermite插值满足在节点上等于给定函数值,而且在节点上的导数值也等于给定的导数值.对于高阶导数的情况,Hermite插值多项式比较复杂,在实际情况中,常常遇到的是函数值与一阶导数给定的情况.在此情 ...
数学建模准备插值（拉格朗日多项式插值，牛顿多项式插值，分段线性插值，分段三次样条插值，分段三次Hermite插值）
文章目录摘要(必看) 0 基础概念什么是插值插值用途什么是拟合插值和拟合的相同点插值和拟合的不同点 1 常用的基本插值方法 1.1 多项式插值法 1.1.1 拉格朗日多项式插值法多项式插 ...
埃尔米特(Hermite)插值及其MATLAB程序
%hermite.m %求埃尔米特多项式和误差估计的MATLAB主程序 %输入的量:X是n+1个节点(x_i,y_i)(i = 1,2, ... , n+1)横坐标向量,Y是纵坐标向量, %以f'(x ...
hermit插值 matlab,埃尔米特(Hermite)插值及其MATLAB程序
%hermite.m %求埃尔米特多项式和误差估计的MATLAB主程序 %输入的量:X是n+1个节点(x_i,y_i)(i = 1,2, ... , n+1)横坐标向量,Y是纵坐标向量, %以f'(x ...
埃尔米特插值c语言编程,131张艳-埃尔米特(Hermite)插值逼近的C语言程序.doc
PAGE PAGE 1 论文题目:埃尔米特(Hermite) 插值逼近的C语言程序院系: 数学科学学院专业: 数学与应用数学姓名: 张艳学号: 指导教师: 侯存贵完成时间: 2 ...
2021-07-21 分段三次Hermite插值及其与三次样条的比较
Numpy数组自身加和与乘积运算 http://www.weixueyuan.net/a/794.html python-scipy -hermite插值,分段三次Hermite插值 import n ...
【经典】纪念米加艾尔·塔利维尔吉耶夫《春天的17个瞬间》音乐会[部分]
米加艾尔·塔利维尔吉耶夫他是俄罗斯音乐的代言人,他用烂漫的童心谱写出恒久的感动.他是影视音乐的大师,他用诗意的旋律勾勒出不凡的人生轨迹.音乐故事,精彩呈现<春天的瞬间>系列节目-- 有作 ...

艾尔米特插值的MATLAB实现,埃尔米特(Hermite)插值

艾尔米特插值的MATLAB实现,埃尔米特(Hermite)插值相关推荐

最新文章

热门文章