版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者問天賣餅2將贈送您 100 個金幣

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時
蟲號: 24557690
注冊: 2020-11-12
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] 振蕩的波函數(shù)如何求積分已有1人參與

我用薛定諤方程分別求出199個phi和199個psi，現(xiàn)在的問題是phi*psi后求積分求不出來，代碼如圖所示。我嘗試了很多種method都不行。求助廣大強(qiáng)大的網(wǎng)友，該怎么求積分呢？

@月只藍(lán) @月只藍(lán) @beefly 發(fā)自小木蟲Android客戶端

回復(fù)此樓

» 猜你喜歡

德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生已經(jīng)有0人回復(fù)
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試已經(jīng)有11人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有65人回復(fù)
2026-CJ開始申報了已經(jīng)有1人回復(fù)
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士已經(jīng)有13人回復(fù)
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有1人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有2人回復(fù)

1樓 2020-11-12 10:03:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時
蟲號: 24557690
注冊: 2020-11-12
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

2樓2020-11-12 10:09:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時
蟲號: 24557690
注冊: 2020-11-12
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

1D的高度振蕩的波函數(shù)的積分

贊一下

回復(fù)此樓

3樓2020-11-12 10:13:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時
蟲號: 24557690
注冊: 2020-11-12
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

mathematica

發(fā)自小木蟲Android客戶端

回復(fù)此樓

4樓2020-11-12 10:17:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時
蟲號: 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}];

贊一下

回復(fù)此樓

5樓2020-11-12 11:06:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時
蟲號: 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}];

贊一下

回復(fù)此樓

6樓2020-11-12 11:47:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時
蟲號: 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}}}}];

贊一下

回復(fù)此樓

7樓2020-11-12 14:29:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機(jī)械動力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(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ì)算積分，這樣可以提高精度

贊一下(1人)

回復(fù)此樓

數(shù)值計(jì)算

8樓2020-11-13 14:34:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主問天賣餅2 的主題更新

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎勵，但不允許灌水，必須填寫15個字符以上)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[論文投稿] EST拒稿重投 5+3	15102603076 2026-03-02	3/150	2026-03-04 00:51 by bobvan
[考研] 化學(xué) 0703求調(diào)劑總分293 一志愿211 +3	土土小蟲 2026-03-03	3/150	2026-03-03 23:17 by zhukairuo
[考研] 化工專碩調(diào)劑 +4	利好利好. 2026-03-03	7/350	2026-03-03 21:30 by L135790
[考研] 理學(xué)，工學(xué)，農(nóng)學(xué)調(diào)劑，少走彎路，這里歡迎您！ +8	likeihood 2026-03-02	11/550	2026-03-03 19:39 by 一葉凡塵?
[考研] 290求調(diào)劑 +9	ErMiao1020 2026-03-02	9/450	2026-03-03 18:03 by linlonghao
[考研] 0805總分292，求調(diào)劑 +12	幻想之殤 2026-03-01	12/600	2026-03-03 15:58 by tgxtgxtgx9
[考研] 材料工程求調(diào)劑 +3	1431251 2026-03-03	3/150	2026-03-03 11:58 by EBSD
[考研] 考研復(fù)試調(diào)劑，過國家線的同學(xué)都可報名 +7	黑！在干嘛 2026-02-28	8/400	2026-03-03 08:39 by ww哈
[考研] 288求調(diào)劑 +3	少71.8 2026-03-02	5/250	2026-03-03 06:01 by tgxtgxtgx9
[考研] 【2026 碩士調(diào)劑】課題組招收調(diào)劑生 +3	考研版棒棒 2026-03-02	5/250	2026-03-03 01:45 by kkky.
[考研] 11408，學(xué)碩276求調(diào)劑 +3	崔wj 2026-03-02	5/250	2026-03-02 22:14 by 255671
[考研] 261求調(diào)劑 +3	陸lh 2026-03-01	3/150	2026-03-02 19:32 by zhukairuo
[考博] 博士自薦 +4	kkluvs 2026-02-28	5/250	2026-03-02 19:19 by 輕松不少隨
[考研] 材料085601調(diào)劑 +5	多多子. 2026-03-02	5/250	2026-03-02 19:15 by zhukairuo
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3	shs1083 2026-03-02	3/150	2026-03-02 17:27 by houyaoxu
[考研] 化工專碩342，一志愿大連理工大學(xué)，求調(diào)劑 +6	kyf化工 2026-02-28	7/350	2026-03-02 10:56 by 無際的草原
[基金申請] 成果系統(tǒng)訪問量大，請一小時后再嘗試。---NSFC啥時候好哦，已經(jīng)兩天這樣了 +4	NSFC2026我來了 2026-02-28	4/200	2026-03-01 22:37 by 鐵門栓
[考研] 272求調(diào)劑 +6	田智友 2026-02-28	6/300	2026-03-01 21:40 by 公瑾逍遙
[考研] 313求調(diào)劑 +3	水流年lc 2026-02-28	3/150	2026-03-01 16:01 by 新能源達(dá)人
[考研] 307求調(diào)劑 +4	73372112 2026-02-28	6/300	2026-03-01 00:04 by ll247

24小時熱門版塊排行榜

[求助] 振蕩的波函數(shù)如何求積分 已有1人參與

» 猜你喜歡

【答案】應(yīng)助回帖

[求助] 振蕩的波函數(shù)如何求積分已有1人參與