当前位置:首页 >> 工学 >>

函数平方逼近多项式的均方误差计算


函数平方逼近多项式的均方误差计算 函数平方逼近多项式的均方误差计算
1 1 + 25 x 2 (1) 求连续函数 f ( x) 在区间[-1,1]上的 3 次最佳平方逼近多项式,计算均方误 差 δc ;

实验要求:设 f ( x) =

(2) 在区间[-1,1]上取 5 个等距结点, f ( x) 的离散 3 次最佳平方逼近多项式, 求 计算均方误差 δ 5 ; (3) 在区间[-1,1]上取 9 个等距结点, f ( x) 的离散 3 次最佳平方逼近多项式, 求 计算均方误差 δ 9 ; (4) 比较 δ 5 和 δ 9 ,应如何合理地定义离散情况下的均方误差?该定义(1)中 的 δ c 有何关系? 实验步骤: (1) 利用 legendre 正交多项式作 f ( x) 在 [-1,1] 上的最佳平方逼近,先计算 ( f ( x), Pk ( x)) , k =0,1,2,3。 1 1 2 ( f ( x), P0 ( x)) = ∫ dx = arctan(5) 2 ?1 1 + 25 x 5 1 1 ( f ( x), P ( x)) = ∫ x dx = 0 1 ?1 1 + 25 x 2 1 3 1 1 2 2 ( f ( x), P2 ( x)) = ∫ ( x 2 ? ) dx = ? × arctan(5) 2 ?1 2 2 1 + 25 x 25 125 1 5 3 1 ( f ( x), P3 ( x)) = ∫ ( x 3 ? ) dx = 0 ?1 2 2 1 + 25 x 2 由方程组计算出系数

?2 ?0 ? 2 ?3 ? ?0

0
2 3

2 3

0
2 5

0
2 5

0

2 ? a rc ta n ( 5 ) 0? ? ? 5 ? ? 2 ? 0 5 ? ? a = ? ( 2 ? a rc ta n ( 5 )) 2 ? 25 ? ? 0? 125 ? ? 2 ? 0 ? ? 7 ?

解得系数为

a = [0.50923 0 ?0.70366 0] 由以上系数和 legendre 正交函数簇可得 f ( x) 在 [-1,1]上的三次最佳平方逼近多 项式 S ( x) 为:
? * S ( x) = a0 P0 ( x) + a2 P2 ( x) = 0.50923 ? 0.70366 x 2

均方误差:

δ c ( x) 2 = f ( x) ? S ( x) = 0.27243 ?
f ( x) 和多项式 S ( x) 在区间[-1,1]上的图形如下所示:

1.2

1

0.8

0.6

0.4

0.2

0

-0.2 -1

-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1

(2) 在区间[-1,1]上等距的取 5 个点 x -1.0000 -0.5000 0 y 0.0385 0.1379 1.0000 由法方程 Ha = d ,得到如下方程组: 0 2.5000 0 ? ?5.0000 ?1.3528 ? ? 0 ? ? 0 ? 2.5000 0 2.1250 ? ? ? a=? ? 2.5000 ? 0.1459 ? 0 2.1250 0 ? ? ? ? ? 2.1250 0 2.0313? ? 0 ? 0 ? 解得系数为 a = [ 0.4855 0 ?0.5609 0] 由此得到 均方误差
S ( x) = 0.5737 ? 0.6063x 2

0.5000 0.1379

1.0000 0.0385

δ5 2 =



4 k =0

( f ( xk ) ? S ( xk )) 2 =0.5945

f ( x) 和多项式 S ( x) 在区间[-1,1]上的图形如下所示:

1.2

1

0.8

0.6

0.4

0.2

0

-0.2 -1

-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1

(3) 在区间[-1,1]上等距的取 9 个点 x -1.000 -0.750 -0.500 -0.250 0 0 0 0 y 0.0385 0.0664 0.1379 0.3902

0 1.000 0

0.250 0 0.390 2

0.500 0 0.137 9

0.750 0 0.066 4

1.000 0 0.038 5

由法方程 Ha = d ,得到如下方程组: 0 3.7500 0 ? ?9.0000 ? 2.2661? ? 0 ? ? 0 ? 3.7500 0 2.7656 ? ? ? a=? ?3.7500 ?0.2694 ? 0 2.7656 0 ? ? ? ? ? 2.7656 0 2.3877 ? ? 0 ? 0 ? 解得系数为 a = [ 0.4855 0 ?0.5609 0] 由此得到 均方误差 S ( x) = 0.4855 ? 0.5609 x 2

δ9 2 =



8 k =0

( f ( xk ) ? S ( xk ))2 = 0.6367

f ( x) 和多项式 S ( x) 在区间[-1,1]上的图形如下所示:

1.2

1

0.8

0.6

0.4

0.2

0

-0.2 -1

-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1

(4) 由 (2) 、 可以看到 δ 5 < δ 9 , (3) 这意味着在该均方误差定义下, 得到的 f ( x) 的信息越多,反而拟合出的曲线误差越大。 可以将离散情况下均方差定义为

δn =



n ?1 k =0

hk ( f ( xk ) ? S ( xk )) 2

其中 hk 是 xk +1 和 xk 之间的距离,,采用均匀分布,在区间[a,b]上 n 个采样点, b?a hk = h = , 在此定义下, , [a b]=[-1,1]时, 5 =0.4203, 9 =0.3183, δ 9 < δ 5 , δ δ 有 n ?1 并且有:当 n 趋近于无穷时, δ n = δ c


