sun1993 发表于 2009-6-3 13:33

积分计算，不明白，请指教

被积函数 r*exp(-6*r/a0)*sin(k＊r) ，变量 r 从 0 到正无穷，a0,k都是常数
直接手算结果 12*a0^3*k/(36+a0^2＊k^2),
用int 命令，结果一大窜，且有 r 的一次项，根本不对
不知道问题出在哪，MATLAB 应该比我聪明的呀

xiezhh 发表于 2009-6-3 14:06

首先得说你打错了一个符号k*r 的*打错了，另外这个积分matlab判断不了是否收敛，出来的结果是一个极限
ans =
limit(-(a0^3*k^3*exp(-6*r/a0)*cos(k*r)*r-exp(-6*r/a0)*sin(k*r)*a0^3*k^2+6*exp(-6*r/a0)*sin(k*r)*r*k^2*a0^2+12*a0^2*k*exp(-6*r/a0)*cos(k*r)-12*k*a0^2+36*a0*k*exp(-6*r/a0)*cos(k*r)*r+36*exp(-6*r/a0)*sin(k*r)*a0+216*r*exp(-6*r/a0)*sin(k*r))*a0/(1296+72*k^2*a0^2+k^4*a0^4),r = Inf)

friendchj 发表于 2009-6-3 14:25

把limit化简后，应该和你的结果是一致的。

sun1993 发表于 2009-6-3 15:02

我以为它会给我最终结果呢

用limit化简不了结果，还是得手动算

[ 本帖最后由 ChaChing 于 2010-6-15 00:22 编辑 ]

无水1324 发表于 2009-6-3 21:21

回复 楼主 sun1993 的帖子

matlab有些东西就是不能处理，如果不是很麻烦还是受动算了！
页: [1]
查看完整版本: 积分计算，不明白，请指教