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文件即可。