 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1709 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
被積函數(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ù)此樓» 猜你喜歡
11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè)
已經(jīng)有6人回復(fù)
-
求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有8人回復(fù)
-
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料工程專碩求調(diào)劑
已經(jīng)有8人回復(fù)
-
297求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有7人回復(fù)
-
281求調(diào)劑
已經(jīng)有7人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有8人回復(fù)
-
求收留
已經(jīng)有6人回復(fù)
-
085600 295分求調(diào)劑
已經(jīng)有9人回復(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ù)
4樓2012-12-31 10:12:19
2樓2012-12-28 23:06:44
3樓2012-12-28 23:08:20
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩333求調(diào)劑 +14 | 北道巷 2026-03-24 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|