计算机仿真技术第三章数值积分法在系统仿真中的应用.pptVIP

? parabolic 调用格式为：u1=parabolic(u0,tlist,b,p,e,t,c,a,f,d) 该命令用来求解抛物线型偏微分方程(3.4.33)。 ? pdeeig 调用格式为：[v,l]=pdeeig(b,p,e,t,c,a,d,r) 该命令用来求解特征值型偏微分方程(3.4.37)。 ? pdenonlin 调用格式为：[u,res]=pdenonlin(b,p,e,t,c,a,f) 该命令使用具有阻尼的Newton迭代法，在由参量p、e、t确定的网格上求解非线性椭圆型偏微分方程(3.4.31)。 ? poisolv 该命令在一个矩形网格上求解Poisson方程。 ? 结果处理：如Matlab的主要特色一样，在PDE工 具箱中提供了丰富的图形显示，因此用户不但可以对产生的网格进行图形显示和处理，对求解的数据也可以选择多种的图形显示和处理方法，甚至包括对计算结果的动画显示。用户可以参考相关资料来使用。 (3) 应用实例 在这里我们给出一个简单的例子，来说明利用PDE工具箱求解偏微分方程的方法。在Matlab的PDE帮助文件中和在线演示中提供了8个计算实例，可供读者仔细参考。 例3.4.2：最小表面问题求解。 最小表面问题方程可以表示为如下形式： 边界条件为： 。显见这是一个非线性问题，我们用命令pdenonlin来求解。 clc % Lets solve the minimal surface problem % -div( 1/sqrt(1+grad|u|^2) * grad(u) ) = 0 % with u=x^2 on the boundary g=circleg; % The unit circle b=circleb2; % x^2 on the boundary c=1./sqrt(1+ux.^2+uy.^2); a=0; f=0; rtol=1e-3; % Tolerance for nonlinear solver pause % Strike any key to continue. clc % Generate mesh [p,e,t]=initmesh(g); [p,e,t]=refinemesh(g,p,e,t); % Solve the nonlinear problem u=pdenonlin(b,p,e,t,c,a,f,tol,rtol); % Solution pdesurf(p,t,u); pause % Strike any key to end. 计算结果如图3.4.4所示。 三、偏微分方程求解界面 在Matlab中的PDE Toolbox 包括一个图形用户界面(GUI)，在Matlab窗口运行pdetool就进入PDE Toolbox，如图3.4.5所示。GUI主要部分是菜单、对话框和工具条。考虑到本书的篇幅和主要内容，读者可以参考PDE Toolbox的在线帮助和其它资料，在此不再详述。 图3.4.4 最小表面问题求解结果 图3.4.5 PDE Toolbox用户界面 第三章 数值积分法在系统仿真中的应用3-5 面向微分方程的仿真程序设计 一 数字仿真程序的构成 一般来讲，系统设计人员所面临的物理系统的数学模型，可以是微分方程、传递函数或其他形式。如果数学模型是一阶微分方程组的形式，设计人员就可以选择前述方法中的一种进行仿真运算，否则还必须做模型变换。为了使系统设计人员摆脱复杂的程序设计工作，而将精力集中于系统性能的研究和分析上，国内外的仿真学者研制成了多种专门用于系统分析的程序包和各种仿真语言。系统设计人员只要熟悉所选的仿真程序的使用方法就可以。但是由于种种原因，如用户的模型特殊性、通用程序包的计算速度等，也常常需要设计人员来设 计适用于自己的仿真程序，或者根据自己问题的需要，修改已有的仿真程序包或仿真语言。组成程序包的一般原则是：仿真程序应使设计人员使用方便，便于输入及改变参数，便于观察仿真结果，甚至进行自动设计，选择最佳参数；同时要求所组成的仿真程序应在该专业领域内具有一定的通用性和灵活性，并充分利用计算机的各种外部设备。 一般的仿真程序的组成可以用图3.5.1来表示。图中每个方块表示它应具有的功能及相互关系。每个方块的功能大致是：  主程序：实现对整个仿真 计算的逻辑控制。 输入或预制参数块：输入 系统的参数初值、 计算步长、计算时 间等参数。 运行管理块