2012深圳杯数学建模夏令营C题3D机房建模（天津农学院） - 图文

2012年 3Ｄ仿真机房建模

摘要：本文主要研究根据机房的基础设施状态，按照行业规范要求合理地布置机柜，分布任务，尽量避免局部地区过热问题，我们根据题目所给定的一些数据和信息分析建立了扩散传播模型、采用了线性回归方法、并基于k~?两方程紊流模型建立数学模型。

首先，本文将绿色机房设计聚焦于机房内热环境分析，然后我们根据附件1给出的数据，通过采用MATLAB软件绘制出机房温度和流场分布，之后求出温度的极值，从而求出温度最高点。

其次，在确定该机房热分布时，考虑到机柜的对称性，本文将每个机柜作为一个热源，建立扩散传播模型，将热量扩散模拟为不同浓度空气在机房进行流通，通过冷通道散发的冷气流与机柜散发的热气流接触，使得机房不同方位的温度不

同，根据最后推导的空气浓度表达式，基于任务量，得出整体机房的热分布情况。

再次，经过对题目的分析，基于上一部分所得模型和附件2中的数据，在充分考虑机房对称性的基础上，选取任务量对称的温度，采用matlab软件计算出温度极值点，再整合机柜任务量对温度极值点的权向量，得出机柜任务量与温度的相关度，再通过线性回归，最后得到最低温度的分配方案。

最后，模型四以机房空气流组织为研究对象，采用k~?两方程紊流模型建立数学模型，对机房的冷通道出口风速与机房室内的温度的关系进行了详细的动态描述，使得气流动态化，从而总结出通过控制冷通道出风口风速来调节室内的温度，达到《电子信息系统机房设计规范》C级要求。

在本次模型建立中，我们得出以下结论： 1.热量在温度极值点会发生集中现象

2.在机房对称的情况下，温度分布还是会收到各种干扰，从而导致温度分布不均

3.当机房开机时应该加大冷通道的通风速度，关机的时候降低冷通道的通风速度

关键词：热分布 机房 气流 扩散传播

1

一、问题重述

1.1问题背景

由于高密度计算、多任务计算的需要，越来越多的高性能数据中心或互联网中心（DC、IDC）正逐渐建成。在现代的数据中心内，由于刀片服务器成本与性价比高，体积小而被广泛使用。由于自身能源与冷却条件限制，这类大规模的数据中心或许每年需要花费数百万美元，主要用于计算设备及系统冷却所需的能源费用。因此有必要提高数据中心设备的能效，极大化数据中心的能源利用率及计算能力。大约在上世纪90年代后期，IBM、HP等公司首先提出绿色数据中心的概念，并受到世界各国的广泛重视。

1.2提出问题

（1）根据附件1的数据，绘出冷、热通道的热分布及流场分布及室内最高温度位置。

（2）建立描述该问题热分布的数学模型及算法，并与测试案例进行比较。 （3）如果定义该机房的总体任务量为1，根据你的模型及附件1的流场数据，确定服务器实际任务量为0.8及0.5的最优任务分配方案，并给出室内最高温度。

（4）如果按照《电子信息系统机房设计规范》（附件3）C级要求控制机房温度，讨论服务器设计任务量一定条件下，如何控制空调的送风速度或送风温度（可以通过送风槽的出口风速与温度来描述）。

二、问题分析

（1）由于整个机房内不同方位的温度是有差异的，根据附件1的数据，采用MATLAB软件绘制出冷热通道的热分布图与流场分布图，并根据热分布图得出了机房室内的最高温度。

（2）由于机柜散发热量及冷通道的气体传播的途径有很多，本题忽略其他因素，只考虑在空气中传播，建立大气扩散模型，通过与测试案例得出机房室内的热分布情况。

（3）如果定义该机房的总体任务量为1，根据你的模型及附件1的流场数据，确定服务器实际任务量为0.8及0.5的最优任务分配方案，并给出室内最高温度。

（4）根据《电子信息系统机房设计规范》C级要求，在一定任务量的基础上，以机房内的气流组织为研究对象，采用k~?两方程紊流模型建立数学模型，并得出了空调如何控制风速来调节机房室内的温度。

三、问题假设

（1）假设每个机柜散发的热量是相通的，且热量分布稳定。

2

（2）假设机柜散发热量是以相同溶度的热空气在机房传播。

