| 4 | 1/1 | 返回列表 |
| 查看: 1105 | 回復(fù): 3 | |||
wbjhpp金蟲 (初入文壇)
|
[交流]
【求助】用mathematic編個(gè)二分法的程序 已有2人參與
|
| 請(qǐng)用mathematic編個(gè)二分法的程序,謝謝嘍!如果方程類似反函數(shù)的話,應(yīng)該怎么排除這種情況? |
至尊木蟲 (正式寫手)
金蟲 (初入文壇)
|
h[x]是薛定諤方程的波函數(shù)(不知道具體形式),x表示能量,下面是我編的找出h[x]==0的解的二分法程序,v1、v2是已知的勢(shì)能;但如果此波函數(shù)是類似反函數(shù)的情況,也滿足h[p1]*h[p3]<0,但它與x軸沒交點(diǎn),也就是說并沒有解,那怎么用程序排除這種情況!謝謝 total = 1000; y1 = v1;y2 = v2; y0 = (y2 - y1)/total; rt[x_, y_] := Module[{p1 = x, p2 = y, p3 = (x + y)/2}, While[Abs[p1 - p2] > 0.0001, If[h[p1]*h[p2] == 0, If[h[p1] == 0, p3 = p1, p3 = p2]; Break[]]; If[Negative[h[p1]*h[p3]], p2 = p3; p3 = (p1 + p2)/2, p1 = p3; p3 = (p1 + p2)/2]]; p3]; last = Table[ If[h[y1 + w*y0]*h[y1 + (w + 1) y0] < 0, rt[y1 + w*y0, y1 + (w + 1) y0], 0], {w, 1, total - 2}]; sz = Module[{i = 0}, Do[If[Part[last, j] != 0, i = i + 1, i], {j, 1, Length[last]}]; i]; netlast = Table[0, {i, 1, sz}]; Module[{j = 1}, Do[If[Part[last, i] != 0, Part[netlast, j] = Part[last, i]; j = j + 1], {i, 1, Length[last]}]]; netlast |
金蟲 (初入文壇)
|
h[x]是薛定諤方程的波函數(shù)(不知道具體形式),x表示能量,下面是我編的找出h[x]==0的解的二分法程序,v1、v2是已知的勢(shì)能;但如果此波函數(shù)是類似反函數(shù)的情況,也滿足h[p1]*h[p3]<0,但它與x軸沒交點(diǎn),也就是說并沒有解,那怎么用程序排除這種情況!謝謝 total = 1000; y1 = v1;y2 = v2; y0 = (y2 - y1)/total; rt[x_, y_] := Module[{p1 = x, p2 = y, p3 = (x + y)/2}, While[Abs[p1 - p2] > 0.0001, If[h[p1]*h[p2] == 0, If[h[p1] == 0, p3 = p1, p3 = p2]; Break[]]; If[Negative[h[p1]*h[p3]], p2 = p3; p3 = (p1 + p2)/2, p1 = p3; p3 = (p1 + p2)/2]]; p3]; last = Table[ If[h[y1 + w*y0]*h[y1 + (w + 1) y0] < 0, rt[y1 + w*y0, y1 + (w + 1) y0], 0], {w, 1, total - 2}]; sz = Module[{i = 0}, Do[If[Part[last, j] != 0, i = i + 1, i], {j, 1, Length[last]}]; i]; netlast = Table[0, {i, 1, sz}]; Module[{j = 1}, Do[If[Part[last, i] != 0, Part[netlast, j] = Part[last, i]; j = j + 1], {i, 1, Length[last]}]]; netlast |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +13 | 稻妻小編 2026-03-26 | 13/650 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 0703化學(xué) +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 【雙一流院校新能源、環(huán)境材料,材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|