1. 练习
作为练习函数和循环的简单途径,用牛顿法实现开方函数。在这个例子中,牛顿法通过选择一个初始点z,然后重复这一选择过程求Sqrt(x)的近似值,这里用图片描述牛顿法的数学公式:
为了做到这个,只需要重复计算 10 次,并且观察在不同的循环次数(1,2,3,……)下z是如何逐步逼近结果的。 然后,修改循环条件,使得当z值停止改变(或改变非常小)的时候退出循环。观察迭代次数是否变化,结果与math.Sqrt接近。
提示:定义并初始化一个浮点值,向其提供一个浮点语法或使用转换:
z := float64(1)
或
z := 1.0
2. 解题思路
定义并初始化一个浮点值z,并套用计算方式: z = z - (z*z-x)/(2*z),循环计算并查看变化(迭代10次,到结果没有变化),目的是得到类似math.Sqrt()的效果。
3. 代码实现
// 循环与函数练习:用牛顿法实现开方
package main
import (
"fmt"
"math"
)
func Sqrt(x float64) float64 {
z := 1.00 // z用来套用公式z = z - (z*z-x)/(2*z)
temp := 0.00 // 临时变量,用来记录z上次的值
for i := 0; i < 10; i++ {
z = z - (z*z-x)/(2*z) // 计算最新的z值
fmt.Println(z) // 观察不同循环序数的z值,如何逐步逼近结果
if math.Abs(z-temp) < 0.000000000000001 {
break // 当值停止改变(或改变非常小)的时候退出循环
} else {
temp = z // 把每次循环的z值赋给temp,做条件比较
}
}
return z // 当z值停止改变(或改变非常小)时跳出循环,返回z值
}
// 比较牛顿法与普通开方函数的输出
func main() {
fmt.Println("牛顿法:", Sqrt(2))
fmt.Println("math.Sqrt(2):", math.Sqrt(2))
}
输出结果为
1.5
1.4166666666666667
1.4142156862745099
1.4142135623746899
1.4142135623730951
1.414213562373095
牛顿法: 1.414213562373095
math.Sqrt(2): 1.4142135623730951