计算方法 6.插值法
需求的由来:
插值法定义:
唯一性:
其实就是一个多元方程组嘛。
这个完全就是数字游戏吧,后面那一项永远是0.
拉格朗日求法:
Ln(x)就是p(x)。里面有yi作为系数,求这样的li应该是什么样。
可以看到,x不等于xi的时候,x总会与某一个xj相等,就会变成零,等于xi的时候就是1。
记的时候可以这样想:每带入一个xi,我只能让一个yi系数为1,那肯定是(x-x0)(x-x1)……(x-xn)(中间没有x-xi),这样其他的就都是零了,自己还要是1,那就还要再除以(xi-x0)(xi-x1)……,就行了。
那插值的余项怎么说?:
(t-xi)连乘的式子求导到最后是(n+1)!,这个值只和最高次数的t有关,也就是从(n+1)一直乘到1.最后的余项就是拉格朗日余项。
牛顿插值法:
拉格朗日算法增加一个结点的时候每个系数都要改,非常麻烦。
公式看起来很复杂,其实用上面的原式推一下就可以了。
导数的那个推论不晓得怎么整,但你可以通过n=1的情况(就是拉格朗日公式)记忆。我突然想到了一个证明方法,f[x0,……xn]就是Nn(x)最后一项的系数,也就是原式最高项的系数,除以n!就是要把n次求导中产生的乘积去掉,留下的就是最高项的真正系数了。
可以这样计算:
注意:虽然算法不同,但是我们之前讲过大家都是n阶的话那就是唯一的,所以这两种方法算出来的式子是一样的。:
那个wn+1(x)就是(x-x0)……(x-xn),这个x一般都是要在x0和xn之间的。
等距节点公式:
计算很简单,但暂时看不出是拿来干嘛的。
从最开始的式子我们可以算出三角fi=fi+1-fi,三角^2fi=三角fi+1-三角fi=fi+2-2*fi+1+fi,以此类推可以推出上面的公式。(x+h)^n-x^n结果当然是n-1阶多项式,第二行的性质也很简单。f[x0,x1……]公式里的那个k!是没有问题的,每一次分母都会变得大一份,一开始是h,然后是2h,3h……所以有阶乘。
插值多项式的收敛与稳定:
稳定性:
分段插值:
最后那个结论使用Rn(x)=f(n+1)(c)/(n+1)! *w(x)推出来的,就是一开始的那个余项公式。
函数逼近:
最佳平方逼近:
递推公式确定,那正交多项式也就唯一了。这里解释一下什么叫正交:
<>是内积(也就是向量直接用坐标相乘)的意思。
这个方法就是把多项式看作向量。这里手写推一下:
接下来就是以此类推,(xi,xj)不是内积神似内积。
接下来介绍几种常见的正交多项式:
这个其实就是刚刚的那个例子,但是有了公式,而且没有要求首项系数为1。
切比雪夫多项式:
判断函数族线性无关:
最佳平方逼近:
这个平方逼近就是要差值最小,平方是为了体现绝对值。
先来考虑连续函数:
对a(也就是那些系数)求导,然后就求出了这个罪孽深重的东西。
那之前讲的正交基地是干嘛的?就是这里用的,正交基地只有对角线上的元素是有数字的,其他的都是0,方程就很好解:
-1到1正好和拉格朗日相符,可以直接用他的基地,这样就有现成的正交方程,不用解复杂的方程组了。
曲线拟合的最小二乘逼近:
注意注意,之前是求积分,现在是求和,()的含义也从积分变成了求和,m个和的原因是点给了m个。(应该是从1加起,可能是打错了)
计算方法 6.插值法相关推荐
- 【数值分析/计算方法】插值法及其余项MATLAB仿真实验
1.Lagrange插值法 方法:对给定的n个插值节点x1,x2,⋯,xn以及它们对应的函数值y1,y2,⋯,yn,利用构造的n-1次Lagrange插值多项式,可以通过下面式子求得插值区间内任意x的 ...
- 【计算方法】插值法多项式的求法--利用Lagrange插值和Newton插值
[说明]插值与拟合 在科学研究和工程中,许多问题都可以用y = f(x)来表示其某种内在规律的数量关系,不仅函数f(x)是各种各样的,而且有的函数很复杂,甚至没有明显的解析表达式 因此可以采用两种方法 ...
- 数值计算方法-拉格朗日插值法
要求: 编写拉格朗日插值函数,输入为(x,f,order1) 其中x代表节点数据,为n*2矩阵 f代表是否分段 oreder1:若分段,oreder1代表函数阶次,只能为1或者2 若不分段,则函数阶次 ...
- SAP Portfolio Analyzer 中 索提诺比率(Sortino Ratio)
索提诺比率(Sortino Ratio)与夏普比率或特雷诺比率类似,但是这个收益是和下行风险进行比较:意味着在计算风险波动率时它所采用的不是标准差,而是下行标准差.也就是说只是考虑不利的波动.因此会更 ...
- 信用评级模型_不寻常的方式来提高您的信用评分
信用评级模型 There are various widely known ways to boost your credit score, from making loan repayments o ...
- SAP Portfolio Analyzer信息比率(Information Ratio)
信息比率(Information Ratio)和Sharpe夏普比率类似,不过和超额收益比较的是跟踪误差,或者可以说是相对风险. 跟踪误差是指组合超额收益率与平均超额收益率之间的差异的标准差,学术认为 ...
- 分段线性插值法实验报告_计算方法实验报告习题1(浙大版)
计算方法实验报告 实验名称: 实验1 从函数表出发进行插值 1 引言 某个实际问题中,函数f (x)在区间[a,b]上存在且连续,但难以找到其表达式,只能通过实验和观测得到有限点上的函数表.有些情况虽 ...
- 数值计算方法第四章—插值法
插值法 本文参考书为马东升著<数值计算方法> 代数插值 插值 已知某些(有限)点的函数值求其余点的函数值 代数多项式插值 n+1n+1n+1 个互异节点处满足插值条件 P(xi)=yi,i ...
- python插值程序_计算方法(2)——插值法(附Python程序)
给定一些数据,生成函数的方式有两种:插值,回归. 插值而得到的函数通过数据点,回归得到的函数不一定通过数据点. 下面给出拉格朗日插值,牛顿插值和Hermite插值的程序, 具体原理可参考课本,不再赘述 ...
- 【计算方法】拉格朗日插值法
概念 设f(x)在N+1个点(x0,y0)-(xn,yn)处的值已知,其中值xk在区间[a,b]上,xk互不相同,满足a≤x0<x1<-<xn≤b,yk=f(xk).求任一插值点x对 ...
最新文章
- 聊一聊 Python 安装中的 --enable-shared
- 存储安全 系统的最后一道防线
- brew安装mysql 卸载_Ubuntu环境下MySQL卸载安装配置远程访问三步曲
- 【Linux】一步一步学Linux——unzip命令(68)
- 每周一算法之六——KMP字符串匹配算法
- HashMap jdk1.7和1.8概述
- 深度剖析Java数据结构之表(三)——ArrayList泛型类的实现
- 剑指Offer替换空格问题
- AD+邮件服务器搭建方案
- Qt第六课 QDateTime使用
- oracle ora-27544,数据库启动后,大量报错ORA-27300
- C#编写的打印排版控件拖拽软件
- AWWWB.COM网站克隆器
- JButton的使用
- routing-controllers工作原理解析
- 华为首款台式机计算机发布,华为首款商用台式机发布,或再成商用PC市场破局者...
- bootstrap——模态框
- 从坚果3的发布来看,锤子未来的发展将依然艰难
- 计算机基础应用在线免费答题,计算机应用基础简答题附答案.doc
- C++学习的一些经历