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

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

24小時(shí)熱門版塊排行榜

返回列表

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

問天賣餅2

捐助貴賓 (初入文壇)

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

[求助] 振蕩的波函數(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ù)此樓

» 猜你喜歡

幫我的英語(yǔ)口語(yǔ)老師找學(xué)生已經(jīng)有4人回復(fù)
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有54人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有1人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有2人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)

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

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

問天賣餅2

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 287
帖子: 10
在線: 2.9小時(shí)
蟲號(hào): 24557690
注冊(cè): 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小時(shí)
蟲號(hào): 24557690
注冊(cè): 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小時(shí)
蟲號(hào): 24557690
注冊(cè): 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小時(shí)
蟲號(hào): 24557690
注冊(cè): 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小時(shí)
蟲號(hào): 24557690
注冊(cè): 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小時(shí)
蟲號(hào): 24557690
注冊(cè): 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小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(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)助才可能被獎(jiǎng)勵(lì)，但不允許灌水，必須填寫15個(gè)字符以上)

普通表情龍兔虎貓

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

[考研] 一志愿華科071000生物學(xué) 338求調(diào)劑 +4	九月九里 2026-03-05	4/200	2026-03-11 20:22 by 求調(diào)劑zz
[考研] 0856化工原理 +5	z2839474511 2026-03-10	5/250	2026-03-11 19:46 by peike
[考研] 289求調(diào)劑 +4	步川酷紫123 2026-03-11	4/200	2026-03-11 19:41 by peike
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪？ +5	MayUxw1 2026-03-08	7/350	2026-03-11 18:09 by 17862566385
[考研] 標(biāo)題:撿漏預(yù)警\|08工科/09農(nóng)學(xué)調(diào)劑!英語(yǔ)要求低，過線即有機(jī)會(huì)! +7	馬超放煙花 2026-03-07	12/600	2026-03-10 23:16 by Equinoxhua
[考研] 求調(diào)劑！ +3	朔朔話 2026-03-09	3/150	2026-03-10 15:48 by houyaoxu
[考研] 337一志愿華南理工材料求調(diào)劑 +5	mysdl 2026-03-07	5/250	2026-03-10 14:54 by ruiyingmiao
[考研] 085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211 +8	等春來， 2026-03-06	9/450	2026-03-10 10:38 by 沙漠之狐994
[考研] 070300化學(xué)求調(diào)劑 +5	撲風(fēng)鈴的貓 2026-03-08	10/500	2026-03-09 13:58 by macy2011
[考研] 0703化學(xué)求調(diào)劑 +4	很老實(shí)人 2026-03-09	4/200	2026-03-09 13:38 by macy2011
[考研] 求調(diào)劑，數(shù)一英一274分 +4	小菲會(huì)努力 2026-03-08	4/200	2026-03-09 12:40 by 一定上岸喲_
[考研] 337求調(diào)劑 +3	睡醒，。 2026-03-09	3/150	2026-03-09 10:02 by 求調(diào)劑zz
[考研] 考研求調(diào)劑，工科，總分304 +4	矮子不想長(zhǎng)高了 2026-03-08	7/350	2026-03-08 22:21 by 小木蟲tim
[考研] 理學(xué)287分求調(diào)劑 +4	看看我. 2026-03-05	4/200	2026-03-08 22:11 by qingfeng258
[考研] 調(diào)劑 +3	13853210211 2026-03-08	3/150	2026-03-08 19:34 by ms629
[考研] 303求調(diào)劑 +8	forgman95 2026-03-05	10/500	2026-03-08 12:41 by 藍(lán)莓都給你
[考研] 化學(xué)工程與技術(shù)316求調(diào)劑 +7	Sixuan wang 2026-03-06	7/350	2026-03-08 11:23 by 我的船我的海
[考研] 材料化工求調(diào)劑 +7	博斯特525 2026-03-06	8/400	2026-03-08 03:51 by lfhuang
[考研] 材料與化工354調(diào)劑 +4	Lucy-xiao 2026-03-06	7/350	2026-03-06 09:57 by hxuq8429
[考研] 紡織、生物、化學(xué)、材料等專業(yè) +3	Eember. 2026-03-05	7/350	2026-03-05 16:35 by Eember.

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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

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