Fortran Coder

查看: 16605|回复: 6
打印 上一主题 下一主题

[求助] fortran怎么把角度值规范到-pi到pi范围内?

[复制链接]

33

帖子

12

主题

0

精华

熟手

F 币
138 元
贡献
93 点
跳转到指定楼层
楼主
发表于 2020-9-11 11:30:21 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
通过微分方程数值解求到的角度值,怎么把角度值通过取余运算,规范到-pi到pi。我试了取余运算mod,当角度是负数的时候,mod(-3.15,6.28)=-3.15并不是我想要的3.13
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

213

帖子

2

主题

0

精华

宗师

F 币
2142 元
贡献
875 点

规矩勋章

沙发
发表于 2020-9-11 21:14:03 | 只看该作者
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 就是正确的余数啊。

835

帖子

2

主题

0

精华

大宗师

F 币
3926 元
贡献
2334 点
板凳
发表于 2020-9-11 21:27:40 | 只看该作者
本帖最后由 li913 于 2020-9-11 21:52 编辑

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

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

QQ截图20200911214951.png

33

帖子

12

主题

0

精华

熟手

F 币
138 元
贡献
93 点
地板
 楼主| 发表于 2020-9-12 10:38:38 | 只看该作者
感谢以上各位铁友的回答,我想了一个简单的方法与大家共享,就是先求角度值得正弦和余弦值sin(A),cos(A),然后利用atan2计算角度值,自然就能规范到-pi到pi,即A=atan2(sin(A),cos(A))

33

帖子

12

主题

0

精华

熟手

F 币
138 元
贡献
93 点
5#
 楼主| 发表于 2020-9-12 10:39:37 | 只看该作者
风平老涡 发表于 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 = - ...

感谢,请看我的回答

835

帖子

2

主题

0

精华

大宗师

F 币
3926 元
贡献
2334 点
6#
发表于 2020-9-12 10:40:14 | 只看该作者
本帖最后由 li913 于 2020-9-12 10:41 编辑
mmSCI 发表于 2020-9-12 10:38
感谢以上各位铁友的回答,我想了一个简单的方法与大家共享,就是先求角度值得正弦和余弦值sin(A),cos(A), ...

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

33

帖子

12

主题

0

精华

熟手

F 币
138 元
贡献
93 点
7#
 楼主| 发表于 2020-9-12 10:40:23 | 只看该作者
li913 发表于 2020-9-11 21:27
b = modulo(a,2*PI); if(b>PI) b = b - 2*PI

感谢,请看我下面的回答
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-21 21:53

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表