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

合集下载

小波分析在函数优化中的应用

小波分析在函数优化中的应用

小波分析在函数优化中的应用摘 要: 本文介绍了小波变换以及它的多尺度性,结合局部优化来解决函数的全局优化问题,提出了一种新的全局优化方法。

首先给出算法,再用matlab 编程并给出若干个多峰值函数求解全局最优值的算例,得到了与理论值十分接近的优化结果。

结果表明该算法具有一定的可行性。

关键词: 小波分析;全局优化;多尺度分析;多峰函数 中图分类号: 文献标识码:Application of Wavelet Analysis in Function OptimizationAbstract :This paper introduces wavelet transform and its property of multi-scale . It combines with local optimization to solve the global optimization, and puts forward a new method of global optimization. First giving the algorithm, giving a number of multimodal functions to solve the global optimal values with matlab programming and obtain the optimal result which is very close to the theoretical value. It shows that the algorithm has certain feasibility.Key words :Global Optimization; Multimodal Function; Wavelet Analysis; Multiscale Analysis0引言“小波”(Wavelet )是当前应用数学中一个迅速发展起来的新领域,经过近十几年的探索研究,已逐步建立并形成了比较完善的数学形式化体系。

(整理)数值分析课件 第3章 函数逼近与曲线拟合

(整理)数值分析课件 第3章 函数逼近与曲线拟合

第三章 函数逼近与曲线拟合1 函数的逼近与基本概念1.1问题的提出多数计算机的硬件系统只提供加、减、乘、除四种算术运算指令,因此为了计算大多数有解析表达式的函数的值,必须产生可用四则运算进行计算的近似式,一般为多项式和有理分式函数.实际上,我们已经接触到两种逼近多项式,一种是泰乐多项式,一种是插值多项式.泰乐多项式是一种局部方法,误差分布不均匀,满足一定精度要求的泰乐多项式次数太高,不宜在计算机上直接使用.例如,设()f x 是[1,1]-上的光滑函数,它的Taylor 级数0()k k k f x a x ∞==∑,()(0)!k k f a k =在[1,1]-上收敛。

当此级数收敛比较快时,11()()()n n n n e x f x s x a x ++=-≈。

这个误差分布是不均匀的。

当0x =时,(0)0n e =,而x 离开零点增加时,()n e x 单调增加,在1x =±误差最大。

为了使[1,1]-的所有x 满足()()n f x s x ε-<,必须选取足够大的n ,这显然是不经济的。

插值函数出现的龙格现象表明,非节点处函数和它的插值多项式相差太大。

更重要的是,实际中通过观测得到的节点数据往往有各种误差,此时如果要求逼近函数过全部节点,相当于保留全部数据误差,这是不适宜的。

如图1所示,给出五个点上的实验测量数据,理论上的结果应该满足线性关系,即图1中的实线。

由于实验数据的误差太大,不能用过任意两点的直线逼近函数。

如果用过5个点的4次多项式逼近线性函数,显然误差会很大。

实验数据真函数插值多项式逼近精确的线性逼近图11.2范数与逼近一、线性空间及赋范线性空间要深入研究客观事物,不得不研究事物间的内在联系,给集合的元素之间赋予某种“确定关系”也正是这样的道理.数学上常把在各种集合中引入某些不同的确定关系称为赋予集合以某种空间结构,并将这样的集合称为空间.最常用的给集合赋予一种“加法”和“数乘”运算,使其构成线性空间.例如将所有实n 维数对组成的集合,按照“加法”和“数乘”运算构成实数域上的线性空间,记作n R ,称为n 维向量空间.类似地,对次数不超过n 的实系数多项式全体,按通常多项式与多项式加法及数与多项式乘法也构成数域R 上一个线性空间,用n H 表示,称为多项式空间.所有定义在[,]a b 上的连续函数集合,按函数加法和数与函数乘法构成数域R 上的线性空间,记作[,]C a b .类似地,记[,]p C a b 为具有p 阶连续导数的函数空间.在实数的计算问题中,对实数的大小、距离及误差界等是通过绝对值来度量的.实践中,我们常常会遇到对一般线性空间中的向量大小和向量之间的距离进行度量的问题,因此有必要在一般线性空间上,赋予“长度”结构,使线性空间成为赋范线性空间.定义1 设X 是数域K 上一个线性空间,在其上定义一个实值函数,即对于任意,x y X ∈及K α∈,有对应的实数x 和y ,满足下列条件(1) 正定性:0x ≥,而且0x =当且仅当0x =;(2) 齐次性:x x αα=;(3) 三角不等式:x y x y +≤+;称为X 上的范数,定义了范数的线性空间就称为赋范线性空间.以上三个条件刻划了“长度”、“大小”及“距离”的本质,因此称为范数公理.对n X 上的任一种范数,n X ∀∈x,y ,显然有±≥-x y x y .n R 上常用的几种范数有:(1) 向量的∞-范数:1max i i nx ∞≤≤=x(2) 向量的1-范数:11n i i x ==∑x(3) 向量的2-范数:12221()n i i x ==∑x (4) 向量的p -范数:11()n p pi p i x ==∑x其中[1,)p ∈∞,可以证明向量函数()p N x x ≡是nR 上向量的范数. 前三种范数是p -范数的特殊情况(lim p p ∞→∞=x x ).我们只需表明(1).事实上1111111max max max n n p pp p i i i i i n i n i n i i x x x x ≤≤≤≤≤≤==⎛⎫⎛⎫≤≤≤ ⎪ ⎪⎝⎭⎝⎭∑∑及max 1p →∞=,故由数学分析的夹逼定理有1l i m ma x i p p i nx ∞→∞≤≤==x x 。

