| 查看: 430 | 回复: 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 |
回复此楼» 猜你喜欢
08工学调剂
已经有5人回复
-
281求调剂（0805）
已经有25人回复
-
085600材料与化工
已经有6人回复
-
265求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有21人回复
-
295材料求调剂，一志愿武汉理工085601专硕
已经有4人回复
-
0856调剂，是学校就去
已经有9人回复
-
298-一志愿中国农业大学-求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
学校已经提交到NSFC，还能修改吗？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 碳纳米管MD模拟中积分时间步长选取时存在的问题
已经有15人回复
2楼2015-05-02 21:08:36
