Fortran Coder

标题: fortran怎么把角度值规范到-pi到pi范围内？ [打印本页]

作者: mmSCI 时间: 2020-9-11 11:30
标题: fortran怎么把角度值规范到-pi到pi范围内？
通过微分方程数值解求到的角度值，怎么把角度值通过取余运算，规范到-pi到pi。我试了取余运算mod，当角度是负数的时候，mod(-3.15,6.28)=-3.15并不是我想要的3.13

作者: 风平老涡 时间: 2020-9-11 21:14
mod(-3.15,6.28) = -3.15 - int(-3.15 / 6.28) * 6.28 = -3.15 - int(-0.5) * 6.28 = -3.15 - 0 * 6.28 = -3.15
-3.15 就是正确的余数啊。

作者: li913 时间: 2020-9-11 21:27
本帖最后由 li913 于 2020-9-11 21:52 编辑

b = modulo(a,2*PI); if(b>PI) b = b - 2*PI

QQ截图20200911214951.png (10.71 KB, 下载次数: 391)

作者: mmSCI 时间: 2020-9-12 10:38
感谢以上各位铁友的回答，我想了一个简单的方法与大家共享，就是先求角度值得正弦和余弦值sin(A)，cos(A)，然后利用atan2计算角度值，自然就能规范到-pi到pi，即A=atan2(sin(A),cos(A))

作者: mmSCI 时间: 2020-9-12 10:39

风平老涡发表于 2020-9-11 21:14
mod(-3.15,6.28) = -3.15 - int(-3.15 / 6.28) * 6.28 = -3.15 - int(-0.5) * 6.28 = -3.15 - 0 * 6.28 = - ...

感谢，请看我的回答

作者: li913 时间: 2020-9-12 10:40
本帖最后由 li913 于 2020-9-12 10:41 编辑

mmSCI 发表于 2020-9-12 10:38
感谢以上各位铁友的回答，我想了一个简单的方法与大家共享，就是先求角度值得正弦和余弦值sin(A)，cos(A)， ...

用三角函数太慢，还会降低精度。

作者: mmSCI 时间: 2020-9-12 10:40

li913 发表于 2020-9-11 21:27
b = modulo(a,2*PI); if(b>PI) b = b - 2*PI

感谢，请看我下面的回答

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)