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
2008信息与计算科学专业计算方法习题参考解答 江世宏编
第四章 数值积分
姓名 学号 班级
习题主要考察点：代数精度的计算，构造插值型求积公式（梯形，辛甫生公式），复化求积的计算，高斯公式的构造。 1给定求积公式
?h?hf(x)dx?af(?h)?bf(0)?cf(h)试确定a,b,c使它的代数精度尽可能
高。（代数精度的应用和计算）
解：分别取f(x)?1,x,x2，使上述数值积分公式准确成立，有；
?a?b?c?2h? ?a(?h)?c(h)?0?a(?h)2?c(h)2?2h3/3?h4hh,c?。 解得：a?,b?333hh4hhf(0)?f(h)。 故求积公式为?f(x)dx?f(?h)??h333hh4hh333?0?(h)3?0 再取f(x)?x，左边=?xdx?0，右边=(?h)??h333再取f(x)?x42h5h4hh42h54?0?(h)?，左边=?xdx?，右边=(?h)?
?h53333h4此求积公式的最高代数精度为3。 2 求积公式
?10f(x)dx?A0f(0)?A1f(1)?B0f?(0)，试确定系数A0,A1及B0，使该求积
公式具有尽可能高的代数精确度，并给出代数精确度的次数。（代数精度的应用和计算） 解：分别取f(x)?1,x,x，使求积公式准确成立，有
2?A0?A1?1??A1?B0?1/2 ?A?1/3?1211,A1?,B0?。 3361211求积公式为?f(x)dx?f(0)?f(1)?f?(0)。
03361121133再取f(x)?x，左边=?xdx???0??1??0?右边
04336解得：A0?故该求积公式的最高代数精度为2。 3数值积分公式
?303f(x)dx?[f(1)?f(2)]，是否为插值型求积公式，为什么？又该公式
213
2008信息与计算科学专业计算方法习题参考解答 江世宏编
的代数精确度为多少？（插值型求积公式特征）
3解：令f(x)?1，dx?3?03?33[1?1]?[f(1)?f(2)] 22f(x)?x，?xdx?03933?[1?2]?[f(1)?f(2)] 2221533?[1?22]?[f(1)?f(2)] 222bf(x)?x，?x2dx?9?20故代数精度为1。由于求积节点个数为2，代数精度达到1次，故它是插值型的求积公式。 4如果f??(x)?0，证明用梯形公式计算积分几何意义。（梯形求积） 解：梯形求积公式
?af(x)dx所得到的结果比准确值大，并说明其
T?b?a[f(a)?f(b)] 2是由过点(a,f(a))，(b,f(b))的线性插值函数
L(x)?x?bx?af(a)?f(b) a?bb?a在[a,b]上的定积分。
注意到：在区间[a,b]上，f??(x)?0，而(x?a)(x?b)?0，有
bbbbI?T??f(x)dx??L(x)dx??[f(x)?L(x)]dx??aaaaf??(?)(x?a)(x?b)dx?0 2!从而I?T。
其几何意义可作以下解释：
在区间[a,b]上，f??(x)?0，故曲线y?f(x)下凹，直线y?L(x)位于曲线之上，因
bb此，曲边梯形的面积I??f(x)dx小于梯形面积T??L(x)dx。
aa
14
2008信息与计算科学专业计算方法习题参考解答 江世宏编
1?1xdx，并估计误差。（复化梯形求积）
2?111?，取求积节点为xi?1??i(i?0,1,?,4) 解：h?4445用n?4的复化梯形公式计算积分
2?2131dx??xi?0xi?1xi?31h11dx??[f(xi)?f(xi?1)]?h[f(x0)?f(x1)?f(x2)?f(x3)?f(x4)]x22i?02114444141171?[??????]??0.6970 424567281680因
?211dx?ln2，则误差大约为：ln2?0.6970?0.0039。 x6设f(?1)?1,f(?0.5)?4,f(0)?6,f(0.5)?9,f(1)?2,则用复化辛甫生公式计算
?1?1f(x)dx，若有常数M使 |f(4)|?M，则估计复化辛甫生公式的整体截断误差限。（复
化辛甫生公式） 解：
?101?1f(x)dx??1?f(x)dx??f(x)dx
0?[141141f(?1)?f(?0.5)?f(0)]?[f(0)?f(0.5)?f(1)] 666666167?[1?4?4?6?6?4?9?2]??11.1667 660I?S2??1?f(4)(?1)(x?1)(x?0.5)2(x?0)dx?4!011?0f(4)(?2)(x?0)(x?0.5)2(x?1)dx 4!M?[?(x?1)(x?0.5)2(x?0)dx??(x?0)(x?0.5)2(x?1)dx] 24?10M?12M(x?0)(x?0.5)(x?1)dx??60210.522?t(0.25?t)dt?0M?0.0042 6?0.008M
17已知高斯求积公式
?1?f(x)dx?f(0.57735)?f(?0.57735) 将区间[0,1]二等分，用复
1化高斯求积法求定积分
11/2?01xdx的近似值。（高斯公式）
解：
?0xdx??0xdx?1/2?xdx
15