 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1707 | 回復(fù): 3 | ||
[求助]
被積函數(shù)很快收斂到0，為什么積分不收斂。
|
|
函數(shù)f4（qz）在qz=1e3*kf時(shí)已經(jīng)為0，將f4在（1，Inf）積分，按說(shuō)，積分上限取到1e3*kf就可以了，事實(shí)上，積分上限在10*kf到1e3*kf之間時(shí)，積分結(jié)果是一致的，但是如果把上限放的非常大，（1,1e8*kf）這時(shí)候，積分結(jié)果變化又非常大。 后進(jìn)末學(xué)，求大賢指點(diǎn)一二，感激涕零，頓首再拜！ 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; %真空介電常數(shù) 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;%壓電系數(shù)間接引用 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);%費(fèi)米 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 |
回復(fù)此樓» 猜你喜歡
070300化學(xué)354求調(diào)劑
已經(jīng)有7人回復(fù)
-
286分調(diào)劑
已經(jīng)有6人回復(fù)
-
287求調(diào)劑
已經(jīng)有11人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿鄭大材料工程290求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有4人回復(fù)
-
化學(xué)308分調(diào)劑
已經(jīng)有8人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有16人回復(fù)
-
085701求調(diào)劑初試286分
已經(jīng)有5人回復(fù)
-
289求調(diào)劑
已經(jīng)有14人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【新人求助】構(gòu)建一PVAc分子鏈，用discover進(jìn)行minimize，怎么也不收斂
已經(jīng)有8人回復(fù)
- siesta-3.1 計(jì)算石墨烯條帶不收斂，不知是不是我的輸入有問(wèn)題。
已經(jīng)有24人回復(fù)
- Lyapunov函數(shù)的選取和控制系統(tǒng)收斂速度快慢是怎么一種關(guān)系?
已經(jīng)有10人回復(fù)
- 計(jì)算時(shí)候總是不能收斂如何解決？
已經(jīng)有5人回復(fù)
- 靜態(tài)計(jì)算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？
已經(jīng)有4人回復(fù)
- 電子步不收斂，跑了60多電子步第一個(gè)離子步還沒(méi)有完成，這是怎么回事？
已經(jīng)有28人回復(fù)
- cl-neb結(jié)果不收斂
已經(jīng)有19人回復(fù)
- 就解方程組，未知數(shù)不收斂，要怎么辦？誠(chéng)求！
已經(jīng)有12人回復(fù)
- 計(jì)算表面時(shí)：Vacuum真空層和slab層厚度收斂性測(cè)試，保持超晶胞體積不變嗎？？
已經(jīng)有10人回復(fù)
- 函數(shù)收斂于零的疑問(wèn)��！
已經(jīng)有3人回復(fù)
- siesta 表面不收斂
已經(jīng)有3人回復(fù)
- 【求助】K點(diǎn)取得越大越不收斂，什么原因！
已經(jīng)有20人回復(fù)
- 【求助】用Orca計(jì)算一個(gè)片簇的單點(diǎn)能，不收斂。
已經(jīng)有7人回復(fù)
- 【交流】為什么不收斂？
已經(jīng)有3人回復(fù)
- 【求助】計(jì)算不收斂怎么辦？
已經(jīng)有4人回復(fù)
- 【求助】ATK優(yōu)化計(jì)算不收斂怎么辦
已經(jīng)有15人回復(fù)
2樓2012-12-28 23:06:44
3樓2012-12-28 23:08:20
4樓2012-12-31 10:12:19
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 li_qun 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +11 | land xuxu 2026-03-26 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7 | nhj_ 2026-03-25 | 7/350 |
|
|
[考研] 283求調(diào)劑（080500） +14 | A child 2026-03-27 | 14/700 |
|
|
|
[考研] 085601材料工程找調(diào)劑 +16 | oatmealR 2026-03-29 | 17/850 |
|
|
|
[考研] 一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 5/250 |
|
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
|
[考研] 311（085601）求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +4 | 美味可樂(lè)雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛(ài)有機(jī) 2026-03-25 | 4/200 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|