作者Gwaewluin (神无月 孝臣)
看板MATLAB
标题Re: [问题] Ei?
时间Thu Sep 14 13:10:25 2006
※ 引述《lwhsyh (KINKI)》之铭言：
: syms x real
: f=((sin(2*x)))/x)*exp(6*x);
: F=int(f,x,-10000,10000)
: F =
: 1/2*i*Ei(1,-60000-20000*i)-1/2*i*Ei(1,-60000+20000*i)
: -1/2*i*Ei(1,60000+20000*i)+1/2*i*Ei(1,60000-20000*i)
: -1/2*i*log(3+i)+1/2*i*log(3-i)+1/2*i*log(-3-i)-1/2*i*log(-3+i)
: 请问各位高手先进
: 我这积分出来的答案表示什麽
: Ei??
: 有无办法可以表示数字
: 感谢大家帮忙
这边的Ei是Symbolic Math Toolbox里面的mfun的东西
使用方法为mfun( 'Ei' , n , z )即可算出数字
其中n必须大於等於0
real(z)必须大於等於0
可是你的答案里面有-60000-20000*i这种东西
real(z)是小於0的
这样子跑下去mfun( 'Ei' , 1,-60000-20000*i )的答案是NaN - Infi
而就算real(z)是大於0的拿去跑mfun( 'Ei' , 1,60000-2000*i )答案也是0而已
而这个积分我用quad和quadl去处理也是跑不出来（答案是NaN）
使用高斯基分法也是跑出来NaN
（不过我没有帮你做变数变换转换区间，只有稍为算一下而已）
你这个可能真的没办法积分
Symbolic Toolbox常常会跑出来mfun的东西
要查询各mfun细节的话可以help mfunlist
--
Deserves death! I daresay he does. Many that live deserve death. And some die
that deserve life. Can you give that to them? Then be not too eager to deal out
death in the name of justice, fearing for your own safty. Even the wise cannot
see all ends.
Gandalf to Frodo
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.120.25.235
1F：推 lwhsyh:感谢帮忙~~我会再试试~谢谢 09/14 14:00