数值分析(6)-函数逼近的基本概念

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/kevin_zhao_zl/article/details/90739508

整理一下数值分析的笔记~
目录:

1. 误差
2. 多项式插值与样条插值
3. 函数逼近(THIS)
4. 数值积分与数值微分
5. 线性方程组的直接解法
6. 线性方程组的迭代解法
7. 非线性方程求根
8. 特征值和特征向量的计算
9. 常微分方程初值问题的数值解

1. 函数逼近与函数空间

1.1 常见的一次近似式 (|x|很小时,泰勒多项式)

(1) n 1 + x z + 1 n x ^n\sqrt{1+x} \approx z+\frac{1}{n}x

(2) s i n x x ( x ) sinx \approx x(x为弧度)

(3) t a n x x ( x ) tanx \approx x(x为弧度)

(4) e x 1 + x e^x \approx 1+x

(5) l n ( 1 + x ) x ln(1+x) \approx x

1.2 函数逼近

  对函数类A中给定的 f ( x ) f(x) ,记作 f ( x ) A f(x) \in A ,要在另一类简单的便于计算的函数类B中求函数 p ( x ) B p(x) \in B ,使 f ( x ) f(x) p ( x ) p(x) 误差在某种度量意义下最小。

   函数类A通常是区间[a,b]上的连续函数,记作 C [ a , b ] C[a,b] ,称为来连续函数空间。函数类B通常为n次多项式,有理函数或分段低次多项式等。

   数学上常把各种集合中引入某些不同的确定关系称为赋予集合以某种空间结构,并将这样的集合称为空间。

   对次数不超过 n n 的实系数多项式全体,按通常多项式与多项式假发及数与多项式乘法构成数域R上的一个线性空间,称为多项式空间。

   所有定义在 [ a , b ] [a,b] 上的连续函数集合,按函数加法和数与函数乘法构成数域R上的线性空间记作 C [ a , b ] C[a,b] C p [ a , b ] C^p[a,b] 称为具有p阶连续导数的函数空间。

扫描二维码关注公众号,回复: 6417368 查看本文章

定义1:设集合 S S 是数域 P P 上的线性空间,元素 x 1 , . . . , x n S x_1,...,x_n \in S ,如果存在不全为零的数 α 1 , . . . , α n P \alpha_1,...,\alpha_n \in P 使得 α 1 x 1 + . . . + α n x n = 0 \alpha_1x_1+...+\alpha_nx_n=0 则称 x 1 , . . . , x n x_1,...,x_n 线性相关,若只对 α 1 = . . . = α n = 0 \alpha_1=...=\alpha_n=0 则称为线性无关。

  若线性空间 S S 是由 n n 个线性无关的元素生成的,即对所有的 x S x\in S ,都有 x = α 1 x 1 + . . . + α n x n x=\alpha_1x_1+...+\alpha_nx_n ,则 x 1 , . . . x n x_1,...x_n 为空间 S S 的一组基,记为 S = s p a n S=span { x 1 , . . . , x n x_1,...,x_n },并称 S S n n 维空间,系数 α 1 , . . . , α n \alpha_1,...,\alpha_n 称为x在基 x 1 , . . . , x n x_1,...,x_n 下的坐标,记作 ( α 1 , . . . , α n ) (\alpha_1,...,\alpha_n) ,如果 S S 中有无限个线性无关元素 x 1 , . . . , x n , . . . x_1,...,x_n,... 则称S为无限维线性空间。

维尔斯特拉斯定理:设 f ( x ) C [ a , b ] f(x)\in C[a,b] ,则对任何 ε > 0 \varepsilon>0 ,总存在一个代数多项式 p ( x ) p(x) ,使得 f ( x ) p ( x ) < ε ||f(x)-p(x)||_{\infty}<\varepsilon [ a , b ] [a,b] 上成立。

Bernstein逼近函数:伯恩斯坦根据函数整体逼近得特性构造出伯恩斯坦多项式:

