DIMENSION REDUCTION OF STATE SPACE BUILDING MODEL
   

    1 引言
    状态空间建筑建筑热模型[1]采有现代控制论中"状态空间"的概念，对多个房间的建筑物的热过程列写动态平衡方程，其中包括各围护结构内部的导热，各表面与空气之间的对流换热，各表面之间的长波辐射，各房间之间的空气流动以及室内外遮阳等过程的细致描述。对于一个建筑物的动态热过程，此模型表达为
    （1）
    式中，T为建筑物各围护结构表面及其内部节点和空气节点的温度构成的向量，W为室外气象参数（空气温度、太阳辐射等）和室内热源等扰量构成的向量，C是热容阵，A是热导阵，B是边界阵，它们都取决于建筑结构热物性， 表示T对时间τ的导数。（参考《建筑中文网》）
    某个房间的空气温度可以写成：
    (2)
    式中m为热扰量数，n为状态空间维数，{φij}和{λi}分别是由矩阵A、B和C导出的系数向量序列和系数序列。利用（2）式可以很好地模拟室温的动态过程，这已被子BTP程序[3]有用，并在实用中得到验证。
    在描述建筑物动态热过程的诸多方法中，状态空间能灵活方便地处理多个房间多个表面之间的耦合关系，各种热边界条件和各种热扰量。而反应系数法不能考虑各个表面之间的长波辐射，谐波反应法在计算时必须先对各种热扰量进行Fourier 分解。另外，状态空间法可以取任意的模拟时间步长，大到1小时，小到几秒钟表，而有限差分法由于受算法稳定性的限制只能限较小的时间步长。但是状态空间法为保证模拟精度，单个房间的热模型要求的维数一般不低于30，房间多时，维数不更大，因此要花不少CPU时间去完成（2）式的计算。为此，笔者在参考Cools等从的降维理论[2]的基础上，提出一种实用可行的状态空间建筑热模型的降维方法，并用BTP程序模拟实际气候条件下的典型结构房间在建筑热模型降维前后的室温过程。结果表明，7维建筑热模型的模拟结果与39维的模拟结果相当一致，室温最大误差不大于0.1℃。因此，降维的状态空间建筑热模型可以广泛应用于建筑热过程的研究。
   
    2 状态空间建筑热模型的降维过程
    状态空间建筑热模型的降维实际上是先找出状态空间中的主要节点，然后把其它节点集结到这些主要节点上，得到降维后的状态空间建筑热模型。
    （3）
    另一方面，各种热扰量对各节点的作用强度不同，如取Wj（τ）为Dirac函数，则（2）式变成
    （4）
    于是，各节点对第j种热扰量的响应为
    （5）
    经归功一化处理，得到各种热扰量对各节点的综合作用强度为
    （6）
    其中，
    （7）
    因此，在模拟建筑热过程时，只要选取综合响应强度大于δ（控制精度的常数，一般可取百分之几）的节点作为主要节点。即选j1 ，j2，…，jq，使
    Ei≥δ，i= j1 ，j2，…，jq (8)
    综合考虑以上两方面，取{ j1 ，j2，…，jq}与{j1 ，j2，…，jq}的交集作为主要节点，记它们相应的为λi为λ*1 ，λ*2，…，λ*k,下面的工作是把其它节点集结到这些节点，也就是寻找{φij}，使
    （9）
    最小，经类似最小二乘法的推导，从（9）式可得到
    （10）
    式中，
   
   
   
    … … … … … … … … … …
   
   
    通过以上分析，可以把n 维建筑热模型降到k维，且k＜＜n。
   

    i

    τi/min

    Ei/%

    i

    τi/min

    Ei/%

    * 1

    3.9

    4.5

    21

    10.3

    0.1

    2

    965.4

    0.3

    22

    34.6

    0.6

    3

    8

    0

    *23

    301.3

    4

    *4

    5176.2

    31.6

    24

    14.6

    0

    5

    7.9

    0

    25

    11.9

    0.2

    6

    20.2

    0

    26

    3.6

    0

    7

    11.6

    0

    *27

    4.4

    2

    8

    101.2

    0.4

    28

    9.9

    0

    9

    8.4

    0

    29

    63.9

    0.2

    10

    24.1

    0.8

    30

    18.1

    0.3

    11

    9

    0

    31

    58.9

    0

    12

    54.6

    0.1

    32

    32.3

    0

    *13

    14.4

    1.1

    33

    11.7

    0

    *14

    227.2

    1.5

    34

    150.9

    0

    15

    7.1

    0

    35

    20.9

    0.2

    16

    3.6

    0

    36

    13.6

    0.5

    17

    9

    0

    37

    1.2

    0

    18

    715.5

    0

    38

    3.3

    16

    19

    31.4

    0

    39

    2.9

    35.3

    20

    9.2

    0.1