函数逼近试验-Read.DOCVIP

函数逼近实验 学号：P0725005 姓名：乔海亮 1. 问题描述 2 2. 解决方法 2 2.1. 最佳平方逼近二次多项式 2 2.2. Chebyshev截断级数 3 2.3. 插值余项极小化 4 3. 数值结果 4 4. 结果分析对比 6 5. 代码和程序说明 6 5.1. 开发环境 6 5.2. 主要变量和函数说明 6 5.3. 程序说明 7 问题描述 设，试求出权函数的最佳平方逼近二次多项式。另外请用Chebyshev截断级数的办法和插值余项极小化方法分别给出近似最佳一致逼近二次多项式。并画出所有曲线图。 解决方法 最佳平方逼近二次多项式 设所求多项式为 其中 则 即有法方程 求解得到 所求最佳平方逼近二次多项式为 最大误差 Chebyshev截断级数 利用Chebyshev多项式进行函数值的计算，主要是利用Chebyshev的系数趋于零的速度比较快，从而级数收敛速度快。 由于讨论区间为，故不能直接用Chebyshev截断，作变换 则有 函数变为 用Chebyshev多项式展开为 其中，系数 Chebyshev多项式为 利用数值积分，得 此时Chebyshev截断级数为 最大误差 插值余项极小化 根据Chebyshev多项式的性质，在区间上的所有最高项系数为1的n次多项式中，与零的偏差最小，其偏差为。我们以的n个零点作为插值节点，这里取三个插值节点（），构造Lagrange多项式。 其中 经数值计算，得最大误差 数值结果 x f(x)原函数 最佳平方逼近 Chebyshev截断 插值余项极小化 1.0 0.000000 -0.021609 -0.013370 -0.015203 1.1 0.104841 0.095702 0.102909 0.101685 1.2 0.2187858 0.2182072 0.2244184 0.2237418 1.3 0.3410735 0.3459083 0.3511588 0.3509683 1.4 0.4710611 0.4788048 0.4831300 0.4833641 1.5 0.6081976 0.6168967 0.6203319 0.6209292 1.6 0.7520058 0.7601840 0.7627645 0.7636636 1.7 0.9020680 0.9086668 0.9104279 0.9115673 1.8 1.058015 1.062344 1.063321 1.064640 1.9 1.219522 1.221218 1.221446 1.222882 2.0 1.386294 1.385287 1.384802 1.386294 2.1 1.558068 1.554551 1.553388 1.554875 2.2 1.734606 1.729011 1.727205 1.728625 2.3 1.915690 1.908666 1.906253 1.907545 2.4 2.101124 2.093517 2.090531 2.091633 2.5 2.290726 2.283563 2.280041 2.280892 2.6 2.484329 2.478804 2.474780 2.475319 2.7 2.681779 2.679241 2.674751 2.674916 2.8 2.882934 2.884873 2.879953 2.879682 2.9 3.087661 3.095701 3.090385 3.089617 3.0 3.295836 3.311724 3.306048 3.304722 Error（误差） 0.0 0.02161 0.01337 0.01520 结果分析对比 由上图可以看出，尽管只做到二次多项式逼近，但三种逼近方法给出的结果，与原函数基本一致，在图中看来，四条曲线已经重叠在一起。之所以看起来全是蓝色，是由于作图顺序的先后顺序造成的。 画图的顺序为：原函数曲线（黑色）、最佳平方逼近二次多项式（红色）、Chebyshev截断级数（绿色）、插值余项极小化（蓝色）。 代码和程序说明 开发环境 Windows XP SP2 + Visual Studio 2005 主要变量和函数说明 变量： double* m_dblX; // x坐标 double* m_dblF; // 原函数 double* m_dblBAS; // 最佳平方逼近 Best Approximation Square double* m_dblCheTr;// Chebyshev截断 Chebyshev Truncation double* m_dblCheZI;//