求解偏微分方程 -- 椭圆方程实例:求解在域Ω上的泊松方程-Δ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文件即可。
上一篇:怎么使用学生专用计算器?