函数逼近与曲线拟合

合集下载

常用函数的逼近和曲线拟合

常用函数的逼近和曲线拟合

常用函数的逼近和曲线拟合在数学中,函数逼近和曲线拟合都是常见的问题。

函数逼近是指找到一个已知函数,尽可能地接近另一个函数。

而曲线拟合则是给定一组数据点,找到一条曲线来描述这些数据点的分布。

本文将讨论常用的函数逼近和曲线拟合方法。

一、函数逼近1. 插值法插值法是最简单的函数逼近方法之一。

它的基本思想是:给定一组已知点,通过构造一个多项式,使得该多项式在这些点处的函数值与已知函数值相等。

插值法的优点是精度高,缺点是易产生龙格现象。

常用的插值多项式有拉格朗日插值多项式和牛顿插值多项式。

拉格朗日插值多项式的形式为:$f(x)=\sum_{i=0}^{n}y_{i}\prod_{j=i,j\neq i}^{n}\frac{x-x_{j}}{x_{i}-x_{j}}$其中,$x_{i}$是已知点的横坐标,$y_{i}$是已知点的纵坐标,$n$是已知点的数量。

牛顿插值多项式的形式为:$f(x)=\sum_{i=0}^{n}f[x_{0},x_{1},...,x_{i}]\prod_{j=0}^{i-1}(x-x_{j})$其中,$f[x_{0},x_{1},...,x_{i}]$是已知点$(x_{0},y_{0}),(x_{1},y_{1}),...,(x_{i},y_{i})$的差商。

2. 最小二乘法最小二乘法是一种常用的函数逼近方法。

它的基本思想是:给定一组数据点,找到一个函数,在这些数据点上的误差平方和最小。

通常采用线性模型,例如多项式模型、指数模型等。

最小二乘法的优点是适用性广泛,缺点是对于非线性模型要求比较高。

最小二乘法的一般形式为:$F(x)=\sum_{i=0}^{n}a_{i}\varphi_{i}(x)$其中,$a_{i}$是待求的系数,$\varphi_{i}(x)$是一组已知的基函数,$n$是基函数的数量。

最小二乘法的目标是使得$\sum_{i=1}^{m}[f(x_{i})-F(x_{i})]^{2}$最小,其中$m$是数据点的数量。

曲线拟合实验报告[优秀范文5篇]

曲线拟合实验报告[优秀范文5篇]

曲线拟合实验报告[优秀范文5篇]第一篇:曲线拟合实验报告数值分析课程设计报告学生姓名学生学号所在班级指导教师一、课程设计名称函数逼近与曲线拟合二、课程设计目的及要求实验目的: ⑴学会用最小二乘法求拟合数据的多项式,并应用算法于实际问题。

⑵学会基本的矩阵运算,注意点乘与叉乘的区别。

实验要求: ⑴编写程序用最小二乘法求拟合数据的多项式,并求平方误差,做出离散函数与拟合函数的图形;⑵用MATLAB 的内部函数polyfit 求解上面最小二乘法曲线拟合多项式的系数及平方误差,并用MATLAB的内部函数plot作出其图形,并与(1)结果进行比较。

三、课程设计中的算法描述用最小二乘法多项式曲线拟合,根据给定的数据点,并不要求这条曲线精确的经过这些点,而就是拟合曲线无限逼近离散点所形成的数据曲线。

