[曲面拟合高大平房仓20仓的粮温情况.doc

应用曲面拟合粮温的方法分析研究 高大平房仓20仓的粮温情况 小草 摘要：应用曲面拟合粮温的方法分析研究高大平房仓20仓的粮温情况，分析研究了粮温与平面几何空间，与空间、时间，与空间、时间和气温的曲面拟合，得出了曲面拟合粮温函数，作了等温线图和细分各层粮温图，并分析研究。 关键词：曲面拟合粮温 粮温与空间、时间、气温的曲面拟合函数 等温线图 细分各层粮温图 分析研究 前 言 应用曲面拟合粮温的方法分析研究粮温情况是一个新方法。我应用该方法对高大平房仓20仓这个实例进行分析研究，以便打开分析研究粮温的新方法做点有益的探索。 主要对高大平房仓20仓的宽度方向上的垂直面、长度上的垂直面和水平面的第3层的粮温T（X，Y）进行了曲面拟合，并作了其等温线图，另外对测温线28号的粮温（中央点）进行了曲面拟合粮温T（X，t时间）和T（X，t气温，t时间），并作了其等温线图和细分各层粮温变化图，而后进行了分析讨论。 处理数据主要方法：曲面拟合法； 主要工具：电子表格中的函数LINEST（）和曲面俯视图功能。 计算结果 1.1、20仓的宽度方向上的垂直面的粮温T（X，Y）曲面拟合函数和其等温线图： 其基函数为双3次多项式：XiYj (i=0,1,2,3 j=0,1,2,3) 曲面拟合粮温: i j T（X，Y）=∑∑Aij XiYj i=0 j=0 表1：高大平房仓20仓第6列垂直剖面拟合粮温系数Aij j 0 1 2 3 i 0 30.4257 0.393892340 -0.009789360 -0.000342456 1 -6.695658263 -0.368433605 -0.047225885 0.003684763 2 0.474295601 -0.065863941 0.039947036 -0.002029261 3 0.048944914 0.023581470 -0.006057786 0.000267199 复相关系数R^2=0.9985，剩余标准差Se=0.7127，自由度f1=15，f2=4，F值=141.7724。 1.2、20仓的长度方向上的垂直面的粮温T（X，Y）曲面拟合函数和其等温线图： 其基函数为双3次多项式：XiYj (i=0,1,2,3 j=0,1,2,3) 曲面拟合粮温: 3 3 T（X，Y）=∑∑Aij XiYj i=0 j=0 表2：高大平房仓20仓第3行垂直剖面拟合粮温系数Aij j 0 1 2 3 i 0 32.3818 0-0.00498169248 0.00003902579 1 -3.87261757850 -1.35040043695 0.06804188530 -0.00099131542 2 -0.77503811068 0.57715134007 -0.03444309649 0.00056406825 3 0-0.06576487342 0.00430606564 -0.00007413815 复相关系数R^2=0.9756，剩余标准差Se=1.4960，自由度f1=15，f2=28，F值=74.8012。 1.3、20仓的水平面的第3层粮温T（X，Y）曲面拟合函数和其等温线图： 其基函数为对称园锥面二次曲面函数： Uk(X,Y)=Q（X，Y，Xj ,Yj）=((X-Xj)^2+(Y-Yj)^2)^0.5 式中：Xj ,Yj为基函数的中心点坐标。 曲面拟合粮温: 了 m m T（X，Y）=∑Aj Q（X，Y，Xj ,Yj ）+A0= ∑Aj((X-Xj)^2+(Y-Yj)^2)^0.5 +A0 j=1 j=1 表3：高大平房仓20仓第3层水平剖面拟合粮温的中心点方位和系数AjA0表 中心点坐标 （0，0） (0，9.08) （0，18.16） （0，27.24） （0，36.32） （0，45.4） 对应系数 1.2749 -1.4030 -0.4562 -1.0686 -0.8403 0.5214 中心点坐标 （4.6，0） (4.6，4.54) (4.6，13.62) (4.6，22.7) (4.6，31.78) (