搜搜做题作业网请大家帮我看看这个mathematica的积分问题
来源：学生作业帮 编辑：搜搜做题作业网作业帮 分类：综合作业 时间：2024/04/28 09:22:00
请大家帮我看看这个mathematica的积分问题
我想用mathematica计算一个积分,这个积分其实是从特征函数求得分布函数,但是用mathematica积分总是计算不出来,苦恼中……,程序如下,诚恳地请各位高手帮我看看问题之所在.
ClearAll[ sigma,psec,lambda,istar,temp,z,b,poutrayl];
sigma = 1;
psec = 1;
lambda = 0.1;
istar = 1*psec;
b = 2;
temp[w_,d_] =
Integrate[(1 - Exp[(I*w*psec*z)/r^(2*b)])*(z*Exp[-z^2/(2*sigma^2)])/
sigma^2*r,{z,0,Infinity},{r,d,Infinity},
GenerateConditions -> False]
poutrayl =
Plot[1 + I/(2*Pi)*
NIntegrate[
Exp[-2*lambda*\[Pi]*temp[w,d]]*(1 - Exp[-I*w*istar])/
w,{w,-Infinity,Infinity},MaxRecursion -> 20,
Exclusions -> w == 0],{d,0.01,5},
PlotStyle -> {Dashing[{0.001}],Black},
PlotRange -> {{0,5},{0,1}}]
我想用mathematica计算一个积分,这个积分其实是从特征函数求得分布函数,但是用mathematica积分总是计算不出来,苦恼中……,程序如下,诚恳地请各位高手帮我看看问题之所在.
ClearAll[ sigma,psec,lambda,istar,temp,z,b,poutrayl];
sigma = 1;
psec = 1;
lambda = 0.1;
istar = 1*psec;
b = 2;
temp[w_,d_] =
Integrate[(1 - Exp[(I*w*psec*z)/r^(2*b)])*(z*Exp[-z^2/(2*sigma^2)])/
sigma^2*r,{z,0,Infinity},{r,d,Infinity},
GenerateConditions -> False]
poutrayl =
Plot[1 + I/(2*Pi)*
NIntegrate[
Exp[-2*lambda*\[Pi]*temp[w,d]]*(1 - Exp[-I*w*istar])/
w,{w,-Infinity,Infinity},MaxRecursion -> 20,
Exclusions -> w == 0],{d,0.01,5},
PlotStyle -> {Dashing[{0.001}],Black},
PlotRange -> {{0,5},{0,1}}]
不要把mathematica看作计算的教父,列出来个式子就让它计算是不靠谱的.不仅工作在复数域而且有双重积分,基本没戏.
既然可以求出来 temp[w_,d_],建议你直接用切割求和的的方法求积分近似就不要使用NIntegrate了.
既然可以求出来 temp[w_,d_],建议你直接用切割求和的的方法求积分近似就不要使用NIntegrate了.