最新matlab绘制温度场

通过在室内的某些位置布置适当的节点，采集回来室内的温湿度以及空气质量等实际参数。

首先对室内空间建模，用一个无限细化的三维矩阵来模拟出室内的温度分布情况，针对采集回来的数据，采用插值法和适当次数的拟合函数的拟合，得出三维矩阵的实际值的分布，最后结合matlab软件绘制出计算出的温度场的三维图像。

一．数据的采集与处理
因为影响人的舒适感的温度层只是室内的某一高度范围内的温度，而温度传感器虽然是布置在一个平面内，但是采用插值法和拟合函数法是可以大致再现出影响人的舒适感的温度层的温度变化的。

同时，在构建出的三维模型中，用第三维表示传感器层面的温度。

在传感器层面，传感器分布矩阵如下：
X=【7.5 36.5 65.5】（模型内单位为cm）
Y=【5.5 32.5 59.5】
Z=【z1 z2 z3；
z4 z5 z6；
z7 z8 z9；】（传感器采集到的实时参数）
采用meshgrid（xi，yi，zi，…）产生网格矩阵；
首先按照人的最小温度分辨值，将室内的分布矩阵按照同样的比例细化，均分，使取值点在坐标一定程度上也是接近于连续变化的，从而才能最大程度上使处理数据得来的分布值按最小分辨值连续变化！
根据人体散热量计算公式：C=hc（tb-Ta）
其中hc为对流交换系数；
结合Gagge教授提出的TSENS热感觉指标可以计算出不同环境下人的对环境温度变化时人体温度感知分辨率，作为插值法的一个参考量，能使绘制出的温度场更加的符合人体的温度变化模式。

例如按照10cm的均差产生网格矩阵（实际上人对温度的分辨率是远远10cm大于这个值的，但是那样产生的网格矩阵也是异常庞大的，例如以0.5cm为例，那么就可以获得116*108=12528个元素，为方便说明现已10cm为例）：
[xi yi]=meshgrid(7.5:10:65.5,5.5:10:59.5)
xi =
7.5000 17.5000 27.5000 37.5000 47.5000 57.5000
7.5000 17.5000 27.5000 37.5000 47.5000 57.5000
7.5000 17.5000 27.5000 37.5000 47.5000 57.5000
7.5000 17.5000 27.5000 37.5000 47.5000 57.5000
7.5000 17.5000 27.5000 37.5000 47.5000 57.5000
7.5000 17.5000 27.5000 37.5000 47.5000 57.5000
yi =
5.5000 5.5000 5.5000 5.5000 5.5000 5.5000
15.5000 15.5000 15.5000 15.5000 15.5000 15.5000
25.5000 25.5000 25.5000 25.5000 25.5000 25.5000
35.5000 35.5000 35.5000 35.5000 35.5000 35.5000
45.5000 45.5000 45.5000 45.5000 45.5000 45.5000
55.5000 55.5000 55.5000 55.5000 55.5000 55.5000
产生网格矩阵之后，就可以在测得的实时数据的基础上，通过相关的温度场的专业的估算函数，以及相关的数值处理函数来估计整个分布面（有最小的分辨率）上的温度了。

即在这些函数的基础之上，对参数进行一些必要的处理。

然后，加上室内的预设温度，扩大温度数值矩阵，例如预设温度为21℃，采用一组试验中测得的温度值，如下：
z =
21 21 21 21 21
21 19 18 17 21
21 18 16 19 21
21 16 15 14 21
21 21 21 21 21
此时再结合matlab中定义的spline（三次样条差值法）和interp2（二维插值法）这两个函数来估算预定分辨率上的温度场的分布。

因为在实际中温度的变化必然是连续的，故而得到的温度场中通过任意点的截面截出的曲线必然是连续可导的，故而必须使用spline（三次样条插值法）得到光滑的插值分布曲线。

同时又因为该温度场现实的是传感器分布面上的温度分布，故而必须采用interp2（二维插值法），同时用第三维的高度值的变化和连续变化的颜色来显示温度的连续变化。

Zi=interp2（x，y，z，xi，yi，‘spline’）
zi =
19.0000 4.5173 15.2048 16.7822 -6.7326 -16.9509
7.1027 -91.1872 -15.4389 -11.6435 -191.3871 -266.8890
18.9402 6.1945 18.1887 11.8409 -29.5509 -44.3035
11.4135 -68.5455 -25.7034 8.1347 -40.1338 -69.1854
-19.6972 -344.0942 -148.6804 -37.1340 -325.8118 -477.3336
-12.7095 -278.7121 -114.6136 -28.6166 -283.7157 -414.3010 由此得出网格矩阵的元素值之后，采用matlab自定义的三维曲面绘制函数，绘制出得到的温度场。

Surf（xi，yi，zi）
经过一定的图像处理之后可得到温度场的分布如下：
从图中坐标也可以看出，分辨率设置的过大之后，插值法会出现较大的误差甚至是错误。

然而当分辨率设置的较为合理之后，二维插值结合三次线条插值能较好地吻合实际的温度场
的分布。

在上图中可以看出温度场变化很突兀，设置较高的人体温度感知分辨率之后图像如下：
下图为温度场分布曲面在传感器分布面上的投影的等高线图，可以比较直观的反映室内
的温度变化情况。

【c，h】=contour（x，y，z）；
Clabel（c，h）；
参考文献：
敢于正视压力
压力很可怕，是个坏东西。

然而，人们又说，井无压力不喷油，人无压力轻飘飘。

这样看来，压力又不可少。

我们应该怎样认识和对待压力呢？
生在世间，存于宇宙。

同顶炎炎烈日，共沐皎皎月辉。

压力对每一人来说是避免不了有的，只是或小或大，或轻或重罢了。

作为工人有失业待岗的压力，作为干部有成绩大小而影响升降的压力，作为我
们学生则有升学的压力……种种压力赋于人们的头上。

怎样看待这一个个压力呢？我要说每个人都要敢于正视压力。

为人做事没有压力是万万不行的。

俗语说，井无压力不喷油，人无压力轻飘飘。

当今世界五百强之一的中国海尔集团，多年前乃一家濒临的企业，面对外债累累，一个将近衰败的企业，张瑞敏走马上任，可以想象那时的张瑞敏心理上没有压力的话，中国就不可能有现在的海尔集团名誉中国大地和世界各国。

是压力给海尔集团带来了经济复燃，是压力给海尔集团带来了生机勃勃。

这就是敢于正视压力，把沉重的压力化做前进的动力，乘风破浪在高浪大风的海上，搏击高浪巨风，巍然屹立于经济技场中。

史铁生年纪青青就残了双腿，曾经无奈过，曾经失落过，也曾一个人摇着轮椅在“那座废墟的荒圆”中虚度。

做人的压力如同大青石般的压在心口，压的史铁生喘不过气来。

也曾因为“双腿残废，无法做人”的压力，他对生命不报太大的希望。

然而，后来史铁生认识到：“死是一件无需着急的事，是一件无论如何耽搁也避免不了的事。

”自此之后，史铁生用坚强的毅力去正视这已成为现实的人生压力，给自己折了一只名为“写作”号的船，将自己从压力的深渊里摆渡出来，终于实现了自己的人生美丽。

这也正是史铁生敢于正视压力，坚强的从压力的深潭里挣脱出自己解救了自己，让自己驾着“写作”号去实现自己的人生价值。

作为年轻人也要敢于正视自己的压力，化压力为动力，为自己也折只能摆渡压力这个深渊的船。

正视自己的压力，努力，尽力，便是
美丽。