vpa( )问题

superhe · 发表于 2009-9-24 13:29

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

为什么vpa( )命令后计算出的数据都带一个虚数的符号i 啊？谢谢
vpa(a4)

a4 =
4634.5731964731538902675737918278+0.*i

ChaChing · 发表于 2009-9-24 15:51

应该不会有这现象!
LZ的a4是什么? 复数吗?

sogooda · 发表于 2009-9-24 16:16

可能是由于精度的问题造成的吧，如果楼主确认计算没有错误的话，可以把上面得到的数字直接取real就可以了。

superhe · 发表于 2009-9-24 17:35

原数据是这样的
125/8589934592*(766791424236068935830*pi^2+30223145490365730000*pi^2*sum(1/12*(-62062884+57857087401*pi^2*_R^2-358749530000*pi^4*_R^4+1249860937500*pi^6*_R^6)*log(283/100-_R)/pi^2/_R/(5859375000*pi^6*_R^6+12382812500*pi^4*_R^4+4628191875*pi^2*_R^2-93713644),_R = RootOf(8789062500*pi^8*_Z^8+24765625000*pi^6*_Z^6+13884575625*pi^4*_Z^4-562281864*pi^2*_Z^2+6895876))-30223145490365730000*pi^2*sum(1/12*(-62062884+57857087401*pi^2*_R^2-358749530000*pi^4*_R^4+1249860937500*pi^6*_R^6)*log(11/1000-_R)/pi^2/_R/(5859375000*pi^6*_R^6+12382812500*pi^4*_R^4+4628191875*pi^2*_R^2-93713644),_R = RootOf(8789062500*pi^8*_Z^8+24765625000*pi^6*_Z^6+13884575625*pi^4*_Z^4-562281864*pi^2*_Z^2+6895876)))^(1/2)
请高手指教啊！

friendchj · 发表于 2009-9-26 04:01

不知道_R和_Z代表什么意思，虚数应该是由RootOf引起的。

allenfieldin · 发表于 2009-9-26 09:26

我也遇到这样的问题，什么意思？

sogooda · 发表于 2009-9-26 10:14

这样看来，就是直接用
real(vpa(...))就可以了。
如果本来不应该有虚数的话，那其中的0*i应该是由于计算精度误差引入的，直接无视就可以了。

ChaChing · 发表于 2009-9-26 12:25

AM4:01这麽晚, 人在国外!?
熬夜基本上对身体不太好!

friendchj · 发表于 2009-9-26 12:37

谢谢关心，呵呵
在国外呆着，时差13小时，呵呵

账号		自动登录	找回密码
密码			我要加入

[综合讨论] vpa( )问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复楼主 superhe 的帖子

回复地板 superhe 的帖子

评分

回复地板 superhe 的帖子

评分

回复 5楼 friendchj 的帖子

回复 8楼 ChaChing 的帖子

[综合讨论] vpa( )问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 楼主 superhe 的帖子

回复 地板 superhe 的帖子

评分

回复 地板 superhe 的帖子

评分

回复 5楼 friendchj 的帖子

回复 8楼 ChaChing 的帖子

回复楼主 superhe 的帖子

回复地板 superhe 的帖子

回复地板 superhe 的帖子