![计算机如何表达函数_第1页](http://file4.renrendoc.com/view11/M02/3D/00/wKhkGWWjpAWAL9iAAAGU8Qko9Ao589.jpg)
![计算机如何表达函数_第2页](http://file4.renrendoc.com/view11/M02/3D/00/wKhkGWWjpAWAL9iAAAGU8Qko9Ao5892.jpg)
![计算机如何表达函数_第3页](http://file4.renrendoc.com/view11/M02/3D/00/wKhkGWWjpAWAL9iAAAGU8Qko9Ao5893.jpg)
![计算机如何表达函数_第4页](http://file4.renrendoc.com/view11/M02/3D/00/wKhkGWWjpAWAL9iAAAGU8Qko9Ao5894.jpg)
![计算机如何表达函数_第5页](http://file4.renrendoc.com/view11/M02/3D/00/wKhkGWWjpAWAL9iAAAGU8Qko9Ao5895.jpg)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第1章插值
实际中,f(x)多样,复杂,通常只能观测到一些离散数据;或者f(x)过于复杂而难以运算。这时我们要用近似函数g(x)来逼近f(x)。
自然地,希望g(x)通过所有的离散点x0x1x2x3x4xg(x)
f(x)定义:为定义在区间上的函数,为区间上n+1个互不相同的点,为给定的某一函数类。求上的函数满足问题是否存在唯一如何构造误差估计设则所以有唯一解,当且仅当m=n,且系数行列式不为0存在唯一定理定理1.1:为n+1个节点,
n+1维空间,则插值函数存在唯一,当且仅当与基函数无关与原函数f(x)无关基函数个数与点个数相同特点:对应于则Vandermonde行列式病态多项式插值的Lagrange型如何找?在基函数上下功夫,取基函数为要求则求,易知:记线性插值二次插值例:算法:fx=0.0for(i=0;i<=n;i++){
tmp=1.0;
for(j=0;j<i;j++)
tmp=tmp*(x-x[j])/(x[i]-x[j]);
for(j=i+1;j<=n;j++)
tmp=tmp*(x-x[j])/(x[i]-x[j]);
fx=fx+tmp*y[i];}returnfx;Lab02
Lagrange插值对函数构造插值,并求为近似的误差。插值节点取为:(1)(2)对N=5,10,20,40比较以上两组节点的结果。Chebyshev点SampleOutput(
representsaspace)第1组节点,误差为n=5
,
0.244934066848e+001n=10
,
0.534607244904e+001...第2组节点,误差为n=5
,
0.244934066848e+001n=10
,
0.534607244904e+001...误差解:求设易知有n+2个零点由a的任意性例:已知分别利用sinx的1次、2次Lagrange插值计算sin50
并估计误差。解:n=1分别利用x0,x1
以及x1,x2
计算
利用这里而
sin50=0.7660444…)185(50sin10
pL0.77614外推
/*extrapolation*/
的实际误差
0.01001
利用sin50
0.76008,内插
/*interpolation*/
的实际误差
0.00596内插通常优于外推。选择要计算的x
所在的区间的端点,插值效果较好。n=2)185(50sin20
pL0.76543
sin50=0.7660444…2次插值的实际误差
0.00061高次插值通常优于低次插值事后误差估计给定任取n+1个构造如:另取则近似则HomeWork1.证明线性无关2.若,记为节点上的k次插值多项式试证明:Lagrange
插值的缺点无承袭性。增加一个节点,所有的基函数都要重新计算为实数Newton型多项式插值易知同样承袭性:而且有:这样:称为k阶差商称为1阶差商定义:差商差商的一个性质:(用归纳法易证)
对称性:定义关键:找不同的元素相减作分母由归纳:Newton插值构造1、先构造差商表例子2点Newton型插值2、利用差商表的最外一行,构造插值多项式差商表求值算法:for(i=1;i<=n;i++){
for(j=n;j>=i;j--)
y[j]=(y[j]-y[j-1])/(x[j]-x[j-i]);}fx=y[n];for(i=n;i<=1;i--){
fx=y[i-1]+(x-x[i-1])fx;}问题:如果要做到增加一个点,而尽可能减少重复计算,要如何改进前面的算法?一些性质性质2误差性质3差商性质总结证明作为作业1.4Hermite插值
有时候,构造插值函数除了函数值的条件以外,还需要一定的连续性条件,如一阶导数值等,这种插值称为Hermite插值。称为二重密切Hermite插值例:设x0
x1
x2,已知f(x0)、f(x1)、f(x2)
和f’(x1),求多项式P(x)
满足P(xi)=f(xi),i=0,1,2,且P’(x1)=f’(x1),并估计误差。模仿Lagrange多项式的思想,设解:首先,P
的阶数=3
+=213)()()()()(=0iiixhx1f’xhxfxP
h0(x)有根x1,x2,且h0’(x1)=0
x1是重根。)()()(22100xxxxCxh--=又:h0(x0)=1C0h2(x)与h0(x)完全类似。其中hi(xj)=
ij,hi’(x1)=0,
(xi)=0,
’(x1)=1
h1
h1h1(x)有根x0,x2
))()(()(201xxxxBAxxh--+=由余下条件h1(x1)=1和
h1’(x1)=0可解。
(x)
h1有根x0,x1,x2
h1))()(()(2101xxxxxxCx---=
h1又:’(x1)=1C1
可解。与Lagrange分析完全类似
仿照Lagrange插值的做法,首先确定多项式插值空间的维数,注意到,我们的条件共有2(n+1)个条件,所以,最高次数为2n+1对二重密切Hermite插值
整个构造步骤如下:1、确定多项式的最高项次数,就是函数空间的维数2、假设一组基函数,列出插值多项式3、列出基函数满足的公式(画表),求基函数称为构造基函数方法误差分析类似Lagrange插值的分析方法二重密切Hermite插值误差例:在[5,5]上考察的Ln(x)。取
-5
-4
-3
-2
-1
0
1
2
3
4
5
-0.5
0
0.5
1
1.5
2
2.5
n
越大,端点附近抖动越大Ln(x)
f(x)
是否次数越高越好呢?分段低阶插值Runge现象例:等距节点构造10次Lagrange插值多项式-0.900.047061.57872-0.700.07547-0.22620-0.500.137930.25376-0.300.307690.235351901年,Runge等距高次插值,数值稳定性差,本身是病态的。分段线性插值每个小区间上,作线性插值(1)(2)在每个小区间上为一个不高于1次的多项式特性误差可以看出收敛,可惜只一阶精度,不够光滑。类似,可以作二重密切Hermite插值关键:
分段、低阶插值三次样条插值
分段低阶插值,收敛性好,但光滑性不够理想。在工业设计中,对曲线光滑性要求高,如:流线型设想这样一曲线:插值,次数不高于3次,整个曲线2阶连续导数,称为三次样条函数插值。每个小区间不高于3次,有4n个未知数,我们的已知条件如下:共3n-3+n+1
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 肺动脉高压患者的麻醉和围术期管理
- 门岗保安人员出入管理制度
- 2024年智能配网产品项目发展计划
- 老年人法律常识-合同篇
- 仙人球花的药用价值
- 基于“两网融合”的数据全生命周期安全建设需求
- 2024年地下智能输电系统项目发展计划
- 2024年压磁应力计项目建议书
- 2024年分纸机项目合作计划书
- 辽宁省葫芦岛市重点中学2022年中考数学模拟预测题含解析
- 监理抽检表 - 04路基土石方工程
- 音乐学科家长会
- 外研版(2019) 选择性必修第二册 Unit 2 Improving Yourself Developing ideas课件
- 2023十八项医疗核心制度考试题库及参考答案
- 2023光伏工程施工队合同书正规范本(通用版)
- 医院师带徒协议书
- 火灾的应急处置流程
- 静电的防止与利用(课件)(人教版2019必修第三册)
- EPC项目提供后续设计服务的保证措施
- 作物“源、流、库”理论及其应用
- 粤教版综合实践活动三年级上册第四单元第一课时《玩具总动员》课件
评论
0/150
提交评论