1. 基于 matlab-PDE Toolbox 的泊松(拉普拉斯)方程求解
在二维电磁场的有限元法计算中,用矩阵方程编制的计算程序长、大,而又复杂,且输入数据要化费很
大的劳动。而 MATLAB 是一种以矩阵运算为基础的交互式语言,它是采用有限元法来求解偏微分方程的。
因此在计算中,我们选用 MATLAB 提供的一个用户图形界面(GUI)的偏微分方程工具箱(PDE Toolbox)
进行数值求解,采用的是三角形网格自动剖分。下面举例说明。
【例 1-1】 横截面为矩形的无限长槽由 3 块接地导体板构成,如图 3-3 所示,槽的盖板接直流电压
100V,求矩形槽的电位分布。
解:这是二维平面场问题。由于电位函数 和电场强度 之间关系为
利用麦克斯韦方程 和关系式 ,得到泊松方程
式中, 为介电常数, 为体电荷密度。由于所求区域内体电荷密度 ,得到拉普拉斯方程:
其边界满足狄里赫利(Didchlet)条件:
,
,
本题运用MATLAB的偏微分方程工具箱(PDE Toolbox)进行数值求解。在命令窗口中输入命令pdetool,
打开 PDE 图形用户界面,计算步骤为:
(1)网格设置:选择菜单 Options 下的 Grid和 Grid Spacing…,将 X-axis linear Spacings 设
置为[-1.5:0.2:1.5],Y-axis linear Spacings取