实验要求
···误差的来源?误差的类型?(模型误差?截断误差?舍入误差﹖浮点运算舍入误差?)
···误差的度量方法:相对误差、绝对误差。
···理解迭代序列的收敛性?误差的收敛阶(定义与表达),以及阶的估计表达。
···误差的传播途径、误差的累积、局部误差、总体误差等。
实验内容
分析讨论题
题目一:
求方程 x 2 + ( α + β ) x + 1 0 9 = 0 x^2+(\alpha+\beta)x+10^9=0 x2+(α+β)x+109=0的根,其中 α = − 1 0 9 , β = − 1 \alpha=-10^9 ,\beta=-1 α=−109,β=−1,讨论如何设计计算格式才能有效地减少误差,提高计算精度。
解答:
对于求解一元二次方程,不妨使用求根公式求解:
x
=
−
b
±
b
2
±
4
a
c
2
a
x=\frac{-b \pm \sqrt{b^2 \pm 4ac}}{2a}
x=2a−b±b2±4ac
因为 α = − 1 0 9 , β = − 1 \alpha=-10^9,\beta=-1 α=−109,β=−1
所以 ∣ b ∣ ≈ b 2 − 4 a c |b|≈\sqrt{b^2-4ac} ∣b∣≈b2−4ac 且b<0,从而计算 x 2 x_2 x2应用原公式的改进形式为:
x 2 = − 2 c b − b 2 − 4 a c x_2=\frac{-2c}{b - \sqrt{b^2-4ac}} x2=b−b2−4ac−2c
所以解得方程的解为:
x 1 = 109 x_1=109 x1=109
x
2
=
1
x_2=1
x2=1
对于求根公式求解方程的根,要根据实际情况来讨论是否修改其形式,对于
∣
b
∣
=
b
2
−
4
a
c
|b|=\sqrt{b^2-4ac}
∣b∣=b2−4ac时,若要修改其形式,当b>0时,
x
1
x_1
x1的求解公式为:
x 1 = − 2 c b + b 2 − 4 a c x_1=\frac{-2c}{b + \sqrt{b^2-4ac}} x1=b+b2−4ac−2c
当b<0是, x 2 x_2 x2的求根公式为:
x 2 = − 2 c b − b 2 − 4 a c x_2=\frac{-2c}{b - \sqrt{b^2-4ac}} x2=b−b2−4ac−2c
题目二
以计算 x 31 x^{31} x31为例,讨论如何设计计算格式才能减少计算次数。
解答:
因为
31
=
2
0
+
2
1
+
2
2
+
2
3
+
2
4
31=2^0+2^1+2^2+2^3+2^4
31=20+21+22+23+24
所以,原式的计算过程可写为:
β 1 = 1 , α 1 = x β_1=1 , α_1=x β1=1,α1=x
β 2 = β 1 α 1 , α 2 = α 1 ∗ α 1 β_2=β_1α_1 , α_2=α_1*α_1 β2=β1α1,α2=α1∗α1
β 3 = β 2 α 2 , α 3 = α 2 ∗ α 2 β_3=β_2α_2 , α_3=α_2*α_2 β3=β2α2,α3=α2∗α2
β 4 = β 3 α 3 , α 4 = α 3 ∗ α 3 β_4=β_3α_3 , α_4=α_3*α_3 β4=β3α3,α4=α3∗α3
β 5 = β 4 α 4 , α 5 = α 4 ∗ α 4 β_5=β_4α_4 , α_5=α_4*α_4 β5=β4α4,α5=α4∗α4
β 6 = β 5 α 5 , α 6 = α 5 ∗ α 5 β_6=β_5α_5 , α_6=α_5*α_5 β6=β5α5,α6=α5∗α5
所以一共需要9次乘法就可以完成,这样就可以减少计算次数
完整实验报告
数值计算实验报告完整版1-9章:
链接:数值计算实验报告