相关文章:
函数平方逼近多项式的均方误差计算.doc
函数平方逼近多项式的均方误差计算 - 函数平方逼近多项式的均方误差计算 函数平方逼近多项式的均方误差计算 1 1 + 25 x 2 (1) 求连续函数 f ( x) 在区间[...
数值分析第一次作业(2).doc
2 上求一函数,使其为 x2 ? C[0,1]的最佳...的经验公式,使它与下列数据拟合,并 计算均方误差。...sin ? x 在[0,1]上的二次 最佳平方逼近多项式...
第三章_函数逼近与计算.ppt
? P ? x ? 2 2 作为度量误差“大小”标准的函数逼近称 为平方逼近均方逼近。 §2 正交多项式一、正交多项式 1.正交 若 ? x, y ? ? 0,则称 x与 ...
数值分析2最佳逼近和最小二乘法.pdf
|| 的计算取决于该内积空间中范数(即内积)的定义。 >>例如连续函数构成的空间...中的最佳平方逼近多项式均方误差。 解: 记? 0 ? 1, ?1 ? x 2 , ?...
第3章+函数逼近与计算1_图文.ppt
2. n ? 1次多项式的n次最佳逼近。 一般情形求近似最佳逼近。 15 2n ? 4个未知量的非线性方程组(很难求解)。 16 §3 函数平方逼近均方误差最小作为...
第3章 函数逼近与计算2_图文.ppt
第3章 函数逼近计算2 - §3 函数平方逼近均方误差最小作为度量标准, 用均方误差最小作为度量标准,研究函数 f ( x) ∈ C[a, b] 的逼近多项式,就...
最佳平方逼近算例.pdf
最佳平方逼近例 - 例:求函数 f ( x) = e x 在[0,1]上的二次最佳平方逼近多项式,并估计平 ,小数点后保留 5 位. 方逼近误差 δ 2 2 解: (解法 ...
计算方法总结.doc
? ai* ( f , ?i ) i ?0 n 求区间 [?1,1] 上函数 f ( x) ? x 在 M ? span ?1, x2 , x4 ? 中的最佳平方逼近 多项式均方误差。 解:...
矩阵理论作业6:两种算法求三次最佳平方逼近多项式.doc
两种算法求 f ( x) 的三次最佳平方逼近多项式摘 要 对于一个较复杂的函数,往往需要求一个简单多项式来逼近。本文选取两种基函数,用两种算 法计算一个函数 f ...
数值分析 第三章 函数逼近与计算_图文.ppt
? P ? x ? 2 2 作为度量误差“大小”标准的函数逼近称 为平方逼近均方逼近。 五、一致逼近多项式的存在性定理 1(维尔斯特拉斯定理) 若f (x)是区间[a...
计算方法 第五章第二节最佳平方逼近_图文.ppt
平方逼近问题,即求拟合 函数,使之与函数 f ( x) 在[a, b] 上的均方差...( x) 正交多项式可由递推公式递推计算。 推论 设{g n ( x)} 为正交...
数值分析第3章 函数逼近与计算.ppt
a b 2 2 作为度量误差“大小”标准的函数逼近称为最佳平 方逼近或均方逼近...x ? n 次最佳平方逼近多项式。 * * 求最佳平方逼近函数 S ( x) ? ? a...
数值分析(计算方法)第七章_图文.ppt
(二) 平方逼近: 采用 ? b [ f ( x) ? p( x)] dx 2 a 作为度量误差的“大小”的标准的函数逼近称为平方逼近均方逼近。 §1 正交多项式 一、...
最佳平方逼近_图文.ppt
类 B 通常是代数多项式,有理多项式,三角多 项式,分段多项式等容易计算的函数。...2、平方逼近(均方逼近) 以 ? ? f (x) ? a b p ( x ) ? dx 2 ...
最佳平方逼近多项式.ppt
最佳平方逼近多项式 - 最佳平方逼近多项式 §5.2 最佳平方逼近多项式 ? 本节内容 1.内积空间 2.两类特殊的函数族 3.函数的最佳平方逼近 4.举例 5.MATLAB程序...
数值计算方法_最佳平方逼近.doc
数值计算方法_最佳平方逼近 - 数值分析最佳逼近 ━基于MATLAB的实现与分析 §1引言 所谓函数最佳逼近就是从指定的一类简单的函数中寻找一个和给定的函 数“最...
第3章+函数逼近与计算1_图文.ppt
要求所得 用插值的方法对这一函数进行近似 要求所得 到的插值多项式经过已知的...§2 函数平方逼近均方误差最小作为度量标准, 用均方误差最小作为度量标准,...
数值计算方法第七章习题 2013.doc
(x) = 1 的最佳平方逼近二次多项式。 20.求函数 f (x)在指定区间上对于?...97.8 用最小二乘法求形如 y = a + bx 2 的经验公式,并计算均方程误差...
研究生数值分析(20)函数的最佳平方逼近_图文.ppt
k k ? 0,1,?, n Legendre多项式的应用 对给定的函数 f (x)∈ C [a,b] ,在区间 [a,b]上求f (x)的n次最佳平方逼近多项式 pn (x) 这个问题,相当...
2011数值分析第一次作业及参考答案.doc
已测得函数 y = f ( x) 的三对数据: (0,1...的经验公式,使它与下列数据拟合,并计算均方 误差。...sin π x 在[0,1]上的二次最佳平方 逼近多项式...
更多相关标签: