Python模拟Lagrange插值结果并分析误差

Python模拟Lagrange插值代码及结果

题目:exp(x)(函数可变)在[0,4]的n阶多项式插值及误差

思路

主要思路分为两部分:

  1. 如何计算n阶多项式插值?
  2. 误差如何计算?

首先第一个部分:
对于lagrange插值函数的计算,首先我们观察公式:
L n ( x ) = ∑ k = 0 n y k l k ( x ) L_n(x)=\sum\limits_{k=0}^n y_k l_k(x) Ln(x)=k=0nyklk(x)
其中, l k ( x ) = ( x − x 0 ) … ( x − x k − 1 ) ( x − x k + 1 ) … ( x k − x n ) ( x k − x 0 ) … ( x k − x k − 1 ) ( x k − x k + 1 ) … ( x k − x n ) l_k(x)=\frac{(x-x_0)\dots(x-x_{k-1})(x-x_{k+1})\dots(x_k-x_n)}{(x_k-x_0)\dots(x_k-x_{k-1})(x_k-x_{k+1})\dots(x_k-x_n)} lk(x)=(xkx0)(xkxk1

实验一 误差分析 一、实验目的及要求 1.了解误差分析对数值计算的重要性。 2.掌握避免或减小误差的基本方法。 二、实验设备 安装有C、C++或MATLAB的计算机。 三、实验原理 误差是指观测值与真值之差,偏差是指观测值与平均值之差。根据不同的算法,得到的结果的精度是不一样的。 四、实验内容及步骤 求方程ax2+bx+c=0的根,其中a=1,b= -(5×108+1),c=5×108 采用如下两种计算方案,在计算机上编程计算,将计算结果记录下来,并分析产生误差的原因。 ////////////////////////////// 实验二 Lagrange插值 一、实验目的及要求 1.掌握利用Lagrange插值法及Newton插值法求函数值并编程实现。 2.程序具有一定的通用性,程序运行时先输入节点的个数n,然后输入各节点的值( ),最后输入要求的自变量x的值,输出对应的函数值。 二、实验设备和实验环境 安装有C、C++或MATLAB的计算机。 三、算法描述 1. 插值的基本原理(求解插值问题的基本思路) 构造一个函数y=f(x)通过全部节点,即 (i=0、1、… n) 再用f(x)计算插值,即 2. 拉格朗日(Lagrange多项式插值 Lagrange插值多项式: 3.牛顿(Newton)插值公式 //////////////////////////////////// 实验三 高斯消去法解方程组 一、实验目的及要求 1.掌握求解线性方程组的高斯消去法---列选主元在计算机上的算法实现。 2.程序具有一定的通用性,程序运行时先输入一个数n表示方程含有的未知数个数,然后输入每个线性方程的系数和常数,求出线性方程组的解。 二、实验设备和实验环境 安装有C、C++或MATLAB的计算机。 三、算法描述 1.高斯消去法基本思路 设有方程组 ,设 是可逆矩阵。高斯消去法的基本思想就是将矩阵的初等行变换作用于方程组的增广矩阵 ,将其中的 变换成一个上三角矩阵,然后求解这个三角形方程组。 2. 利用列选主元高斯消去法求解线性方程组
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值