数值分析课件 第3章5-7节.ppt

3.5 曲线拟合的最小二乘法 3.5.2 用正交多项式做最小二乘拟合 3.6 最佳平方三角逼近与快速傅里叶变换 3.6.1 最佳平方三角逼近与三角插值 3.6.2 快速傅氏变换（FFT） 3.7 有 理 逼 近 3.7.2 帕德逼近 另一方面,若对（7.2）式用辗转相除可得到 的 一种连分式展开 （7.3） （7.2） （7.4） （7.3）右端为 的无穷连分式的前5项，最后式子 若取（7.3）的前2，4，6，8项，则可分别得到 的以下有理逼近: 是它的紧凑形式. 并逐步把 累加到 中去，最后就可得到所求的 用这种方法编程序不用解方程组，只用递推公式，并 且当逼近次数增加一次时，只要把程序中循环数加1，其余 不用改变. 这里 可事先给定或在计算过程中根据误差确定. 拟合曲线 当 是周期函数时，显然用三角多项式逼近 比用代数多项式更合适，本节主要讨论用三角多项式做最 小平方逼近及快速傅里叶变换，简称FFT算法. 设 是以 为周期的平方可积函数，用三角多 项式 （6.1） 作为最佳平方逼近函数. 由于三角函数族 在 上是正交函数族，于是 在 上的最小 平方三角逼近多项式 的系数是 称为傅里叶系数. 函数 按傅里叶系数展开得到的级数 （6.3） 就称为傅里叶级数. （6.2） 只要 在 上分段连续，则级数（6.3）一致 收敛到 . 对于最佳平方逼近多项式（6.1）有 由此可以得到相应于（4.11）的贝塞尔不等式 因为右边不依赖于 ，左边单调有界，所以级数 （6.3） （6.1） （4.11） 当 只在给定的离散点集 上已知时，则可类似得到离散点集正交性与相应的离散傅 里叶系数. 下面只给出奇数个点的情形. 收敛，并有 可以证明对任何 成立 令 这表明函数族 在点集 上正交. 若令 其中 则 的最小二 乘三角逼近为 当 时 于是 （6.4） 就是三角插值多项式，系数仍由（6.4）表示. 由于 所以函数族 在区间 上是正交的. 一般情形，假定 是以 为周期的复函数，给定 在 个等分点 上的值 函数 在等距点集 上的值 组成的向量记作 当 时， 个复向量 具有 如下正交性： （6.5） 事实上，令 于是 即 若 若 则有 则 从而 于是 若 这就证明了（6.5）成立. 即 是正交的. 则 于是 因此， 在 个点 上的 最小二乘傅里叶逼近为 （6.5） （6.6） 其中 （6.7） 在（6.6）中， 若 ， 则 为 在点 上的插值函数， 于是由（6.6）得 即 （6.8） 而（6.8）是由 求 的过程，称为反变换. （6.7）是由 求 的过程， 称为 的离散 傅里叶变换. 简称DFT， （6.7） （6.8） 不论是按（6.7）式由 求 ， 由 求 ， （6.9） 其中 (正变换) 或 (反变换)， 还是由（6.4）计算傅里叶逼近系数 都可归结为计算 是已知复数序列. 或是按（6.8） （6.4） 当 较大且处理数据很多时，就是用高速的电子计算 机，很多实际问题仍然无法计算， 如直接用（6.9）计算 ，需要 次复数乘法和 次 复数加法，称为 个操作，计算全部 共要 个操作. 直到20世纪60年代中期产生了FFT算法，大大提高了 运算速度，从而使傅氏变换得以广泛应用. FFT算法的基本思想就是尽量减少乘法次数. （6.9） 用（6.9）计算全部 ， 表面看要做 个乘法， 实际上所有