20138 静态仿真31[优质文档].pptVIP

系统动态的仿真模型;一. 系统仿真（Simulation） 1. 系统仿真： 使用计算机对一个系统的结构和行为进行动态模拟 为决策提供必要的参考信息。 2. 仿真模型： 由计算机程序控制模型的运行 从数值上模仿实际系统的动态行为。;3. 仿真模型的特点： 对象真实、复杂，进行模仿。 4. 关于仿真技术 仿真技术已经成为最重要的和最流行的分析动态系统模型的方法。 ?微分方程的精确解方法具有局限性。 非常多的微分方程我们不知道如何求解。 ?定性分析方法可用于讨论系统的动态行为， 但是对某些问题我们需要定量的答案。 ?仿真技术非常灵活。 可以不困难地将诸如时滞或随机因素等复杂的属性引入模型。 这些是难以用解析的方法处理的。;4. 仿真过程 ? 现实系统的分析： 了解背景，明确目的，提出总体方案。 ?组建模型： 确定变量, 明确关系, 设计流程, 编制程序 ?运行检验： 确定初始状态，参量数值， 运行程序，检验结果，改进模型。 ?输出结果： 清单、记录、重要的中间结果等。; ;回顾——机理模型; ;变量、参量 时间 t，体积 V(t), 盐量 S(t), 浓度 p(t); 流入流速 rI=6, 流入浓度 pI=0.5, 流出流速 rO=4, 改变时候的盐水浓度 p*=0.2，改变的时刻t*， 终止时候的盐水浓度 p*=0.1，终止的时刻t**。 时间步长 Δt , 打印步长 T. 关系: 在 [t, t+Δt] 内有 ; 动态系统仿真的伪代码 运算 池水含盐动态系统模拟 变量 V(n)=时刻 n 池中盐水体积 p(n)=时刻 n 池中盐水浓度 S(n)=时刻 n 池中盐水含盐量 Δt = 时间单位 N = 仿真时间长度 输入 Δt，V(0), p(0), S(0), N ;过程 Begin for n=0 to N do Begin V(n+1)←V(n)+(rI-r0) Δt S(n+1)←S(n)+[ripi-r0p(n)]Δt p(n+1)←S(n+1)/V(n+1) End End 输出 V(1), V(2), … , V(n) S(1), S(2), … , S(n) p(1), p(2), … , p(n) ;系统仿真流程图;考虑到系统运行过程中输入参数的变化，设置指标ind来标志输入条件的转换， 令ind=1表示输入池中盐水的浓度为pI（≠0） ind=0表示流入池中的水为清水（pI=0） 在仿真的过程，可以每运行m个时间步长打印输出一次计算的结果。 ;输入条件改变的系统流程图;Matlab 程序 t=1; v=[2000];s=[2];p=[1/1000];% 初始状态 ri=6;ro=4;po=0.5;p1=0.2;p2=0.1;% 参数 V=[v(end)]; S=[s(end)]; P=[p(end)]; x=[0];% 打印记录 while p(end)p1%调整输入前动态 t=t+1; %时钟步进 v=[v, 1]; s=[s, 1]; p=[p, 0]; %变量步进 v(t)=v(t-1)+ri-ro; s(t)=s(t-1)+ri*po-ro*p(t-1); p(end)=s(end)/v(end); %仿真计算 end; t1=t-1; %调整输入时间 for k=1:floor(t1/20);% 打印结果 x=[x,20*k];V=[V,v(20*k+1)];S=[S,s(20*k+1)];P=[P,p(20*k+1)]; end if floor(t1/20)t1/20 x=[x,t1];V=[V,v(end)];S=[S,s(end)];P=[P,p(end)]; end;po=0;% 调整输入 while p(end)=p2%调整后动态 t=t+1; %时钟步进 v=[v, 1]; s=[s, 1]; p=[p, 0]; %变量步进 v(t)=v(t-1)+ri-ro; s(t)=s(t-1)+ri*po-ro*p(t-1); p(end)=s(end)/v(end); %仿真计算 end;