| 查看: 413 | 回复: 1 | ||
[求助]
在求无穷积分时，上限及步长的选取
|
|
虫友们，我算无穷积分时，发现积分上限和步长会影响结果，甚至出现NoN值，公式已在附件里，敬请多多指教！谢谢！如果能帮忙化简下，更好，运行太慢了 clear all lambda=1; tend=599.99;%积分上限 gamma=0.2; theta=0:0.001:pi/2; nn=1000;%步长决定积分的精确度 s= (tend-0)/nn; t=0:s:tend; m=length(t); d=sqrt(lambda.^2-2*lambda*gamma); for j=1:m k(j)=exp(-lambda*t(j)/2).*(cosh(d*t(j)/2)+lambda./d*sinh(d*t(j)/2)); kdot(j)=-0.5*lambda*k(j)+0.5*d*exp(-lambda*t(j)/2).*(sinh(d*t(j)/2)+lambda./d*cosh(d*t(j)/2)); p(j)=k(j)^(2); pdot(j)=2*k(j).*kdot(j); for q=1:length(theta) sigma1(j,q)=(cos(2*theta(q))^2*p(j)+0.5*sin(2*theta(q))^2)*pdot(j)./sqrt(cos(2*theta(q))^2*p(j)^2+sin(2*theta(q))^2*p(j)); sigma(j,q)=0.5*(sigma1(j,q)+abs(sigma1(j,q))); end end for q=1:length(theta) trap(q)=0; for i=1:nn-1 trap(q)=trap(q)+0.5*s*(sigma(i,q)+sigma(i+1,q)); end end [x,y]=max(trap)%积分后对theta求最大值  fig.png |
回复此楼» 猜你喜欢
实验室接单子
已经有7人回复
-
假如你的研究生提出不合理要求
已经有11人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 碳纳米管MD模拟中积分时间步长选取时存在的问题
已经有15人回复
2楼2015-05-02 21:08:36