3函数逼近和曲线拟合

3函数逼近和曲线拟合
n
其他性质:
Pk ( x)

n
0;
k0
Pk ( x)

n

n xk (1
k0 k
x)nk
1.
若 f ( x) ,x [0,1],则
Bn (
f
, x)

n

k0
f

k n

Pk
(
x)

max
0 x1
f
n
( x) Pk ( x)
k0
( x) Pn
( 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
(5)设{ pn( x)}0 为以( x)为权函数的[a,b]上的正交多项式
序列. 则pn( x)(n 1)的n个根都是在(a,b)内的单重实根;并有
pn1( x) (n x n ) pn( x) n pn1( x), n 1,2,, (2.4)
其中
p1( x)
)
(u2, u2 )

(un , u2 )
(u1, un ) (u2, un ) (un , un )
称为Gram矩阵,则G非奇异的充要条件是u1, u2,, un线性

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

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

(ϕ 0 , ϕ 0 ) (ϕ 0 , ϕ1 ) ⋯ (ϕ 0 , ϕ n ) (ϕ1 , ϕ 0 ) (ϕ1 , ϕ1 ) ⋯ (ϕ1 , ϕ n ) G = ⋮ ⋮ ⋮ (ϕ , ϕ ) (ϕ , ϕ ) ⋯ (ϕ , ϕ ) n 1 n n n 0
∑α
j =1
n
j
u j = 0.
Th3的证明
G非奇异当且仅当齐次方程组
∑ (u
j =1
n
j
, u k ) α j = 0, k = 1, 2, ⋯ , n.
只有零解,即 ( ∑ α j u j , u k ) = 0, k = 1, 2, ⋯ , n.
j =1
n
只有零解,即
∑α
j =1
n
j
u j = 0 只有零解,即
C[a.b]带权内积的定义
• [a,b]上的非负函数 ρ (x) 称为[a,b]上的权 函数,若满足: b ① ∫ x k ρ ( x)dx 存在且为有限值(k=0,1,…). a ② 对[a,b]上的非负连续函数g (x) ,若

b
a
ρ ( x ) g ( x ) dx = 0 ,
则 g ( x) ≡ 0
2
= =

n
x
2 i
i = 1

m a x
1 ≤ i ≤ n
x
i
n与C[a,b]上范数的扩充关系 R
= = =
• 向量范数: 范数: 范数
x x x
1

n
x
i = 1
i
2

n∞
max
x
1 ≤ i ≤ n

数值分析 第3章 函数逼近与曲线拟合)

数值分析 第3章 函数逼近与曲线拟合)

