空调制冷装置与系统仿真第13章常见热工及其它系统的建模与仿真算例.pptVIP

状态空间问题 状态空间问题（集总参数）： Sfun_statespace.m ..\..\..\..\..\ProgramFiles\matlab13\work\sfun_statespacemdl.mdl 第五课 s-function的高级应用 为s-function封装，可以较容易地输入变量参数。 example_sfun_bengjibianliang.mdl bengjitiao_sfun_bianliang.m derivative 1点 sys(1)=k/(h*h)*(x(2)-2*x(1)+t0) 2点sys(2)=k/(h*h)*(x(3)-2*x(2)+x(1)) 3点sys(3)=k/(h*h)*(x(4)-2*x(3)+x(2)) 4点sys(4)=k/(h*h)*(x(5)-2*x(4)+x(3)) 5点sys(5)=k/(h*h)*(2*x(4)-2*x(5)) output sys(1)=x(1); sys(2)=x(2); sys(3)=x(3); sys(4)=x(4); sys(5)=x(5); 输出计算结果图形 10个点 % Return the derivatives for the continuous states. %============================================ %function sys=mdlDerivatives(t,x,u,t0,k,h) sys(1) = [k/(h*h)*(x(2)-2*x(1)+t0)]; sys(2) = [k/(h*h)*(x(3)-2*x(2)+x(1))]; sys(3) = [k/(h*h)*(x(4)-2*x(3)+x(2))]; sys(4) = [k/(h*h)*(x(5)-2*x(4)+x(3))]; sys(5) = [k/(h*h)*(x(6)-2*x(5)+x(4))]; sys(6) = [k/(h*h)*(x(7)-2*x(6)+x(5))]; sys(7) = [k/(h*h)*(x(8)-2*x(7)+x(6))]; sys(8) = [k/(h*h)*(x(9)-2*x(8)+x(7))]; sys(9) = [k/(h*h)*(x(10)-2*x(9)+x(8))]; sys(10) = [k/(h*h)*(2*x(9)-2*x(10))]; % end mdlDerivatives 第三课蓄热加热炉的数字仿真(分布参数)问题描述 蓄热砖填充高度H,初始时刻 t=t0时砖温为Ts0,砖间气流温度为TF0, 通入气流率为G，进口处气流温度随时间变化的规律为：TF(0,t)=F(t),现计算炉内任何高度h处，任何时间t时，气流及炉砖的温度TF(h,t),TS(h,t)。 H h G TF(0,t)=F(t) TF(H,t) TF(h,t) Ts(h,t) 建立数学模型 描述该系统的 偏微分方程组为： 其边界条件为TF(0,t)=F(t) 初始条件为: TF(h,0)=TF0 Ts(h,0)=Ts0 其中ABC均由系统的物理性质决定 建立仿真模型 TF和TS均为h和t的函数，但在任何固定位置h处， TF和TS将仅是t的函数。现将炉膛沿高度方向分为n段，hi为其中任意一点，用中心差分公式代换右式 仿真模型 可得常微分方程组 利用S-function进行计算机仿真 建立框图模型：examtemperature.mdl 较为复杂 X(0)=[20,20,20,20,20,20,20,20,20,20,20,20] TF(0)=F(t)=50*1(t) A=-2,b=-0.1,deth=1,c=0.04 Derative(状态方程) 1点气流sys(1)=[a/(2*deth)*(x(3)-u)+b*(x(1)-x(2))]; 1点砖sys(2)=[c*(x(1)-x(2))]; 2点气流sys(3)=[a/(2*deth)*(x(5)-x(1))+b*(x(3)-x(4))]; 2点砖sys(4)=[c*(x(3)-x(4))]; 3点气流sys(5)=[a/(2*deth)*(x(7)-x(3))+b*(x(5)-x(6))]; 3点砖sys(6)=[c*(x(5)-x(6))]; Derative(状态方程) 4点气流sys(7)=[a/(2*deth)*(x(9)-x(5))+b*(x(7)-x(8))]; 4点砖sys(8)=[c*(x(7)-x(8))]; 5点气流sys(9)=[a/(2*deth)*(x(11)-x(7))+b*(x(9)-x(10))]; 5点砖sys(10)=[c*(x(9)-x(10))]; 6点气流sys(11)=