实验报告
课程名称：过程工程原理实验（乙） 指导老师： 杨国成 成绩：__________________ 实验名称：传热综合实验 实验类型：工程实验 同组学生姓名： 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得
横管对流传热系数的测定
1 实验目的：
1.1 掌握空气在普通和强化传热管内的对流传热系数的测定方法，了解影响传热系数的因素和强化传热的途径。
1.2 把测得的数据整理成n
B N Re u 形式的准数方程式，并与教材中相应公式进行比较。 1.3 了解温度、加热功率、空气流量的自动控制原理和使用方法。
2 装置与流程：
2.1 实验装置如图1所示：
图1.装置示意图
专业：
姓名： 学号：
日期： 2015.12.04
地点： 教十1206
2.2 流程介绍：
实验装置由蒸汽发生器、孔板流量变送器、变频器、套管换热器及温度传感器、智能显示仪等构成。 空气—水蒸气换热流程：来自蒸汽发生器的水蒸气进入套管换热器，与被风机抽进的空气进行热交换，冷凝水经排出阀排入盛水装置。空气经孔板流量计进入套管换热管内（紫铜管），流量通过变频器调节电机转速达到自动控制，热交换后从风机出口排出。
本实验中，普通管和强化管实验通过管路上的切换阀门进行切换。
2.3 横管对流传热系数测定实验数据符号说明表：
名称
符号 单位 备注
冷流体流量 V
紫铜管规格： Φ19mm ×1.5mm ， 即内径为16mm ， 有效长度为1020mm ， 冷流体流量范围： 3~18 m^3/h
冷流体进口温度 t 1 ℃ 普通管冷流体出口温度 t 2 ℃ 强化管冷流体出口温度 t 2’
℃ 蒸汽发生器内蒸气温度 T 1 ℃ 普通管热流体进口端壁温 T W1 ℃ 普通管热流体出口端壁温 T W2 ℃ 普通管外蒸气温度 T ℃ 强化管热流体进口端壁温 T W1 ‘ ℃ 强化管热流体出口端壁温 T W2 ’ ℃ 强化管外蒸气温度
T ’
℃
3 基本原理：
间壁式换热器：冷流体之间有一固体壁面，两流体分别在固体壁面的两侧流动，两流体不直接接触，通过固体壁面进行热量交换。
本装置主要研究汽—气综合换热，包括普通管和强化管。其中，水蒸气空气通过紫铜管间接换热，空
气走紫铜管内，水蒸气走紫铜管外，采用逆流换热。
所谓强化管，是在紫铜管内加了弹簧，增大了绝对粗糙度，进而增大了空气流动的湍流程度，使换热效果更明显。
3.1 空气在传热管内对流传热系数的测定
间壁式传热过程由热流体对固体壁面的对流传热，固体壁面的热传导和固体壁面对冷流体的对流传热所组成，如图2：
T
b
TW tw
t
图2.间壁式传热过程示意图
间壁式传热元件在传热稳定后，有公式：
（1） 式中： Q — 传热量，W ；
m1 — 热流体质量流量，kg/s ；
Cp1 — 定性温度下热流体比热，J/(kg ·℃)； T1 — 热流体进口温度，℃； T2 — 热流体出口温度，℃； m2 — 冷流体质量流量，kg/s ；
Cp2 — 定性温度下冷流体比热，J/(kg ·℃)； t1 — 冷流体进口温度，℃； t2 — 冷流体出口温度，℃；
α1 — 热流体与固体壁面的对流传热系数，W/()； A1 — 热流体测的传热面积，；
(T-Tw)m — 热流体与固体壁面的对数平均温度，℃； α2 — 冷流体与固体壁面的对流传热系数，W/()； A2 — 冷流体测的传热面积，；
(t-tw)m — 冷流体与固体壁面的对数平均温度，℃； K — 以传热面积A 为基准的总传热数，W/()； A — 传热面积，；
Δt — 冷、热流体的对数平均温差，℃。
其中：
'ρV
V =
2'm ρV =
空气质量流量算法：
— 空气之际体积流量，；
— 空气在孔板处的密度，kg/。 本实验中即为空气在进口温度下对应的密度。
空气在无纸记录仪上记录的体积流量与空气流过孔板时候的密度有关。考虑到实际过程中，空气的进口温度不是定值，为了处理上的方便，无纸记录仪上显示的体积流量是孔板处的空气密度为1 kg/时候读数，因此，如果空气的密度不等于该值，则空气的实际体积流量应按照下式进行校正。
（5）
V — 无纸记录仪上显示的空气的体积流量，。
在本装置的套管换热器中，环隙内通水蒸气，紫铜管内通空气，水蒸气在紫铜管表面冷凝放热而加热空气。当内管材料导热性能很好，即λ值很大，且管壁厚度较薄时，可以认为同一截面处换热管两侧壁温近似相等，即：
， （6） 在传热达到稳定后，可得：
（7）
即：
（8）
实验中测定出紫铜管的壁温t W1和t W2 ；冷流体的进出口温度t 1和t 2 ，并查得t 平均=0.5(t 1+t 2)下冷流体对应的C p2以及实验用铜管的长度l 、内径d 2，A 2=πd 2l 和冷流体的质量流量m 2，即可计算出α2。 3.2 准数方程式
对于流体在圆形直管内作强制湍流对流传热时，传热准数经验式为：
（9）
式中： —努赛尔数；
—雷诺数； —普兰特数；
对于流体是被空气冷却的情况，Pr ≈常数，n=0.3，式(9)可简化为：
在双对数坐标中作图，找出直线斜率，即为方程的指数n 。在直线上任取一点的函数值代入方程中，则可得到系数B 。
4 实验步骤：
4.1 检查仪表、风机、蒸汽发生器及测温点是否正常； 4.2 蒸汽发生器灌水至接近最高点处；
4.3 打开总电源开关、仪表电源开关，启动控制软件。开启蒸汽发生器开始加热； 4.4 温度上升到水的沸点后，开启风机；
4.5 选择普通管管路，与仪表显示一致后，通过控制软件上的“流量设定”，按12 m 3/h 、10 m 3/h 、8 m 3
/h 、6 m 3
/h 的数据改变冷流体的流量；
4.6在每个实验点时，观察温度变化，特别注意冷流体出口温度的变化，等待系统稳定。系统稳定后，
冷流体出口温度不再有明显的变化时，点击读取数据按钮，保存得到该冷流体流量下的系统数据；
4.7 完成普通管的测定后，切换到强化管，按第4.5、4.6两步相同的操作进行测定和数据的采集；
4.8完成实验后，按序切断电加热器电源、水源，停风机，最后导出实验数据，关闭计算机。
5 实验数据记录和处理：
5.1 基本参数：紫铜管规格：Φ19mm×1.5mm，即内径为16mm，A=0.050m2，有效长度为1020mm。
5.2 实验数据：
冷流体进口温度
/℃
fAirTemp1 fAirTemp1 19.186 19.75 20 20 20.374 20.436 20.25
冷流体出口温度
/℃
fAirTemp2 fAirTemp2 65.874 68.436 70.934 72.434 78.124 75.436 75.684
冷流体流量(m3/h)
fV fV 20.0055 14.7686 10.9605 7.99963 8.03193 12.6724 16.4429 蒸汽进口处壁温
/℃
fVaporTemp1 fVaporTemp1 102.124 102 102 101.812 101.934 102.812 102.812 蒸汽出口处壁温
/℃
fVaporTemp2 fVaporTemp2 102.374 102.434 102.25 102.124 101.374 102.124 102.124
管外蒸汽温度/℃
ai04 ai04 101.75 101.686 101.686 101.624 99.624 100.812 100.874 加热器汽相温度
/℃
fTemp fTemp 102.06 101.934 101.934 101.874 101.124 102 102
管道名称
iPipeName iPipeName 普通管1 普通管2 普通管3 普通管4 强化管1 强化管2 强化管3
5.3 数据处理：
计算示例：以强化管第一组数据为例：
t1=20.374 ℃ t2’=78.124℃ V=8.03193m 3/h T 1=101.124 ℃ TW1’=101.934℃ TW2’=101.374℃ T ’=99.624℃
t1=20.374℃，空气的密度为1.2035 kg/m 3，空气平均温度40~50℃下空气热容为1005J/(kg ·K)，导热
系数为0.0280/[W/(m ·K)]，空气粘度为1.94×10-5
/Pa ·s.
'ρV
V =
=
2035
.103193
.8=7.3214m 3/h
02'm ρV ==7.3214/3600*1.2035=0.002448 kg/s
Q=)(m 122*2t t c p -=0.002448*1005*（78.124-20.374）=142.06 J/kg
2
w21w12w21w1m w t -t t -t ln
)
t -(t -)t -(t t -t =
）（= [(20.374-101.934)-(78.124-101.374)]/ln[(20.374-101.934)/
(78.124-101.374)]= 46.46℃ α2=
m
w *2122*2t)-(t )(m A t t c p - =142.06/0.051/46.46=59.95W/m2?K
Nu='*22λαd = 59.95*0.016/0.0280=34.26
Re=
μ
ρ
**d 2u =0.016*0.002448/（0.008*0.008*3.14*1.94）*105=1.005*104
分别计算，得到如下表格
在双对数坐标内做Nu~Re 图像为：
普通管y = 0.0249x 0.7637
10100
1000
10000100000
Re
N u
系列1
乘幂 (系列1)
T W1/℃
T W2/℃
t 1/℃
t 2/℃
V ‘/(m 3
/s) m 2/(Kg/s )
Q/(J/k g)
(t W -t )m /K
α/(W /m 2
·K )
Nu Re
（104）
管道类
型
102.124
102.374 19.186 65.874 0.00506 0.00611 286.0
3
56.57 99.14 56.66 2.509 普通管
102 102.434 19.75
68.436 0.00373 0.00450 220.10 54.61
6 79.01
45.12 1.846 普通管 102 102.25 20 70.934 0.00277 0.00333 170.6
7 52.65 63.56 36.29 1.368 普通管 101.812 102.124 20 72.434 0.00202 0.00244
128.6
6 46.84 49.06 28.01 1.002 普通管 101.934 101.374 20.374 78.124 0.00203 0.00245 142.0
6 46.46 59.95 34.26 1.005 强化管 102.812 102.124 20.436 75.436 0.00321 0.00386 213.3
6 49.41 84.6
7 48.35 1.584 强化管 102.812
102.124
20.25
75.684
0.00416
0.00501
279.11 49.29 111.1
63.44
2.056
强化管
加强管y = 0.0135x 0.8491
10100
1000
10000100000
Re
N u
系列1
乘幂 (系列1)
6 实验结果分析： 6.1准数方程式：
根据图像及双对数坐标下直线方程，可得： 对于普通管，有Nu=7637
.0e 0249.0R ； 对于强化管，有Nu=8491
.0e
0135.0R 。
6.2 结果分析：
6.2.1 从实验结果可以发现，热流体流速越大，横管对流传热系数越大，热量交换越迅速，越明显，换热效果也越好。另外，随着热流体温度的升高，平均温差变化很小，可以认为基本上没有发生变化，故在生产中在负载允许的情况下，应适当增大流速，以更好地换热任务。
6.2.2 在管内加装弹簧，增大空气流动时的湍流程度后，空气在换热过程中所获得的热量增大，所以强化管的α较普通管有明显的增大，说明传热效果更好。因此在工程上，使用列管式换热器时，常通过增加壳程挡板来增大湍流程度，以获得更好的换热效果。
6.2.3 通过加弹簧方法来提高传热效果，由于污垢不断积累所产生的热阻和流体运行阻力的增加，使得效果的提高并没有想象中那么理想，所以实际生产中很少用这样的方法。 6.3 误差分析：
无论是普通管还是强化管，实验得到的Nu 与Re 的关系式与公认的经验式（Nu=0.023Re 0.8
）有一定的偏差。分析起来，产生偏差的主要原因可能有：
（1）测量仪器本身的系统误差和外界因素的干扰。比如由于长期有人做实验，导致管内污垢层积，使管内热阻增大，导致测得的数据发生偏离。
（2）实验数据处理过程中，假设内管材料导热性能很好，即λ值很大，且管壁厚度较薄时，所以可以认为同一截面处换热管两侧壁温近似相等，即：， 而实际上内管材料的性能究竟有没有这么好，有待确定。
（3）测定实验数据时，要求在稳定传热状态下，要求稳定时间在8min 以上，而事实上测量数据的时候，在温度稳定5min 左右就已经开始采集数据，而有可能此时传热尚未达到真正的稳定。即使每次数据都在度数稳定8min 采集，也不能保证此时传热已经稳定，因为过程工程原理实验是工程实验，实验结果受实际复杂情况影响较大。
（4）数据处理的时候也会引入一定的误差。
（5）加强管的测量误差较大，可能是因为测量的数据较少，公式未能完全拟合.
6.4思考题：
6.4.1 实验中冷流体和蒸汽的流向，对传热效果有何影响？
在不考虑错流和折流的情况下，实验中冷流体和蒸汽的流向有两种：一种是两者流动方向相同，即并流；另一种是两者流动方向相反，即逆流。
这两种情况的传热效果是不相同的，因为在同样的进、出口温度下，逆流的传热温差比并流的大，所以当传热管的传热面积和传热系数不变时，逆流时传热速率更大，传热效果也更好。
6.4.2 在计算冷流体质量流量时所用到的密度值与求雷诺数时的密度值是否一致？它们分别表示什么位置的密度，应在什么条件下进行计算？
不一致。计算冷流体质量流量时的密度是孔板流量计测压的地方的温度对应的密度，应该将压力换算到常压下进行计算。本实验测压处的压力近似认为是常压。
求雷诺数时，所也用的密度是传热管中的定性温度所对应的气体密度，它所表示的是管内空气的平均温度。在处理数据的时候，将雷诺数Re的公式进行转化：Re=
6.4.3 实验过程中，冷凝水不及时排走，会产生什么影响？如何及时排走冷凝水？
冷凝水不及时排走将会附着在管壁上，由于水的导热性能远差于紫铜，所以未及时排走的冷凝水将降低传热的速率，使传热效果变差。要及时排走冷凝水，可以使管束有一定的倾斜，同时也可以设置冷凝液排泄挡板和改善冷凝表面状况等。
7 讨论心得：
7.1 对间壁式传热器的认识：间壁式换热器是一种传统的、应用最广泛的热交换设备。由于它结构坚固，且能选用多种材料制造，故适应性极强，尤其在高温、高压和大型装置中得到普遍应用。按照结构区分，可分为管式和板式，实验中所用为管式，且是横管。
7.2 对强化传热的认识：强化传热的目的是用较少的传热面积或较小的设备完成同样的传热任务（设计），或力求使换热设备在单位时间、单位面积传递的热量尽可能地大。从传热速率方程可以看出，增大传热系数、扩展传热面积和增大传热平均温差均可提高传热效率。
传热膜系数测定实验（第四组） 一、实验目的 1、了解套管换热器的结构和壁温的测量方法 2、了解影响给热系数的因素和强化传热的途径 3、体会计算机采集与控制软件对提高实验效率的作用 4、学会给热系数的实验测定和数据处理方法 二、实验内容 1、测定空气在圆管内作强制湍流时的给热系数α1 2、测定加入静态混合器后空气的强制湍流给热系数α1’ 3、回归α1和α1’联式4.0Pr Re ??=a A Nu 中的参数A 、a * 4、测定两个条件下铜管内空气的能量损失 二、实验原理 间壁式传热过程是由热流体对固体壁面的对流传热，固体壁面的热传导和固体壁面对冷流体的对流传热三个传热过程所组成。由于过程复杂，影响因素多，机理不清楚，所以采用量纲分析法来确定给热系数。 1）寻找影响因素 物性：ρ，μ ，λ，c p 设备特征尺寸：l 操作：u ，βg ΔT 则：α=f （ρ，μ，λ，c p ，l ，u ，βg ΔT ） 2）量纲分析 ρ[ML -3]，μ[ML -1 T -1]，λ[ML T -3 Q -1]，c p [L 2 T -2 Q -1]，l [L] ，u [LT -1]， βg ΔT [L T -2]， α[MT -3 Q -1]] 3）选基本变量（独立，含M ，L ，T ，Q-热力学温度） ρ，l ，μ, λ 4）无量纲化非基本变量 α：Nu ＝αl/λ u: Re ＝ρlu/μ c p : Pr ＝c p μ/λ βg ΔT : Gr ＝βg ΔT l 3ρ2/μ2 5）原函数无量纲化 6）实验 Nu ＝ARe a Pr b Gr c 强制对流圆管内表面加热：Nu ＝ARe a Pr 0.4 圆管传热基本方程： 热量衡算方程： 圆管传热牛顿冷却定律： 圆筒壁传导热流量：)]/()ln[)()()/ln(11221122121 2w w w w w w w w t T t T t T t T A A A A Q -----?-?=δλ 空气流量由孔板流量测量：54.02.26P q v ??= [m 3h -1，kPa] 空气的定性温度：t=(t 1+t 2)/2 [℃]
实验报告 ( 2016-2017年度第二学期) 名称：《现代控制理论基础》 题目：状态空间模型分析 院系：控制科学与工程学院 班级： ___ 学号： __ 学生姓名： ______ 指导教师： _______ 成绩： 日期： 2017年 4月 15日
线控实验报告 一、实验目的： l.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验内容 1 第一题：已知某系统的传递函数为G (s) S23S2 求解下列问题： （1）用 matlab 表示系统传递函数 num=[1]; den=[1 3 2]; sys=tf(num,den); sys1=zpk([],[-1 -2],1); 结果： sys = 1 ------------- s^2 + 3 s + 2 sys1 = 1 ----------- (s+1) (s+2) （2）求该系统状态空间表达式： [A1,B1,C1,D1]=tf2ss(num,den); A = -3-2 10 B = 1 C = 0 1
第二题：已知某系统的状态空间表达式为： 321 A ,B,C 01：10 求解下列问题： （1）求该系统的传递函数矩阵： （2）该系统的能观性和能空性： （3）求该系统的对角标准型： （4）求该系统能控标准型： （5）求该系统能观标准型： （6）求该系统的单位阶跃状态响应以及零输入响应：解题过程： 程序： A=[-3 -2;1 0];B=[1 0]';C=[0 1];D=0; [num,den]=ss2tf(A,B,C,D); co=ctrb(A,B); t1=rank(co); ob=obsv(A,C); t2=rank(ob); [At,Bt,Ct,Dt,T]=canon(A,B,C,D, 'modal' ); [Ac,Bc,Cc,Dc,Tc]=canon(A,B,C,D, 'companion' ); Ao=Ac'; Bo=Cc'; Co=Bc'; 结果： （1） num = 0 01 den = 1 32 （2）能控判别矩阵为： co = 1-3 0 1 能控判别矩阵的秩为： t1 = 2 故系统能控。 （3）能观判别矩阵为： ob = 0 1
实验报告 课程名称：控制理论（乙） 指导老师：林峰 成绩：__________________ 实验名称：MATLAB 仿真实验 实验类型：________________同组学生姓名：__________ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 实验九 控制系统的时域分析 一、 实验目的： 1．用计算机辅助分析的办法，掌握系统的时域分析方法。 2．熟悉Simulink 仿真环境。 二、实验原理及方法： 系统仿真实质上就是对系统模型的求解，对控制系统来说，一般模型可转化成某个微分方程或差分方程表示，因此在仿真过程中，一般以某种数值算法从初态出发，逐步计算系统的响应，最后绘制出系统的响应曲线，进而可分析系统的性能。控制系统最常用的时域分析方法是，当输入信号为单位阶跃和单位冲激函数时，求出系统的输出响应，分别称为单位阶跃响应和单位冲激响应。在MATLAB 中，提供了求取连续系统的单位阶跃响应函数step ，单位冲激响应函数impulse ，零输入响应函数initial 等等。 二、实验内容： 二阶系统，其状态方程模型为 ? 1x -0.5572 -0.7814 1x 1 = + u ? 2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 四、实验要求： 1．编制MATLAB 程序，画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应； （1）画出系统的单位阶跃响应曲线； A=[-0.5572 -0.7814;0.7814 0 ]; B=[1;0];
华北电力大学 实验报告| | 实验名称状态空间模型分析 课程名称现代控制理论 | | 专业班级：自动化1201 学生姓名：马铭远 学号：2 成绩： 指导教师：刘鑫屏实验日期：4月25日
状态空间模型分析 一、实验目的 1.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验仪器与软件 1. MATLAB7.6 环境 三、实验内容 1 、模型转换 图 1、模型转换示意图及所用命令 传递函数一般形式： MATLAB 表示为： G=tf(num,den)，，其中 num,den 分别是上式中分子，分母系数矩阵。 零极点形式： MATLAB 表示为：G=zpk(Z,P,K) ，其中 Z，P ，K 分别表示上式中的零点矩阵，极点矩阵和增益。 传递函数向状态空间转换：[A,B,C,D] = TF2SS(NUM,DEN)； 状态空间转换向传递函数：[NUM,DEN] = SS2TF(A,B,C,D,iu)---iu 表示对系统的第 iu 个输入量求传递函数；对单输入 iu 为 1。
例1：已知系统的传递函数为G（S）= 2 2 3 24 11611 s s s s s ++ +++ ,利用matlab将传递函数 和状态空间相互转换。 解：1.传递函数转换为状态空间模型： NUM=[1 2 4];DEN=[1 11 6 11]; [A,B,C,D] = tf2ss(NUM,DEN) 2.状态空间模型转换为传递函数： A=[-11 -6 -11;1 0 0;0 1 0];B=[1;0;0];C=[1 2 4];D=[0];iu=1; [NUM,DEN] = ss2tf(A,B,C,D,iu); G=tf(NUM,DEN) 2 、状态方程状态解和输出解 单位阶跃输入作用下的状态响应： G=ss(A,B,C,D);[y,t,x]=step(G);plot(t,x). 零输入响应 [y,t,x]=initial(G,x0)其中，x0 为状态初值。
现代控制理论实验报告
实验一系统能控性与能观性分析 一、实验目的 1.理解系统的能控和可观性。 二、实验设备 1.THBCC-1型信号与系统·控制理论及计算机控制技术实验平台； 三、实验内容 二阶系统能控性和能观性的分析 四、实验原理 系统的能控性是指输入信号u对各状态变量x的控制能力，如果对于系统任意的初始状态，可以找到一个容许的输入量，在有限的时间内把系统所有的状态引向状态空间的坐标原点，则称系统是能控的。 对于图21-1所示的电路系统，设iL和uc分别为系统的两个状态变量，如果电桥中 则输入电压ur能控制iL和uc状态变量的变化，此时，状态是能控的。反之，当 时，电桥中的A点和B点的电位始终相等，因而uc不受输入ur的控制，ur只能改变iL的大小，故系统不能控。 系统的能观性是指由系统的输出量确定所有初始状态的能力，如果在有限的时间内根据系统的输出能唯一地确定系统的初始状态，则称系统能观。为了说明图21-1所示电路的能观性，分别列出电桥不平衡和平衡时的状态空间表达式： 平衡时：
由式（2）可知，状态变量iL和uc没有耦合关系，外施信号u只能控制iL的变化，不会改变uc的大小，所以uc不能控。基于输出是uc，而uc与iL无关连，即输出uc中不含有iL的信息，因此对uc的检测不能确定iL。反之式（1）中iL与uc有耦合关系，即ur的改变将同时控制iL和uc的大小。由于iL与uc的耦合关系，因而输出uc的检测，能得到iL的信息，即根据uc的观测能确定iL（ω） 五、实验步骤 1．用2号导线将该单元中的一端接到阶跃信号发生器中输出2上，另一端接到地上。将阶跃信号发生器选择负输出。 2．将短路帽接到2K处，调节RP2，将Uab和Ucd的数据填在下面的表格中。然后将阶跃信号发生器选择正输出使调节RP1，记录Uab和Ucd。此时为非能控系统，Uab和Ucd没有关系（Ucd始终为0）。 3．将短路帽分别接到1K、3K处，重复上面的实验。 六、实验结果 表20－1Uab与Ucd的关系
实验报告 课程名称：过程工程原理实验（乙） 指导老师： 叶向群 成绩：__________________ 实验名称：吸收实验 实验类型：工程实验 同组学生姓名： 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 填料塔吸收操作及体积吸收系数测定 1 实验目的： 1.1 了解填料吸收塔的构造并熟悉吸收塔的操作； 1.2 观察填料塔的液泛现象，测定泛点空气塔气速； 1.3 测定填料层压降ΔP 与空塔气速u 的关系曲线； 1.4 测定含氨空气—水系统的体积吸收系数K y a 。 2 实验装置： 2.1 本实验的装置流程图如图1： 专业： 姓名： 学号： 日期：2015.12.26 地点：教十2109
2.2物系：水—空气—氨气。惰性气体由漩涡气泵提供,氨气由液氮钢瓶提供，吸收剂水采用自来水，他们的流量分别通过转子流量计。水从塔顶喷淋至调料层与自下而上的含氮空气进行吸收过程，溶液由塔底经过液封管流出塔外，塔底有液相取样口，经吸收后的尾气由塔顶排至室外，自塔顶引出适量尾气，用化学分析法对其进行组成分析。 3 基本原理： 实验中气体流量由转子流量计测量。但由于实验测量条件与转子流量计标定条件不一定相同，故转子流量计的读数值必须进行校正。校正方法如下：
3.2 体积吸收系数的测定 3.2.1相平衡常数m 对相平衡关系遵循亨利定律的物系（一般指低浓度气体），气液平衡关系为： 相平衡常数m与系统总压P和亨利系数E的关系如下： 式中：E—亨利系数，Pa P—系统总压（实验中取塔内平均压力），Pa 亨利系数E与温度T的关系为： lg E= 11.468-1922 / T 式中：T—液相温度（实验中取塔底液相温度），K。 根据实验中所测的塔顶表压及塔顶塔底压差△p，即可求得塔内平均压力P。根据实验中所测的塔底液相温度T，利用式（4）、（5）便可求得相平衡常数m。 3.2.2 体积吸收常数 体积吸收常数是反映填料塔性能的主要参数之一，其值也是设计填料塔的重要依据。本实验属于低浓气体吸收，近似取Y≈y、X≈x。 3.2.3被吸收的氨气量,可由物料衡算 (X1-X2) 式中：V—惰性气体空气的流量，kmol/h；
现代控制理论 实 验 报 告 学院：机电学院 学号:XXXXX 姓名：XXXXX 班级：XXXX
实验一系统的传递函数阵和状态空间表达式的转换 一、实验目的 1.熟悉线性系统的数学模型、模型转换。 2.了解MATLAB 中相应的函数 二、实验内容及步骤 1.给定系统的传递函数为 150 3913.4036 18)(23++++= s s s s s G 要求（1）将其用Matlab 表达；（2）生成状态空间模型。 2.在Matlab 中建立如下离散系统的传递函数模型 y (k + 2) +5y (k +1) +6y (k ) = u (k + 2) + 2u (k +1) +u (k ) 3.在Matlab 中建立如下传递函数阵的Matlab 模型 ?????? ??????+++++++++++=7266 11632256 51 2)(2 32 2s s s s s s s s s s s s G 4.给定系统的模型为 ) 4.0)(25)(15() 2(18)(++++= s s s s s G 求（1）将其用Matlab 表达；（2）生成状态空间模型。 5.给定系统的状态方程系数矩阵如下： []0 , 360180,001,010001 1601384.40==???? ? ?????=????? ?????---=D C B A 用Matlab 将其以状态空间模型表示出来。 6.输入零极点函数模型，零点z=1,-2;极点p=-1,2,-3 增益k=1；求相应的传递函数模型、状态空间模型。 三、实验结果及分析 1.程序代码如下： num = [18 36]; den = [1 40.3 391 150]; tf(num,den) ss(tf(num,den))
太原科技大学全国硕士研究生招生考试 业务课考试大纲（初试） 科目代码：837 科目名称：化工原理 1.前言 化工原理课程研究生入学考试主要测试考生化工单元操作的掌握情况。测试分两个方面：一是化工单元过程原理，测试考生基本概念，过程计算和熟悉程度；二是综合应用化工单元过程原理能力，从而对考生有较全面的评价。 2．题型说明 化工原理考试采用闭卷考试，试卷由以下三部分构成： （1）基本概念题：由选择题、填空题和解答题构成。 （2）计算题：包括过程计算、公式推导。 （3）实验题：包括实验设计、实验原理和实验现象解释。 3．考试内容 3.1绪论 （1）化学工程及其发展。 （2）化工原理课程的性质、内容和任务。 （3）四个基本关系:物料衡算、热量衡算、平衡关系及速率关系。 3.2流体流动 （1）流体静力学方程及其应用。 （2）流量与流速、定态与非定态流动、连续性方程式、能量衡算式、柏努利方程式的应用。 （3）牛顿粘性定律与流体的粘度、非牛顿型流体的概念、流动类型与雷诺准数、滞流与湍流、边界层的概念。 （4）流体在直管中的流动阻力、摩擦系数、因次分析、管路上的局部阻力、管路系统中的总能量损失。 （5）并联管路与分支管路。 （6）测速管、孔板与文丘里流量计和转子流量计。 3.3流体输送设备 （1）离心泵的工作原理和主要部件、离心泵的基本方程式、离心泵的性能参数与特性曲线、离心泵的性能改变和换算、离心泵的气蚀现象与允许吸上高度、离心泵的工作点与调节、离心泵的联用、离心泵的类型与选用。其它类型泵,如往复泵、旋转泵、漩涡泵的工作原理和适用范围。 （2）离心通风机的结构、性能参数和选择,离心鼓风机和压缩机、旋转鼓风机、真空泵。 3.4非均相物系的分离 （1）沉降速度、降沉室、沉降槽。 （2）过滤操作的基本概念、过滤基本方程式、恒压过滤、恒速过滤与先恒速后恒压过滤、过滤常数的测定、过滤设备、滤饼的洗涤、过滤机的生产能力。
聊城大学实验报告 课题名称：化工原理实验 实验名称：总传热系数的测定 姓名：元险成绩： 学号：1989 班级： 实验日期：2011-9-18 实验内容：测定套管换热器中水—水物系在常用流速范围内的总传热系数K，分析强化传热效果的途径。
总传热系数的测定 一、实验目的 1.了解换热器的结构，掌握换热器的操作方法。 2.掌握换热器总传热系数K 的测定方法。 3.了解流体的流量和流向不同对总传热系数的影响 二、基本原理 在工业生产中，要完成加热或冷却任务，一般是通过换热器来实现的，即换热器必须在单位时间内完成传送一定的热量以满足工艺要求。换热器性能指标之一是传热系数K 。通过对这一指标的实际测定，可对换热器操作、选用、及改进提供依据。 传热系数K 值的测定可根据热量恒算式及传热速率方程式联立求解。 传热速率方程式： Q =kS ?t m （1） 通过换热器所传递的热量可由热量恒算式计算，即 Q =W h C ph （T 1－T 2）=W c C pc （t 2－t 1）＋Q 损 （2） 若实验设备保温良好，Q 损可忽略不计，所以 Q =W h C ph （T 1－T 2）=W c C pc （t 2－t 1） （3） 式中，Q 为单位时间的传热量，W ；K 为总传热系数，W/(m 2·℃)；?t m 为传热对数平均温度差，℃；S 为传热面积（这里基于外表面积），m 2；W h ,W c 为热、冷流体的质量流量，kg/s ；C ph ,C pc 为热、冷流体的平均定压比热，J/(kg ·℃)；T 1,T 2为热流体的进出口温度，℃；t 1,t 2为冷流体的进出口温度，℃。 ?tm 为换热器两端温度差的对数平均值，即 12 1 2ln t t t t t m ???-?=? （4） 当212≤??t t 时，可以用算术平均温度差(2 12t t ?+?)代替对数平均温度差。由上式所计算出口的传热系数K 为测量值K 测。 传热系数的计算值K 计可用下式进行计算： ∑+++=S i R K λδαα11 10计 （5） 式中，α0为换热器管外侧流体对流传热系数，W/(m 2·℃)；αi 为换热器管内侧流体对流传热系数，W/(m 2·℃)；δ为管壁厚度，m ；λ——管壁的导热系数，W/(m 2·℃)；R S 为污垢热阻，m 2·℃/W 。 当管壁和垢层的热阻可以忽略不计时，上式可简化成：
河南工业大学 现代控制理论实验报告姓名：朱建勇 班级：自动1306 学号：201323020601
现代控制理论 实验报告 专业: 自动化 班级: 自动1306 姓名: 朱建勇 学号: 201323020601 成绩评定: 一、实验题目： 线性系统状态空间表达式的建立以及线性变换 二、实验目的 1. 掌握线性定常系统的状态空间表达式。学会在MATLAB 中建立状态空间模型的方法。 2. 掌握传递函数与状态空间表达式之间相互转换的方法。学会用MATLAB 实现不同模型之 间的相互转换。 3. 熟悉系统的连接。学会用MATLAB 确定整个系统的状态空间表达式和传递函数。 4. 掌握状态空间表达式的相似变换。掌握将状态空间表达式转换为对角标准型、约当标准 型、能控标准型和能观测标准型的方法。学会用MATLAB 进行线性变换。 三、实验仪器 个人笔记本电脑 Matlab R2014a 软件 四、实验内容 1. 已知系统的传递函数 (a) ) 3()1(4)(2++=s s s s G
(b) 3486)(22++++=s s s s s G
(c) 6 1161)(232+++++=z z z z z z G （1）建立系统的TF 或ZPK 模型。 （2）将给定传递函数用函数ss( )转换为状态空间表达式。再将得到的状态空间表达式用函 数tf( )转换为传递函数，并与原传递函数进行比较。 （3）将给定传递函数用函数jordants( )转换为对角标准型或约当标准型。再将得到的对角 标准型或约当标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。 （4）将给定传递函数用函数ctrlts( )转换为能控标准型和能观测标准型。再将得到的能控标 准型和能观测标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。
化工原理实验仿真软件简介 在教育领域中，计算机不仅是一门学科，而且正逐渐成为有效的教学媒体和教育管理的有力工具。计算机辅助教学是以计算机为媒介，通过学生——计算机之间的交互活动达到教学目的的一种手段。 1. 化工原理实验模拟的发展 实验模拟(Experiment Imitation)是利用计算机的高级图形功能模拟真实的实验环境，通过计算机与操作者之间的交互活动，达到辅助实验教学的目的。实验模拟既是计算机辅助教学的一个重要组成部分，也可以自成体系，这种现代化的新方法，有助于培养学生分析问题、处理问题、解决问题的能力。 化工原理实验模拟系统为辅助化工原理实验教学而设计的软件包。近年来，国内许多高校在化工原理实验模拟方面做了大量的工作，因为化工原理实验模拟必须依托实际的实验装置，而各高校的化工原理实验装置不尽相同，再加上实验模拟投资小、运行费用低、安全、高效等特点，因而受到了高度的重视。北京化工大学早在1985年就开发了一套多功能的单元操作实验模拟软件系统，该系统具有动态画面、音响效果、启发教学、错误处理、自动评分等功能特点。由于开发时间较早，其最大的缺陷是不能独立于西文DOS系统运行，而需要CCDOS 中文汉字系统支撑。华南理工大学开发的化工原理实验模拟系统则较先进，该系统自带中文字库，可以脱离中文汉字系统运行，并且具有窗口式中文界面提示、画面清晰、动画与声响结合。此外，浙江大学开发的化工原理实验模拟系统软件的特点是以该校的实际装置为依托，图像具有3D立体效果。从以上的开发成果可以看出，化工原理实验模拟软件从最初的非中文界面，发展到依托中文操作系统，再发展到自带中文字库脱离汉字操作系统，最后发展到充分利用多媒体技术和3Ｄ图像技术，而且界面日趋友好，功能日渐增多。 2. 化工原理实验模拟的特点 化工原理实验模拟通过计算机模拟真实的实验操作，使学生能快速地掌握如何操作化工单元过程，熟练地测定、整理实验数据，而且可以提高学生对化工原理理论课程的学习兴趣。它具有如下特点： (1) 实验模拟可以模拟传统实验过程，形象生动、简明易懂，既有科学性又富有趣味性，有利于增强教学效果，可在较短的时间内使学生了解化工原理实验单元操作的方法和技巧。 (2) 实验模拟可以快速完成耗费时间很长的实验，并可不断地重复各个实验过程，有利于提高实验教学效果，降低实验运行费用。 (3) 实验模拟可以按实验者的意图任意改变“实验条件”，模拟许多非正常的操作，有利于改善学生在实验装置上操作的安全性。 (4) 实验模拟可以清晰地观察实验的变化规律，使学生获得更多的感性认识，有利于培养学生理论联系实际的能力。 3. 化工原理实验模拟系统的组成 开发化工原理实验模拟系统的目的在于将先进的模拟技术与传统的实验教学相结合，改进实验教学的效果，提高实验教学水平。该模拟系统以基于Windows
实验三 对流传热实验 一、实验目的 1.掌握套管对流传热系数i α的测定方法，加深对其概念和影响因素的理解，应用线性回归法，确定关联式4.0Pr Re m A Nu =中常数A 、m 的值； 2.掌握对流传热系数i α随雷诺准数的变化规律； 3．掌握列管传热系数Ko 的测定方法。 二、实验原理 ㈠ 套管换热器传热系数及其准数关联式的测定 ⒈ 对流传热系数i α的测定 在该传热实验中，冷水走内管，热水走外管。 对流传热系数i α可以根据牛顿冷却定律，用实验来测定 i i i S t Q ??= α （1） 式中：i α—管内流体对流传热系数，W/(m 2?℃)； Q i —管内传热速率，W ； S i —管内换热面积，m 2； t ?—内壁面与流体间的温差，℃。 t ?由下式确定： 2 2 1t t T t w +- =? （2） 式中：t 1，t 2 —冷流体的入口、出口温度，℃； T w —壁面平均温度，℃； 因为换热器内管为紫铜管，其导热系数很大，且管壁很薄，故认为内壁温度、外壁温度和壁面平均温度近似相等，用t w 来表示。 管内换热面积： i i i L d S π= （3） 式中：d i —内管管内径，m ； L i —传热管测量段的实际长度，m 。
由热量衡算式： )(12t t Cp W Q m m i -= （4） 其中质量流量由下式求得： 3600 m m m V W ρ= （5） 式中：m V —冷流体在套管内的平均体积流量，m 3 / h ； m Cp —冷流体的定压比热，kJ / (kg ·℃)； m ρ—冷流体的密度，kg /m 3。 m Cp 和m ρ可根据定性温度t m 查得，2 2 1t t t m +=为冷流体进出口平均温度。t 1，t 2， T w ， m V 可采取一定的测量手段得到。 ⒉ 对流传热系数准数关联式的实验确定 流体在管内作强制湍流，被加热状态，准数关联式的形式为 n m A Nu Pr Re =. （6） 其中： i i i d Nu λα= ， m m i m d u μρ=Re ， m m m Cp λμ=Pr 物性数据m λ、m Cp 、m ρ、m μ可根据定性温度t m 查得。经过计算可知，对于管内被加热的空气，普兰特准数Pr 变化不大，可以认为是常数，则关联式的形式简化为： 4.0Pr Re m A Nu = （7） 这样通过实验确定不同流量下的Re 与Nu ，然后用线性回归方法确定A 和m 的值。 ㈡ 列管换热器传热系数的测定 管壳式换热器又称列管式换热器。是以封闭在壳体中管束的壁面作为传热面的间壁式换热器。这种换热器结构较简单，操作可靠，可用各种结构材料（主要是金属材料）制造，能在高温、高压下使用，是目前应用最广的类型。由壳体、传热管束、管板、折流板（挡板）和管箱等部件组成。壳体多为圆筒形，
本科实验报告 实验名称：控制理论基础（实验）
实验一：控制系统的模型建立 一、实验目的 1.掌握利用MATLAB 建立控制系统模型的方法。 2.掌握系统的各种模型表述及相互之间的转换关系。 3. 学习和掌握系统模型连接的等效变换。 二、实验原理 1、系统模型的 MATLAB描述 系统的模型描述了系统的输入、输出变量以及内部各变量之间的关系，表征一个系统的模型有很多种，如微分方程、传递函数模型、状态空间模型等。这里主要介绍系统传递函数（TF）模型、零极点增益（ZPK）模型和状态空间（SS）模型的MATLAB 描述方法。 1）传递函数（TF）模型 传递函数是描述线性定常系统输入-输出关系的一种最常用的数学模型，其表达式一般为 在MATLAB 中，直接使用分子分母多项式的行向量表示系统，即 num = [bm, bm-1, … b1, b0] den = [an, an-1, … a1, a0] 调用tf 函数可以建立传递函数TF对象模型，调用格式如下： Gtf = tf(num,den) Tfdata 函数可以从TF对象模型中提取分子分母多项式，调用格式如下： [num,den] = tfdata(Gtf) 返回cell 类型的分子分母多项式系数 [num,den] = tfdata(Gtf,'v') 返回向量形式的分子分母多项式系数 2）零极点增益（ZPK）模型 传递函数因式分解后可以写成
式中, z1 , z2 , …,zm 称为传递函数的零点， p1,p2,…,pn称为传递函数的极点，k 为传递系数（系统增益）。 在MATLAB 中，直接用[z,p,k]矢量组表示系统，其中z，p，k 分别表示系统的零极点及其增益，即： z=[z1,z2,…,zm]; p=[p1,p2,…,pn]; k=[k]; 调用zpk 函数可以创建ZPK 对象模型，调用格式如下： Gzpk = zpk(z,p,k) 同样，MATLAB 提供了zpkdata 命令用来提取系统的零极点及其增益，调用格式如下： [z,p,k] = zpkdata(Gzpk) 返回cell 类型的零极点及增益 [z,p,k] = zpkdata (Gzpk,’v’) 返回向量形式的零极点及增益 函数pzmap 可用于求取系统的零极点或绘制系统得零极点图，调用格式如下： pzmap(G) 在复平面内绘出系统模型的零极点图。 [p,z] = pzmap(G) 返回的系统零极点，不作图。 3）状态空间（SS）模型 由状态变量描述的系统模型称为状态空间模型，由状态方程和输出方程组成： 其中：x 为n 维状态向量；u 为r 维输入向量； y 为m 维输出向量； A 为n×n 方阵，称为系统矩阵； B 为n×r 矩阵，称为输入矩阵或控制矩阵；C 为m×n 矩阵，称为输出矩阵； D为m×r 矩阵，称为直接传输矩阵。 在MATLAB 中，直接用矩阵组[A,B,C,D]表示系统，调用ss 函数可以创建ZPK 对象模型，调用格式如下： Gss = ss(A,B,C,D) 同样，MATLAB 提供了ssdata 命令用来提取系统的A、B、C、D 矩阵，调用格式如下： [A,B,C,D] = ssdata (Gss) 返回系统模型的A、B、C、D 矩阵 4）三种模型之间的转换
实验二：总传热系数的测定 一、实验目的 1、了解换热器的结构与用途; 2、学习换热器的操作方法； 3、掌握传热系数k计算方法； 4、测定所给换热器的逆流传热系数k。 二、实验原理 在工业生产过程中冷热流体通过固体壁面（传热元件）进行热量传递，称为间壁式换热。间壁式换热过程由热流体对固体壁面的对流传热，固体壁面的热传导和固体壁面对冷流体的对流传热三部分组成。本实验热流体采用饱和蒸汽走壳程，冷流体为空气走管程。 当传热达到稳定时，总传热速率与冷流体的传热速率相等时， 而即为， 综上可得，其中。 T --- 热流体; t --- 冷流体; V --- 冷流体进口处流量计读数; ---冷流体平均温度下的对应的定压比热容; ρ --- 冷流体进出口平均温度下对应的密度. 三、实验设备及流程 1、实验设备
传热单元实验装置(换热器、风机、蒸汽发生器) ，整套实验装置的核心是一个套管式换热器，它的外管是一根不锈钢管，内管是一根紫铜管。根据紫铜管形状的不同，我们的实验装置配有两组换热器，一种是普通传热管换热器，另一种是强化传热管换热器，本实验以普通传热管换热器为例，介绍总传热系数的测定。 2、实验流程 来自蒸汽发生器的水蒸气从换热器的右侧进入换热器的不锈钢管。而来自风机的冷空气从换热器的左侧进入换热器的紫铜管，冷热流体通过紫铜管的壁面进行传热。冷空气温度升高而水蒸汽温度降低，不凝气体和冷凝水通过疏水阀排出系统，而冷空气通过风机的右侧排出装置。 四、实验步骤 需测量水蒸气进口温度，出口温度，冷空气进口温度，出口温度，冷空气的体积流量以及紫铜管的长度及管径。前四项通过仪表读数可获得，冷空气进口温度可以由另外一块仪表盘读数计算可获得。紫铜
现代控制理论课程总结 学习心得 从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。现代控制论是用状态空间方法表示,概念抽象,不易掌握。对于《现代控制理论》这门课程，在刚拿到课本的时候，没上张老师的课之前，咋一看，会认为开课的内容会是上学期学的控制理论基础的累赘或者简单的重复，更甚至我还以为是线性代数的复现呢！根本没有和现代控制论联系到一起。但后面随着老师讲课的风格的深入浅出，循循善诱，发现和自己想象的恰恰相反，张老师以她特有的讲课风格，精心准备的ppt 课件，向我们展示了现代控制理论发展过程，以及该掌握内容的方方面面，个人觉得，我们不仅掌握了现代控制理论的理论知识，更重要的是学会了掌握这门知识的严谨的逻辑思维和科学的学习方法，对以后学习其他知识及在工作上的需要大有裨益，总之学习了这门课让我受益匪浅。 由于我们学习这门课的课时不是很多，并结合我们学生学习的需求及所要掌握的课程深入程度，张老师根据我们教学安排需要，我们这学期学习的内容主要有：1.绪论；2.控制系统的状态表达式；3.控制系统状态表达式的解；4.线性系统的能空性和能观性；5.线性定常系统的综合。而状态变量和状态空间表达式、状态转移矩阵、系统的能控性与能观性以及线性定常系统的综合是本门课程的主要学习内容。当然学习的内容还包括老师根据多年教学经验及对该学科的研究的一些深入见解。 在现代科学技术飞速发展中,伴随着学科的高度分化和高度综合,各学科之间相互交叉、相互渗透,出现了横向科学。作为跨接于自然科学和社会科学的具有横向科学特点的现代控制理论已成为我国理工科大学高年级的必修课。 经典控制理论的特点 经典控制理论以拉氏变换为数学工具，以单输入－单输出的线性定常系统为主要的研究对象。将描述系统的微分方程或差分方程变换到复数域中，得到系统的传递函数，并以此作为基础在频率域中对系统进行分析和设计，确定控制器的结构和参数。通常是采用反馈控制，构成所谓闭环控制系统。经典控制理论具有明显的局限性，突出的是难以有效地应用于时变系统、多变量系统，也难以揭示系统更为深刻的特性。当把这种理论推广到更为复杂的系统时，经典控制理论就显得无能为力了，这是因为它的以下几个特点所决定。 1．经典控制理论只限于研究线性定常系统，即使对最简单的非线性系统也是无法处理的；这就从本质上忽略了系统结构的内在特性，也不能处理输入和输出皆大于1的系统。实际上，大多数工程对象都是多输入－多输出系统，尽管人们做了很多尝试，但是，用经典控制理论设计这类系统都没有得到满意的结果；2．经典控制理论采用试探法设计系统。即根据经验选用合适的、简单的、工程上易于实现的控制器，然后对系统进行分析，直至找到满意的结果为止。虽然这种设计方法具有实用等很多完整，从而促使现代控制理论的发展：对经典理论的精确化、数学化及理论化。优点，但是，在推理上却是不能令人满意的，效果也
现代控制理论实验报告 组员： 院系：信息工程学院 专业： 指导老师： 年月日
实验1 系统的传递函数阵和状态空间表达式的转换 [实验要求] 应用MATLAB 对系统仿照[例]编程，求系统的A 、B 、C 、阵；然后再仿照[例]进行验证。并写出实验报告。 [实验目的] 1、学习多变量系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数相互转换的方法； 2、通过编程、上机调试，掌握多变量系统状态空间表达式与传递函数相互转换方法。 [实验内容] 1 设系统的模型如式示。 p m n R y R u R x D Cx y Bu Ax x ∈∈∈?? ?+=+=& 其中A 为n ×n 维系数矩阵、B 为n ×m 维输入矩阵 C 为p ×n 维输出矩阵，D 为传递阵，一般情况下为0，只有n 和m 维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式示。 D B A SI C s den s num s G +-== -1)() () (()( 式中，)(s num 表示传递函数阵的分子阵，其维数是p ×m ；)(s den 表示传递函数阵的按s 降幂排列的分母。 2 实验步骤 ① 根据所给系统的传递函数或（A 、B 、C 阵），依据系统的传递函数阵和状态空间表达式之间的关系如式，采用MATLA 的编程。注意：ss2tf 和tf2ss 是互为逆转换的指令； ② 在MATLA 界面下调试程序，并检查是否运行正确。 ③ [] 已知SISO 系统的状态空间表达式为，求系统的传递函数。
, 2010050010000100001 0432143 21u x x x x x x x x ? ? ??? ? ??????-+????????????????????????-=????????????&&&&[]??? ? ? ???????=43210001x x x x y 程序: A=[0 1 0 0;0 0 -1 0;0 0 0 1;0 0 5 0]; B=[0;1;0;-2]; C=[1 0 0 0]; D=0; [num,den]=ss2tf(A,B,C,D,1) 程序运行结果: num = 0 den = 0 0 0 从程序运行结果得到:系统的传递函数为: 2 4253 )(s s s S G --= ④ [] 从系统的传递函数式求状态空间表达式。 程序: num =[0 0 1 0 -3]; den =[1 0 -5 0 0]; [A,B,C,D]=tf2ss(num,den) 程序运行结果: A = 0 5 0 0 1 0 0 0 0 1 0 0
总传热系数测定实验 一、实验目的 1． 观察水蒸气在换热管外壁上的冷凝现象，并判断冷凝类型； 2． 测定饱和水蒸气在圆形管外壁上的冷凝给热系数； 二、基本原理 在套管换热器中，环隙通以水蒸气，内管管内通以空气，水蒸气冷凝放热以加热空气，在传热过程达到稳定后，有如下公式： V ρC P (t 2-t 1)=K A m t ? 其中： V ：空气体积流量，m 3/s A ：内管的外壁的传热面积，m 2 ρ：空气密度，kg/m 3 C P ：空气平均比热，J/(kg ℃) t 1、t 2：空气进、出口温度，℃ T 1、T 2：蒸汽进、出口温度，℃ m t ?：对数平均温差，℃ 1 2211221ln ) ()(t T t T t T t T t m -----= ? 若能测得被加热流体的V 、t 1、t 2，内管的换热面积A 以及水蒸气温度T 1、T 2，即可计算实测的水蒸气（平均）冷凝给热系数。 三、实验装置与流程 实验装置如下图
水蒸气～空气换热流程图 来自蒸汽发生器的水蒸气进入玻璃套管换热器，与来自风机的风进行热交换，冷凝水经疏水器排入地沟。冷空气经孔板（转子）流量计进入套管换热器内管（紫铜管），热交换后排出装置外。 2．设备与仪表规格 （1）紫铜管规格：直径φ21×2.8mm，长度L=1000mm （2）外套玻璃管规格：直径φ100×5mm，长度L=1000mm （3）压力表规格：0～0.1MPa 四、实验步骤与注意事项 1．打开总电源空气开关，打开仪表及巡检仪电源开关，给仪表上电。 2．打开仪表台上的风机电源开关，让风机工作，同时打开冷流体入口阀门，让套管换热器里冲有一定量的空气。 3．打开冷凝水出口阀，注意只开一定的开度，开的太大会让换热桶里的蒸汽跑掉，关的太小会使换热玻璃管里的蒸汽压力集聚而产生玻璃管炸裂。 4．在做实验前，应将蒸汽发生器到实验装置之间管道中的冷凝水排除，否则夹带冷凝水的蒸汽会损坏压力表及压力变送器。关闭蒸汽进口阀门，打开装置下面的排冷凝水阀门，让蒸汽压力把管道中的冷凝水带走，当听到蒸汽响时关闭冷凝水排除阀。 5．刚开始通入蒸汽时，要仔细调节蒸汽进口阀门的开度，让蒸汽徐徐通入换热器中，
实验三 利用MATLAB 导出连续状态空间模型的离散化模型 实验目的： 1、基于对象的一个连续时间状态空间模型，导出其相应的离散化状态空间模型； 2、通过编程、上机调试，掌握离散系统运动分析方法。 实验原理： 给定一个连续时间系统的状态空间模型： ()()()()()() x t Ax t Bu t y t Cx t Du t =+=+ （3.1） 状态空间模型（3.1）的输入信号()u t 具有以下特性： ()(),u t u kT kT t kT T =≤≤+ （3.2） 已知第k 个采样时刻的状态()x kT 和第k 个采样时刻到第1k +个采样时刻间的输入()()u t u kT =，可得第1k +个采样时刻(1)k T +处的状态 (1)((1))((1))()((1))()k T kT x k T k T kT x kT k T Bu d τττ++=Φ+-+Φ+-? （3.3） 其中： ((1))((1))A k T kT AT k T kT e e +-Φ+-== ((1))((1))A k T k T e ττ+-Φ+-= 由于输入信号在两个采样时刻之间都取常值，故对式（3.3）中的积分式进行一个时间变量替换(1)k T στ=+-后，可得 0((1))()()()AT A x k T e x kT e d Bu kT τ σσ+=+? （3.4） 另一方面，以周期T 对输出方程进行采样，得到 ()()()y kT Cx kT Du kT =+ 在周期采样的情况下，用k 来表示第k 个采样时刻kT 。因此，连续时间状态空间模型
（3.1）的离散化方程可以写成 (1)()()()()()()() x k G T x k H T u k y k Cx k Du k +=+=+ （3.5） 其中： 0()()()AT A G T e H T e d B τσσ==? （3.6） 已知系统的连续时间状态空间模型，MATLAB 提供了计算离散化状态空间模型中状态矩阵和输入矩阵的函数： [G ,H]=c2d(A,B,T) 其中的T 是离散化模型的采样周期。 实验步骤 1、导出连续状态空间模型的离散化模型，采用MA TLAB 的m-文件编程； 2、在MA TLAB 界面下调试程序，并检查是否运行正确。 例3.1 已知一个连续系统的状态方程是 010()()()2541x t x t u t ????=+????--???? 若取采样周期0.05T =秒，试求相应的离散化状态空间模型。 编写和执行以下的m-文件： A=[0 1;-25 –4]; B=[0;1]; [G ,H]=c2d(A,B,0.05) 得到 G= 0.9709 0.0448 -1.1212 0.7915 H= 0.0012 0.0448 因此，所求的离散化状态空间模型是 0.97090.04480.0012(1)()()1.12120.79150.0448x k x k u k ????+=+????-????
-150-100 -50 50 实验一 典型环节的模拟研究及阶跃响应分析 1、比例环节 可知比例环节的传递函数为一个常数： 当Kp 分别为0.5，1，2时，输入幅值为1.84的正向阶跃信号，理论上依次输出幅值为0.92，1.84，3.68的反向阶跃信号。实验中，输出信号依次为幅值为0.94，1.88，3.70的反向阶跃信号， 相对误差分别为1.8%,2.2%,0.2%. 在误差允许范围内可认为实际输出满足理论值。 2、 积分环节 积分环节传递函数为： （1）T=0.1(0.033)时，C=1μf (0.33μf )，利用MATLAB ，模拟阶跃信号输入下的输出信号如图： T=0.1 T=0.033 与实验测得波形比较可知，实际与理论值较为吻合，理论上T=0.033时的波形斜率近似为T=0.1时的三倍，实际上为8/2.6=3.08，在误差允许范围内可认为满足理论条件。 3、 惯性环节 i f i o R R U U -=TS 1 CS R 1Z Z U U i i f i 0-=-=-=15 20
惯性环节传递函数为： K = R f /R 1,T = R f C, （1） 保持K = R f /R 1 = 1不变，观测T = 0.1秒，0.01秒（既R 1 = 100K,C = 1μf ， 0.1μf ）时的输出波形。利用matlab 仿真得到理论波形如下： T=0.1时 t s （5%）理论值为300ms,实际测得t s =400ms 相对误差为：（400-300）/300=33.3%,读数误差较大。 K 理论值为1，实验值2.12/2.28， 相对误差为（2.28-2.12）/2.28=7%与理论值 较为接近。 T=0.01时 t s （5%）理论值为30ms,实际测得t s =40ms 相对误差为：（40-30）/30=33.3% 由于ts 较小，所以读数时误差较大。 K 理论值为1，实验值2.12/2.28， 相对误差为（2.28-2.12）/2.28=7%与理论值较为接近 （2） 保持T = R f C = 0.1s 不变，分别观测K = 1，2时的输出波形。 K=1时波形即为（1）中T0.1时波形 K=2时，利用matlab 仿真得到如下结果： t s （5%）理论值为300ms,实际测得t s =400ms 相对误差为：（400-300）/300=33.3% 读数误差较大 K 理论值为2，实验值4.30/2.28， 1 TS K )s (R )s (C +-=