第二章,数值数组及其运算.ppt

* * MATLAB可以求解一般的偏微分方程，也可以利用偏微分方程工具箱（PDE　Toolbox）中给出的相当函数求解一些偏微分方程。 1、偏微分方程组的求解 考虑如下的偏微分方程： (四) 偏微分方程 该偏微分方程可以编写下面的函数表示，其入口为： 其中pdefun为函数名。由给定的输入变量即可计算出c,f,s这三个函数。 * * 边界条件可以用以下函数描述： 边值函数可以编写为一个MATLAB函数： 初始条件函数可以表示为： 即： MATLAB提供了pdepe()函数来求解该问题的数值解。其基本调用格式为： sol = pdepe(m,@pdefun,@pdeic,@pdebc,x,t) 其中有t0ttf ，axb,，m=0、1或2。如果m0，则必须a0。 * * 【例】试求解下面的偏微分方程。 其中：　　　　　　　 初始条件： 边界条件： * * 原方程改写为： 由此可见m=0，且： 描述偏微分方程的函数： function [c,f,s]=c7mpde(x,t,u,du) c=[1;1];y=u(1)-u(2);F=exp(5.73*y)-exp(-11.46*y); s=F*[-1;1]; f=[0.024*du(1);0.17*du(2)]; * * 写出边值方程： 左边界：　 　　　　　右边界： 描述边界条件的MATLAB函数： function [pa,qa,pb,qb]=c7mpbc(xa,ua,xb,ub,t) pa=[0;ua(2)];qa=[1;0];pb=[ub(1)-1;0];qb=[0;1]; 描述初值的MATLAB函数： u0=@(x)[1;0]; 求解此微分方程，并得出解. x=0:.05:1;t=0:.05:2;m=0; sol=pdepe(m,@c7mpde,u0,@c7mpbc,x,t);surf(x,t,sol(:,:,1)) figure;surf(x,t,sol(:,:,2)) * * U1(t,x) * * U2(t,x) * * MATLAB除提供了函数求解偏微分方程之外，还有自己的偏微分工具箱PDE Toolbox，计算机仿真求解的偏微分方程类型可分为： 椭圆型方程： 抛物型方程： 双曲型方程： 特征值问题： 　　特征值偏微分方程中不含参数f。 2、二阶偏微分方程的数学描述 * * 在MATLAB提示符下键入pdetool即可以启动偏微分方程求解界面。包括以下几部分： 菜单系统 工具栏 集合编辑 求解区域 绘制矩形 绘制椭圆 绘多边形 区域边界 输入方程 稀疏网格 加密网格 求解方程 图形设置 局部放大 可解方程类型列表 集合编辑 (Set Formula) * * 偏微分求解区域的绘制 首先：绘制出以下一些区域，然后用集合编辑栏目中将原来的内容修改为(R1+E1+E2)-E3，表示从矩形R1，椭圆E1和E2的并集中剔除掉E3。单击区域边界按钮即可得到求解区域。 * * 偏微分求解区域的绘制 其次，选择Boundary—Remove All Subdomain Borders,即可消除相邻区域中间的分隔线，得到下图。 * * 偏微分求解区域的绘制 然后，单击稀疏网格按钮划分成若干网格，得到下图。 * * 偏微分求解区域的绘制 如果感觉网格不够密，可以单击右侧的加密网格，得到更密的网格图，如下图。一般来说，网格越密，计算结果越精确，但代价是计算时间越长。 * * 【例25】求解热传导方程　　　　　，边界条件是齐次类型，定解区域自定。 求解过程： （1）启动MATLAB，键入pdetool，并回车进入GUI。在Options菜单下选择Grid命令，打开栅格，栅格可以使用户容易确定所绘图形的大小。 （2）选定定解区域。 由于本题要求定解区域自定，故定解区域为如下图所示：E1－E2＋R1－E3，然后在Set formula中进行编辑并用算术运算符对图形对象名称连接起来。 或删去默认的表达式，直接键入： E1－E2＋R1－E3 * * 所讨论定解问题的区域 * * （3）选取边界 　　选择Boundary菜单中Boundary Mode命令，进入边界模式。然后单击Boundary菜单中Remove All Subdomain Borders选项，从而去掉子域边界，如下图。单击Boundary菜单中Specify Boundary Conditions选项，打开Boundary Conditons对话框，输入边界条件。本例取默认条件，即将全部边界设为齐次Dirichlet条件，边界显示为红色。如果想将几何与边界信息存储，可选择Boundary菜单中的Export Decomposed Geometry，Bound