（3）假设机柜散发热与冷通道的冷气在机房的分布是不受机房风速影响。 （4）假设机房室内气体为不可压缩流体，且满足Boussinesq假设：认为流体密度的变化仅对浮升力产生影响。

（5）假定流场具有高的紊流Re数，流体的紊流粘性具有各向同性。 （6）气流为低速流动，可忽略由流体粘性力做功所引起的耗散热。 （7）假设一点的温度只与一个机柜相关。

（8）假设温度与机柜任务量满足指数函数关系。 （9）假设温度最高点只发生在温度极值点上

四、符号说明

ki?i?x，y，z?：机房内任意一点的扩散系数 C?x，y，z，t?：机房内任意一点热空气的浓度 Ux ，Uy，Uz：为x，y，z方向风速 Kx，Ky，Kz：为x，y，z方向上的扩散系数 p0：单位时间内机柜散发热空气的速率 Q0：从机柜散发出的热流总量 Q1：在?t，t??t?内通过某空间域的热流量 Q2：在某空间域热流量的增量 ?x，?y，?z：分别为用浓度标准差表示的x，y，z轴上热量的扩散参数 vs：热空气与冷空气接触的系数； S：机柜散发热的扩散的系数 K、t：分别为风速与任意扩散时刻 H：机柜散发出的热源距地面的高度 G：为重力加速度； K：为流体紊流动能； P：为时均压力； C：为热源强度； Pr：为充分紊流时的普朗特数； T：为流体温度； T0：为参考温度； Ui：为速度分量，当i=1、2、3时分别代表X、Y、Z方向；

3

ε：为紊流能量耗散率； ?、?：分别为层流和紊流的动力粘性系数； ρ：为流体密度； β：为流体体积膨胀系数。 五、模型建立及求解

5.1模型一

经过对附件1的数据和题目的分析，直接使用MATLAB软件将该机房的二维热分布及二维的流场分布图绘制如下图 1 由附件1得到的通道二、三热分布图和流场分布图：

图1.1 通道二热分布图 图1.2 通道二流场分布图

图1.4 通道三热分布图 图1.4 通道三流场分布图

图 1 由附件1得到的通道二、三热分布图和流场分布图

同时，通过MATLAB求出极值（MATLAB计算程序见附录Ⅰ）： 通道二：

Imregionalmax = 41.3707 Points =

3.3697 2.7000 通道三：

Imregionalmax = 47.8571

4

Points =

3.3697 2.4091

5.2模型二

机房的热分布模型

记热量开始扩散的时刻t?0，并且以热源作为扩散中心，则可将t时刻机房内任意一点（x，y，z）的热空气浓度记为C?x，y，z，t?。假设单位时间内通过单位法向量面积的热流量与浓度梯度成正比，则有

（1）

q??ki?gradC其中k?i?x，y，z?是扩散系数，grad表示浓度梯度。

i假设机房空间?的体积为V，包括机房空间?的曲面为一规则的球面，设其表面面积S表，外法线向量为

?xy?，n???，?，1?zz??则在?t，t??t?内流通过空间域?的热流量可以表示为：

Q1??t??tt??q?nd?dts （2）

在空间?所包围的区域内空气中的热流量增量可表示为：

Q2?????C?x，y，z，t??t?-C?x，y，z，t??dVV （3）

由机柜散发的的总热流量为：

Q0??t??tt??pdVdt0? （4）

由质量守恒可得出：Q?Q?Q，即，

012????C?x,y,z,t??t?-C?x,y,z,t??dV??Vt??tt??q?nd?dt??st??tt???pdVdt0? （5）

根据曲面积分Gauss公式得：

??q?nd?????divqdVsV （6）

则式子（5）可以转换成

?C?x，y，z，t??t?-C?x，y，z，t????tdV??????t?V?

5

??t??tt???divqdVdt??Vt??tt???pdVdt0?t??t （7）

由于：

?CC?x，y，z，t??t?-C?x，y，z，t???lim?limt?t?0?t?t?0?tkdiv?gradC?dt

?t（8）

故式子（7）即可转换成：

??C?dV??t????divqdV??t????p0dV??t??????t?V?V?化简得：

（9）

??C?dV????divqdV?p0??????t?V?V （10）

根据A.Fick扩散微分方程式中：

?C?C?C?C?2V?2V?2V

?Ux?Uy?Uz?KX?KY?KZ22?t?X?Y?Z?X?Y?Z2（11）

其中：C为常温气体浓度；t为时间；Ux ，Uy，Uz为x，y，z方向风速；

Kx，Ky，Kz为x，y，z方向上的扩散系数。

假设机房热量是在无风条件下扩散，此时有，

?C?2V?2V?2V

?Kx?Ky?Kz22?t?X?Y?Z2（12）

结合式子（11）可解出（10）的结果为：

（13）

?x?p0??yz-HC?x，y，z，t??exp---??4Kt4Kt4Kt?4??t?1.5?KxKyKz?0.5?xyz???

由此得到机房里通过机柜散发热后的空气浓度：

C1?x，y，z，t??vsC?x，y，z，t?

?x2?y2?z2?s2t2?

（14）

由（13）可知机柜的热分布是一个向四周散发且在冷通道里的温度是降低的。

进而可知在测试案例中机房的热分布呈一个向四周发散的状态，特别是在冷通道的温度明显的可以看出是在降低（相对于其他通道）。

5.3模型三

机房设计是一个对称结构，也就是机柜位置对称，空调位置对称，冷热通道对称。所以我们选取任务量相同的温度组（即附件二中第E、M、P列）通过MATLAB作图并寻求温度极值点（见图 2 机柜不同任务量不同高度的热分布图）：

6

图2.1.1 高度为0.2 图2.1.2 高度为1

图2.1.3 高度为1.8

图2.1.4 高度为2.6

图2.1 第E列即每个机柜任务量为0.5

图2.2.1 高度为0.2 图2.2.2 高度为1

7

图2.2.3 高度为1.8

图2.2.4 高度为2.6

图2.2 第M列即每个机柜任务量为0.2

图2.3.1 高度为0.2 图2.3.2 高度为1

图2.3.3 高度为1.8

第P列即每个机柜任务量为0.3

图2.3.4 高度为2.6

图 2 机柜不同任务量不同高度的热分布图

通过MATLAB计算出温度最高的几个极值点（MATLAB计算程序见附录

8

Ⅱ）：

每个机柜的任务量为0.5： Imregionalmax = 32.7584 32.8246 32.7247 Points =

2.4162 4.3727 4.0889 4.3727 7.2485 4.3222 每个机柜任务量为0.2： Imregionalmax = 28.6021 28.6308 28.5428 Points =

2.4162 4.3727 4.1818 4.3727 7.2485 4.3222 每个机柜任务量为0.3： Imregionalmax = 30.6982 30.7181 30.6863 Points =

2.4162 4.3727 4.1818 4.3727

7.2485 4.3222

根据数据可以发现三个任务量下极值点相同，由此，可以得出极值点即为?2.4162,4.3727?

?4.0889,4.3727??7.2485? ,4.3222但由于高度上的分布是离散的，所以必须根据不同位置做出温度分布图，这样可以减小误差。因为得到的极值点x值即通道位置就为给出数据中通道二、通道三、通道四的数据，所以，仅需要做出这三个通道的温度分布，见图 3 不同通道的热分布图：

图3.1 通道二 图3.2 通道三 图3.3 通道四

9

图 3 不同通道的热分布图

通过MATLAB计算出温度最高的几个极值点（MATLAB计算程序见附录Ⅲ）：

每个机柜的任务量为0.5： Imregionalmax = 32.5008 Points =

4.2717 1.7758

每个机柜的任务量为0.2： Imregionalmax = 32.4353 Points =

4.2717 1.7758

每个机柜的任务量为0.3： Imregionalmax = 32.5858 Points =

4.2717 1.7758

根据数据可以发现三个任务量下极值点相同，由此，可以得出极值点即为

?2.4162,4.3?。由于我们没有温度极值点的数据，所以通过点727,1.7758?4.0889,4.3?727,1.7758?7.2485 ,4.3222,1.7758??2.4,4.1,1.8??4.1,4.1,1.8?进行模拟，从而得到不同任务量下温度极值点的温度分布表（表 1 ?7.2 ,4.1,1.8?机柜不同任务量下温度极值点的温度分布表）：

表 1 机柜不同任务量下温度极值点的温度分布表

（注其中x、y、z分别为通道位置，距离空调位置，高度）

机柜一 0.5 0.8 0.5 0.5 0.5 0.3 0.5 0.5 0.2 0.3 0.3

机柜二 0.5 0.5 0.8 0.5 0.5 0.5 0.5 0.2 0.2 0.2 0.2 机柜三 0.5 0.5 0.5 0.8 0.5 0.5 0.2 0.5 0.2 0.2 0.2 机柜四 0.5 0.5 0.5 0.5 0.8 0.5 0.5 0.5 0.2 0.2 0.3 x y z 2.4 4.1 1.8 32.1 35.9 32.4 32.2 32.1 30.2 28.4 31.9 28.1 30.0 30.1 4.1 4.1 1.8 32.0 32.0 35.9 32.0 32.0 32.0 32.0 28.1 28.0 28.0 30.0 7.2 4.1 1.8 32.2 32.2 32.2 32.6 35.8 32.2 32.2 32.2 28.1 28.1 30.2 10

0.3 0.3 0.3 0.3 30.1 30.0 30.2 分析表格可知当机柜一工作量提高到0.8时点二，点三温度没有变化，同样当机柜二工作量提高到0.8时点一，点四温度没有变化，当机柜四工作量提高到0.8时点一，点二温度没有变化，而当机柜三工作量提高到0.8时所有点的温度都几乎没有变化。同样，当温度降低时也是如此。根据上述分析，我们做出假设：点一的最高温度只与机柜一有关，点二的最高温度只与机柜二有关，点三的最高温度只与机柜四有关。做出散点图（图 4 温度极值点与相关机柜任务量的拟合），由图像确定使用指数函数拟合。

40383634

3230282600.10.20.30.40.50.60.70.80.918? 图4.1 点?2.4,4.1,1.

11

40383634

3230282600.10.20.30.40.50.60.70.80.918? 图4.2 点?4.1,4.1,1.40383634

3230282600.10.20.30.40.50.60.70.80.918? 图4.3 点?7.2,4.1,1.图 4 温度极值点与相关机柜任务量的拟合

分析数据，选取温度极值点对应的机柜号拟合得到方程（MATLAB计算程

12

序见附录Ⅳ）：

T173 4268.3?1.T2;0368 5206.1.T3;426.?x1??x2??856 4760.1x4

其中Ti为第i个温度极值点的温度

xi为第i个相关机柜的任务量

运用MATLAB分配任务量，并求出机房最高温度的最低值（MATLAB计算程序见附录Ⅴ）：

任务量为0.8时： maxint =35.1919 x11 =0.7400 x12 =0.7300 x13 =1.0000 x14 =0.7300

任务量为0.5时： maxint =30.1168 x11 =0.3200 x12 =0.3500 x13 =1.0000 x14 =0.3300

5.4模型四

k~?两方程紊流模型建立数学模型

对于不可压缩、定常的流动，根据以上的假设，该模型的控制方程可用如下的通用形式：

(1) 连续性方程：

??Ui?Uj?????UiUj??????i??i?????X?Xi?Xi??Xij???????UT??(2) 能量方程：i?Xi?Xi???u????v????w????0

?X1?X2?X3（15）

???????T0?T??g???? （16）

???i?i???Pr??T????T?q???X??CP ?i?（17）

(3) 紊流动能方程(k方程)：

???Uik?????Xi?Xi??Ui?Uj??i?k???????X???i??X??Xi?ji?K????Uj??C1?i??Ui????X?Xij????Ui??T??????gi（18） ??XPr?Xji? (4) 紊流动能耗散率方程(ε方程)：

???Ui??????Xi?Xj??i?k????Xj??2??Ui??T???C2?C3?gi(19) ??XkPr?Xij?紊流粘性系数的表达式：

13

?i?C?CD?k2? （20）

5.4.1计算结果及分析

点三的相关机柜为机柜二，根据题三拟合的公式：。由于点（7.2,4.1,1.8）就在机柜和冷通道交界处。根据题设，冷通道的出风口为通道1/3。我们假设风速恒定，那么

，又因为

由此可得

即

。根据上面所设立的

模型，可以计算出机柜中的温度是……根据附件2中的数据，可以得到任务量和

温度之间的散点图。把散点图和上面个所得的公式进行比较，算出R2值。同理，可以得到点2的散点图和拟合图。如图（matlab程序见附录Ⅵ）：

