我们求解的问题是
上述问题的精确解为 u=sin(3*pi*x+pi/4)*sin(2*pi*y+pi/4);
接下来用Matlab工具箱求解
首先打开matlab,命令窗口处输入:pdetool
将看到如下界面:
显示栅格:Options->Grid
选择求解区域:Draw (这个地方我们的求解区域为[0 1]X[0 1])->Rectange/square 在面板上随意拖动,双击出现的区域,输入如下数值
进入边界条件设置:
上边界与下边界 sin(3*pi*x+pi/4)*sin(pi/4)
(注意:此处不能为\'x\',后同)
左边界 3*pi*sin(2*pi*y+pi/4)*cos(pi/4)
右边界 -3*pi*sin(2*pi*y+pi/4)*cos(pi/4)
方程设置:13*pi*pi.*sin(3*pi*x+pi/4).*sin(2*pi*y+pi/4)
网格
求解