内容发布更新时间 : 2024/12/23 0:20:14星期一 下面是文章的全部内容请认真阅读。
2008信息与计算科学专业计算方法习题参考解答 江世宏编
第三章 函数逼近
姓名 学号 班级
习题主要考察点:最小二乘法,最佳平方逼近,正交多项式的构造。
1 设f(x)?sin?x,求f(x)于[0,1]上的线性最佳平方逼近多项式。(最佳平方逼近) 解:??span{1,x}
11(?1,?1)??dx?1,(?1,?2)??xdx?,(?2,?2)??x2dx?
230001111(f,?1)??sin?xdx?021?,(f,?2)??xsin?xdx??0x?cos?x?11?2sin?x?01?
法方程组为
1??2?2??a1?????
??1?1??a??2???3????2解得:a1?,a2?0
???1?1??2线性最佳平方逼近多项式为:??*2?。
x2 令f(x)?e,?1?x?1,且设p(x)?a0?a1x,求a0,a1使得p(x)为f(x)于[?1,1]
上的最佳平方逼近多项式。(最佳平方逼近) 解:??span{1,x}
2(?1,?1)??dx?2,(?1,?2)??xdx?0,(?2,?2)??x2dx?
3?1?1?11x?11111(f,?1)??edx?e?e,(f,?2)??xexdx?2e?1
?1?1法方程组为
?20??a1??e?e?1?2??????1? ?0a?3????2??2e?1e?1?1解得:a1?(e?e),a2?
23e?e?1e?1?x。 线性最佳平方逼近多项式为:p(x)?23
9
2008信息与计算科学专业计算方法习题参考解答 江世宏编
3证明:切比雪夫多项式序列
Tk(x)?cos(karccosx)
在区间??1,1?上带权?(x)?1/1?x2正交。(正交多项式的证明) 解:对于l?k,有
1(Tl,Tk)??10?11?x2cos(larccosx)cos(karccosx)dx
?2???11?cost?cos(lt)cos(kt)(?sint)dt??cos(lt)cos(kt)dt
01??[cos(l?k)t?cos(l?k)t]dt 20111?[sin(l?k)t?sin(l?k)t]?0?0 2l?kl?k对于l?k,有
1(Tk,Tk)??10??11?x21cos2(karccosx)dx
????1?cos2tcos(kt)(?sint)dt??cos2(kt)dt
02111???[1?cos(2k)t]dt?[t?sin(2k)t]?? 02022k2故,序列{Tk(x)}在[-1,1]上带权?(x)?11?x2正交。
?x1?x2?3?4求矛盾方程组:?x1?2x2?4的最小二乘解。(最小二乘法)
?x?x?22?1解法一:求x1与x2,使得
f(x1,x2)?(x1?x2?3)2?(x1?2x2?4)2?(x1?x2?2)2
达到最小。于是,令
?f?2(x1?x2?3)?2(x1?2x2?4)?2(x1?x2?2)?0 ?x1 10
2008信息与计算科学专业计算方法习题参考解答 江世宏编
?f?2(x1?x2?3)?2(x1?2x2?4)?2?2(x1?x2?2)(?1)?0 ?x2即:??3x1?2x2?9?x1?2.5714,其最小二乘解为:?。
?2x1?6x2?9?x2?0.6429解法二:
?11??3??12??x1???4?,记作AX?b,该矛盾方程组的最小二乘解,应满足以下方程组 ???x????2??2???1?1?????32??x1??9?AAX?Ab,即???x???9?
26???2???TT解之,得??x1?2.5714。
?x2?0.64295 已知一组试验数据
xk 2 4 2.5 4.5 3 6 4 8 5 8.5 5.5 9 yk 试用直线拟合这组数据. (计算过程保留3位小数)。(最小二乘线性逼近) 解:作矩阵
?12??4??12.5??4.5??????13??6?A???,y???
14???8??15??8.5?????15.5??????9??法方程为
(ATA)X?(ATy)
即
22??a??40??6?2290.5???b???161.25? ??????解得:a?1.2288,b?1.4831。
x。 其直线拟合函数为y?1.2288?1.48316 用最小二乘原理求一个形如y?a?bx的经验公式,使与下列数据相拟合.
2 11
2008信息与计算科学专业计算方法习题参考解答 江世宏编
xk 19 19 25 32.3 31 49 38 44 73.3 97.8 yk (最小二乘二次逼近) 解:等价于对数据表
2 xk361 19 625 32.3 961 49 1444 1936 73.3 97.8 yk 作线性拟合。其法方程组为:
5327??a??271.4??5?53277277699???b???369321?
.5??????解得:a?0.9726,b?0.0500 故经验公式为 y?0.9726?0.05x2。
12