图 5 点2散点图和拟合线

14

图 6 点3散点图和拟合线

根据图可知，假设符合实际，温度控制时应该满足

六、模型的检验

1. 本文中大多数图形都用MATLAB 软件绘制，而MATLAB软件具有较强的仿真性，功能也十分强大，精确度也很高，从这个角度可证明我们结果的可靠性与方法的合理性。

2.对于机房热分布采用模型二，将机柜散发热看成一个热源可以逼真的模拟出机柜散发热的传播途径，使得机房的热分布描绘形象逼真。

3.充分考虑了机柜的对称性，模型三分析不同任务量的条件下，选取最高温度来确定机柜的任务量，从而来控制机房整体温度，这样更能够确保到机房的安全性。

4. 从模型与实际情况的对比来看，虽然我们建立模型有一些假设，但这些假设都是合乎常理的，故我们所建的模型的可靠性与真实性也是值得认可的。

七、模型的推广及发展前景

本文第一个模型根据题干给出的数据，采用MATLAB软件编程绘制的热分布及流场分布图可以应用于各个领域，如学校可以推广到消息的传播、文化的传播、疾病的传播、或者流言蜚语的传播，从而对现今这个信息时代做出一个正确的把握。同样的方法可以运用到气流运动、地质勘探和测量等方面。

本文只是阐述创建绿色机房一个因素——温度，其实还可以把它推广到温度对机房设备的寿命的影响，以及将机房机柜产生的热量转换成其他有效能源，从而有利益节约资源，更好的创造绿色机房，绿色社会。

此模型符合当今社会的发展状况，目前，人类都在关注者社会的有限资源，创建友好型、环保型的绿色社会，应用本文的模型及分析方法，通过对机房的温度叙述，可以避免因温度过高而引发的不良事故，同时也为创建绿色机房提供了可靠的理论和实践依据。绿色数据中心的建立，实现了最佳数据中心基础设施能源利用效率（PUE）和动态智能制冷，精确了送配风系统，同时采用区域化和模块化设计－－高热区和低热区不同的散热方式，实现对不同负载的有效支持。对

15

大型数据中心进行模块化设计理念整合的智能的机房监控系统 (动力设施,环境与IT设施,平台统一)实现自动化管理。

八、模型的优缺点

一、优点：

1.本文中模型的建立充分考虑到实际情况，从热源出发，推广到全机房，从一定程度上减小了误差。

2.在分析问题的时候，减少了一些次要因素的影响，虽在一定程度上降低了精度，但是简化了过程，增加了可读性。

3.本题的解答过程中，充分使用数学软件，根据图像，求出极值点，更加准确

二、缺点：

1.思考问题时过于简单，忽略了次要因素的影响，可能导致精度不高。 2.模型的假设，未能定量分析，只能做到定性分析

16

九、参考文献

（1）Robinson D F，Hassan Hassan H A. Further Development of the

(Enstrophy) Turbulence Closure Model. Journal of AIAA，1998，36：1825~1833

（2）Pratap，V.S.and Spalding，D.B，Numerical Computions of the Flow in

Curved Ducts, The Aeronautical Quarterly, Aug.(1975)

（3）Eiseman，P.Retal，Development of a Three-Dimensional Turbulent Duct

Flow Analysis，NASA CR3029，（1978）

（4）Chen Y S，Kim S W. Com putation of Turbulent Flows Using an Extended

K-? Turbulent Closure Model. ASA - C R-179204，1987

（5）Tumber M G，Jennions IK. An Investigation of Turbulent Modeling in

Transonic Fan Including a Novel Implementation of an Implicit K-? Turbulence Model. ASME92-GT-308 （6）姜启源、谢金星、叶俊，《数学模型》（第三版），高等教育出版社 （7）尹贞勤，程控交换机房空调设计探讨[ J]，安徽建筑，1999( 4)

（8）杨世铭，自由运动放热的基本准则式[ J]，西安交通大学学报，1962( 2)

17

十、附录

附录Ⅰ：

clear%根据附件一做出的通道二热分布图 clc

A=[2.4 0.3 13;2.4 0.9 13;2.4 1.5 17;2.4 2.1 30;2.4 2.7 30;5 0.3 13;5 0.9 13;5 1.5 25; 5 2.1 30;5 2.7 30;7.2 0.3 13;7.2 0.9 13;7.2 1.5 19;7.2 2.1 30;7.2 2.7 30;];

