数值分析--第3章函数逼近与快速傅里叶变换

例如将所有实 n维向量组成的集合，按向量加法及向量 与数的乘法构成实数域上的线性空间， 记作 R n，称为 n维 向量空对次数不超过 n( n为正整数)的实系数多项式全体， 按通常多项式与多项式加法及数与多项式乘法也构成数域 R上一个线性空间，用 H n 表示，称为多项式空间.
(u1,un)
(u2,u1) (u2,u2)
(u2,un)
(un,u1) (un,u2)
(un,un)
（1.7）
称为格拉姆(Gram) 则 非G奇异的充分必要条件是
u1,u线2,性无,u关n .
2020/6/16
课件
24
证明 方程组
G非奇异等价于deG t，其0充要条件是齐次
n
n
( juj,uk) (uj,uk)j 0, k1 ,2 ,,n（1.8）
（1.15）
1
f(x)2(f(x),f(x))2
b a
(x)
f
1
2(x)dx2
（1.16）
称（1.15）和（1.16）为带权 (x) 的内积和范数.
2020/6/16
课件
32
常用的是 (x) 1的情形，即
b
(f(x)g ,(x))af(x)g(x)d.x
1
f
(x)
2020/6/16
课件
20
定义了内积的线性空间称为内积空间. 定义中（1）的右端 ( u称, v为) 的(u共, 轭v)， 当K为实数域R时 (u,v).(v,u) 如果 (u,v，)则0称 与 正u 交，v 这是向量相互垂 直概念的推广.
2020/6/16
课件
21
定理2
设X为一个内积空间， 对u,v有X,
（1.3）
ln i m B n (m )(f,x)f(m )(x).
这个结果不但证明了定理1，而且由(1.3)给出了 f (x) 的一个逼近多项式，但它收敛太慢，实际中很少使用.
2020/6/16
课件
11
更一般地，可用一组在 C[a,b] 上线性无关的函数集合
i(x)in0 来逼近 f(x)C[a,b]，此时元素
使得
1x1 nxn0 ,
（1.1）
则称 x1, , xn 线性相关.
否则，若等式(1.1)只对 1 2 n0 成立，
则称 x1, , xn线性无关.
2020/6/16
课件
6
若线性空间 S是由 n个线性无关元素 x1, , xn生成的， xS 都有
x1x1nxn
则 x1, , xn称为空间 S的一组基，记为
所有定义在 [a, b] 上的连续函数集合，按函数加法和 数与函数乘法构成数域 R上的线性空间，记作 C[a,b].
类似地, 记 C p[a,b]为具有 p阶连续导数的函数空间.
2020/6/16
课件
5
定义1 设集合 S是数域 P上的线性空间，元素
x1,,xnS, 如果存在不全为零的数 1,,nP，
即ju11,u2,线,性un无j关1 .
n
在内积空ju间j X1u12u2nun 0
j1
uX,
记
u (u,u)
（1.10）
2020/6/16
课件
26
利用
(uv)2u2 2 uvv2 ( u ,u ) 2 ( u ,v ) ( v ,v ) (uv,uv)uv2
两端开方即得三角不等式
uv uv
取 (u,v，)/v (,v) 代入上式右端，得
2020/6/16
课件
22
(u,v)2 (u,v)2
(u,u)2
0
(v,v) (v,v)
即得 v 0时
(u,v)2(u,u)(v,v).
2020/6/16
课件
23
定理3 矩阵
设X为一个内积空间，u1,u2,,unX ,
(u1,u1) G (u1,u2)
g(x)0.
则称 (x) 为[a, b]上的一个权函数.
2020/6/16
课件
31
例2 C[a,b] 上的内积.
设 f(x)g ,(x) C [a,b],(x) 是 [a, b] 上给定的权函数,
则可定义内积
b
(f(x )g ,(x ) )a (x )f(x )g (x )d.x
由此内积导出的范数为
n
Bn(f,x)
k0
f(k n)Pk(x)
（1.3）
2020/6/16
课件
10
其中
Bn(P fk,(xx)) kn 0 knfx(k n k()1P k(xx))nk,
k nn(n1)k为(!n二k项式1)展开系数，并证明了
ln i m B n(f,x)f(x)在 [ 0 ,1上] 一致成立；
若 f在( x) 上[ 0 ,阶1]导数m 连续，则
R n 上的加权内积为
n
(x, y) ixi yi i1
（1.13）
相应的范数为
n
1
x ( 2
i xi2)2.
i1
当 i 1 (i1 ,2,,n)时，（1.13）就是前面定义的内积.
2020/6/16
课件
29
如果 x, yCn，带权内积定义为
n
(x, y) ixi yi i1
这里 { i } 仍为正实数序列，y i 为y i 的共轭.
而满足‖·‖1 =1 的向量 x1 x2 1则为对角线长 度为1的菱形.
2020/6/16
课件
16
所以说，范数是对向量长度的度量，度量方式不同， 结果也不一样，但不同范数之间是存在等价关系的.
2020/6/16
课件
17
类似地，对连续函数空间 C[a,b] ，若 f(x)C[a,b],
f mafx(x), axb
第3章 函数逼近与快速傅里叶变换
3.1 函数逼近的基本概念 3.2 正交多项式 3.3 最佳平方逼近 3.4 曲线拟合的最小二乘法 3.5 有理逼近 3.6 三角多项式与快速傅里叶变换
2020/6/16
课件
1
3.1 函数逼近的基本概念
3.1.1 函数逼近与函数空间
问题 1、数值计算中经常要计算函数值，如计算机中计算 基本初等函数及其他特殊函数；
在 C[a,b]上也可以类似定义带权内积.
（1.14）
2020/6/16
课件
30
定义4 设 [a, b] 是有限或无限区间，在 [a, b] 上的非负 函数 (x)满足条件：
（1）
b
a
xk
(x)dx存在且为有限值
(k0,1,),
（2） 对 [a, b] 上的非负连续函数 g (x) ，如果
bg(x)(x)dx0, a
xi ,
i 1
n
1
x ( 2
xi2)2 ,
i1
称为 范数或最大范数，
称为 1-范数， 称为 2-范数.
2020/6/16
课件
15
实际上任何向量的实值函数，只要满足上述三个条件，
就可以定义成一种向量范数. 在 R 2 中，满足‖·‖2 =1 ，即
为单位圆，
x1 2 的x2向2量1
满足‖·‖∞ =1 ，即 max1 x,x{ 2}1 的向量为单位正 方形，
b
f f(x)dx,
1
a
称为 范数，
称为 1-范数，
f
(
b
1
f 2(x)d x )2,
2
a
称为 2-范数.
可以验证这样定义的范数均满足定义2中的三个条件.
2020/6/16
课件
18
3.1.3 内积与内积空间
在线性代数中， 中R n两个向量 y(y1,,yn)T 的内积定义为
x(x1及,,xn)T
（1.11）
2020/6/16
课件
27
例1 R n 与 C n 的内积.
设 x, yRn, x(x1,,xn)T, y(y1,,yn)T,
n
(x, y) xi yi i1
（1.12）
向量2-范数为
1
n
1
x (x,x)2 ( 2
xi2)2
i1
2020/6/16
课件
28
若给定实数 i 0(i1 ,2,,n),称{ i }为权系数，
（1.2）
它由 n 1个系数 (a0,a1,,an)唯一确定.
1, x, , xn是线性无关的，它是 H n 的一组基，故
H nsp{ 1 a ,x,n ,xn},
且 (a0,a1,,an)是 p(x) 的坐标向量，H n 是 n 1维的.
2020/6/16
课件
8
对连续函数 f(x)C[a,b]，它不能用有限个线性无关的 函数表示，故 C[a,b] 是无限维的，但它的任一元素 f ( x) 均可用有限维的 p(x)Hn逼近，使误差
2020/6/16
课件
12
3.1.2 范数与赋范线性空间
为了对线性空间中元素大小进行衡量，需要引进范数 定义，它是 R空n 间中向量长度概念的直接推广.
2020/6/16
课件
13
定义2 设S为线性空间，xS，若存在唯一实数‖·‖，
满足条件：
（1） x 0, 当且仅当 x0时， x 0; （正定性）
(2 ( u ) , v ) ( u , v ), K , u , v X ;
(3 ( u v ) , w ) ( u , w ) ( v , w ) , u , v , w X ; (4) ( u ,u ) 0 , 当且 u 0 时 仅 ( u ,u ) ， 当 0 .
则称 (u,为v)X上 与u的内v积.