我怎么能把一个非常小的符号数(vpa)双没有四舍五入到零?
6视图(30天)
显示旧的评论
我需要用修改后的第一类贝塞尔函数的计算。在我的情况下,
ν
参数的函数,
I = besseli(ν,Z)
~ 1800。为了得到正确的近似,我需要使用精度高、128位的
vpa
,
数字(128);
besselResult = besseli (vpa (abs(ν)),2 * c1_sqrt);% abs(ν)~ 1800
这个数,
besselResult
涉及到的订单大约10 ^ -2400。对于速度的考虑,我不想使用符号计算,所以想把
besselResult
马上到一个双(这说得通吗?)。但这结果四舍五入为0。我当然可以减少精密数字很少,工作时保留指数的符号形式:
besselResult = vpa (besselResult 5);
我怎样才能保留尾数(尽管在低精度)和指数仍然能够将它转换成一个双吗?
