 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1667 | 回復(fù): 2 | ||
cumtsyf01鐵蟲 (初入文壇)
|
[求助]
請(qǐng)教高手：PDE工具箱能解兩個(gè)偏微分方程相互迭代求解嗎
|
|
請(qǐng)教各位高手：PDE工具箱能解兩個(gè)偏微分方程相互迭代求解嗎？請(qǐng)求圓域上二維薛定諤方程泊松方程迭代自洽的MATLAB解法 眼看越來(lái)越快臨近畢業(yè)了，老板催交論文時(shí)間很急！只是苦于遇到個(gè)技術(shù)性障礙一直沒(méi)有解決，故在此懇請(qǐng)各位高手不吝賜教！小弟感激不盡�。�！請(qǐng)教高手： 問(wèn)題【1】—— PDE工具箱能解這兩個(gè)偏微分方程相互迭代求解嗎？(其中一個(gè)是本征值方程）能的話怎么實(shí)現(xiàn)？求圓域上二維薛定諤方程泊松方程迭代自洽的MATLAB解法， 這兩個(gè)偏微分方程如下： ① 倒三角laplace算子*f(x,y)+V(x,y)*f(x,y)=E*f(x,y)E為本征值，這是薛定諤方程 ② 倒三角laplace算子*V(x,y)=a+B(x,y)+c(E)*|f(x,y)|^2這是泊松方程 以上兩個(gè)方程相互迭代。由①求出f(x,y)和E來(lái)代入②中。再由②求出V(x,y)代入①中，再由①重新求出f(x,y)和E來(lái)代入②中。。。。如此相互迭代直到相鄰兩次迭代的數(shù)值相差足夠小則認(rèn)為收斂。 求解區(qū)域是個(gè)圓域邊界條件是圓周上f(x,y)和V(x,y)為0.其中a是常數(shù)，函數(shù)B(x,y)c(E)形式已知。倒三角laplace算子 是 偏比偏x平方+偏比偏y平方 現(xiàn)在遇到的最大的技術(shù)障礙是： 引申到 問(wèn)題【2】—— 末學(xué)發(fā)現(xiàn)PDE工具箱雖然很強(qiáng)大可視化也很好，但是不是只能解已知了“解析表達(dá)式”的系數(shù)的方程？換句話說(shuō)只能求解一次已知系數(shù)的，不能連續(xù)解多次由上一次得出系數(shù)的？只能解單個(gè)偏微分方程，不能解兩個(gè)相互迭代的偏微分方程？但這里由于要相互迭代，上一步PDE解出來(lái)的不可能是“解析”表達(dá)式的V(x,y),而是一堆離散的xi，yi數(shù)據(jù)點(diǎn)上的Vi“數(shù)值”的值，若這些離散數(shù)值能輸出，請(qǐng)問(wèn)這“數(shù)值”的系數(shù)【xi，yi，Vi】能代入作為下一個(gè)偏微分方程的系數(shù)用PDEtool來(lái)解嗎？能的話怎么實(shí)現(xiàn)？ 再引申到 問(wèn)題【3】—— 退一步想，如果能想辦法找到二維曲面擬合函數(shù)把這些離散的xi，yi數(shù)據(jù)點(diǎn)上的Vi“數(shù)值”擬合出一個(gè)V(x,y)的表達(dá)式，但由于這個(gè)表達(dá)式預(yù)先并不知道，沒(méi)法像這樣一次性手動(dòng)添加事先已知的表達(dá)式放在單引號(hào)里： ‘x.^2+y.^2’ 用這樣的pdetool代碼來(lái)手動(dòng)添加系數(shù)‘x.^2+y.^2’： %PDEcoefficients: pdeseteq(4,... '2.0',... 'x.^2+y.^2',...%此處就是對(duì)應(yīng)的系數(shù)，但問(wèn)題是這個(gè)表達(dá)式事先不知道啊，怎么自動(dòng)添加到pdeseteq里？ '10.0',... '1.0',... '0:10',... '0.0',... '0.0',... '[0100]') setappdata(pde_fig,'currparam',... ['2.0';... 'x.^2+y.^2';...%此處也是對(duì)應(yīng)的系數(shù)，但問(wèn)題是這個(gè)表達(dá)式事先不知道啊，怎么自動(dòng)添加到setappdata里？ '10.0';... '1.0']) 請(qǐng)問(wèn)假如能擬合成功，能把此表達(dá)式自動(dòng)讓它生成函數(shù)表達(dá)式句柄，而自動(dòng)添加到PDEtool的系數(shù)里【pdeseteq和setappdata】嗎？ 問(wèn)題【1】【2】【3】中的任何一個(gè)技術(shù)障礙解決掉，整個(gè)問(wèn)題都能解決掉 另外（最好用MATLAB的PDE函數(shù)解決，而不是差分法，因?yàn)橐痪S的薛定諤方程泊松方程自洽用差分法的我是算過(guò)的，但二維極坐標(biāo)下的用差分法各種細(xì)節(jié)實(shí)在太過(guò)繁瑣麻煩也容易出錯(cuò)。。。。所以目前寄希望于MATLAB的PDE函數(shù)。而PDE工具箱用指令實(shí)現(xiàn)倒是大有希望，只可惜遇到以上的技術(shù)障礙，眼看畢業(yè)越來(lái)越近了，老板也催論文催的挺緊，可是在這個(gè)問(wèn)題一直【困擾】耽擱了末學(xué)很久很久！遲遲得不到解決！懇請(qǐng)各位高手不吝賜教！末學(xué)在此感激不盡！！！末學(xué)郵箱cumtsyf03@163.com |
回復(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é)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有235人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有29人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有53人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 新人請(qǐng)教兩個(gè)問(wèn)題：1.聯(lián)合培養(yǎng)2年是否可行？2。。。
已經(jīng)有3人回復(fù)
- 2012版申請(qǐng)書二級(jí)標(biāo)題內(nèi)容的格式問(wèn)題請(qǐng)教
已經(jīng)有7人回復(fù)
- 請(qǐng)教晶體材料的禁帶寬度和激發(fā)波長(zhǎng)的關(guān)系！
已經(jīng)有4人回復(fù)
- 求解譜，有質(zhì)譜、氫譜、碳譜，多謝！
已經(jīng)有14人回復(fù)
- 請(qǐng)教origin三維作圖問(wèn)題
已經(jīng)有4人回復(fù)
- 請(qǐng)教有什么方法可以確定高分子聚合物的粘結(jié)性能呢？
已經(jīng)有15人回復(fù)
- 【求助】請(qǐng)教牛人--索氏抽提粉末樣怎么包？？？
已經(jīng)有9人回復(fù)
- 【求助】請(qǐng)教一個(gè)基層工作最低年限的問(wèn)題
已經(jīng)有5人回復(fù)
- 請(qǐng)教“2區(qū)文章”
已經(jīng)有29人回復(fù)
- 求張文生《科學(xué)計(jì)算中的偏微分方程有限差分法》和胡建偉，湯懷民《微分方程數(shù)值解法》
已經(jīng)有3人回復(fù)
- 【求助】偏微分方程解法
已經(jīng)有19人回復(fù)
q68
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 3 (幼兒園)
- 貴賓: 1.615
- 金幣: 14714.2
- 散金: 10
- 紅花: 3
- 帖子: 2309
- 在線: 149.8小時(shí)
- 蟲號(hào): 325340
- 注冊(cè): 2007-03-17
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
2樓2011-07-14 15:07:21
cumtsyf01
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 22.5
- 帖子: 42
- 在線: 31.3小時(shí)
- 蟲號(hào): 304570
- 注冊(cè): 2006-12-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2011-07-18 14:11:01
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cumtsyf01 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +9 | 秋有木北 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|