| 查看: 1863 | 回复: 2 | ||
cumtsyf01铁虫 (初入文坛)
|
[求助]
请教高手：PDE工具箱能解两个偏微分方程相互迭代求解吗
|
|
请教各位高手：PDE工具箱能解两个偏微分方程相互迭代求解吗？请求圆域上二维薛定谔方程泊松方程迭代自洽的MATLAB解法 眼看越来越快临近毕业了，老板催交论文时间很急！只是苦于遇到个技术性障碍一直没有解决，故在此恳请各位高手不吝赐教！小弟感激不尽！！！请教高手： 问题【1】—— PDE工具箱能解这两个偏微分方程相互迭代求解吗？(其中一个是本征值方程）能的话怎么实现？求圆域上二维薛定谔方程泊松方程迭代自洽的MATLAB解法， 这两个偏微分方程如下： ① 倒三角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这是泊松方程 以上两个方程相互迭代。由①求出f(x,y)和E来代入②中。再由②求出V(x,y)代入①中，再由①重新求出f(x,y)和E来代入②中。。。。如此相互迭代直到相邻两次迭代的数值相差足够小则认为收敛。 求解区域是个圆域边界条件是圆周上f(x,y)和V(x,y)为0.其中a是常数，函数B(x,y)c(E)形式已知。倒三角laplace算子 是 偏比偏x平方+偏比偏y平方 现在遇到的最大的技术障碍是： 引申到 问题【2】—— 末学发现PDE工具箱虽然很强大可视化也很好，但是不是只能解已知了“解析表达式”的系数的方程？换句话说只能求解一次已知系数的，不能连续解多次由上一次得出系数的？只能解单个偏微分方程，不能解两个相互迭代的偏微分方程？但这里由于要相互迭代，上一步PDE解出来的不可能是“解析”表达式的V(x,y),而是一堆离散的xi，yi数据点上的Vi“数值”的值，若这些离散数值能输出，请问这“数值”的系数【xi，yi，Vi】能代入作为下一个偏微分方程的系数用PDEtool来解吗？能的话怎么实现？ 再引申到 问题【3】—— 退一步想，如果能想办法找到二维曲面拟合函数把这些离散的xi，yi数据点上的Vi“数值”拟合出一个V(x,y)的表达式，但由于这个表达式预先并不知道，没法像这样一次性手动添加事先已知的表达式放在单引号里： ‘x.^2+y.^2’ 用这样的pdetool代码来手动添加系数‘x.^2+y.^2’： %PDEcoefficients: pdeseteq(4,... '2.0',... 'x.^2+y.^2',...%此处就是对应的系数，但问题是这个表达式事先不知道啊，怎么自动添加到pdeseteq里？ '10.0',... '1.0',... '0:10',... '0.0',... '0.0',... '[0100]') setappdata(pde_fig,'currparam',... ['2.0';... 'x.^2+y.^2';...%此处也是对应的系数，但问题是这个表达式事先不知道啊，怎么自动添加到setappdata里？ '10.0';... '1.0']) 请问假如能拟合成功，能把此表达式自动让它生成函数表达式句柄，而自动添加到PDEtool的系数里【pdeseteq和setappdata】吗？ 问题【1】【2】【3】中的任何一个技术障碍解决掉，整个问题都能解决掉 另外（最好用MATLAB的PDE函数解决，而不是差分法，因为一维的薛定谔方程泊松方程自洽用差分法的我是算过的，但二维极坐标下的用差分法各种细节实在太过繁琐麻烦也容易出错。。。。所以目前寄希望于MATLAB的PDE函数。而PDE工具箱用指令实现倒是大有希望，只可惜遇到以上的技术障碍，眼看毕业越来越近了，老板也催论文催的挺紧，可是在这个问题一直【困扰】耽搁了末学很久很久！迟迟得不到解决！恳请各位高手不吝赐教！末学在此感激不尽！！！末学邮箱cumtsyf03@163.com |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有287人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有6人回复
-
散金币祈福
已经有33人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 新人请教两个问题：1.联合培养2年是否可行？2。。。
已经有3人回复
- 2012版申请书二级标题内容的格式问题请教
已经有7人回复
- 请教晶体材料的禁带宽度和激发波长的关系！
已经有4人回复
- 求解谱，有质谱、氢谱、碳谱，多谢！
已经有14人回复
- 请教origin三维作图问题
已经有4人回复
- 请教有什么方法可以确定高分子聚合物的粘结性能呢？
已经有15人回复
- 【求助】请教牛人--索氏抽提粉末样怎么包？？？
已经有9人回复
- 【求助】请教一个基层工作最低年限的问题
已经有5人回复
- 请教“2区文章”
已经有29人回复
- 求张文生《科学计算中的偏微分方程有限差分法》和胡建伟，汤怀民《微分方程数值解法》
已经有3人回复
- 【求助】偏微分方程解法
已经有19人回复
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 15672.2
- 散金: 10
- 红花: 3
- 帖子: 2457
- 在线: 149.9小时
- 虫号: 325340
- 注册: 2007-03-17
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2011-07-14 15:07:21
cumtsyf01
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 22.5
- 帖子: 42
- 在线: 31.3小时
- 虫号: 304570
- 注册: 2006-12-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-07-18 14:11:01
