计算方法参考答案2011

四川工业学院研究生课程考试试题

课程名称： 计算方法 考试类型（考试或考查）： 考试 年 级： 2000 学时： 80 考试时间： 120 专 业： 学生姓名: 学号:

一、(8分)计算f??2?1，取2?1.4，利用下列等式计算，哪一个得到的结果最好？

?61?

解: 方法一

2?1?6， 3?2?2，

?31?3?22?3， 99?70。2

设y?(x?1)，若x?若通过

612，p?1.4，则|e(p)|??10?1。

21?2?1?6计算y值，即计算函数f(x)?(x?1)?6在x?2处的值，由于

f?(x)??6?x?1?，故e(f(p))?|f?(p)|e(p)，得

|e(f(1.4))|??6(1.4?1)?7|e(p)|?0.013|e(p)|；（2分）

若3?22?7??计算y值，即计算函数f(x)?(3?2x)323在x?2处的值，由于

f?(x)??6?3?2x?，故得

|e(f(1.4))|??6(3?2?1.4)2|e(p)|?0.24|e(p)|；（2分）

若3?22???3计算y值，即计算函数f(x)?(3?2x)?3在x??42处的值，由于

f?(x)??6?3?2x?，故得

|e(f(1.4))|??6(3?2?1.4)?4|e(p)|?0.005|e(p)|；（2分）

若99?702计算y值，即计算函数f(x)?99?70x在x?2处的值，由于

f?(x)??70，故得

（2分） |e(f(1.4))|??70|e(p)|?70|e(p)|；

比较4个结果得，通过

1?3?22?3计算得到的结果最好。

第 1 页(共 页)

方法二

根据数值计算原则：（1）避免两个相近的数相减；（3分）（2）简化计算步骤，减少运算次数。（3分） 可以判断得出：通过

1?3?22?3计算得到的结果最好。（2分）

二、(10分) 设

?1?1?1???A???11?1?

??1?11???计算||A||1，||A||2，||A||?

解： ||A||1?3，||A||??3（6分）。由于A是对称矩阵，所以||A||2?max(|?(A)|)（2分），而A的特征值为?1,2,2，故||A||2?2（2分）

三、(10分)为求方程x?x?1?0在1.5附近的根，设将方程改写成下列等价形式，并建立相应的迭代公式。 （1）x?1?3211，迭代公式； （2）x?x?1?k?122xxk1，迭代公式xk?1?x?11。 xk?1试分析每种迭代公式的收敛性。 解：考虑区间[1.3,1.6]（2分）。

3[1,.6.] （1）当x?1g(x)?1?时，

122??[1.3,1.6]，|g(x)|????0.910?1，233xx1.3故迭代xk?1?1?1在[1.3,1.6]上收敛。（4分） 2xk1， x?1 （2）当x?[1.3,1.6]时，令g(x)?|g?(x)|??11??1.076?1。

2(x?1)322(1.6?1)32故迭代xk?1?

1发散。（4分） xk?1第 2 页(共 页)

四、(10分)用LU分解法求解方程组：

?x1?x2?x3?1??x1?2x2?2x3?0。 ??2x?x?x?1123??11?1???解：设A??12?2?，X?(x1,x2,x3)T，b?(1,0,1)T，Y?(y1,y2,y3)T。

??211??? （1）对A进行LU分解，

?100??11?1?????A?LU??110??01?1?。（4分）

??231??002????? （2）解方程LY?b得：Y?(1,?1,6)T；（3分） （3）解方程UX?Y得：X?(2,2,3)T（3分）

五、(10分)线性方程组：

?x1?2x2?2x3?1??x1?x2?x3?3 ?2x?2x?x?5123?考察用Jacobi迭代和Gauss-Seidel迭代解此方程组的收敛性。

?0?22???解：1、BJ???10?1?，BJ的特征值为?1??2??3?0，?(BJ)?0?1。故Jacobi

??2?20???迭代收敛。（5分） 2、BGS?0?22?????02?3?，BGS的特征值为?1??2?2，?3?0，?(BGS)?1。故?002???Gauss-Seidel迭代发散。（5分）

六、(10分)现在你没有计算器，也没有计算机。请你用线性插值和二次插值计算出43的近似值，并估计误差。

解：问题可化为：已知函数f(x)?插值法求函数f(x)?（1）线性插值：

第 3 页(共 页)

x在xi?0,1,4,9,16,?的值yi?0,1,2,3,4,?，通过

x在x?43的值。

43介于36与49之间，以x0?36，x1?49为节点，进行线性插值得

f(x)?L1(x)?x?42，于是43?6.5385， 131(43?36)(49?43)?0.0243。误差为：|R1(42)|?（5分）

2?4?3632（2）二次插值：

离43最近的三个节点为：25，36 49，用这三个节点进行二次失插值得

f(x)?L2(x)，于是43?6.5629，

误差为：|R1(43)|?

七、(8分)已知y?f(x)的数据表如下：

3(43?25)(43?36)(49?43)?0.0061。（5分） 526?8?36xi yi -2 0 -1 0.2 0 0.5 1 0.8 2 1 求一次式p(x)?ax?b，使得p(x)为f(x)的最小二乘一次近似；

??2???1解：记A??0??1?2?1??0????1??0.2??a?1?，y??0.5?，c???，要使p(x)?ax?b为f(x)的最小二乘一次????b?1?0.8???1?1????近似，a,b为超定方程Ac?y的最小二乘解。（4分）

TT 超定方程Ac?y的法方程AAc?Ay为??100??a??2.6???????,解之

?05??b??2.5?得:a?0.26,b?0.5 （2分），最后得p(x)?0.26x?0.5为f(x)的最小二乘一次近似。（2分）

八、八、(8分)如果f??(x)?0，证明用梯形公式计算积分I??f(x)dx所得结果比准确值

abI大，并说明其几何意义。

证明：由梯形公式的余项

(b?a)3RT(f)??f??(?)，??(a,b)，（3分）

12知若f??(x)?0，则RT(f)?0，因而

第 4 页(共 页)