用正交多项式做最小二乘拟合.pptxVIP

1 接下来根据给定节点 及权函数 构造带权 正交的多项式 . 注意 ，用递推公式表示 ，即（5.10）根据 的这里 是首项系数为1的 次多项式，正交性，得第1页/共57页第一页，共58页。 2（5.11） 下面用归纳法证明这样给出的 是正交的. 第2页/共57页第二页，共58页。 3 假定 对 及要证 对 均成立. 由（5.10）有 由（5.10）第二式及（5.11）中 的表达式，有 均成立，（5.12）第3页/共57页第三页，共58页。 4 而 ，于是由（5.12），当 时， 另外， 是首项系数为1的 次多项式，它可由由归纳法假定，当 时的线性组合表示.由归纳法假定又有第4页/共57页第四页，共58页。 5由假定有 再考虑 （5.13）利用（5.11）中 表达式及以上结果，得 第5页/共57页第五页，共58页。 6至此已证明了由（5.10）及（5.11）确定的多项式  组成一个关于点集 的正交系. 用正交多项式 的线性组合作最小二乘曲线拟合，只要根据公式（5.10）及（5.11）逐步求 的同时，相应计算出系数最后，由 和 的表达式（5.11）有 第6页/共57页第六页，共58页。 7并逐步把 累加到 中去，最后就可得到所求的 用这种方法编程序不用解方程组，只用递推公式，并且当逼近次数增加一次时，只要把程序中循环数加1，其余不用改变. 这里 可事先给定或在计算过程中根据误差确定. 拟合曲线第7页/共57页第七页，共58页。 83.6 最佳平方三角逼近与快速傅里叶变换 当 是周期函数时，显然用三角多项式逼近 比用代数多项式更合适，本节主要讨论用三角多项式做最小平方逼近及快速傅里叶变换，简称FFT算法. 第8页/共57页第八页，共58页。 9最佳平方三角逼近与三角插值 设 是以 为周期的平方可积函数，用三角多项式 （6.1）做最佳平方逼近函数. 由于三角函数族 在 上是正交函数族，于是 在 上的最小平方三角逼近多项式 的系数是 第9页/共57页第九页，共58页。 10 称为傅里叶系数. 函数 按傅里叶系数展开得到的级数 （6.3）就称为傅里叶级数.（6.2）第10页/共57页第十页，共58页。 11 只要 在 上分段连续，则级数（6.3）一致收敛到 . 对于最佳平方逼近多项式（6.1）有 由此可以得到相应于（4.11）的贝塞尔不等式 因为右边不依赖于 ，左边单调有界，所以级数 第11页/共57页第十一页，共58页。 12 当 只在给定的离散点集 上已知时，则可类似得到离散点集正交性与相应的离散傅里叶系数. 下面只给出奇数个点的情形. 收敛，并有 第12页/共57页第十二页，共58页。 13可以证明对任何 成立 令第13页/共57页第十三页，共58页。 14这表明函数族 在点集上正交. 若令则 的最小二其中 乘三角逼近为第14页/共57页第十四页，共58页。 15当 时 于是 （6.4）就是三角插值多项式，系数仍由（6.4）表示. 第15页/共57页第十五页，共58页。 16由于 所以函数族 在区间 上是正交的. 一般情形，假定 是以 为周期的复函数，给定 在 个等分点 上的值函数 在等距点集 上的值组成的向量记作第16页/共57页第十六页，共58页。 17当 时， 个复向量 具有如下正交性： （6.5）第17页/共57页第十七页，共58页。 18事实上，令于是 即 若若则有则从而第18页/共57页第十八页，共58页。 19于是 若这就证明了（6.5）成立. 即 是正交的. 则于是 因此， 在 个点