多语言展示
当前在线:1911今日阅读:145今日分享:43

如何求解偏微分方程 -- 椭圆方程实例:

求解偏微分方程 -- 椭圆方程实例:求解在域Ω上的泊松方程-ΔU=1,边界条件Ω上U=0的数值解,其中Ω是一个单位圆。
方法/步骤
1

步骤1:启动pdetool界面。在matlab命令窗口中输入pdetool,然后画一个单位圆,并单击图标,如下图所示:

2

步骤2:在Boundary下选择Specify boundary conditions选项,将边界条件选中为Dirichlet条件并设置h=1,r=0,如图2所示:

3

步骤3:单击PDE按钮,会弹出图3的对话框,选中椭圆按钮并设置c=1,a=0,f=1.

4

步骤4:划分单元。单价三角按钮,弹出如图4所示的图像。继续单击双三角按钮,弹出如图5所示的图像。

5

步骤5 求解方程。单击等号按钮,显示求解方程值的分布。

6

步骤6:对比精确解的绝对误差值。选择Plot下的Paremeters选项,弹出选择框,如图7所示,在Preperty下选择user entry选项,并在其中输入方程的精确解u-(1-x.^2-y.^2)/4,并按Plot键,弹出如图8所示的绝对误差值。

7

步骤7:选择File下的Save as选项,选择一个文件存放路径。最后将结果保存为M文件即可。

推荐信息