 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1708 | 回復(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ù)此樓» 猜你喜歡
0703化學(xué)321分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有19人回復(fù)
-
298求調(diào)劑
已經(jīng)有5人回復(fù)
-
346求調(diào)劑 一志愿070303有機(jī)化學(xué)
已經(jīng)有6人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有3人回復(fù)
-
324求調(diào)劑
已經(jīng)有11人回復(fù)
-
307求調(diào)劑
已經(jīng)有13人回復(fù)
-
071010 323 分求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有10人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)321分求調(diào)劑 +3 | 三dd. 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +3 | lllllos 2026-03-30 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 研研，接電話 2026-03-24 | 8/400 |
|
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 6/300 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒(méi)煩惱 2026-03-23 | 4/200 |
|