在[a, b]上一致成立 。
定理:设X为一个内积空间,u1,u2,…,un∈X,矩阵
(u1, u1) (u2 , u1)
G
(u1, u2
(u1, un
) )
(u2 , u2 )
(u2 , un )
(un , u1)
(un , u2 )
(un
, un
)
称为格拉姆矩阵,则G非奇异的充分必要条件是 u1,u2,…,un线性无关 。
n1(x) (x an )n (x) n n1(x)
(n 0,1,...)
其中 0 (x) 1, -1(x) 0, n (xn (x),n (x)) /(n (x),n (x)), n (n (x),n (x)) /(n1(x),(n1(x))
(n 1,2,.....)
并且(
中找一个元素 * (x) 使 f (x) *(x) 在某种意义下
最小.
3、 范数的定义
设S为线性空间,x∈S,若存在唯一实数 || || 满足条件:
(1)‖x‖≥0;当且仅当x=0时,‖x‖=0; (正定性)
(2)‖αx‖=|α|‖x‖,α∈R; (齐次性)
(3)‖x+y‖≤‖x‖+‖y‖,x,y∈S. (三角不等式)
类较简单的便于计算的函数类B中,求函数 P(x) B , 使P(x)与f(x)
之差在某种度量意义下最小” . 函数类A通常是区间[a,b]上的连续 函数,记作C[a,b];函数类B通常是代数多项式,分式有理函数或 三角多项式.
2、函数空间 数学上常把在各种集合中引入某些不同的确定关系称为赋予
集合以某种空间结构,并将这样的集合称为空间.
1 2n n!
dn dxn
{(

基于小波模糊网络的非线性函数逼近方法的研究

基于小波模糊网络的非线性函数逼近方法的研究

・322・ 计算机测量与控制.2006.14(3) Computer Measurement &Control 工业控制中华测控网收稿日期:2005-07-16; 修回日期:2005-08-27。

作者简介:李春鑫(1980-),男,辽宁省大连市人,硕士研究生,主要从事军事航海信息及控制关键技术的研究。

文章编号:1671-4598(2006)03-0322-02 中图分类号:TP183 文献标识码:A基于小波模糊网络的非线性函数逼近方法的研究李春鑫1,李天伟2,王孝通2(1.海军大连舰艇学院研究生2队,辽宁大连 116018; 2.海军大连舰艇学院航海系,辽宁大连 116018)摘要:针对非线性函数逼近问题,提出了一种新的融合策略———小波模糊网络;该网络将模糊模型引入小波网络,采用正交最小二乘法筛选小波,利用推广卡尔曼滤波算法调整网络非线性参数,避免陷入局部最优,提高学习速度,并采用最小二乘法修正权值,在不增加小波基函数的基础上提高网络的逼近精度;通过仿真,该网络的准确性和泛化能力都优于传统的小波神经网络,具有广泛的应用前景。

关键词:小波模糊网络;小波神经网络;推广卡尔曼滤波;函数逼近R esearch on Approximating Non -linear FunctionsB ased on W avelet Fuzzy N et w orkLi Chunxin 1,Li Tianwei 2,Wang Xiaotong 2(1.Postgraduate Team 2,Dalian Naval Academy ,Dalian 116018,China ;2.Dept.of Navigation ,Dalian Naval Academy ,Dalian 116018,China )Abstract :Aiming at t he problem of approximating non -linear functions ,a new syncretized strategy -wavelet fuzzy network is presen 2ted.The network brings fuzzy model in wavelet network ,uses Ort honormal Least Square to select wavelet ,and to avoid into part optimiza 2tion and improve learning speed ,Extended Kalman Filter is adopted to correct weight s ,t herefore ,it can improve approach precision wit hout increasing wavelet bases functions.The simulation result illustrates it is superior to traditional wavelet neural network in veracity and general 2ization ability ,and it s application foreground is broad.K ey w ords :wavelet fuzzy network ;wavelet neural network ;Extended Kalman Filter ;function approach0 引言小波神经网络是用小波基作为神经元函数而生成的一种网络。

曲线拟合

曲线拟合

根据待检信号分 析结果,与待检 信号拟合系数距 离最小的标准故 障模式,即为待 检信号所反映的 故障。

3、在故障诊断中的应用
现以柴油机的故障诊断为例,通过曲线拟合法进行故障诊断。 1)建立标准故障模式 对于每一种故障状态,提取所采集的多组信号的多个特征参数,求每组特征 参数的平均值,然后分别将不同的特征参数的平均值作为拟合曲线的纵坐标,即:
有理论模型的曲线拟合
无理论模型的曲线拟合
2、基本原理及实现方法
2.2 曲线拟合的方法——最小二乘法
已知试验数据点(xi.yi)(i=1,2,…,n),假设实验数据点可以用线性模型 拟合,解析式为: 其中,β0,β1是待求参数,误差ε服从N(1,σ2) 将n个实验点分别带入表达式(1)得到:
y=β0+β1x+ε
y=[u1,u2,…,un]
同时取自然数横坐标,即:
x=[1,2,…,n]
然后,运用最小二乘法进行多项式曲线拟合,求出拟合系数,这样便可以得到 不同故障状态下的多项式拟合系数模式,即:
Mk=[ank,an-1k,…,a2k,a1k]
式中,n 为拟合多项式的拟合阶数。
3、在故障诊断中的应用
对柴油机取六种工作状态,每种工作状态取五个时域特征参数,然后根据 大量的实验数据求出特征参数平均值,如下表所示:
那么它们在相应的权重下生成的网络函数就可看作 是泛化了的曲线模型。在曲线拟合中最常用的激活 函数是Sigmoid 函数
3、在故障诊断中的应用
建立标准 故障模式
分析待检 信号
故障判断
对每种故障状态 ,根据试验测试 数据进行曲线拟 合,得到相应的 拟合系数,形成 对应的标准故障 模式。
对待检信号的参 数进行拟合,并 求出待检信号与 各种标准故障模 式拟合系数之间 的距离。

DNA测序信号去噪分析的一种新方法

DNA测序信号去噪分析的一种新方法
方法。
可得 到一系列频率通带空 间,ef 波空 间。目的是获得信号 P, l , £的离散小波分解在不同尺度下 的带通 项 。Maa 提 出了 () lt l

种快速分解算法L , 4 利用小波变换将信号 分解成不 同的频 ]
, £ c () c () (): A1 £+ D1£ () 2
的部分 ,而
c () D1£ 为频率介 于 2 和 2 之间的部分 。 。
() £ 还可以继续
由于小 波分析 中用到的小波函数 具有多样性 , 用不 同的 小波基 函数分析 同一个问题会得到不 同的结果 ,因此 ,如何
选择最优 的小波基 函数非常重要 ; 外 , 此 如何 选择合适 的小 波分解层数 与去噪 阈值也 直接关 系到信号去 噪处 理 的质量 。
D A测序 信 号去 噪分 析 的一种 新方 法 N
郑 华 王立 强 , 岩 , 洁 , 一, 石 汪 陆祖康
1 江 大 学 现代 光 学 仪 器 国 家重 点 实 验 室 , 江 大 学 国家 光 学 仪 器 工程 技 术 研 究 中心 , 江 杭 州 3 02 .浙 浙 浙 10 7
维普资讯
第5 期
光谱学与光谱分析
12 17
量 , 以可以从 细节部分去除噪声 ,而不会 影响信 号 中的有 所
用成分 。
理电路的本底噪声 等。因此 我们 以实 际实验 系 统 ( 1 中, 图 ) D NA样 品电泳实验未出荧光 信号 时采集到 的噪声 信号 替代 () 8 式产生的随机噪声来构建 DN A测序仿真信号 ,对去噪效 果进 行更 加客观和准确地评估 。 由于 目前毛细管电泳激 光诱导荧 光 ( E LF 信号 的信 C -I)
2 .福建师范大学物理与光电信息科技学院 , 福建 福州 3 0 0 507

小波变换、神经网络和小波网络的函数逼近能力分析与比较

小波变换、神经网络和小波网络的函数逼近能力分析与比较
为 简 单 起 见, 这 里 以 人 工 神 经 网 络 中 最 常 用 的
Span{ 2 D( 2 t - k ) } ; V j+ 1 = V j G W j ; D( t) 是相 应 小波的尺度函数 , 由 上 面 的 讨 论 , 可 以 把 空 间 L2 ( R) 中 的 函 数 f ( t) 进行分解 [8] , f ( t) =
式中 , S 称为框架算子 , 在小 波 变 换 的 实 际 计 算 中 , 式 ( 6) 中 多 项 式 的 有 限项之和就可以对函数 f E L2 ( R) 进行有效的逼 近, 即,
N N
f=
j, k
zj, k Pj, k ( t) =
j, k
t - bj zj, k P aj
元 数分别为 n, n 1 , 1 的 BP 子网络 , 来逼近 Fourier 级
馈神经网络之间的关系 O 在此基础上讨论了小波网 络 的函数 逼近能 力 以 及 三 者 之 间 的 对 比 9 并 利 用 不 同的方法对一典型函数进行了仿真 O
1
小波变换的函数逼近能力
设 V n C 2 ( 92T) 9V n = span{ e iM 9.9e i 919 ei 9.9ei } 9 其中 : 2 ( 92T) 是指在 ( 92T) 内的平方 可 积空间 S n = + M + 1S 9M 为正整数 O 若函数 f ( I) E
i
k z xz
=
k
6k e
ik x
( 8) 当 N - O 时 , 满足 ,
N- O
f ( t) Pj, k dt, 在介
lim
[O , 1]
n
l f ( x ) - f F ( x , N, f ) l d x = O

小波讲义

小波讲义
16
5.2.4.尺度函数和小波函数的性质 在讨论MRA理论下的尺度函数 (t ) 和小波函数 (t ) 的性质之前,先介绍著名的Poisson公式。该公式用于描 述正数平移系列的正交归一性在频域的表示。 (1)设 f (t k ), k Z 是一组正交归一的函数集合:
f (t k ) f (t k )dt (k
小波讲义
由于 V0 V1 W1 ,对任意函数 f (t ) V0 , 可以分解为两部分:细节部分 W1 和大尺度 平滑逼近部分 V1 。 而大尺度平滑逼近部分可以近一步分解为细 节部分和大尺度平滑逼近部分。如此重复就 可以得到任意尺度(或分辨率)上的细节部 分和大尺度平滑逼近部分。这就是多分辨率 分析的思想。
18
1 ik1 i t 证明 f (t k1 ) F ( ) e e d 2 R 1 ik2 i t f (t k2 ) F ( ) e e d 2 R 将上面两式代入(5.33):
左边= 1 4 2 R R R 1 = 2 F ( ) F ( )[ ei ( )t dt ]e ik1 eik2 dd 4 R R R 1 = F ( ) F ( ) ( )]e ik1 eik2 dd 2 R R 1 2 i ( k k ) = F ( ) e d 2 R
2 i ( k k ) 2 1
d
(5.39)
(5.40)
=右面= (k1 k2 )

M ( ) F ( 2k )
kZ
2
,则(5.37)式可以看
成周期为 2 的函数 M ( ) 的Fourier级数系数的 表示式。(5.40)说明 M ( ) 只有常数项,而没有 谐波分量。故 M ( )的值恒定。所以下式几乎处处 成立: 2 (5.41) M ( ) F ( 2k ) 1

函数逼近与数据拟合

函数逼近与数据拟合

第三章 函数逼近与数据拟合整理:赖志柱本章继续讨论用简单函数近似代替较复杂函数的问题。

上章提到的插值就是近似代替的方法之一,插值的近似标准是在插值点处误差为零。

但在实际应用中,有时不要求具体某些点误差为零,而要求考虑整体的误差限制,这就引出了逼近和拟合的概念。

最佳平方逼近和最小二乘拟合有密切的关系,当后者的节点数无限加密时,拟合就变成了逼近。

§3.1 正交多项式正交多项式是函数逼近的重要工具,在数值积分中也有广泛的应用。

3.1.1 正交多项式概念 定义3.1 设函数()x ρ满足 (1)()0x ρ≥,(,)x a b ∀∈;(2)()bk a x x dx ρ⎰存在且有限的(0,1,k =);(3)若对于[,]a b 上的非负连续函数()g x ,有()()0bax g x dx ρ=⎰,则()0g x ≡;则()x ρ称为[,]a b 上的一个权函数。

定义3.2 设(),()[,]f x g x C a b ∈,()x ρ为[,]a b 上的权函数,称(,)()()()baf g x f x g x dx ρ=⎰为函数(),()f x g x 在[,]a b 上带权函数()x ρ的内积。

定义3.3 ()f x 在[,]a b 上的范数定义为1()|()|bafx f x dx ρ=⎰ 1222()[()]ba fx f x dx ρ⎡⎤=⎢⎥⎣⎦⎰[,]max ()x a b ff x ∞∈=定义3.4 若(,)0f g =,则称(),()f x g x 在[,]a b 上带权()x ρ正交,记为f g ⊥。

性质(1)Cauchy-Schwartz 不等式:22|(,)|f g fg ≤⋅;(2)三角不等式:222f g fg +≤+;(3)平行四边形定理:222222222()f g f g fg ++-=+; (4)若f 与g正交,则222222f gfg +=+; 定义3.5 若函数系{}{}01():0(),(),,(),n n x n x x x ϕϕϕϕ≥=满足22(,)()()()0bi j i j ai i j x x x dx i j ϕϕρϕϕϕ≠⎧⎪==⎨>=⎪⎩⎰则称{}{}01():0(),(),,(),n n x n x x x ϕϕϕϕ≥=是区间[,]a b 上的正交函数系。

数值分析函数逼近与曲线拟合

数值分析函数逼近与曲线拟合
y
f (x)
P1 ( x)
E1
a
x2
bx
最佳一次逼近多项式例题1(继续)
最佳一次逼近多项式例题2(返回)
切比雪夫定理图示(定理)
E2
P2 (x) f (x)
E3
E4
P4 (x) f (x)
P3 (x) f (x)
最佳平方逼近问题(返回)
法方程的建立(特例)
C[0,1]上的最佳平方逼近(例题)
C[0,1]上的最佳平方逼近例题(返回)
用正交函数做最佳平方逼近(返回)
最佳平方逼近多项式(例题)
最佳平方逼近多项式例题(返回)
线性模型例题(返回)
线性模型图例(返回)
指数模型例题(返回)
指数模型图例(返回)
双曲模型图例(返回)
S-曲线模型图例(返回)
§3.6最佳平方三角逼近与FFT(返回)
§3.2 正交多项式(返回)
正交函数族与正交多项式 正交多项式的性质 勒让德(Legendre)多项式 切比雪夫(Chebyshev)多项式 其他正交多项式
§3.3 最佳一致逼近多项式(返回)
偏差与偏差点 最佳一致逼近多项式 切比雪夫定理 最佳一致逼近多项式的构造 最佳一次逼近多项式
T0
T0
T3
T2 T3
TT11
T2
偏差与偏差点(返回)
最佳一致逼近多项式(返回)
切比雪夫定理(返回)
最佳一致逼近 多项式的构造(例题)
切比雪夫多项式 与零的偏差(定理)
最佳一致逼近例题(继续)
最佳一致逼近例题(返回)
最佳一次逼近多项式(例题)
最佳一次逼近多项式图示(返回)
哈尔(Haar)条件(法方程)

梁振动响应曲线滑动拟合法及在移动荷载识别中的应用

梁振动响应曲线滑动拟合法及在移动荷载识别中的应用

c uei i a ao o swi vn v rg t o a s s n lg u t mo ig a ea e meh d.Mo e e trs l s o h tt e mo ig f t g t h d lts eut h ws t a h vn i i tn
dvdn h e p n e it mals cin . Th o tn iy o h ep n tt es cin n d sc n b iiigt ers s n o s l e t s o o ec n iut ft e r s s a h e t o e a e o e o g aa te y t ef t g d t v r p e h e t n n d s u rn e d b h i i aao e l p d t es ci o e .Thsmeh d i c l d mo ig f t g b — tn a o i t o s al vn i i e e tn
好, 尤其是移动荷载 引起 的结构响应是各个进 出桥 时间不同的移动荷载引起结构 响应 的叠加 , 每个移 动荷载进或出桥时都造成响应 曲线上一个折点 , 在
此点虽然曲线连续 , 曲线的斜率和曲率都不连续。 但 因此 , 应该对响应 曲线进行分段拟合 , 如果对曲线的 细部拟合有较高的要求 , 分段 曲线应该较短 , 分段数 目较多。这时如果采用分段点曲线直到二阶导数连
meh sefciei vn a sie t ia in. t o i fet nmo ig l d d n i c t d v o f o Ke r s vb aina d wa e la sie t i t n; uv vn it g;ain l rcinf n t n ywo d : ir t n v ; d d ni c i c r emo ig f i rt a a t ci o o fa o tn o f o u o

函数逼近与曲线拟合

函数逼近与曲线拟合

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

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

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

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

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

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

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

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

这样,自然地会有下面两种不同的度量标准:⎰-=-ba nn dx x px f x p x f 22)]()([)()( (5-1)使用这个度量标准的函数逼近称为均方逼近或平方逼近;)()(max )()(x p x f x p x f n bx a n -=-≤≤∞(5-2)使用这个度量标准的函数逼近称为一致逼近或均匀逼近。

关于一致逼近的问题,在数学分析中有以下结论。

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

对于函数插值而言,如果插值余项也能满足对任意的0>ε,ε<-=)()()(x p x f x R n n 都成立的话,则插值多项式()n p x 是()f x 的一致逼近多项式。

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

第五章曲线拟合和函数逼近
0
注:多项式拟合实际上是选取
0 ( x), 1 (x ), , ( ) 为 1, x,, x 。 n x
n
11
非线性最小二乘拟合
可化为线性拟合问题的常见函数类型
1 y ae (a 0) 设y ln y, x x y 1 /(a be x )(a 0) 设x e x , y 1 / y
中求函数 P( x) akk ( x) ,使得平方和
k 0
n
Q i P( xi ) yi
i 0
2
n i akk ( xi ) yi min 。 k 0 i 0
m
2
10
问题如何求解系数 a0 , a1,, an ?
b/ x
y a bx (a ln a ) y a bx y ax b y a bx y ax2 bx c y ax2 bx c y a bx cx 2
12
y 1 /(ax b) y x /(ax b) y 2 ax2 bx c y x /(ax bx c)
i 0 i 0
3
3
i
2, xi 6, yi 8.2, xi yi 14.1,
2 i 0 i 0 i 0
3
3
3
4 2 a0 8.2 对应的法方程组为 a 2 6 14.1 1
得到 a0 1.05, a1 2 ,
Байду номын сангаас
( x ) ( x ) a ( x ) ( x ) a
i n i 0 i i n i 1 i

函数逼近与曲线拟合

函数逼近与曲线拟合

8
x p x ( p 时),
所以
x

也是
x
的特例
p
例1.求下列向量的各种常用范数
x (1,4,3,1)T
解: x 1 x1 x2 x4 9 x 2 ( x1 2 x2 2 x4 2 )12 27 3 3
x
max 1i4

a0 a1
an



m i0 m i0
m i0
yi xi yi
xin yi

称之为正规方程组或者法方程组。如何简化表示呢?
17
引入符号
1
x0
A=

1
x1
1 xm
x0n x1n



xmn
a0


=

a1



an
则正规方程组可简化为
y0

Y
=
y1



ym
AT A ATY
一个有趣的问题是AT能不能消去?
18
若强行消除,则正规方程组变化为
xi
4
9
设f (x) C[a,b],可定义范数如下:
f
(
b
1
f p (x) dx ) p ,
p
a
f (x)的p 范数, p 1
b
f f (x) dx
1
a
f (x)的1 范数
f
(
b
f
2 (x)
dx
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

{u1 , u2 ,⋯ , un } 线性无关.
Rn带权内积的定义
• 满足 ωi > 0 (i = 1,2,⋯ , n.)的实数组{ω1 , ω2 ,⋯ , ωn } 称为权系数。 Rn上加权内积定义为:
( x, y ) = ∑ ωi xi yi
i =1
n
相应的范数为
n 2 x 2 = ∑ ωi xi i =1
n
(ϕ n , ϕ j ) = 0 ( j = 0,1, ⋯ n − 1).
a j ( j = 0,1,⋯ n − 1.)
确定出
勒让德(Legendre)多项式
• [-1,1]上权函数 ρ ( x) ≡ 1 时,由 1, x, x 2 , ⋯, x n ,⋯ 正交化得到的多项式 P0 ( x), P ( x),⋯ 1 称为勒让德多项式。其解析表达式为:
∑α
j =1
n
j
u j = 0.
Th3的证明
G非奇异当且仅当齐次方程组
∑ (u
j =1
n
j
, u k ) α j = 0, k = 1, 2, ⋯ , n.
只有零解,即 ( ∑ α j u j , u k ) = 0, k = 1, 2, ⋯ , n.
j =1
n
只有零解,即
∑α
j =1
n
j
u j = 0 只有零解,即
所有定义在[ a,b ]上的连续函数做成的集合 在函数加法、数乘下做成无穷维线性空间, 记为 C[a,b] . C[a,b] 上内积定义为:
( f , g ) = ∫ f ( x ) g ( x ) dx
a
b
Rn上的范数举例
Ex. Rn上范数: 对X=(x1,…, xn )∈Rn , n 1. (1-范数): x = ∑ x
R = { (a1 , a2 ,⋯ , an ) ai ∈ R, i = 1,2,⋯ , n.}
n
全体n 维实向量的集合
x = (a1 , a2 ,⋯ , an ), y = (b1 , b2 ,⋯ , bn ), ( x, y ) = a1b1 + a2b2 + ⋯ + anbn .
内积空间C[a,b]
2
= =

n
x
2 i
i = 1

m a x
1 ≤ i ≤ n
x
in与C[a,b]Fra bibliotek范数的扩充关系 R
= = =
• 向量范数: 范数: 范数
x x x
1

n
x
i = 1
i
2

n
x
2 i
i = 1

max
x
1 ≤ i ≤ n
i
• 函数范数:
f f f
1
= = =

b a
f ( x ) dx
b 2 a
2 ∞
内积空间
设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正交.

=
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
0 ≤ (u + λv , u + λv) = (u, u) + 2λ(u, v) + λ 2 (v, v).
特别地,取 λ = − (u , v) (v, v) 代入上式得:
内积空间的性质
(u, v) (u, v) (u, u ) − 2 + ≥ 0, (v, v ) (v, v )
( u , v ) ≤ ( u , u )( v , v ) .
1
(n + 1) Pn +1 ( x) = (2n + 1) Pn ( x) − nPn −1 ( x)
各类最佳逼近的一般提法
• 设 Φ ⊆ C[a, b] 是某个子空间,对任 * f ( x) ∈ C[a, b],寻找 S ( x) ∈ Φ, 满足: f ( x ) − S * ( x ) = min f ( x ) − S ( x ) .
∑α
j =1
n
j
u j = 0.
Th3的证明

n
∑α
j =1
j
uj = 0 ⇒
n ∑ α j u j , u k = 0, k = 1, 2, ⋯ , n. j =1
Th3的证明
故,
n ∑ α ju j , j =1 ⇔
u k = 0, k = 1, 2, ⋯ , n.
• 在一般线性空间上定义了线性空 间的范数概念,然后推广到函数 空间C[a,b]上得到度量两个函数 间距离(差异)的多个范数。相应 得到多个逼近概念和每一逼近概 念下的最佳(误差最小)多项式的 特征或求法。
线性空间

设L是一个非空集合,K是实(或复)数域,并 可 在其上定义“加法”,“数乘”运算,而且满足以 下公理 • ฀ 加法交换律:x+y= y+x • ฀ 加法结合律:(x+y)+z= x+(y+z) • ฀ 存在零元:x+0=x • ฀ 存在逆元:x+(-x)=0 • ฀ 数乘:1x=x • ฀ a(bx)= (ab)x • ฀ (a+b)x=ax+bx • ฀ a(x+y)=ax+ay 则称L是数域K上的线性空间
S ( x )∈Φ
在不同的范数,不同的子空间 Φ 下形成了不同的最佳逼近.
最佳一致逼近多项式
对 f ( x) ∈ C[a, b], H n = span {1, x,⋯ , x n } ,满足:
f ( x) − P ( x) = max f ( x) − P ( x) =
* n ∞ a≤ x≤b * n P ( x)∈Hn n
(ϕ 0 , ϕ 0 ) (ϕ 0 , ϕ1 ) ⋯ (ϕ 0 , ϕ n ) (ϕ1 , ϕ 0 ) (ϕ1 , ϕ1 ) ⋯ (ϕ1 , ϕ n ) G = ⋮ ⋮ ⋮ (ϕ , ϕ ) (ϕ , ϕ ) ⋯ (ϕ , ϕ ) n 1 n n n 0
内积、范数、距离之间的关系
• 由每一内积可以导出一范数: x = ( x, x). • 由每一范数可以导出一距离: ρ ( x, y ) = x − y . 注:有的范数并没有导出它的内积。
n 维实(或复)Euclid空间Rn

在向量加法、数乘下为n维线性空间. Rn且 为距离空间,赋范线性空间,内积空间. Rn的内积为:
距离空间定义
• ฀ 设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, ρ).
C[a.b]带权内积的定义
• C[a,b]上带权 ρ (x ) 的内积和范数定义为:
( f ( x), g ( x)) = ∫ ρ ( x) f ( x) g ( x)dx,
a
b
f ( x) 2 =

b
a
ρ ( x) f ( x)dx .
2
C[a,b]上线性无关函数簇的格拉姆矩 阵判定定理
• C[a,b]上函数簇 {ϕ 0 , ϕ1 , ⋯ , ϕ n } 线性无关的充要条件是其格拉姆矩阵G非奇 异。
非奇异当且仅当 {u1 , u2 ,⋯ , un } 线性无关.
Th3的证明
[证明] 首先,
n ∑ α ju j , j =1 n ∑ α ju j , j =1
u k = 0, k = 1, 2, ⋯ , n. ⇒ n ∑1 α j u j = 0 ⇒ j=
范数与赋范线性空间
฀ 设X是实(或复)线性空间,如果对于X中 每 x ||x|| 个元素x,按照一定的法则对应于实数||x||, 且满足: • ฀ ||x||≥0,||x||=0当且仅当x=0; • ฀ ||ax||=|a|||x||,a是实(或复)数; • ฀ ||x+y||≤||x|| + ||y||. 则称X是实(或复)赋范线性空间,||x||称为x的 范数.
C[a.b]带权内积的定义
• [a,b]上的非负函数 ρ (x) 称为[a,b]上的权 函数,若满足: b ① ∫ x k ρ ( x)dx 存在且为有限值(k=0,1,…). a ② 对[a,b]上的非负连续函数g (x) ,若

b
a
ρ ( x ) g ( x ) dx = 0 ,
则 g ( x) ≡ 0
“函数逼近” 问题提法
• 对[a,b]上的连续函数类C[a,b]中的给 定函数f(x),在另一“简单”函数类----n次多项式函数类Hn中求函数p(x), 使p(x)与f(x)的误差在某种度量意义下 最小。 • p(x)与f(x)的误差的度量 的误差的度量是逼近论的核 的误差的度量 心问题。
相关文档
最新文档