如何用Go的math/big库实现高精度公式Ao=Bo*(1-((Bi/(Bi+Ai))**(Wi/Wo)))

你的核心问题是将big.Float转换为float64调用math.Pow导致精度丢失——float64有效位数有限，无法容纳极大数的全部精度。要全程基于big.Float完成计算，可利用指数与对数的数学关系：x^y = e^(y * ln(x))，结合big.Float的Log和Exp方法实现高精度幂运算。

以下是修正后的代码：

package main

import (
	"math/big"
)

func calculateAo(Bo, Bi, Ai, Wi, Wo *big.Float) *big.Float {
	// 根据实际需求设置足够精度，数值越大精度需越高，示例用256位
	precision := uint(256)
	
	// 初始化所有临时变量并统一设置精度
	one := new(big.Float).SetPrec(precision).SetInt64(1)
	temp0 := new(big.Float).SetPrec(precision)
	wiDivWo := new(big.Float).SetPrec(precision)
	logTemp0 := new(big.Float).SetPrec(precision)
	expVal := new(big.Float).SetPrec(precision)
	subResult := new(big.Float).SetPrec(precision)
	Ao := new(big.Float).SetPrec(precision)

	// 计算 Bi/(Bi+Ai)
	temp0.Add(Bi, Ai)
	temp0.Quo(Bi, temp0)

	// 计算 Wi/Wo
	wiDivWo.Quo(Wi, Wo)

	// 计算 ln(temp0) → 再乘以 Wi/Wo → 最后取指数得到 temp0^(Wi/Wo)
	logTemp0.Log(temp0)
	logTemp0.Mul(logTemp0, wiDivWo)
	expVal.Exp(logTemp0, nil)

	// 计算 1 - temp0^(Wi/Wo)
	subResult.Sub(one, expVal)

	// 最终计算 Ao = Bo * (1 - ...)
	Ao.Mul(Bo, subResult)

	return Ao
}

关键注意事项： #

• 精度设置：必须为每个big.Float实例调用SetPrec，默认64位精度无法满足极大数计算需求，需根据数值规模调整精度值（如512位、1024位）。
• 无类型转换：所有计算全程基于big.Float完成，彻底避免了向float64转换的精度损失。
• 幂运算替代方案：由于big.Float未内置幂运算方法，通过Log+Exp的数学转换是实现高精度幂运算的标准方式。

内容的提问来源于stack exchange，提问作者Oluwatobi Femi