思路分析 : 从整体上考虑近似函数)(x p 同所给数据点)(i iy x , 误差i i iy x p r -=)(的大小,常用的方法有三种:一就是误差i i iy x p r -=)(绝对值的最大值im ir≤≤ 0max ,即误差向量的无穷范数;二就是误差绝对值的与∑=miir0,即误差向量的 1成绩评定范数;三就是误差平方与∑=miir02的算术平方根,即类似于误差向量的 2 范数。

前两种方法简单、自然,但不便于微分运算,后一种方法相当于考虑 2 范数的平方,此次采用第三种误差分析方案。

算法的具体推导过程: 1、设拟合多项式为:2、给点到这条曲线的距离之与,即偏差平方与:3、为了求得到符合条件的 a 的值,对等式右边求偏导数,因而我们得到了:4、将等式左边进行一次简化,然后应该可以得到下面的等式5、把这些等式表示成矩阵的形式,就可以得到下面的矩阵:⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∑∑∑∑∑∑∑∑∑∑∑=====+==+====niininiiknikinikinikinikiniiniinikiniiyyyaax x xx x xx x11i11012111111211 1an MMΛM O M MΛΛ 6.将这个范德蒙得矩阵化简后得到⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡n kkn nkkyyyaaax xx xx x M MΛM O M MΛΛ21102 21 1111 7、因为 Y A X = * ,那么 X Y A / = ,计算得到系数矩阵,同时就得到了拟合曲线。

实验二函数逼近与曲线拟合

实验二函数逼近与曲线拟合

《数值分析》课程设计实验报告实验二函数逼近与曲线拟合一、问题提出从随机的数据中找出其规律性,给出其近似表达式的问题,在生产实践和科学实验中大量存在,通常利用数据的最小二乘法求得拟合曲线。

在某冶炼过程中,根据统计数据的含碳量与时间关系,试求含碳量与时间t 的拟合曲线。

二、实验步骤先写出线性最小二乘法的M文件function c=lspoly(x,y,m)% x是数据点的横坐标组成的向量,y是纵坐标组成的向量% m是要构成的多项式的次数,c是多项式由高到低次的系数所组成的向量n=length(x);b=zeros(1:m+1);f=zeros(n,m+1);for k=1:m+1f(:,k)=x.^(k-1);enda=f'*f;b=f'*y';c=a\b;c=flipud(c);方法一:近似解析表达式为:y(t)=a1t+a2t2+a3t3第二步在命令窗口输入:lspoly([0,5,10,15,20,25,30,35,40,45,50,55],[0,1.27,2.16,2.86,3.44 ,3.87,4.15,4.37,4.51,4.58,4.02,4.64],2)回车得到:ans =0.0000-0.00520.26340.0178即所求的拟合曲线为y=-0.0052t2+0.2634t+0.0178在编辑窗口输入如下命令:>>x=[0,5,10,15,20,25,30,35,40,45,50,55];y=[0,1.27,2.16,2.86,3.44, 3.87,4.15,4.37,4.51,4.58,4.02,4.64];>> t=0:0.1:55;>> z=-0.0052*t.^2+0.2634*t+0.0178;>> plot(x,y,'ro',t,z);grid命令执行得到如下图(图2-1)0102030405060拟合多项式与数据点的关系方法二:假设近似表达式为:y(t)=c0+c1t+c2t2第一步在命令窗口输入:>>lspoly([0,5,10,15,20,25,30,35,40,45,50,55],[0,1.27,2.16,2.86,3. 44,3.87,4.15,4.37,4.51,4.58,4.02,4.64],2)回车得到:ans =-0.00240.20370.2305即所求的拟合曲线为y=-0.0024t2+0.2037t+0.2305在编辑窗口输入如下命令:>>x=[0,5,10,15,20,25,30,35,40,45,50,55];y=[0,1.27,2.16,2.86,3.44,3.87,4.15,4.37,4.51,4.58,4.02,4.64];>> t=0:0.1:55;>> z=-0.0024*t.^2+0.2037*t+0.2305;>> plot(x,y,'ro',t,z);grid命令执行得到如下图(图2-2)拟合多项式与数据点的关系三、实验结论在利用数据的最小二乘法求拟合曲线时,选取合适的近似表达式很重要,应通过不断的试验找出较为合适的近似表达式,这样才能尽可能的提高拟合精度。

5函数逼近与曲线拟合.docx

5函数逼近与曲线拟合.docx

第5章函数逼近与曲线拟合上一章讨论的是函数插值问题,通常都是用一个多项式来代替一个已知的函数,它们在 给定的插值基点上有相同的函数值,是对原函数的一-种近似。

然而,在实际应用中插值问题 仍有明显的缺点:对于有解析式的函数而言,在其它点上误差可能很大,如龙格现象;对于 离散(表)函数而言,给定的数据点本身是有误差的,刚性地让插值函数通过这些点不仅没 有意义,而且会影响对原函数的近似程度。

另外,泰勒展示也是对连续函数的一种低阶近似, 它在展开点附近误差较小,但在展开点远处,误差会很大。

本章讨论在新的函数谋旁度最条件下的函数近似问题,对连续函数称之为函数逼近问题, 对于离散函数称之为dii 线拟合问题。

主要内容有:函数最佳逼近的概念,正交多项式,最佳 均方逼近少最小二乘曲线拟合问题等。

5.1函数最佳逼近的概念希望能有一种方法寻求出一个近似多项式,使它在整个区间上既均匀的逼近/(%),所需 的计算呆又小,这就是函数逼近要解决的问题。

为了刻划“均匀逼近”,设P n (x)是定义在区 间[a,b ]上原函数/(x)的近似多项式。

我们用||/(x) -p n (x)||来度量p n (x)与/(x)近似逼近 程度。

这样,自然地会有下面两种不同的度暈标准:fM- p n (x)使丿IJ 这个度量标准的函数逼近称为均方逼近或平方逼近;/W 一 p n (x) = max f(x) 一 p n (x) 使用这个度量标准的函数逼近称为一致逼近或均匀逼近o关于一致逼近的问题,在数学分析中有以下结论。

设函数/(X )在区间[a,b ]上连续,若£>0,则存在多项式P(x)使|/(x)-P(x)|<£,在区间[a,b ]上一致成立。

对于函数插值而 言,如果插值余项也能满足对任意的£〉0, \R n (x)\ = \f(x)-p n M\<e 都成立的话,贝闹 值多项式P n M 是/(Q 的一致逼近多项式。

计算方法-最佳一致逼近多项式-切比雪夫多项式

计算方法-最佳一致逼近多项式-切比雪夫多项式

xn )
|
要使 max 1 x 1
|
(x
x0 )(x
x1) … (x
xn )
|
取极小值, 只需令:
(x x0 )(x x1) … (x xn)
1 2n
Tn1(x),
最佳一致 逼近0的 多项式
而上式成立的充分必要条件是x0, x1,…xn是切比雪夫 多项式的0点。
将Lagrange插值多项式Ln(x)的节点取为Tn1(x) 的0点 :
最佳一致逼近多项式
§3 最佳一致逼近多项式
一、基本概念及其理论
不超过n次的实系 数多项式的全体
本节讨论f(x) C[a, b], 求多项式pn* (x) Hn , 使得误差
||
f(x)
pn* (x)
||
min
pn Hn
||
f(x)
pn(x)
||
此即所谓最佳一致逼近 或切比雪夫逼近问题 。
Hn
设f(x) Cn1[a, b], 则函数通过变换
x a b b a t, 1 t 1
2
2
化为
f(x)
f(a b 2
b 2
a t)
g(t)
针对g(t) 使用定理7
例如:为将[0, 1] [-1, 1],可以令:

x
0
2
1
1 0 2
t
1 (t 2
1)
f(x) f(1 (t 1)) g(t) , 1 t 1. 2
f(x)
p(x)
|
可以证明存在唯一的(a*0 , a1* , … , an* ), 使得
(a*0 , a1* , … , an* )
min{max

第三章 1 逼近论

第三章 1 逼近论

( x)dx

2mn1m!n!11
dm dxm
[(x2
1)m
]
dn dxn
[(x2
1)n
]dx

1
dm
2mn m!n!dxm
[( x 2

1)m
]ddxnn11 [( x2

1)n
1
]
1

2m

1 n m!
n!
11
dm1 dx m 1
[(
x
2

i 1
1
||
x
||2

n

xi2

2

i1
称为2 范数.
类似地,对C[a,b]上的f ( x),可定义三种常用范数:
|| f || max | f ( x) |, 称为 范数,
a xb
|| f ||1 ab| f ( x) | dx, 称为1 范数,
i 1
定义4 设( x)是区间[a,b]上的非负函数, 如果满足条件
(1)
ab xk ( x)dx存在,
k

0,1,2,; 可以有限或
无限区间
(2) 对于[a,b]上的非负连续函数g( x),若abg( x)( x)dx 0,
则在[a,b]上g( x) 0;
就称( x)为[a,b]上的权函数.
f
,
x)

n

k0
f

k n

Pk
(
x),
(1.3)
其中Pk
(
x)


n k

xk

小波分析之函数逼近与曲线拟合

小波分析之函数逼近与曲线拟合


=
max
f (x)
a≤ x≤b
绝对值与
n上范数的扩充关系 R
• 数a的绝对值(a离开原点0的距离):∣a∣ • 数a与b的差异(距离): ∣a-b∣ • 向量A=( 1, a2,…,an)的范数(A离开0向量 A=(a , 的范数 A=( 的距离) : n • x = ∑ x i
1 i = 1
x x
距离空间定义
• ฀ 设X是非空集合,对于X中的任意两元素x与y ,按某一法则都对应唯一的实数ρ(x, y),并满足 以下三条公理: • 1.非负性:ρ(x, y) ≥0,ρ(x, y) =0当且仅当x=y; • 2.对称性:ρ(x, y) =ρ(y, x); • 3.三角不等式;对任意的x, y, z ρ(x, y) ≤ρ(x, z) + ρ(z, y), 则称ρ(x, y)为x与y间的距离(或度量),并称X是 以ρ为距离的距离空间(或度量空间),记为(X, ρ).
2
2
2

内积空间的性质
定理 设 X 为内积空间,{u1 , u2 ,⋯ , un } ⊆ X , 格拉姆(Gram)矩阵
(u1 , u1 ) (u2 , u1 ) ⋯ (u n , u1 ) (u1 , u2 ) (u2 , u2 ) ⋯ (u n , u2 ) G= ⋮ ⋮ ⋮ (u , u ) (u , u ) ⋯ (u , u ) 2 n n n 1 n
内积空间
设X 是定义在实(或复)数域K上的线性空 间,若对于X中 任意一对有序元素x,y, 恒对应 数域K的值(x, y),且满足: • (x, x) ≥0,且(x, x)=0的充要条件是x=0; • (ax, y) = a(x, y);฀ • (x+y, z) = (x, z) + (x, z).฀ ฀ 则称X为内积空间,(x, y)称为x, y的内积. 正交: 正交 若(x, y)=0,称x与y正交.

第3讲-函数逼近

第3讲-函数逼近

5
函数插值是插值函数 p(x) 与被插函数 f(x) 在节点处函数 值相同,即 P ( x i )
f ( x i ), ( i 0 ,1, , n ) ;而曲线拟合函数 ( x ) 不
要求严格地通过所有数据点 ( x i , y i ) 。 即,拟合函数 (x) 在xi处的偏差(亦称残差)
5
4
实 验 数 据
3
y
2 1
0
1
2
3
4
5
x
17
例 设某实验数据如下: i xi yi 1 0 5 2 1 2 3 2 1 4 3 1 5 4 2 6 5 3
6 6 6 2 a0m a1 xi a2 xi yi i 1 i 1 i 1 6 6 6 6 2 3 a0 xi a1 xi a2 xi xi yi i 1 i 1 i 1 i1 6 6 6 6 2 3 4 2 a0 xi a1 xi a2 xi xi yi i 1 i 1 i 1 i1
13.8434
70.376
i
i 1 4
2 i
i
a0 3.9374 a1 7.4626
4 a 0 7 .32 a1 70 .376 7 .32 a 0 13 .8434 a1 132 .12985
x y
i
132.12985
y 3.9374 7.4626 x
2
2

e

1 2


i0
n
2 i

( x
i0
n
i
) f ( x i )
2
最小。
这种要求误差(偏差)平方和最小的拟合称为曲线拟合的 最小二乘法。

函数逼近与曲线拟合

函数逼近与曲线拟合

函数逼近与曲线拟合3.1函数逼近的基本概念3.1.1 函数逼近与函数空间在数值计算中常要计算函数值,如计算机中计算基本初等函数及其他特殊函数;当函数只在有限点集上给定函数值,要在包含该点集的区间上用公式给出函数的简单表达式,这些都涉及到在区间上用简单函数逼近已知复杂函数的问题,这就是函数逼近问题.上章讨论的插值法就是函数逼近问题的一种.本章讨论的函数逼近,是指“对函数类A中给定的函数,记作,要求在另一类简单的便于计算的函数类B中求函数,使与的误差在某种度量意义下最小”.函数类A通常是区间上的连续函数,记作,称为连续函数空间,而函数类B通常为n次多项式,有理函数或分段低次多项式等.函数逼近是数值分析的基础,为了在数学上描述更精确,先要介绍代数和分析中一些基本概念及预备知识.数学上常把在各种集合中引入某些不同的确定关系称为赋予集合以某种空间结构,并将为样的集合称为空间.例如将所有实n维向量组成集合,按向量加法及向量与数的乘法构成实数域上的线性空间,记作,称为n维向量空间.类似地,对次数不超过n(n为正整数)的实系数多项式全体,按通常多项式与多项式加法及数与多项式乘法也构成数域上的一个线性空间,用表示,称为多项式空间.所有定义在上的连续函数集合,按函数加法和数与函数乘法构成数域上的线性空间,记作.类似地,记为具有p阶的连续导数的函数空间.定义1设集合S是数域P上的线性空间,元素,如果存在不全为零的数,使得, (3.1.1)则称线性相关.否则,若等式(3.1.1)只对成立,则称线性无关.若线性空间S是由n个线性无关元素生成的,即对都有则称为空间S的一组基,记为,并称空间S为n维空间,系数称为x在基下的坐标,记作,如果S中有无限个线性无关元素,…,则称S为无限维线性空间.下面考察次数不超过n次的多项式集合,其元素表示为, (3.1.2)它由个系数唯一确定.线性无关,它是的一组基,故,且是的坐标向量,是维的.对连续函数,它不能用有限个线性无关的函数表示,故是无限维的,但它的任一元素均可用有限维的逼近,使误差(为任给的小正数),这就是著名的Weierstrass定理.定理1(Weierstrass)设,则对任何,总存在一个代数多项式,使在上一致成立.这个定理已在“数学分析”中证明过.这里需要说明的是在许多证明方法中,伯恩斯坦1912年给出的证明是一种构造性证明.他根据函数整体逼近的特性构造出伯恩斯坦多项式, (3.1.3)其中,其中,并证明了在上一致成立;若在上阶导数连续,则.这不但证明了定理1,而且由(3.1.3)给出了的一个逼近多项式.它与拉格朗日插值多项式很相似,对,当=1时也有关系式. (3.1.4)这只要在恒等式中令就可得到.但这里当时还有,于是是有界的,因而只要对任意成立,则有界,故是稳定的.至于拉格朗日多项式,由于无界,因而不能保证高阶插值的稳定性与收敛性.相比之下,多项式有良好的逼近性质,但它收敛太慢,比三次样条插值效果差得多,实际中很少被使用.更一般地,可用一组在上线性无关的函数集合来逼近,元素,表示为. (3.1.5) 函数逼近问题就是对任何,在子空间中找一个元素,使在某种意义下最小.3.1.2 范数与赋范线性空间为了对线性空间中元素大小进行衡量,需要引进范数定义,它是空间中向量长度概念的直接推广.定义2.1.2 设为线性空间,,若存在唯一实数,满足条件:(1)正定性:,(2)当且仅当时,(3);(4)齐次性:,(5);(6)三角不(7)等式:,(8).则称为线性空间上的范数,与一起称为赋范线性空间,记为.例如,在上的向量,三种常用范数为类似地对连续函数空间,若可定义三种常用范数如下:可以验证这样定义的范数均满足定义3.1.2中的三个条件.3.1.3 内积与内积空间在线性代数中,中两个向量及的内积定义为.若将它推广到一般的线性空间,则有下面的定义.定义3.1.3设是数域上的线性空间,对,有中一个数与之对应,记为,它满足以下条件:(1);(2);(3);(4),当且仅当时,.则称为上与的内积.定义了内积的线性空间称为内积空间.定义中(1)的右端称为的共轭,当为实数域时.如果=0,则称与正交,这是向量相互垂直的概念的推广.关于内积空间性质有以下重要定理.定理3.1.2设为一个内积空间,对,有(3.1.6) 称为Cauchy-Schwarz不等式.[证明]当时(3.1.6)式显然成立.现设,则,且对任何数有.取,代入上式右端,得,即得时.定理证毕定理3.1.2设为一个内积空间,,矩阵(3.1.7)称为Gram矩阵,则G非奇异的充分必要条件是线性无关.[证明]G非奇异等价于,其充分必要条件是齐次方程组(3.1.8) 只有零解.而(3.1.9) 从以上的等价关系可知,等价于从(3.1.8)推出.而后者等价于从(3.1.9)推出,即线性无关.定理证毕在内积空间上可以由内积导出一种范数,即对于,记(3.1.10) 容易验证它满足范数定义的三条性质,其中三角不等式(3.1.11)可由定理3.1.2直接得出,即两端开方即得(3.1.11).例1与的内积.设,,,则其内积定义为(3.1.12)由此导出的向量2-范数为.若给定实数,称为权系数,则在上可定义加权内积为(3.1.13)相应的范数为.不难验证(3.1.13)给出的满足内积定义的4条性质,当时,(3.1.13)就是(3.1.12).如果,带权内积定义为(3.1.14) 这里仍为正实数序列,为的共轭.在上也可以类似定义带权内积,为此先给出权函数的定义.定义3.1.4 设是有限或无限区间,在上的非负函数满足条件:(1)存在且为有限值;(2)对上的非负连续函数,如果,则.则称为上的一个权函数.例2上的内积.设,是上给定的权函数,则可定义内积. (3.1.15)容易验证它满足内积定义的4条性质,由此内积导出的范数为. (3.1.16)称(3.1.15)和(3.1.16)为带权的内积和范数.特别常用的是的情形,即若是中的线性无关函数族,记,它的Gram矩阵为(3.1.17)根据定理3.1.3可知线性无关的充分必要条件是.3.2 正交多项式正交多项式是函数逼近的重要工具,在数值积分中也有着重要的应用.3.2.1 正交函数族与正交多项式定义3.2.1 若,为上的权函数且满足, (3.2.1)则称与在上带权正交.若函数族满足关系(3.2.2)则称是上带权的正交函数族;若,则称之为标准正交函数族.例如,三角函数族就是在区间上的正交函数族.因为对有,而对,当时有定义3.2.2 设是上首项系数的次多项式,为上权函数,如果多项式序列满足关系式(3.2.2),则称多项式序列为在上带权正交,称为上带权的次正交多项式.只要给定区间及权函数,均可由一族线性无关的幂函数,利用逐个正交化手续构造出正交多项式序列;,(3.2.3) 这样得到的正交多项式序列有以下性质:(1)是具有最高次项系数为1的次多项式.(2)任何次多项式均可表示为的线性组合.(3)当时,,且与任一次数小于的多项式正交.(4)成立递推关系.其中这里.(5)设是在上带权的正交多项式序列,则的个根都是在区间内的单重实根.3.2.2 勒让德多项式当区间为[-1,1],权函数时,由正交化得到的多项式就称为勒让德(Legendre)多项式,并用表示.这是勒让德于1785年引进的,1814年罗德利克(Rodrigul)给出了简单的表达式由于是2次的多项式,求阶导数后得,于是得首项系数为,显然最高项系数为1的勒让德多项式为.(3.2.6) 勒让德多项式有下述几个性质:性质1正交性(3.2.7) [证明]令,则.设是在区间[-1,1]上的阶连续可微的函数,由分部积分知下面分两种情况讨论:(1)若是次数小于的多项式,则,故得(2)若,则,于是由于,故,于是(3.2.7)得证.性质2奇偶性(3.2.8)[证明]由于是偶次多项式,经过偶次求导仍为偶次多项式,经过奇次求导则为奇次多项式,故为偶数时为偶函数,为奇数时为奇函数,于是(3.2.8)成立.性质3递推关系(3.2.9) [证明]考虑+1次多项式,它可表示为两边乘以,并从-1到1积分,得.当时,的次数小于-1,上式左端积分为0,故得.当时.为奇函数,左端积分仍为0,故.于是.其中,代入上式整理可得(3.2.9).例1由利用性质3可得性质4在区间[-1,1]内有个不同的实零点.3.2.3 切比雪夫多项式当权函数,区间为[-1,1]时,由序列正交化得到的多项式就称为切比雪夫(Chebyshev)多项式,它可表示为(3.2.10)若令,则.切比雪夫多项式有很多重要性质:性质1递推关系(3.2.11) 这只要由三角不等式.令即得.由(3.2.11)就可推出由递推关系(3.2.11)还可得到的最高次项系数是.性质6切比雪夫多项式在区间[-1,1]上带权正交,且(3.2.12) 事实上,令,则,于是性质7只含的偶次幂,只含有的奇次幂.这性质由递推关系直接得到.性质8在区间[-1,1]上的个零点此外,实际计算中时常要求用的线性组合,其公式为. (3.2.13) 例如:结果如下:3.2.4 其他常用的正交多项式一般说,如果区间及权函数不同,则得到的正交多项式也不同.除上述两种最重要的正交多项式外,下面再给出三种较常用的正交多项式.第二类切比雪夫多项式在区间[-1,1]上带权的正交多项式称为第二类切比雪夫多项式,其表达式为. (3.2.14)令,可得即是[-1,1]上带权的正交多项式族.还可得到递推关系式.拉盖尔多项式在区间上带权的正交多项式称为拉盖尔(Laguerre)多项式,其表达式为. (3.2.15)其正交性为和递推关系.3. 埃尔米特多项式在区间上带权的正交多项式称为埃尔米特多项式.其表达式为, (3.2.16)其正交性为递推关系为.3.3 最佳一致逼近多项式3.3.1 基本概念及其理论本节讨论,在中求多项式,使其误差.这就是通常所谓最佳一致逼近或切比雪夫逼近问题.为了说明这一概念,先给出以下定义.定义3.3.1 设,,称. (3.3.1) 为与在上的偏差.显然,的全体组成一个集合,记为{},它有下界0.若记集合的下确界为(3.3.2)则称之为在上的最小偏差.定义3.3.2 假定,若存在,使得, (3.3.3)则称是在上的最佳一致逼近多项式或最小偏差逼近多项式,简称最佳逼近多项式.注意,定义并未说明最佳逼近多项式是否存在,但可证明下面的存在定理.定理4若,则总存在,使.为了研究最佳逼近多项式的特性,先引进偏差点的定义.定义3.3.3设,,若在上有,就称是的偏差点.若,称为“正”偏差点.若,称为“负”偏差点.由于函数在上连续,因此,至少存在一个点,使,也就是说的偏差点总是存在的.下面给出反映最佳逼近多项式特征的切比雪夫定理.定理3.3.2是的最佳逼近多项式的充分必要条件是在上至少有个轮流为“正”、“负”的偏差点,即有个点,使. (3.3.4) 这样的点组称为切比雪夫交错点组.[证明]只证充分性.假定在上有个点使(3.3.4)成立,要证明是在上的最佳逼近多项式.用反证法,若存在,使.由于在点上的符号与一致,故也在个点上轮流取“+”、“-”号.由连续性质,它在内有个零点,但因是不超过次的多项式,它的零点不超过.这矛盾说明假设不对,故就是所求最佳逼近多项式.充分性得证,必要性证明略,可参看[5].定理5说明用逼近的误差曲线是均匀分布的.由这定理还可得以下重要推论.推论1若,则在中存在唯一的最佳逼近多项式.证明略.利用定理5可直接得到切比雪夫多项式的一个重要性质,即定理3.3.3 在区间[-1,1]上所有最高次项系数为1的次多项式中与零的偏差最小,其偏差为.[证明]由于,且点是的切比雪夫交错点组,由定理5可知,区间[-1,1]上在中最佳逼近多项式为,即是与零的偏差最小的多项式.定理证毕例3求在[-1,1]上的最佳2次逼近多项式.解由题意,所求最佳逼近多项式应满足由定理3.3.3可知,当时,多项式与零偏差最小,故就是在[-1,1]上的最佳2次逼近多项式.3.3.2 最佳一次逼近多项式定理3.3.2给出了最佳逼近多项式的特性,但要求出却相当困难.下面讨论的情形.假定,且在内不变号,我们要求最佳一次逼近多项式.根据定理3.3.2可知至少有3个点,使由于在内不变号,故单调,在内只有一个零点,记为,于是,即.另外两个偏差点必是区间端点,即,且满足由此得到(3.3.5) 解出, (3.3.6) 代入(3.3.5)得. (3.3.7)这就得到了最佳一次逼近多项式,其几何意义如图3-3所示.直线与弦MN平行,且通过MQ的中点D,其方程为.图3-3一次最佳一致逼近多项式几何意义例4 求在上的最佳一次逼近多项式。

关于几种曲线拟合基本方法的比较

关于几种曲线拟合基本方法的比较

关于几种曲线拟合基本方法的比较学院:材料科学与工程学院专业:材料学(博)姓名:郑文静学号:1014208040 在实际工作中,变量之间的关系未必都是线性关系,更多时候,它们之间呈现出了曲线关系,在科学实验或社会活动中,通过实验或观测得到一些x和y数据,为了对位置点进行研究,很多时候,我们通过曲线拟合的方式,将这些离散点近似为一条连续的曲线,从而来预测或者得到所需结果。

曲线拟合的方法很多,本文中,主要讨论了曲线拟合的三种基础方法--插值法、磨光法、最小二乘法的特点,并对其在科学实验和生产实践中的应用性进行了比较。

插值法是函数逼近的一种基本方法,插值法就是通过函数在有限个点处的取值情况,估算出函数在其他点处的近似值。

插值法中,选取不同的插值公式,来满足实际或运算需求,得到拟合的函数。

其中,最基础的插值方法是三弯矩法,该方法是利用拉格朗日插值为基础,已知平面中的n+1个不同点,寻找一条n次多项式曲线通过这些点。

该曲线具有唯一性。

另外,还有三转角法,该方法是利用Henmiter插值为基础,其思路与三弯矩法相同,已知条件有所差别,在Henmiter插值中,不仅已知函数在一些点的函数值,而且,还知道它在这些点的导数值,甚至知道其高阶导数值,要求所求函数不仅满足过这些点,同时也要求其导函数,甚至高阶导函数满足条件。

采用Henmiter插值法求得的多项式比拉格朗日法求得的多项式有较高的光滑逼近要求。

此外,还有以分段和B-样条函数为基础的δ-基函数法,其中,样条函数是:对于[a,b]上的划分,称函数S(x)为[a,b]上关于划分△的k次样条函数,记做S k,△[a,b]。

该方法避免了高次插值可能引起的大幅度波动现象,在实际中通常采用分段低次插值来提高近似程度。

插值法常用于填充图像变换时像素之间的空隙。

磨光法是适应保凸性要求的数据拟合方法。

积分可以改变函数的光滑度,而微商是积分的逆运算,对函数进行积分,然后在微商,可以将函数还原。

第三章函数逼近和曲线拟合

第三章函数逼近和曲线拟合
则称 x1, x2 ,..., xn 为空间S的一组基,记为:
S=span{ x1,..., xn}
并称该空间为n维空间。1,2 ,...,n P
称为x在这组基下的坐标。 例:n次多项式
p(x) Hn , p(x)=a0 + a1x ... an xn Hn span{1, x, x2 ,..., xn}
4
11
4.5
12
4.6
强 度 yi 编 号 拉伸倍数 xi
1.4
13
5
1.3
14
5.2
1.8
15
6
2.5
16
6.3
2.8
17
6.5
2.5
18
7.1
3
19
8
2.7
20
8
4
21
8.9
3.5
22
9
4.2
23
9.5
3.5
24
10
强 度 yi
5.5 5
5.5
6.4 6
5.3 6.5
7 8.5
8 8.1 8.1
6
内积与内积空间 定义3:设X为数域K(R或C)上的线性空
间,满足条件:
u, v X , k (u, v) K, st.
(1) (u, v) (v, u)
(2) (u, v) (u, v), for K
(3) (u v, w) (u, w) (v, w), for w X
(4) (u, u) 0, u 0 iff (u, u) 0
存在唯一实数 g ,满足条件:
(1) x 0; x 0 iff x 0
(2) x x , R
(3) x y x y , x, y R

第3章 函数逼近与曲线拟合

第3章  函数逼近与曲线拟合

例如、 三角函数系:1,cosx,sinx,cos2x,sin2x,…是 区间[-π,π]上的正交函数系,因为



s in k xs in jxdx 0,
jxdx 0,

( j k)
( j k)
cos k xcos


s in k xcos


a
(1)
则称f(x)与g(x)在[a,b]上带权正交. 若函数族
0 ( x), 1 ( x), , n ( x),
b
满足关系
j k; 0, ( j , k ) ( x) j ( x)k ( x)dx a Ak 0, j k ;
(2)
则称 k ( x) 是[a, b]上带权 ρ(x)正交函数族 ;
如果范数 取为 2,即 f ( x) P ( x)
2 2
min f ( x) P( x)
PH n PH n
2 2
b min a [ f ( x) P( x)]2 dx ,
则称P ( x)为f ( x)在[a, b]上的最佳平方逼近多项式 。 若f ( x)是 [a, b]上的一个列表函数,在a x0 x1 xm b 上给出f ( xi )(i 0,1,, m), 要求P 使 f P 2 min f P 2 min
PH n
则称P ( x)是f ( x)在[a, b]上的最佳逼近多项式 。 若取 ,即 f ( x) P ( x) min f ( x) P( x)
PH n
min max f ( x) P( x) ,
PH n a x b

Matlab中的曲线拟合与函数逼近

Matlab中的曲线拟合与函数逼近

Matlab中的曲线拟合与函数逼近Matlab是一种功能强大的数学工具,广泛应用于科学计算、工程分析和数据处理等领域。

在Matlab中,曲线拟合和函数逼近是常见的数学问题,它们可以帮助我们从一组离散数据中找到合适的函数形式,从而更好地理解数据的规律和趋势。

本文将介绍Matlab中的曲线拟合和函数逼近的常见方法和技巧,并通过实例来说明其应用。

一、简单线性回归拟合简单线性回归是最基本的曲线拟合方法之一,在Matlab中使用polyfit函数可以实现。

假设我们有一组离散的数据点,分别表示自变量x和因变量y的取值,我们可以通过简单线性回归来寻找一条直线,使得该直线与这些数据点的拟合误差最小。

```matlabx = [1, 2, 3, 4, 5];y = [2.1, 3.9, 6.2, 8.2, 9.8];coefficients = polyfit(x, y, 1);```在上述代码中,x和y分别是自变量和因变量的数据点。

polyfit函数的第三个参数表示我们希望拟合的曲线的阶数,这里是1表示直线拟合。

函数返回的coefficients是拟合曲线的系数,其中第一个元素表示直线的斜率,第二个元素表示直线的截距。

我们可以利用polyval函数来计算拟合直线上的点的函数值,从而与原始数据进行比较。

```matlaby_fit = polyval(coefficients, x);```通过绘制拟合直线和原始数据,我们可以直观地看到拟合效果。

```matlabplot(x, y, 'o', x, y_fit, '-')legend('原始数据', '拟合直线')```二、多项式拟合除了简单线性回归,Matlab还提供了多项式拟合的方法,可以通过增加拟合曲线的阶数来逼近更复杂的数据。

```matlabx = [1, 2, 3, 4, 5];y = [2.1, 3.9, 6.2, 8.2, 9.8];coefficients = polyfit(x, y, 2);```在上述代码中,polyfit函数的第三个参数是2,表示进行二次多项式拟合。

曲线拟合 经验和专业知识 概述及解释说明

曲线拟合 经验和专业知识 概述及解释说明

曲线拟合经验和专业知识概述及解释说明1. 引言1.1 概述在科学研究和工程实践中,曲线拟合是一种常见的分析方法,它用于描述和预测数据集中的趋势。

曲线拟合通过选择适当的数学模型,并使用统计技术对模型参数进行估计,从而找到最佳拟合曲线。

1.2 文章结构本文将对曲线拟合涉及的经验和专业知识进行综述与解释。

首先,在第二部分我们将介绍曲线拟合的基本概念和定义,以及常见的曲线拟合方法。

接着,在第三部分我们将探讨曲线拟合在不同领域中的应用,并提供实例分析。

然后,在第四部分我们将介绍与曲线拟合相关的算法和数值计算技术,并讨论数值稳定性与误差分析方面的考虑。

最后,在第五部分我们将总结文章主要观点和研究成果,同时展望未来发展趋势和可能的研究方向。

1.3 目的本文旨在帮助读者了解曲线拟合所需的经验和专业知识,并为他们在实际问题中正确应用此方法提供指导。

我们希望通过介绍曲线拟合的基本概念、常见方法和实例分析,读者们能够深入理解曲线拟合在不同领域中的应用,并能够正确选择适当的数学模型和参数估计方法。

此外,我们还将讨论与曲线拟合相关的算法和数值计算技术,以及数值稳定性和误差分析方面的问题,帮助读者更好地理解这些技术并掌握其应用。

以上是“1. 引言”部分内容的详细清晰撰写,请参考。

2. 曲线拟合的经验和专业知识2.1 定义和基本概念曲线拟合是一种数学方法,它通过使用已知数据点来构建一个与这些数据最匹配的函数曲线。

在进行曲线拟合时,我们通常选择一个特定的函数形式(例如多项式、指数、对数等)来代表所要拟合的关系。

基本概念包括目标函数、误差函数和参数估计。

目标函数是需要找到的逼近实际数据的理论模型。

这个函数可以是多种形式,我们根据具体问题选择适当的函数类型。

误差函数是用来度量实际数据点与拟合曲线之间的偏离程度。

参数估计则是通过最小化误差函数来确定在所选模型中使用的参数值。

2.2 常见的曲线拟合方法在进行曲线拟合时,有几种常见的方法可供选择:- 最小二乘法:这是最常用且简单直观的方法。

第三章函数逼近1

第三章函数逼近1

6.7941 -5.3475
-653..32457859
5.1084 -49.0086
63.2589 10-0429..50086

-2126...763718262837
Go!
用Gauss列主元消去法,得

a b c


0--.0113..020467113035
y(x) a0 a1x
为拟合函数,其基函数为
0(x) 1 1(x) x
建立法方程组 根据内积公式,可得
(0 ,0 ) 24 (0 , f ) 113.1
法方程组为
(0 ,1 ) 127.5 (1 ,1 ) 829.61 (1 , f ) 731.6
m

*
2 2

(S * ( xi ) yi )2
i0
m
min S ( x)
2 2

min
S ( x)
i0
(S(xi )

yi
)2
m
其中S(x) a j j (x)为中的任意函数 j0
---------(3)
n
称满足条件(3)的求函数S * (x) a*j j (x)的方法为 j0
mn
m
a j j (xi )k (xi ) yik (xi )
i0 j0
i0
mn
m
a j j (xi )k (xi ) yik (xi )
i0 j0
i0
nm
m
[ j (xi )k (xi )]a j yik (xi )
mn
是 (a0 , a1 , , an ) ( a j j ( xi ) yi )2 的最小值 i0 j0

数值计算(分析)实

数值计算(分析)实

数值计算(分析)实验报告2南昌航空大学数学与信息科学学院实验报告课程名称:《数值计算方法》实验名称:曲线拟合实验类型:验证性■综合性□设计性□实验室名称:数学实验室班级学号: 09072113学生姓名:邢宪平任课教师(教师签名):成绩:一、实验目的实验目的:实验目的:了解函数逼近与曲线拟合的基本原理,并且运用MATLAB 软件进行实践操作。

二、实验原理、程序框图、程序代码等 实验题目:题目1:试分别用抛物线2y a bx cx =++和指数曲线bxy ae =拟合下列数据并比较两个拟合函数的优劣。

题目2:已知实验数据如下:试用形如2y a bx =+的抛物线进行最小二乘拟合。

实验原理:1、逼近方式 假设()[,]f x C a b ∈,2{1,,,...,}n nHspan x x x =,()nnP x H ∈,称(,)|||||()()|max n n n a x bf P F P f x P x ≤≤=-=-V 为()f x 与()|nP x 在[,]a b 上的偏差。

若存在*()nnP x H ∈,使得**(,)|||||()()|max inf n nn nn P H a x bf P f Pf x P x ∞∈≤≤=-=-V 则称*()nP x 是()f x 在[,]a b 上的最佳一致逼近多项式。

假设()[,]f x C a b ∈及[,]C a b 的一个子集01{(),(),,...()}nspan x x x ϕ=ϕϕϕ,若存在*()S x ϕ∈,使*22222()()||()()||||()()||()[()()]min min bS x S x af x S x f x S x x f x S x dxϕϕρ∈∈-=-=-⎰则称*()S x 是()f x 在子集[,]C a b ϕ⊂中的最佳平方逼近数。

2、曲线拟合上述函数的最佳平方逼近法中,若()f x 是以一组离散点集的形式给出的,即给出了函数()f x 在一些离散点上的值{(,),0,1,...,}iix y i m =,则该方法就是所说的曲线拟合。

(完整版)高等数值分析48课时教案

(完整版)高等数值分析48课时教案

高等数值分析48课时教案
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
南华大学教案
2010 ~ 2011 学年第 1 学期
课程:高等数值分析授课教师(职称):王礼广(副教授)班级: 2010级理工科研究生
11。

数值分析第三章函数逼近与 曲线拟合习题答案

数值分析第三章函数逼近与    曲线拟合习题答案

6。对,定义 问它们是否构成内积。 解: 令(C为常数,且) 则 而 这与当且仅当时,矛盾 不能构成上的内积。 若,则 ,则 若,则 ,且 即当且仅当时,. 故可以构成上的内积。 7。令,试证是在上带权的正交多项式,并求。 解: 若,则 令,则,且,故 又切比雪夫多项式在区间上带权正交,且 是在上带权的正交多项式。 又 8。对权函数,区间,试求首项系数为1的正交多项式 解: 若,则区间上内积为 定义,则 其中 9。试证明由教材式给出的第二类切比雪夫多项式族是上带权的正交多 项式。 证明: 若 令,可得 当时, 当时, 又,故 得证。 10。证明切比雪夫多项式满足微分方程 证明:
若 且,则 则法方程组为 解得 故关于的最佳平方逼近多项式为 17。求函数在指定区间上对于的最佳逼近多项式: 解: 若 且,则有 则法方程组为 从而解得 故关于的最佳平方逼近多项式为 若 且,则有 则法方程组为 从而解得 故关于的最佳平方逼近多项式为 若 且,则有 则法方程组为 从而解得 故关于的最佳平方逼近多项式为 若 且则有 则法方程组为 从而解得 故关于最佳平方逼近多项式为 18。,在上按勒让德多项式展开求三次最佳平方逼近多项式。 解: 按勒让德多项式展开 则 从而的三次最佳平方逼近多项式为 19。观测物体的直线运动,得出以下数据:
切比雪夫多项式为 从而有 得证。 11。假设在上连续,求的零次最佳一致逼近多项式? 解: 在闭区间上连续 存在,使 取 则和是上的2个轮流为“正”、“负”的偏差点。 由切比雪夫定理知 P为的零次最佳一致逼近多项式。 12。选取常数,使达到极小,又问这个解是否唯一? 解: 令 则在上为奇函数 又的最高次项系数为1,且为3次多项式。 与0的偏差最小。 从而有 13。求在上的最佳一次逼近多项式,并估计误差。 解: 于是得的最佳一次逼近多项式为 即 误差限为 14。求在上的最佳一次逼近多项式。 解: 于是得的最佳一次逼近多项式为 15。求在区间上的三次最佳一致逼近多项式。 解: 令,则 且 令,则 若为区间上的最佳三次逼近多项式应满足 当 时,多项式与零偏差最小,故 进而,的三次最佳一致逼近多项式为,则的三次最佳一致逼近多项式为 16。,在上求关于的最佳平方逼近多项式。 解:

《数值分析》课程实验报告范文

《数值分析》课程实验报告范文

《数值分析》课程实验报告范文《数值分析》课程实验报告姓名:学号:学院:机电学院日期:2022年某月某日目录实验一函数插值方法1实验二函数逼近与曲线拟合5实验三数值积分与数值微分7实验四线方程组的直接解法9实验五解线性方程组的迭代法15实验六非线性方程求根19实验七矩阵特征值问题计算21实验八常微分方程初值问题数值解法24实验一函数插值方法一、问题提出对于给定的一元函数的n+1个节点值。

试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。

实验二函数逼近与曲线拟合一、问题提出从随机的数据中找出其规律性,给出其近似表达式的问题,在生产实践和科学实验中大量存在,通常利用数据的最小二乘法求得拟合曲线。

在某冶炼过程中,根据统计数据的含碳量与时间关系,试求含碳量与时间t的拟合曲线。

t(分)051015202530354045505501.272.162.863.443.874.154.374.51 4.584.024.64二、要求1、用最小二乘法进行曲线拟合;2、近似解析表达式为;3、打印出拟合函数,并打印出与的误差,;4、另外选取一个近似表达式,尝试拟合效果的比较;5、某绘制出曲线拟合图。

三、目的和意义1、掌握曲线拟合的最小二乘法;2、最小二乘法亦可用于解超定线代数方程组;3、探索拟合函数的选择与拟合精度间的关系四、实验步骤:第一步先写出线性最小二乘法的M文件functionc=lpoly(某,y,m)n=length(某);b=zero(1:m+1);f=zero(n,m+1); fork=1:m+1f(:,k)=某.^(k-1);enda=f'某f;b=f'某y';c=a\b;c=flipud(c);第二步在命令窗口输入:>>lpoly([0,5,10,15,20,25,30,35,40,45,50,55],[0,1.27,2.16,2.86,3.44,3.87,4.15,4.37,4.51,4.58,4.02,4.64],2)回车得到:an=-0.00240.20370.2305即所求的拟合曲线为y=-0.0024某2+0.2037某+0.2305在编辑窗口输入如下命令:>>某=[0,5,10,15,20,25,30,35,40,45,50,55];>>y=-0.0024某某.^2+0.2037某某+0.2305;>>plot(某,y)命令执行得到如下图五、实验结论分析复杂实验数据时,常采用分段曲线拟合方法。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。


sin kx cos kx 0,




si n2 kx
Байду номын сангаас


cos2 kx 0,
实际上,这就是付里叶(Fourier)逼近的基函数.
2)如何构造正交多项式 • 只要给定区间[a, b]及权函数,均可由一组线性无 关的幂函数{1,x,…, xn,…},利用逐个正交化手续 构造出正交多项式序列 {n ( x)}0 :
段低次多项式等.
数学上常把在各种集合中引入某一些不同的确 定关系称为集合以某种空间结构赋予,并将这 样的集合称为空间。
例1、按向量的加法和数乘构成实数域 R 上的线性空间--- R n 例2、对次数不超过 n 的实系数多项式,按 加法和数乘构成数域上的多项式线性 空间---- H n
例3、所有定义在 [a,b] 集合上的连续函数
3)几种常用范数
在Rn上的向量 x=(x1,x2,…,xn)T∈Rn, 三种常用 范数为称为:
称为 -范数或最大范数
n ||x|| | x |, 1 i i 1
称为 1-范数
1 2
||x||2= xi2 i 1
n
称为 2-范数
类似的对连续函数空间C[a, b], 若f∈C[a,b] 可定义以下三种常用函数的范数
• 函数的内积
记区间[a,b]上所有连续函数的全体为C[a,b],
可以证明C[a,b]是一个线性空间,把所有次数不超
过n的多项式全体记为Pn,则 Pn是C[a,b]的子空间. 若(x),g(x)C[a,b],
则称
满足

b
a
f ( x) g ( x)dx 为(x)与g(x)的内积,记为(,g),
(2)奇偶性
(3)递推关系
pn ( x) (1) n pn ( x)
(n 1) pn1 ( x) (2n 1) xpn ( x) npn1 ( x) n (1,2)
且有以下常用公式
p 0 ( x) 1 p1 ( x ) x p 2 ( x ) (3 x 2 1) / 2 p 3 ( x ) (5 x 3 3 x ) / 2 p 4 ( x ) (35x 4 30x 2 x ) / 8 p 5 ( x ) (63x 5 70x 3 15x ) / 8 p 6 ( x ) ( 231 6 315x 4 105x 2 5) / 16 x
p0 ( x) 1, 1 dn 2 n Pn ( x) n {( x 1) } n 2 n! dx (n 1,2,)
• 最高项系数为1的勒让德多项式为
n! dn ~ P ( x) [( x 2 1) n ]. n ( 2n)! dx n
• 勒让德多项式的性质
(1)正交性
m n; 0, 1 1 Pn ( x) Pm ( x)dx 2 , m=n. 2n 1
( x, 0 )
1
0
1 x ln xdx 4
( x , 0 )
2
1
0
于是
1
1 1 ( x) x 4
1 x ln xdx 9
2
1 1 2 1 1 7 2 (1 , 1 ) ( ln x)( x ) dx (ln x)( x x )dx 0 0 4 2 16 144
法之一,插值的近似标准是在插值点处误差为
零. 但在实际应用中,有时不要求具体某些点
误差为零,而要求考虑整体的误差限制 ,这就
引出了拟合和逼近的概念.
什么是函数逼近
对函数类A中给定的函数 f(x),记作f(x)∈A, 要求在另一类简单的便于计算的函数类 B 中求函数 p(x)∈B ,使 p(x)与 f(x)的误差在 某种意义下最小.函数类A通常是区间[a, b] 上的连续函数,记作C[a, b],称为函数逼近空 间;而函数B通常为n次多项式,有理函数或分
a
b
f
2


b
a
( x) f 2 ( x)dx
这里函数(x)是非负连续函数,称为[a,b]上的权函数.
它的物理意义可以解释为密度函数.
3.2 什么是正交多项式
1) 正交的定义 若f(x),g(x)∈C[a,b],ρ(x)为[a,b]上的权函数且满 足 b
( f ( x), g ( x)) ( x) f ( x) g ( x)dx 0,
n1 ( x) ( x an ) n ( x) n n1 ( x)
(n 0,1,...)
-1 (x) 0, n ( x n ( x), n ( x)) /( n ( x), n ( x)), n ( n ( x), n ( x)) /( n1 ( x), ( n1 ( x))
1 5 ( x , 1 ) ( ln x) x ( x )dx 0 4 144
2 1 2
于是
1 5 1 5 17 2 2 ( x) x ( x ) x x 9 7 4 7 252
2
3)几种常用的正交多项式
• 勒让德多项式 当区间[-1,1],权函数ρ(x) ≡1时,由 {1,x,…,xn,…}正交化得到的多项式就称为 勒让德多项式,并用P0(x),P1(x),…,Pn(x),… 表示. 其简单的表达式为
b a 2 n
其中 0 ( x) 1,
(n 1,2,.....)
并且( x n ( x), n ( x))= x ( x) ( x)dx.
{ n ( x)} 是在[a, b]上带权ρ(x)的正交多项式 (5)设 0 序列, 则 n (x) (n≥1)的n个根都是在区间(a,
(4)
pn ( x) 在区间[-1,1]内有n个不同的实零点。
切比雪夫多项式
当权函数
( x)
1 1 x
2
,
区间为[-1,1]
时,由序列{1,x,…,xn,…}正交化得到的多项 式就是切比雪夫多项式,它可表示为 Tn(x)=cos(narccosx), | x |≤1 若令x=cosθ,则Tn(x)=cosnθ, 0≤θ≤π.
全体,按函数的加法和数乘构成连续函数 空间---- C[a, b]
3.1 函数逼近的基本概念
1)线性无关
设集合S是数域P上的线性空间,元素
x1,x2,…,xn∈S,如果存在不全为零的数
a1,a2,…,an∈P,使得
a1 x1 a2 x2 ... an xn 0,
则称x1,x2,…,xn线性无关.
2)范数的定义 设S为线性空间,x∈S,若存在唯一实数 || || 满足条件: (1)‖x‖≥0;当且仅当x=0时,‖x‖=0; (正定性) (2)‖αx‖=|α|‖x‖,α∈R; (齐次性) (3)‖x+y‖≤‖x‖+‖y‖,x,y∈S. (三角不等式) 则称 || || 为线性空间S上的范数, S与 || || 一起称为赋范线性空间,记为X.
b)内的单重实根.
例题:利用 Gram-schmidt 方法构造 [0,1] 上带权 1 ( x ) ln 的前3个正交多项式 0 , 1 , 2 x
0 ( x) 1,
n ( x) x
n
n 1
j ( x) j 0 ( j ( x), j ( x))
( x , j ( x))
n
(n 1,2,...).
例题:利用 Gram-schmidt 方法构造 [0,1] 上带权 1 ( x ) ln 的前3个正交多项式 0 , 1 , 2 x
解:利用正交化公式来求
( x ) ln 1 ln x x
0 ( x) 1
(u1 , u1 ) (u2 , u1 ) (u1 , u2 ) (u2 , u2 ) G (u , u ) (u , u ) 2 n 1 n (un , u1 ) (un , u2 ) (un , un )
称为格拉姆矩阵,则G非奇异的充分必要条件是 u1,u2,…,un线性无关 。
0 ( x) 1,
n ( x) x
n
n 1
j ( x) j 0 ( j ( x), j ( x))
( x , j ( x))
n
(n 1,2,...).
如此得到的正交多项式有如下性质: (1) n (x) 是具有最高次项系数为1的n次
(2)任何n次多项式Pn(x)∈Hn均可表示为 0 ( x), 1 ( x),..., n ( x) 的线性组合. ( j ( x), k ( x)) 0, 且 k ( x) (3)当k≠j时, 与任一次数小于k的多项式正交. (4)成立递推关系
a
(1)
则称f(x)与g(x)在[a,b]上带权正交. 若函数族
0 ( x), 1 ( x), , n ( x),
b
满足关系
j k; 0, ( j , k ) ( x) j ( x)k ( x)dx a Ak 0, j k ;
(2)
则称 k ( x) 是[a, b]上带权 ρ(x)正交函数族 ;
|| f || max | f ( x) |,
a x b
称为 范数
|| f ||1 | f ( x) | dx,
a
b
称为 范数 1
||f ||2 ( f 2 ( x) dx) ,
a
b
1 2
称为2 范数
柯西-施瓦次不等式
设X为一个内积空间,对
2
u,v∈X有
(2) c2=|c|2;
(3) +g22+g2 (4) (,g)2 g2
相关文档
最新文档