| 查看: 431 | 回复: 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 |
回复此楼» 猜你喜欢
材料学硕297已过四六级求调剂推荐
已经有9人回复
-
一志愿西南交通 专硕 材料355 本科双非 求调剂
已经有5人回复
-
316求调剂
已经有4人回复
-
一志愿武汉理工材料工程专硕调剂
已经有8人回复
-
329求调剂
已经有8人回复
-
350求调剂
已经有4人回复
-
材料与化工 322求调剂
已经有3人回复
-
一志愿武理材料工程348求调剂
已经有4人回复
-
265求调剂
已经有13人回复
-
290求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 碳纳米管MD模拟中积分时间步长选取时存在的问题
已经有15人回复
2楼2015-05-02 21:08:36
50