{ B n ( f , x ) = k = 0 n f ( k n ) P k ( x ) P k ( x ) = C n k x k ( 1 x ) m k \begin{cases} B_n(f,x)=\sum_{k=0}^nf\left(\frac{k}{n}\right)P_k(x)\\ P_k(x)=C^k_nx^k(1-x)^{m-k} \end{cases}

l i m n B n ( f , x ) = f ( x ) , x [ 0 , 1 ] lim_{n \rarr \infty} B_n(f,x)=f(x),x \in[0,1] 成立,但是收敛速度慢且收敛依赖于多项式次数 n n\rarr\infty

2. 范数与赋范线性空间

定义2: 设S为线性空间, x S x\in S 若存在唯一实数, ||\cdot|| 满足条件:

  • x 0 ||x||\geq 0 ,当且仅当x=0时 x = 0 ||x||=0 (正定性)

  • α x = α x α R ||\alpha x||=|\alpha|||x||,\alpha \in R (齐次性)

  • x + y x + y x , y S ||x+y||\leq||x||+||y||,x,y\in S (三角不等式)

  则称 ||\cdot|| 为线性空间 S S 上得范数, S S ||\cdot|| 一起称为赋范线性空间,记为 X X

  在 R n R^n 上得向量 x ( x 1 , . . . , x n ) T R n x\in (x_1,...,x_n)^T \in R^n ,三种常用范数为 x = m a x 1 i n x i 最大范数,||x||_{\infty}=max_{1\leq i \leq n}|x_i| 1 x 1 = i = 1 n x i 1-范数 ,||x||_1=\sum_{i=1}^n|x_i| , 2 x 2 = ( i = 1 n x i 2 ) 1 2 2-范数, ||x||_2=(\sum_{i=1}^nx_i^2)^{\frac{1}{2}} ,类似地对连续函数空间 C [ a , b ] C[a,b] ,若 f ( x ) C [ a , b ] f(x) \in C[a,b] 可定义三种常用范数 f i n f t y = m a x a x b x i 无穷范数,||f||_{infty}=max_{a \leq x \leq b}|x_i| 1 f 1 = a b f ( x ) d x 1-范数,||f||_1=\int_a^b|f(x)|dx 2 f 2 = ( a b f 2 ( x ) d x ) 1 2 2-范数,||f||_2=(\int_a^bf^2(x)dx)^{\frac{1}{2}}

3. 内积和内积空间

定义3 X X 是数域 K K 上的线性空间,对所有的 u v X , K u,v\in X,K 中一个数与之对应,记为 ( u , v ) (u,v) 它满足以下条件:

  1. ( u , v ) = ( u , v ) ˉ ( ) ( u , v ) X (u,v)=\bar{(u,v)}(共轭),(u,v)\in X

  2. ( α u , v ) = α ( u , v ) , α K , u , v K (\alpha u,v)=\alpha(u,v),\alpha \in K,u,v \in K

  3. ( u + v , w ) = ( u , w ) + ( v , w ) u , v , w X (u+v,w)=(u,w)+(v,w),u,v,w\in X

  4. ( u , u ) 0 u = 0 ( u , u ) = 0 (u,u)\geq 0当且仅当u=0时(u,u)=0

则称 ( u , v ) (u,v) X X u u v v 的内积,内积的线性空间称为内积空间,当 K K 为实数域 R R ( u , v ) = ( v , u ) (u,v)=(v,u) ,若 ( u , v ) = 0 (u,v)=0 则称 u u v v 正交。

柯西-施瓦茨不等式:设 X X 为一个内积空间,对所有的 u , v X u,v\in X ( u , v ) 2 ( u , u ) ( v , v ) |(u,v)|^2 \leq (u,u)(v,v) ,证明:

v = 0 v=0 显然成立;

v ! = 0 v!=0 ( u , v ) > 0 (u,v)>0 且对任何数 λ \lambda 0 ( u + λ v , u + λ v ) = ( u , u ) + 2 λ ( u , v ) + λ 2 ( v , v ) 0\leq(u+\lambda v,u+\lambda v)=(u,u)+2\lambda(u,v)+\lambda^2(v,v) ,取 λ = ( u , v ) / ( v , v ) \lambda=-(u,v)/(v,v) 代入上式右端得 ( u , u ) 2 ( u , v ) 2 ( v , v ) + ( u , v ) 2 ( v , v ) 0 (u,u)-2\frac{|(u,v)|^2}{(v,v)}+\frac{|(u,v)|^2}{(v,v)}\geq0 ( u , v ) 2 ( u , u ) ( v , v ) |(u,v)|^2 \leq(u,u)(v,v)

定理3:设 X X 为一个内积空间 μ 1 , . . . μ n X \mu_1,...\mu_n \in X ,矩阵:

G = [ ( μ 1 , μ 1 ) ( μ 2 , μ 1 ) . . . ( μ n , μ 1 ) ( μ 1 , μ 2 ) ( μ 2 , μ 2 ) . . . ( μ n , μ 2 . . . . . . . . . . . . ( μ 1 , μ n ) ( μ 2 , μ n ) . . . ( μ n , μ n ) ] G=\left[ \begin{matrix} (\mu_1,\mu_1) & (\mu_2,\mu_1) & ... & (\mu_n,\mu_1) \\ (\mu_1,\mu_2) & (\mu_2,\mu_2) & ... &(\mu_n,\mu_2 \\ ...&...&...&...\\ (\mu_1,\mu_n)&(\mu_2,\mu_n) &... & (\mu_n,\mu_n) \end{matrix} \right]

称为格拉姆矩阵,G非奇异的充分必要条件是 μ 1 , . . . , μ n \mu_1,...,\mu_n 线性无关。

定义:设 [ a , b ] [a,b] 是有限或无限区间,在 [ a , b ] [a,b] 上的非负函数 ρ ( x ) \rho(x) 满足条件:

  1. a b x k ρ ( x ) d x \int_a^bx^k\rho(x)dx 存在且为有限值(k=0,1,2,…)

  2. [ a , b ] [a,b] 上的非负连续函数 g ( x ) g(x) ,如果 a b g ( x ) ρ ( x ) d x = 0 , g ( x ) = 0 \int_a^bg(x)\rho(x)dx=0,则g(x)=0

则称 ρ ( x ) \rho(x) [ a , b ] [a,b] 上的一个权函数。

4. 逼近的度量标准

4.1 最佳一致逼近多项式

f C [ a , b ] f\in C[a,b] ,寻求多项式 P n ( x ) P^*_n(x) ,使其误差: f P n = m a x a x b f ( x ) P n ( x ) = m i n P n H n f P n ||f-P^*_n||_{\infty}=max_{a\leq x\leq b}|f(x)-P^*_n(x)|=min_{P_n \in H_n}||f-P_n||

这种意义下的函数逼近称为最佳一致逼近或均匀逼近。

4.2 最佳平方逼近

采用这个标准: f ( x ) S ( x ) 2 2 = m i n S ( x ) f ( x ) S ( x ) 2 2 = m i n S ( x ) a b ρ ( x ) [ f ( x ) S ( x ) ] 2 d x ||f(x)-S^*(x)||_2^2=min_{S(x) \in \empty}||f(x)-S(x)||_2^2\\=min_{S(x)\in \empty}\int_a^b \rho(x)[f(x)-S(x)]^2dx

这种意义下的函数逼近称为最佳平方逼近或均方逼近。

4.3 最小二乘拟合

采用这个标准: f ( x ) P ( x ) 2 2 = m i n P ( x ) e m p t y f ( x ) P ( x ) 2 2 = m i n P ( x ) i = 0 m [ f ( x i ) P ( x i ) ] 2 ||f(x)-P^*(x)||_2^2=min_{P(x)\in empty}||f(x)-P(x)||^2_2\\=min_{P(x)\in \empty}\sum^m_{i=0}[f(x_i)-P(x_i)]^2 ,这种意义下的函数逼近称为最小二乘拟合。


{持续更新}
欢迎扫描二维码关注微信公众号 深度学习与数学   [每天获取免费的大数据、AI等相关的学习资源、经典和最新的深度学习相关的论文研读,算法和其他互联网技能的学习,概率论、线性代数等高等数学知识的回顾]
在这里插入图片描述

猜你喜欢

转载自blog.csdn.net/kevin_zhao_zl/article/details/90739508