x=A(:,1);y=A(:,2);z=A(:,3);

[X,Y,Z]=griddata(x,y,z,linspace(2.4,7.2)',linspace(0.3,2.7),'v4');%插值 pcolor(X,Y,Z);shading interp%伪彩色图 figure,surf(X,Y,Z)%三维曲面 ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)]

clear%根据附件一做出的通道三热分布图 clc

A=[2.4 0.3 27; 2.4 0.9 29; 2.4 1.5 29; 2.4 2.1 30; 2.4 2.7 29; 5 0.3 30; 5 0.9 30; 5 1.5 29; 5 2.1 31; 5 2.7 32; 7.2 0.3 27; 7.2 0.9 31; 7.2 1.5 31; 7.2 2.1 32; 7.2 2.7 31;];

x=A(:,1);y=A(:,2);z=A(:,3);

[X,Y,Z]=griddata(x,y,z,linspace(2.4,7.2)',linspace(0.3,2.7),'v4');%插值 pcolor(X,Y,Z);shading interp%伪彩色图 figure,surf(X,Y,Z)%三维曲面 ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)]

clear%根据附件一做出的通道二流场分布图 clc

A=[2.4 0.3 0.6;

18

2.4 0.9 0.6; 2.4 1.5 0.9; 2.4 2.1 1.1; 2.4 2.7 1.1; 5 0.3 0.4; 5 0.9 0.4; 5 1.5 0.5; 5 2.1 0.6; 5 2.7 0.6; 7.2 0.3 0.4; 7.2 0.9 0.2; 7.2 1.5 0.2; 7.2 2.1 0.2; 7.2 2.7 0.2;];

x=A(:,1);y=A(:,2);z=A(:,3); scatter(x,y,5,z)%散点图 figure

[X,Y,Z]=griddata(x,y,z,linspace(2.4,7.2)',linspace(0.3,2.7),'v4');%插值 pcolor(X,Y,Z);shading interp%伪彩色图 figure,contourf(X,Y,Z) %等高线图 figure,surf(X,Y,Z)%三维曲面 ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)]

clear%根据附件一做出的通道三流场分布图 clc

A=[2.4 0.3 0.4; 2.4 0.9 0.6; 2.4 1.5 0.7; 2.4 2.1 0.8; 2.4 2.7 0.9; 5 0.3 0.4; 5 0.9 0.5; 5 1.5 0.6; 5 2.1 0.7; 5 2.7 0.6; 7.2 0.3 0.4; 7.2 0.9 0.6; 7.2 1.5 0.6; 7.2 2.1 0.6; 7.2 2.7 0.5;];

x=A(:,1);y=A(:,2);z=A(:,3); scatter(x,y,5,z)%散点图

19

figure

[X,Y,Z]=griddata(x,y,z,linspace(2.4,7.2)',linspace(0.3,2.7),'v4');%插值 pcolor(X,Y,Z);shading interp%伪彩色图 figure,contourf(X,Y,Z) %等高线图 figure,surf(X,Y,Z)%三维曲面 ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] 附录Ⅱ：

clf%做出每个机柜任务量为0.5时每个高度温度分布图 clear

A=load('data41.txt');%%读入矩阵 for n=4:1:7

x=A(:,1);y=A(:,2);z=A(:,n);

[X,Y,Z]=griddata(x,y,z,linspace(0,9.2)',linspace(2.1,7.1),'v4');%插值 pcolor(X,Y,Z);

shading interp%伪彩色图 figure,surf(X,Y,Z) figure

ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] end

clf%做出每个机柜任务量为0.2时每个高度温度分布图 clear

A=load('data49.txt');%%读入矩阵 for n=4:1:7

x=A(:,1);y=A(:,2);z=A(:,n);

[X,Y,Z]=griddata(x,y,z,linspace(0,9.2)',linspace(2.1,7.1),'v4');%插值 pcolor(X,Y,Z);

shading interp%伪彩色图 figure,surf(X,Y,Z) figure

ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] end

clf%做出每个机柜任务量为0.3时每个高度温度分布图 clear

A=load('data412.txt');%%读入矩阵 for n=4:1:7

x=A(:,1);y=A(:,2);z=A(:,n);

20

本文来源：https://www.bwwdw.com/article/2m5.html