| 查看: 1929 | 回复: 3 | ||
[求助]
被积函数很快收敛到0，为什么积分不收敛。
|
|
函数f4（qz）在qz=1e3*kf时已经为0，将f4在（1，Inf）积分，按说，积分上限取到1e3*kf就可以了，事实上，积分上限在10*kf到1e3*kf之间时，积分结果是一致的，但是如果把上限放的非常大，（1,1e8*kf）这时候，积分结果变化又非常大。 后进末学，求大贤指点一二，感激涕零，顿首再拜！ global EF A D h2 x Epsilon0 Epsilons Epsilonh e m m2 Ns kb T B DP ul ut rho E0 kf K2 KL KT; kb=1.38e-23; T=3; %温度 h2=1.055e-34; x=0.3; %Mg% Epsilon0=8.8542e-12; %真空介电常数 Epsilons=7.9; % 间接 Epsilonh=3.7; % a1=(3.250+0.05*x)*1e-10; c1=(5.204-0.17*x)*1e-10; a2=3.250e-10;c2=5.204e-10; % e=1.6022e-19; %电子电荷 DP=3.8*e; rho=5760; ul=5e3; % ut=2.1e3; % cl=14.1e10; % ct=2.47e10; % h15=-0.37;%压电系数间接引用 h31=-0.62; h33=0.96; ela2=(8*(2*h15+h31)^2+12*(2*h15+h31)*h33+15*h33^2)/105; % eta2=(6*(h33-h15-h31)^2+16*(h33-h15-h31)*h15+48*h15^2)/105; K2=ela2/(Epsilon0*Epsilons*cl)+eta2/(Epsilon0*Epsilons*ct); % KL=ela2/(Epsilon0*Epsilons*cl); KT=eta2/(Epsilon0*Epsilons*ct); m0=9.11e-31; m=0.26*m0; %自由电子 Vc=0.9*2.145*x*e;% sigmat=-0.054-0.038*x; sigmai=0.038*x;sigmab=0.054; deltaz=0.1e-9;%变分间隔 sigmad=1e17/1e-6*e*deltaz;%一 Ns=1.4e15; B=(12*m*e^2/(Epsilon0*Epsilons*h2^2))^(1/3)*(11*Ns/32)^(1/3); kf=sqrt(2*pi*Ns);%费米 E0=h2^2*B^2/(8*m)+3*e^2/(Epsilon0*Epsilons*B)*(11*Ns/16);% EF=E0+pi*h2^2*Ns/m;% f4=@(qz)B.^6./(B.^2+qz.^2).^3.*DP.^2.*h2.*sqrt(kf.^2+qz.^2)./(2.*rho.*ul).*(1./(exp(h2.*ul.*sqrt(kf.^2+qz.^2)./(kb.*T))-1)).*(1-1./(exp((h2.^2.*kf.^2./(2.*m)+E0+h2.*ul.*sqrt(kf.^2+qz.^2)-EF)./(kb.*T))+1)); >> quadgk(f4,1,10*kf) ans = 2.8824e-063 >> quadgk(f4,1,100*kf) ans = 2.8824e-063 >> quadgk(f4,1,1e8*kf) ans = 4.0394e-094 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【新人求助】构建一PVAc分子链，用discover进行minimize，怎么也不收敛
已经有8人回复
- siesta-3.1 计算石墨烯条带不收敛，不知是不是我的输入有问题。
已经有24人回复
- Lyapunov函数的选取和控制系统收敛速度快慢是怎么一种关系?
已经有10人回复
- 计算时候总是不能收敛如何解决？
已经有5人回复
- 静态计算能收敛，同样的参数弛豫不能收敛，这是为什么？
已经有4人回复
- 电子步不收敛，跑了60多电子步第一个离子步还没有完成，这是怎么回事？
已经有28人回复
- cl-neb结果不收敛
已经有19人回复
- 就解方程组，未知数不收敛，要怎么办？诚求！
已经有12人回复
- 计算表面时：Vacuum真空层和slab层厚度收敛性测试，保持超晶胞体积不变吗？？
已经有10人回复
- 函数收敛于零的疑问！！
已经有3人回复
- siesta 表面不收敛
已经有3人回复
- 【求助】K点取得越大越不收敛，什么原因！
已经有20人回复
- 【求助】用Orca计算一个片簇的单点能，不收敛。
已经有7人回复
- 【交流】为什么不收敛？
已经有3人回复
- 【求助】计算不收敛怎么办？
已经有4人回复
- 【求助】ATK优化计算不收敛怎么办
已经有15人回复
2楼2012-12-28 23:06:44
3楼2012-12-28 23:08:20
4楼2012-12-31 10:12:19
10