插值方法——Lagrange插值公式

问题

已知f(x)=e^x(3x-e^x)，利用插值节点x₀=1.00，x₁=1.02，x₂=1.04，x₃=1.06，构造三次Lagrange插值公式，由此计算f(1.03)的近似值，并给出其实际误差。

原理

根据线性插值和抛物线插值的基函数构造方法，令

其中(i=0,1,..n)为n次多项式，满足

可得：=

则：

根据上面知识可以得到本题的公式

误差：

程序框图

结果比较

误差：

结论

1. 误差为这个单位级别可以忽略。

2. =。所以插值方法算出来基本接近原值。

附件：程序

函数文件fun.m

function y=fun(x)

y=exp(x)*(3*x-exp(x));

主文件main.m

x=[1.00,1.02,1.04,1.06]; %

xt=1.03;

y=[0 0 0 0];

for i=1:4

y(i)=fun1(x(i)); %f(xi)值

end

sums=0;

for i=1:4

prods=1;

for j=1:4

if(i~=j)

prods=prods*(xt-x(j))./(x(i)-x(j)); %联加

end

sums=sums+prods*y(i); %联乘

end

sums

fun1(xt)

转载于:https://www.cnblogs.com/lecone/archive/2011/05/10/2041965.html

插值方法——Lagrange插值公式相关推荐

线性插值、抛物插值、Lagrange插值 | Lagrange拉格朗日插值法（一）
Lagrange(拉格朗日)插值法 Lagrange插值法是一种多项式插值方法. 1. 线性插值(两点插值或一次插值) 线性插值就是通过两个采样点 ( x 0 , y 0 ) (x_0,y_0) (x ...
Lagrange插值法
Lagrange插值公式: P(x)=∑k=0nLn,k(x)f(xk)P(x) = \sum_{k=0}^{n}L_{n,k}(x) f(x_{k})P(x)=∑k=0nLn,k(x)f(xk ...
用C++语言程序实现拉格朗日插值公式
C++程序实现Lagrange插值公式 Lagrange插值公式,是属于数值分析方面的内容,关于其应用,在这里就不多说.此处我想用C语言程序来实现n各插值节点插值公式的求解,并求出在某一个插值节点对应 ...
matlab 算法设计,计算方法——算法设计及其MATLAB实现
引论 0.1 算法重在设计 0.2 直接法的缩减技术 0.3 迭代法的校正技术 0.4 算法优化的松弛技术小结习题0 章插值方法 1.1 插值平均 1.2 Lagrange插值公式 1.3 逐步 ...
matlab用牛顿差值计算三次差值多项式,计算方法用Newton插值多项式求函数的近似值.docx...
计算方法用Newton插值多项式求函数的近似值计算方法课程设计题目: 用Newton插值多项式处理磁化曲线学院: 理学院班级: 学生姓名: 学生学号: 指导教师: 2017 ...
MATLAB两点式方程,§
§4.2Lagrange插值方法 4.2.1拉格朗日插值多项式的基本概念 1线性插值由存在性唯一性定理的证明过程可知,每一个代数插值问题实质上对应一个线性代数方程组的求解问题.对于插值方法的研究实际 ...
计算机数值方法期末考试,《计算机数值方法》测试题二
<计算机数值方法>测试题二上传人:文*** 文档编号:84429114 上传时间:2020-06-05 格式:DOC 页数:3 大小:57KB 下载提示(请认真阅读)1.请仔细阅读文档, ...
matlab的讲稿ppt,Matlab初步(讲稿200508)之四
八．举例(用Matlab作计算) 小题1:随机产生一个满足"标准正态分布"的100维数组,记作a . 再把数组a的这100个数的次序完全颠倒,得到一个新数组,记作b,输出b . ...
matlab中a2=poly(p2),插值与拟合matlab实现
插值与拟合的Matlab实现王正盛编写在科技工程中,除了要进行一定的理论分析外,通过实验.观测数据,做分析.处理也是必不可少的一种途径.由于实验测定实际系统的数据具有一定的代表性,因此在处理时必须 ...

插值方法——Lagrange插值公式

问题

原理

程序框图

结果比较

结论

附件：程序

插值方法——Lagrange插值公式相关推荐

最新文章

热门文章