| 查看: 493 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 碳纳米管MD模拟中积分时间步长选取时存在的问题
已经有15人回复
2楼2015-05-02 21:08:36
20