SS_分布参数系统仿真.ppt

温度场ex06_pde_heat 分布参数系统求解步骤 问题描述 单位矩阵网格化 边界条件和仿真时间 方程求解 动画显示 视频存储 g=squareg; b=squareb1; c=1;a=0;f=1;d=1; [p,e,t]=initmesh(g); u0(ix)=ones(size(ix)); tlist=linspace(0,1,nframes); uu=parabolic(u0,tlist,b,p,e,t,c,a,f,d); u=tri2grid(p,t,u1(:,j),tn,a2,a3); i=find(isnan(u)); u(i)=zeros(size(i)); surf(x,y,u);contour(u,10); colormap(cool) movie2avi(mv,heat.avi,FPS,2); 典型的偏微分方程 高阶系统通常可以转化为一阶微分方程组，典型的偏微分方程包括： 双曲线方程 对流方程 波动方程 抛物线方程 扩散方程 椭圆方程 泊松方程 2.1 差分法原理 基本原理：将时间和空间两个尺度的变量离散化，从而将偏微分方程化为一系列的代数差分方程，利用边界条件和初始条件递推求解。 时空划分如图 示例 设一个长度为L的棒体，棒体表面绝热，t=0时，温度初始分布函数为 ?(x),现在将其两端固定于温度为w1(t)和w2(t)的物体是，试计算不同时间t，棒体上各点的温度演变情况？ 差分法原理 该问题是一个扩散问题 差分法原理 差分法原理 差分法原理 差分法原理 差分法原理 显式和隐式差分法 隐式公式 差分法的稳定性 显式差分法的稳定性： 一阶差分计算公式 从前面的分析可以看出，差分的形式对计算稳定性，计算速度有影响，一阶差分公式如下推导。 一阶差分计算公式 三点中心差分公式 二阶差分计算公式 二阶差分公式可以通过一阶差分公式推导派生出来 2.2 线上积分法 基本原理：将偏微分方程中的空间变量离散化，而时间保持连续，从而将偏微分方程化为一组常微分方程，利用数值积分法求解。 线上积分法原理 在各个离散点出将偏微分转化为常微分 线上积分法原理 对方程可以采用任何一种微分方程数值解法。 比如欧拉法： 线上积分法步骤 将空间变量x等分为M份，共M+1个点 用差分对空间变量求导，f0(u,t)主要利用边界条件， 从时刻t=0开始，采用数值积分方法求出下一个时刻空间各点的函数值。 用差分对空间变量求导fm(u,t) 采用数值积分方法求出t=m+1时刻空间各点的函数值。 重复4,5步骤，直到规定的时刻为之。 线上积分法特点 优点： 方法直观 程序简单，可用现成的数值积分仿真程序 缺点： 误差不易控制，尤其是空间划分及精度难于控制 差分方程很多，不同的差分公式结果的精度不同，较难选择 实例求解（自编程序） 设一个长度为L的棒体，棒体表面绝热，t=0时，温度初始分布函数为 ?(x),现在将其两端固定于温度为w1(t)和w2(t)的物体是，试计算不同时间t，棒体上各点的温度演变情况？ 第三章 有限元方法 有限元（FEM, finite element method ）首出现于Clough R W在1960 年9月发表的一篇讨论弹性力学问题的论文中。在高性能计算机的支持下，得到了蓬勃发展。 有限元的特点： 易于适应区域形状的任意性 易于计算机编程实现 本章节就有限元的基本原理作介绍，Matlab中的PDE工具就是基于有限元的分布参数系统仿真工具。 3.1场论的基础知识 设有一个区域（有限或无限），如果内每一点，都对应着某个物理量的一个确定的值，则称在区域中确定了该物理量的一个场. 若该物理量是数量，则称此场为数量场；温度场、密度场、电位场等为数量场， 若是矢量，则称此场为矢量场。力场、速度场等为矢量场。 此外，若物理量在场中各点处的对应值不随时间而变化，则称该场为稳定场；否则为不稳定场。 当取顶了直角坐标系以后，空间中的点由它的三个坐标(i,j,k)确定 数量场可以用一个数性函数表示： 矢量场可用一个矢性函数表示： 3.1场论的基础知识 标量：只有大小而没有方向的量。如电压U、电荷量Q、电流I、面积S 等。 矢量：具有大小和方向特征的量。如电场强度矢量、磁场强度矢量、作用力矢量、速度矢量等。 标量场：在指定的时刻，空间每一点可以用一个标量唯一地描述，则该标量函数定出标量场。例如物理系统中的温度、压力、密度等可以用标量场来表示。 矢量场：在指定的时刻，空间每一点可以用一个矢量唯一地描述，则该矢量函数定出矢量场。例如流体空间中的流速分布等可以用矢量场来表示。 场论的三度-梯度 数量场u的梯度 场论的三度-梯度 数量场u的梯度