Chapter 1 误差理论

实验要求

···误差的来源?误差的类型?(模型误差?截断误差?舍入误差﹖浮点运算舍入误差?)
···误差的度量方法:相对误差、绝对误差。
···理解迭代序列的收敛性?误差的收敛阶(定义与表达),以及阶的估计表达。
···误差的传播途径、误差的累积、局部误差、总体误差等。

实验内容

figure

分析讨论题

题目一:

求方程 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=2ab±b2±4ac

因为 α = − 1 0 9 , β = − 1 \alpha=-10^9,\beta=-1 α=109,β=1

所以 ∣ b ∣ ≈ b 2 − 4 a c |b|≈\sqrt{b^2-4ac} bb24ac 且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=bb24ac 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=b24ac 时,若要修改其形式,当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+b24ac 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=bb24ac 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章:
链接:数值计算实验报告

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值