第三讲：函数逼近快速傅里叶变换.pptVIP

第三章 函数逼近与快速傅里叶变换 /*Chapter3 Function Approximation and Fast Fourier Transformation */;§3.1.1 函数逼近与函数空间 ; 记作 ，; 定义1; 如果 中有无限个线性无关元素 则称 ;（1.2）;使误差; 更一般地，可用一组在 上线性无关的函数集合 ; §3.1.2 范数与赋范线性空间; 例如，在 上的向量 三种常 用范数为 ; 而满足‖·‖1 =1 的向量 则为对角线长 度为1的菱形.; 类似地，对连续函数空间 ，若; 在线性代数中， 中两个向量 及; 定义3; 定义中（1）的右端 称为 的共轭，; 定理3; 证明; 从以上等价关系知，;两端开方即得三角不等式 ; 例1;相应的范数为 ; 如果 ，; 定义4; 例2; 常用的是 的情形，即 ; 若 是 中的线性无关函数族，;3.2 正交多项式/*Orthogonal Polynomials*/ ;三角函数族 ; 例如, 三角函数族 1, cosx, sinx, cos2x, sin2x, …是在区间[-π,π]上的正交函数系，因为对k=1,2,…有;满足关系式(2.2)，;　　这样得到的正交多项式?n(x)，其最高次项系数为1. 反之，若 {?n(x)}是正交多项式，则 ?0,?1,…,?n在区间[a, b]上是线性无关的.;由于;3.2.2 勒让德(Legendre)多项式 ;勒让德多项式重要性质： ; 性质2;3.2.3 切比雪夫(Chebyshev)多项式 ; 切比雪夫多项式有很多重要性质： ; 性质7; 解 利用T5(x)的零点和区间变换可知节点;利用(2.15)式可得误差估计;3.2.5 其他常用的正交多项式;令x=cosθ，可得;2. 拉盖尔多项式;3. 埃尔米特多项式;3.3 最佳平方逼近;法方程组 /*normal equations */; Hilbert阵！; 解 由公式有;故得所求的一次最佳平方逼近多项式为;3.3.2 用正交函数族作最佳平方逼近;立刻得到法方程的解为;于是f(x)∈C[a, b]在?中的最佳平方逼近函数为;若f(x)∈C[a, b], 按正交函数族; 解 先计算(f(x), Pk(x)) (k=0,1,2,3);得三次最佳平方逼近多项式为; 如果f(x)∈C[a, b]，求[a, b]上的最佳平方逼近多项式，做变换;3.4 曲线拟合的最小二乘法;确定多项式 ，对于一组数据(xi, yi) (i = 1, 2, …, n) 使得 达到极小，这里 n m。;§3.4 L-S Approximating Polynomials;方案二：设; 解 根据所给数据，在坐标纸上标出各点，见图. 从图中看到各点在一条直线附近，故可选择线性函数做拟合曲线，即令;解得a0=2.5648,a1=1.2037. 于是所求拟合曲线为; 例 设数据(xi, yi)(i=0,1,2,3,4)由下表给出，表中第4行为lnyi=zi，可以看出数学模型为y=aebx，用最小二乘法确定a及b.; 对方程y=aebx 两边取对数得lny=lna+bx，如果令z=lny,A=lna，则z=A+bx,?={1,x}. 为确定A和b，先将数据(xi, yi)可转化为(xi, zi) ，见数据表.;故有法方程;解 作散点图如右，;代入法方程(4.6)式有;解之可得;P94 习题：4、5、8、12、17、18 计算实习题：1、2