Ch7-求解常微分方程-3.docVIP

在Plot模式下，单击鼠标左键可以把光标所在位置的三角形编号和该处的图形解的属性显示在信息栏。释放鼠标，光吉林省所在位置显示在信息栏的信息消失。 4）当通过菜单定义边界条件和偏微分方程时，虽然没有边界和子区域被选中，但是输入的值在默认的情况下将应用于所有的边界段和子区域 5）在输入参数表达式时一定要注意按照MATLAB合法的书写形式输入，例如函数4xy应该书写为 4*x.*y。 7-5 不同类型偏微分方程的求解 在7-3-2节中通过一个椭圆型方程（泊松方程）的求解讲述了利用PDE工具箱的图形用户界面求解偏微分方程的过程，本节再给出其他类型偏微分方程的影响。 范例7-1 求解双曲弄分方程，该方程是一个波动方程，描述了分布在单位圆上的薄膜的横向振动，方程为： = 初始条件为t=0时的和： 程序设计： 我们使用PDE工具箱的GUI来解这个波动方程，步骤如下： 打开PDE工具箱的GUI窗口，在绘图区里绘出一个椭圆，然后双击椭圆，在弹出的对话框（见图7-2）里各参数设置如下：X-center为0，Y-center为0，A-semiaxes为1，B-semiaxes为1，另两项使用默认值，单击OK按钮，关闭对话框。 单击工具栏按钮，边界条件彩默认值。 单击工具栏PDE按钮，在弹出的对话框（见图7-20）中，选中Hyperbolic项，系数f高为0。单击OK按钮，关闭对话框。 先对单位圆域进行网格剖分，然后再细分网格一次。 点开下拉菜单Solve，单击Parameters,在弹出的对话框（见图7-26）中time设置为0:5:31，设置为，设置为。 点开下拉菜单Plot，单击Parameters，在弹出的对话框（见图7-28）中选中Height（3D-plot）、Show mesh和Animation，其他项采用默认设置。最后单击Plot按钮观看图形解。 输出结果（如图7-30所示）： 图7-30 波动方程的图形（动画）解 程序说明： 我们看到得到的图形解是个动画，展示了5秒钟簿膜振动的情形。在求解参数设置对话框中time设置为0:5:31的意义为0～5秒中计算31个时间点的解值。 数值解u是一个矩阵，通过菜单项Slove中的Export Solution把解导入到工作窨，然后在命令行输入u就可以看到这个矩阵了。 如果用户的PC配置不是很高的话，求解过程花费的时间可能较长。 范例7-2 求解抛物型偏微分方程，该方程是一个热传导方程，描述了通过一根电阻丝（弯成圆形）加热一种介质时，电阻丝所在平面的固体介质截面温度的变化。假设电阻丝的温度恒为5℃，电阻丝圈产半径为1；介质初始温度为0，大小为半径R=2的圆面，边界温度恒为0，外界温度为0，密度为1，热容量C=5，导热系数K=1，热对流系数为1。 程序设计： 我们使用PDE工具箱的GUI求解方程，具体步骤如下： 道德点开Options菜单，在Application中选择Heat Transfer应用模式。 进入Draw模式，在绘图区绘制一个中心在原点，半径为2的圆。 进入Boundary模式，边界条件采用默认值。 进入PDE模式，打开方程类型及其参数设置对话框，设置如图7-31所示。 图7-31 热传导方程的参数设置 单击OK确认关闭对话框。 进入Mesh,模式，对圆域进行网格剖分，然后再细分网格一次。 打开Plot菜单单击Parameter项，在弹出的对话框（见图7-28）中选中Height(3-Dplot)、Animation和Show mesh项，其他选项采用默认设置，单吉Done按钮确认并关闭对话框。 点开Solve菜单，求解参数采用默认设置，然后单击Solve PDE开始计算。 输出结果（如图7-32所示）： 图7-32 热传导方程的图形（动画）解 程序说明： 介质的各个物理参数都设置在图7-31的对话框里，其中热源Q声明为，意为单位圆上间谍为5，表达式中不要出现等号“=”。 从动画中我们可以看到离单位圆的弧线越近，温度长高的越快、越高。 范例7-3 求解一个正方形区域上的特征值： 正方形区域为-1≤x≤1,-1≤y≤1。 程序设计： 我们使用PDE工具箱的GUI求解方程，应用模式使用Generic Scalar，具体步骤如下： 进入Draw模式，在绘图区绘制一个矩形，然后双击矩形，在弹出的对话框（见图7-3）中Left 设置为-1，Botttom设置为-1，Width设置为2，Height设置为2，确认并关闭对话框。 进