 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1359 | 回復(fù): 7 | ||
| 【懸賞金幣】回答本帖問題，作者問天賣餅2將贈送您 100 個(gè)金幣 | ||
問天賣餅2捐助貴賓 (初入文壇)
|
[求助]
振蕩的波函數(shù)如何求積分 已有1人參與
|
|
|
我用薛定諤方程分別求出199個(gè)phi和199個(gè)psi，現(xiàn)在的問題是phi*psi后求積分求不出來，代碼如圖所示。我嘗試了很多種method都不行。求助廣大強(qiáng)大的網(wǎng)友，該怎么求積分呢？ @月只藍(lán) @月只藍(lán) @beefly 發(fā)自小木蟲Android客戶端 |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有68人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有6人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有3人回復(fù)
問天賣餅2
捐助貴賓 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 287
- 帖子: 10
- 在線: 2.9小時(shí)
- 蟲號: 24557690
- 注冊: 2020-11-12
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2020-11-12 10:09:41
問天賣餅2
捐助貴賓 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 287
- 帖子: 10
- 在線: 2.9小時(shí)
- 蟲號: 24557690
- 注冊: 2020-11-12
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
3樓2020-11-12 10:13:41
問天賣餅2
捐助貴賓 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 287
- 帖子: 10
- 在線: 2.9小時(shí)
- 蟲號: 24557690
- 注冊: 2020-11-12
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
4樓2020-11-12 10:17:43
問天賣餅2
捐助貴賓 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 287
- 帖子: 10
- 在線: 2.9小時(shí)
- 蟲號: 24557690
- 注冊: 2020-11-12
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
L = 1000*10^-10; F = 10^7; m = 0.0665*9.1*10^-31; mhole = 0.34*9.1*10^-31; e = 1.6*10^-19; \[HBar] = 1.05*10^-34; Ne = 199;eqn1 = (-(\[HBar]^2/(2 m))* (\[Phi]^\[Prime]\[Prime])[x] + e*F*x*\[Phi][x]);{phivals, phifuns} = NDEigensystem[{eqn1, DirichletCondition[\[Phi][x] == 0, x \[GreaterSlantEqual] L/2 || x <= -L/2]}, \[Phi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" -> {"FiniteElement", \ {"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}];eqn2 = (-(\[HBar]^2/(2 mhole))* (\[Psi]^\[Prime]\[Prime])[x] - e*F*x*\[Psi][x]);{psivals, psifuns} = NDEigensystem[{eqn2, DirichletCondition[[Psi][x] == 0, x [GreaterSlantEqual] L/2 || x <= -L/2]}, [Psi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" -> {"FiniteElement", \ {"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}];Etable = Table[(psivals[] + phivals[[j]])/(1.6*10^-22), {i, Ne}, {j, Ne}];Itable = Table[ Abs[NIntegrate[psifuns[]*phifuns[[j]], {x, -L/2, L/2}, Method -> {"ExtrapolatingOscillatory"}]], {i, Ne}, {j, Ne}]; |
5樓2020-11-12 11:06:59
問天賣餅2
捐助貴賓 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 287
- 帖子: 10
- 在線: 2.9小時(shí)
- 蟲號: 24557690
- 注冊: 2020-11-12
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
L = 1000*10^-10; F = 10^7; m = 0.0665*9.1*10^-31; mhole = 0.34*9.1*10^-31; e = 1.6*10^-19; \[HBar] = 1.05*10^-34; Ne = 199;eqn1 = (-(\[HBar]^2/(2 m))* (\[Phi]^\[Prime]\[Prime])[x] + e*F*x*\[Phi][x]);{phivals, phifuns} = NDEigensystem[{eqn1, DirichletCondition[\[Phi][x] == 0, x \[GreaterSlantEqual] L/2 || x <= -L/2]}, \[Phi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" -> {"FiniteElement", \ {"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}];eqn2 = (-(\[HBar]^2/(2 mhole))* (\[Psi]^\[Prime]\[Prime])[x] - e*F*x*\[Psi][x]);{psivals, psifuns} = NDEigensystem[{eqn2, DirichletCondition[[Psi][x] == 0, x [GreaterSlantEqual] L/2 || x <= -L/2]}, [Psi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" -> {"FiniteElement", \ {"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}];Etable = Table[(psivals + phivals[[j]])/(1.6*10^-22), {i, Ne}, {j, Ne}];Itable = Table[ Abs[NIntegrate[psifuns[]*phifuns[[j]], {x, -L/2, L/2}, Method -> {"ExtrapolatingOscillatory"}]], {i, Ne}, {j, Ne}]; |
6樓2020-11-12 11:47:37
問天賣餅2
捐助貴賓 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 287
- 帖子: 10
- 在線: 2.9小時(shí)
- 蟲號: 24557690
- 注冊: 2020-11-12
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
L = 1000*10^-10; F = 10^7; m = 0.0665*9.1*10^-31; mhole = 0.34*9.1*10^-31; e = 1.6*10^-19; [HBar] = 1.05*10^-34; Ne = 199; eqn1 = (-(\[HBar]^2/(2 m))* (\[Phi]^\[Prime]\[Prime])[x] + e*F*x*\[Phi][x]); eqn2 = (-(\[HBar]^2/(2 mhole))* (\[Psi]^\[Prime]\[Prime])[x] - e*F*x*\[Psi][x]) {phivals, phifuns} = NDEigensystem[{eqn1, DirichletCondition[\[Phi][x] == 0, x \[GreaterSlantEqual] L/2 || x <= -L/2]}, \[Phi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" ->{"FiniteElement", {"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}]; {psivals, psifuns} = NDEigensystem[{eqn2, DirichletCondition[\[Psi][x] == 0, x \[GreaterSlantEqual] L/2 || x <= -L/2]}, \[Psi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" -> "FiniteElement", \{"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}]; |
7樓2020-11-12 14:29:12
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時(shí)
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動力學(xué)
- 管轄: 計(jì)算模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
問天賣餅2: 金幣+100, ★★★很有幫助 2020-11-13 15:48:19
感謝參與，應(yīng)助指數(shù) +1
問天賣餅2: 金幣+100, ★★★很有幫助 2020-11-13 15:48:19
| 既然phi和psi都已經(jīng)計(jì)算出來（離散數(shù)據(jù)點(diǎn)），那么phi*psi也是確定的離散點(diǎn)，直接用梯形公式求解，當(dāng)然這種精度也較低�？梢圆捎脴訔l插值增加點(diǎn)數(shù)，然后再計(jì)算積分，這樣可以提高精度 |
8樓2020-11-13 14:34:11
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 問天賣餅2 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | danyyyy 2026-03-04 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 328求調(diào)劑 +5 | vuzhdkfjkx 2026-03-04 | 6/300 |
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +10 | pin8023 2026-02-28 | 13/650 |
|
|
|
[考研] 招調(diào)劑 +5 | 帆船哥 2026-03-04 | 7/350 |
|
|
|
[考研] 332材料求調(diào)劑 +4 | zjy101327 2026-03-05 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 曼殊2266 2026-02-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 12/600 |
|
|
|
[基金申請] 沒有青基直接申請面上，感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | 杰尼龜aaa 2026-03-03 | 4/200 |
|
|
|
[考研] 0856材料工程，初試313調(diào)劑 +7 | 賣個(gè)關(guān)子吧 2026-03-03 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +8 | llx0610 2026-03-02 | 9/450 |
|
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 322求調(diào)劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|