数学建模论文-储油罐的变位识别与罐容表标定1.doc

题目： 储油罐的变位识别与罐容表标定 摘 要本文根据题目要求，针对具体问题建立相应的数学模型，并对问题进行精确求解，逐层深入比较题给数据与理论求解值同时分析它们之间的关系，据此得出合理的结论。对于问题一,首先，采用截面法建立数学模型对无变位及有变位的理论油容量与油位高度之间的关系表达式。然后，根据附件1所给的无变位进油时油高求出对应理论油量并与所给实际油量值相比较，发现两者的相对偏差基本成一稳定值3.371%，并根据两者之间的关系求出理论油量的修正函数，同时将偏差函数运用于无变位出油中，通过比较相邻两时刻绝对出油量的相对误差（均值为）检验修正函数的正确性。接着，根据对应倾斜变位进油油高的校正理论油量与实际油量相比较，分析得出纵向变位时液位较小时罐容表读数偏小，液位较高时读数偏大的结论。最后运用matlab编写程序对变位后罐容按油位高度间隔为1cm的进行标定，结果见附表1。对于问题二，由于横向倾斜不影响体积的计算，只影响油面高度实际值与罐容表测量值之间的关系。于是，先根据三重积分与截面法计算只有纵向倾斜的情况下罐容与倾斜角度a和的关系。然后，根据油面高度与横向偏转角度b及关系得罐内储油量与油位高度及纵向倾斜角度a和横向偏转角度b之间的一般关系。取（无变位时），根据附件2中油位高度计算出相应油量，发现与附件所给油量容积几乎完全一致。接着利用附件2中实际出油量与理论出油量的差值的平方和最小的条件，运用遍历法和ndepso算法两种方法分别求得结果如下：（度）（度）误差遍历法2.14.44.9e-5ndepso算法2.1294.2443.7e-6根据误差最小原则选取，求得罐容表标定值部分数据如下（其余数据见附表4）：h(m) 00.10.20.30.40.50.6v(l)45.88355.421068.682228.223709.18 5440.16 7379.69利用第二次注油后的数据，带入已求出的、，和表中的h求出实际出油量与理论出油量近似相等，验证模型的正确性和可靠性。并进行了模型推广，进行了基于matlab的gui设计。关键词：截面法，三重积分，遍历法，ndepso算法，gui设计一 问题的重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为a=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。二 模型的假设1、假设储油罐纵向偏斜时偏斜角较小，符合地基变化的实际情况。2、假设忽略储油罐偏斜时浮子测量所带来的误差。3、假设忽略出油与进油速度对油位高度的影响4、假设图1中的储油罐尺寸上不存在工艺误差三 符号约定 -截椭圆面面积 -截椭圆柱体体积-纵向倾斜角度-截圆面面积-球冠高度-左侧球冠处油体积-中间圆柱体处油体积-右侧球冠处油体积-横向偏转角度-整个球冠体积四 模型的分析与建立4.1对问题一的分析与模型建立4.1.1问题一的分析为了掌握罐体变位后对罐容表的影响，首先对无变位的情况进行分析，通过建立无变位体积求解模型，求得理论值。与附件一中实验测得值进行比较，分析得出两者之间的关系，然后据此对理论值进行修正，同利用无变位出油数据对其进行检验。截椭圆柱体体积与油位高度之间的关系：采用截面法建立数学模型对无变位及有变位的理论油容量进行计算，得出了两端平头的椭圆柱体油罐的罐容与油位高度之间的关系。另外，以所得的偏差为依据对所计算有变位时的理论油容量进行校正，再将其与试验测得数据进行比较，分析罐体变位后对罐容表的影响并对给出罐体变位后油位高度间隔为1cm的罐容表标定值。最后，运用曲线拟合的方法来形象直观地说明理论值、校正值以及实验测试值之间的关系。4.1.2问题一的模型建立4.1.2.1椭圆截面面积公式的推导：通过对问题一的分析，我们建立求解截椭圆面面积的模型：设椭圆的长短轴长分别为、；截面的油位高度为h。图1 截面椭圆示意图设椭圆方程为：；被截椭圆面的液面高为h,其中，图中带阴影部分为储油横截面，先用微元法求解微元面积：；则：采用换元法，令，则与积分区间变换关系如下表：表1 积分区间变化关系则：积分得：带上下限可得 ：既得截面面积与截面油高之间的关系。4.1.2.2储油罐无变位时体积公式推导：储油罐无变位时，当油位高度一定时，储油罐的横截面是一个面积恒定的截椭圆，故此时的体积即为等截面柱体体积的求解。设柱体高为l，其中此时的横截面图形如图1所示：可得无变位时油罐储油量与油位高度之间的关系，如下式所示： 4.1.2.3储油罐纵向变位时体积公式推导：当储油罐纵向变位时，其中储油横截面为变面积截面，故考虑使用截面法采用积分来对体积进行求解，具体推导过程如下图所示：图2 平头储油罐纵切面示意图其中h为油位高度，h为对应每个时所截截面的油面高度；根据上述图形所建立的坐标系，沿y轴对截面s进行积分既得体积，即：其中根据的积分限将积分区间分成三段进行积分，如下所示：当油位高度（即，）时，罐内油可看成一个截面面积一端为0的截椭圆柱体；积分下限，积分上限当油位高度（即，）时。 积分下限，积分上限当油位高度（即，）时，储油罐内的油体积可以看做整个椭圆柱体的体积减去上面一部分一个一端为0的小截椭圆柱体所得的部分；积分下限，积分上限 综上所述,建立椭圆油罐体油容积与油位高度h 的关系式为：4.2问题二的分析与模型建立4.2.1问题二的分析对于问题二，它既有纵向的倾斜又有横向偏转，首先，考虑只有纵向倾斜的情况下得出罐容与倾斜角度a和油位高度的关系。然后，由于在纵向倾斜角一定的情况下，横向偏斜不影响油体积的求取，其实际油位高度与横向偏斜与否无关，根据实际油位高度与油浮子测得的高度及其横向偏转角度b之间的关系，既可得，罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。通过分析附件2的数据我们看出给出的显示油高和显示容积应该为无偏向时的罐容表，我们就可以利用先前得到的一般关系式，令两偏斜角度为零时，即无偏斜的情况下，计算相应油高时的理论容积，再与附件2中的显示容积比较，得出两者偏差，并可达到间接检验设立的模型的正确性的目的。再根据模型，将一定的偏斜角度带入，根据对应的油高，即可得对应的容积，于是我们得到两对应高度下之差的理论容积之差，再与附件2的对应的高度下的实际出油量（容积之差）相比，若两偏斜角度选取合适，则相比之差应很小，于是我们可以采用某种算法，寻找这样的两偏斜角度，使得所有的理论容积之差与实际出油量偏差很小，故我们可以计算多个这样的偏差平方和，将问题转化为求解使得该偏差平方和最小的两偏斜角度。再通过这两偏斜角度带入附件2中数据，如上所示的求解所谓的平方和，即这两偏斜角度对应下的总偏差，检验其值，即可检验模型的正确性和可靠性。4.2.2问题二的模型建立4.2.2.1圆截面面积公式的推导：通过对问题二的分析，我们建立求解截圆面面积的模型：设圆半径为r；截面的油位高度为r。图3 截圆面示意图设圆方程为：根据椭圆与圆之间的对称关系根据椭圆截面面积与油位高度h之间的关系：将椭圆的长短轴、均用圆的半径进行替换，可得圆截面面积与油位高度之间的关系：4.2.2.2体积公式的推导：现要对一个主体为圆柱体，两端为球冠体的储油罐储油体积进行计算。主体圆柱体的半径为r；圆柱体高为l；两端球冠高为；球冠对应球体半径为r；油浮子距圆柱体两端的距离分别为、；罐体纵向倾斜角度为。此时油容积由三部分组成，分别为左侧球冠处油体积，中间圆柱体处油体积，右侧球冠处油体积。现分别对这三部分体积公式进行推导。4.2.2.3罐体纵向变位时体积公式推导先只考虑罐体只有纵向变位时的体积的推导：（1）、左侧球冠处油体积公式推导对左侧球冠处油体积采用三重积分进行计算，由于左侧球冠的球心与与其相交的圆柱体圆的圆心不共心， 现对球冠所在球的半径进行求解：图4 球冠半径求取示意图如图所示可得满足方程：即：，解得：现建立空间直角坐标系如下图所示：以球冠球心为坐标原点，垂直于圆柱高向上方向为z轴，平行于圆柱高方向指向球冠为y轴，垂直圆柱高向里方向为x轴。图5 纵向倾斜储油罐坐标系建立示意图由以球心为坐标原点，得：球面方程为：球冠与圆柱体的交面平行于xoz平面（，），距坐标原点的距离，故交面方程为：现将过球心的纵切面绘制如下：图6 过球心纵切面示意图据图分析可知，油平面平行与x轴，油平面法线为图中n所示，n垂直于平面且与z轴成角，则平面的法向量为（0，-sin,cos）,从图可知平面过点，其中：据平面的点法式方程可知，油平面方程为：积分次序为：先z后x最后y；根据z的积分限将积分分成两部分，如下所述：当油平面与球冠的交线在球冠左端点以下(即b点以下)时，油平面以下，交面以左与球面所围成的球冠面处体积即为所求，积分区域为：其中z轴方向的积分限下限为球面，上限为油平面；由球面方程解出可得z轴方向的积分下限，通过油平面方程解出z的表达式既得z轴方向的积分上限：x轴方向的积分限上、下限为球面与油平面的交线在xoy平面的投影，由球面方程与油平面联立如下所示：由上述联立方程消去z即得在xoy平面的投影方程求解可得x轴方向的积分下限，上限：y轴方向的积分限下限为球冠与圆柱体的交面所对应的y坐标,即，积分限上限为过球心的纵切面、球面、油平面的交点对应的y坐标：图7 y轴方向积分限求取示意图由图可知直线方程斜率为，且过点，为直线方程与圆方程的交点：，联立方程求解可得： 即，解得：，又有综上将三重积分化为一重积分可得： 其中：当油平面与球冠的交线在球冠左端点以上(即b点以上)时，油平面以下，交面以左与球面所围成的球冠面处体积即为整个球冠的体积减去积分区域为部分的体积：整个球冠的体积积分区域（注：与第一种情况相比被积函数区别只在于z轴的积分上、下限分别为、） ，化三重积分为一重积分得：（注：与第一种情况相比被积函数区别只在于前的系数为负）。其中：，其他积分上、下限同第一种情况。由于积分限、为油位高度h跟纵向偏角的函数，故得到了左侧球冠油容积。综上两种情况可得：当油位高度h在（0，2r）之间变化时，左球冠容积表达式为：（2）、中间圆柱体处油体积公式推导由已求得圆截面面积与油位高度之间的关系：图8 积分限分段示意图同样与截椭圆柱体体积公式计算雷同（在此不做详细推导），跟据其对称关系，可得出：其中：（3）、右侧球冠处油体积公式推导：与左侧球冠体积求法类似，将坐标原点建立于球冠对应球的球心。垂直于圆柱高向上方向为z轴，平行于圆柱高方向指向球冠为y轴，垂直圆柱高向外方向为x轴。得球面方程、球冠与圆柱面交面方程、油平面方程如下:其中： ； 。其中： 其中：，其余积分上、下限表达式同左侧球冠体积表达式。由于纵向偏斜时，左侧球冠处油体积，中间圆柱体处油体积，右侧球冠处油体积。三部分体积表达式均为以 、为变量且按照分段的表达式，而与坐标系的建立无关。故可得主体为圆柱体，两端为球冠体的储油罐纵向倾斜时，油容积与油位高度h之间的关系为三部分体积表达式的和，即：=+4.2.2.4罐体纵向变位加上横向偏斜时体积公式推导由于当罐体纵向偏斜角时，在此基础上横向偏斜，此时圆截面所对应的液面高度保持不变，故横向偏斜以后油实际容积不发生改变，改变的只是油浮子与测的的油位高度发生变化，油浮子从原来跟油面垂直变换到与油面成角，从而导致油浮子测得的油位高度与实际油位高度发生了偏差，而在上步罐体纵向偏斜时体积求解是油浮子所测的油位高度的函数，当时油浮子测量高度与实际高度相等。因此，此时，只要将油面实际高度h用罐体横向偏斜时油浮子所测得的油位高度来表示，即，而= 联立 ，可得油面实际高度h用罐体横向偏斜时油浮子所测得的油位高度两者之间的关系用下图表示：图9 横向偏转时液面高度变化关系示意图由图可知：与均过圆心，油面高度高于圆心或者低于圆心，均可得：带入表达式=，便可得，既得罐内储油量与油位高度以变位参数（纵向倾斜角及横向偏转角之间的）一般关系。五 模型的求解5.1模型一的求解 根据问题一可知，两端平头的圆柱体具体参数为：小椭圆油罐横截面椭圆长轴，短轴，,油罐总长。5.1.1确定储油罐实际储油量与理论储油量之间的偏差根据对问题一所建立的模型，罐体无变位时油罐储油量与油位高度之间的关系为：其中h 即为油位高度，将附件1中所给的无变位进油累加进油量与罐内油量初值相加，得到各个油位高度对应的油量。对应附件1中无变位进油的油位高度根据上式运用excell，求出对应进油高度所计算的理论油量，将理论油量与附件中所给的实际油量进行比较分析，得出两者之间的相对偏差，从数据（见附表2）可以看出相对偏差很稳定，偏差的方差很小，故通过取平均得出平均相对偏差。偏差分析：从实测数据可知理论油量大于实际油量，据分析可知，造成这一现象主要是因为题给的是一个简化模型，忽略了壁厚将外径等同于内径，导致带来较大的偏差。同时由于储油罐内的油浮子、管道等也占据一定的空间（这部分只占很少的一部分）造成理论值大于实际值。偏差修正：根据理论油量与实际油量之间的偏差关系，可得油量的修正函数为：。偏差校正检验（见程序2），对应附件附件1中无变位出油的油位高度求出与之对应的理论油容量，并按照上述修正函数进行修正后，将两两相邻的油容量相减即可得出两个时间点之间的绝对出油量与实际所给绝对出油量（相邻两累计出油量之间的差值）取差值，求其相对误差，再对相对误差求和取平均得，故验证了所求得的修正函数关系是正确的。5.1.2罐体变位后对罐容表影响的研究根据对问题一所建立的模型，罐体无变位时油罐储油量与油位高度之间的关系： (详见4.1.2问题一模型建立)根据上述罐体变位后油容积与油位高度h之间的函数关系，对应附件1中倾斜变位进油中各个油位高度所对应的倾斜变位时的理论油量进行求解，同样运用5.1.1中的偏差校正函数进行校正得出各油位高度对应的油容量实际值进行比较分析，计算出两者之间的现对偏差（见附表3）。从数据可以看出随着油位高度从411.29mm到1035.36mm逐渐增加，两者之间的偏差满足从一个正值不断减少到零，再继续递减到一个负值的趋势。变位对罐容表影响分析：从数据分析可知在液位低于一定值时理论值大于实际显示值，并满足随着液位降低两者之间的偏差逐渐增大；当液位高于一定值时理论值小于实际显示值，且满足随着液位高度的增加两者之间的偏差呈增大趋势。这是由于储油罐发生了纵向倾斜而罐容表的读数没有进行相应的修正，这就导致当油位高度h=0时，实际上是存在油容量的而罐容表标注误以为油容量为0；当h达到最大时，实际油并没有满罐，而罐容表误以为油满，故出现h较小时理论值大于显示值，当h较大时理论值小于显示值，符合两者偏差的变化规律。故可得纵向变位时液位较小时读数偏小，液位较高时读数偏大。5.1.3罐体变位后罐容表标定罐体按照倾斜角为a=4.10的纵向变位后，当h从0到1.2之间按1cm间隔递增时，根据模型一所建立的变位后容积与油位高度的函数关系，运用matlab运算出其对应的油容积，并对其进行标注。标注结果见（附表1）。5.2问题二的求解问题要去对一个主体为圆柱体，两端为球冠体的储油罐储油体积进行计算。主体圆柱体的半径为r=1.5；圆柱体高为l=8；两端球冠高为；球冠对应球体半径为r；油浮子距圆柱体两端的距离分别为、；罐体纵向倾斜角度为，横向偏转角为。根据球冠对应球体半径满足关系式，解得r=1.625。根据4.2.2建立的模型，所得主体为圆柱体，两端为球冠体的储油罐储油体积与纵向倾斜角度为，横向偏转角为以及油面高度之间的关系满足：（详见问题二模型建立）。对应于附件2中显示各个油高，运用matlab（见程序4）求出其对应的油量容积，与附件2中的数据进行比较发现与显示油容量容积几乎完全一致，可以准确到小数点后两位，据此可以验证我们所建立的罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系模型是正确的。得到后，我们利用附件二两次注油间的数据，采用用两种方法分别计算出了、。利用第二次注油后的数据检验了、的准确性。计算、：方法一：采用遍历法穷举、。根据jjg 266-1996卧式金属罐容积检定规程，使、分布在在可能出现的极大范围内，故认为其分别在(0，10）间变化,步长均为0.1。对于每组确定的、再由附件2中两个相邻的高h,可求得相邻的之差，和表中给出的相邻v之差比较，越使两者接近的一组、约可靠，故可转化为求最小值问题。附件二中，两次注油间的数据约三百组，从而找出最可靠的、，计算结果见下表。计算程序见（程序6）。方法二：采用改进差分进化算法的基础上与粒子群算法进行融合，在方法一的基础上求解最小值。其融合寻优策略如下： pso算法与de算法都是基于群体的启发式算法，它们的最主要区别是新个体的产生方式。pso算法在优化求解过程中前期收敛速度较快，后期易陷入局部最优。差分进化算法在选择操作中采用的是一种“贪婪”搜索策略，即经过变异和交叉操作后的个体与父代个体进行竞争，只有当其适应度好于父代时才被选作为子代，否则直接进入下一代，该机制可以增加算法的收敛速度，但也容易陷入局部最优点而使算法停滞5。针对de和pso单个算法均异陷入局部最优的缺点，本文考虑将两种算法进行融合。由于两种算法存在差异，故一种算法的局部最优可能对另一种算法不适应，即两者均有能力跳出对方的局部最优解，通过两种算法的结合便可以很好的跳出局部最优去搜索全局最优解。同样两者之间信息的交流，给双方一个扰动，有利于跳出局部最优。下面通过图示来说明其中的关系： 图10 de模型中个体运动方式 图11 ndepso模型中个体运动方式图10给出了在基本nde模型中个体运动方式,随着迭代次数的增加种群的平均适应度方差不断减小，即所有个体都趋同于群体中最优个体（gebest）,若此时的最优个体为非全局最优，同时所给扰动不够，则会陷入局部最优而无法逃逸。若将 nde跟pso进行融合，则个体运动方式如图11所示。若陷入局部最优，但是存在来自于pso的信息交流，此时则会引导个体逃逸原来局部最优点，不断去搜寻全局最优解，从而实现精确求解。nde和pso算法融合策略：将随机产生种群按适应度升序排列后分成大小相等的较小规模种群，一半适应度排序较劣称为nde子种群，另一半适应度排序较优称为pso种群；nde子种群按照nde进行进化，pso子种群按照pso进行进化，两个子种群各自进化形成的新子种群相组合形成我们的下一代种群，然后继续按照上述步骤反复操作，直至寻找到最优解。算法流程如图12所示：图12 ndepso算法流程由于两种算法结合带来计算的时间复杂度的增加，为了克服这一问题，本文在ndepso算法上继续改进。通过改进选择策略，提高收敛速度。改进选择策略为：将父代与试验个体一起入“选择池”按适应度排序，适应度优的留下，称之为串行选择。串行选择的优点在于若父代与对应试验个体在适应度排序中都较优则都被留下来，而在基本算法中父代与试验个体只能选其一，这样加速了进化速度，加快了收敛速率。参数设置：函数维数为2，种群规模设为20，最大迭代次数100，de的缩放因子，交叉概率因子，学习因子均为1.4962，惯性权重为0.7298。种群中每个个体都含有变量（、），经过优化寻优，可得结果如下表所示。计算程序见(程序7、8)，响应曲线(见附图8)误差方法一2.14.44.9e-5方法二2.1294.2443.7e-6由于方法一采用步长较大，显然方法二较有优势。故采用方法二数据，认定=2.219、=4.244。将所得的（、）带入中便较容易得出罐容比，计算程序见（程序5）。结果见（附表4）。模型二可靠性分析：利用附件二第二次注油后的数据对、值的可靠性进行分析。由计算出的、以及标中每个油标高h带入，相邻差值与表中对应相邻v的差值比较，。求解程序见（程序12）=1.5*e-4。六 模型的评价6.1 模型优点本文分别采用了截面法、三重积分法等基本原理，使得计算值更为精确，较好的阐明了罐内储油量、油位高度以及变位参数这一动态平衡数学模型，同时研究了一种混合智能计算方法ndepso，使得计算机能够精确地寻找变位参数，并控制误差在极低的范围，以实现储油罐的准确变位识别，并能根据所得到的变位参数，高精度地进行罐容表标定。同时根据本命题本文进行了基于matlab储油罐罐容表标定gui设计，图形开发界面充分发挥了matlab的强大计算功能，使得用户可以在直观简洁的操作界面上，只需输入相关数据，按按鼠标就可以完成复杂的计算，并可以得到可视化的计算结果。所设计的模型系统具有：1） 适用性强，涉及范围广，针对不同规格的储油罐，可灵活根据参数进行罐容表，且操作方便、简单、明确。2） 能针对不同的进/出油及对应油位高度之差，进行储油罐变位识别，且计算精确、可靠性高。6.2 模型缺点及改进任何模型、系统都受到实际生活中的各种限制，本模型也不例外，为了简化模型，基本假设很多都是理想状态。实际上由于受到外界温度，以及进出油所产生的压力变化和流速等因素的影响，油的密度等自身属性将发生轻微变化，造成系统与实际值产生偏差。同时计算机计算的精度同样会影响到最后的结果值。此外，本模型在运行及gui设计方面还存在缺陷，算法复杂度较高，运行所需时间较长，这些都有待于以后再作进一步研究。七 模型的推广本文中所建立的模型不仅对于解决题目中的问题具有很强的实用性，而且具有很好的推广性。例如，本文提出了一种ndepso算法，该算法结合了改进的de和pso算法的优良性能，在解决局部最优问题上得到了很好的应用，寻优能力较强，并且选择操上采用了串行选择策略，加速了进化速度，同遍历算法进行比较，在寻找最优解的过程中具有一定的优越性。由于本文是在基于大量的matlab实验获得的，程序较多，运行操作较复杂。本文创造性的进行了基于matlab的gui设计，使得用户可以进行直观的操作，并得到准确的罐容表。经设计后的gui设计，只需在matlab命令窗口中输入question即可进入选择系统界面，经选择可进入相应的的系统界面，这时可在文本框中输入储油罐的各参数，按下“求解罐容表”按钮，即可输出相应的结果值，下一步可选择进入另一系统，亦可选择退出。八 参考文献【1】刘卫国 matlab程序设计教程，中国水利水电出版社，2006年3月第三版【2】韩中庚 数学建模方法及其应用，高等教育出版社，2005年6月第一版【3】frank r.giordano等，数学建模 机械工程出版社，2005年1月第三版【4】朱德通，优化模型与实验，统计大学出版社，2003年4月第一版【5】price k，storn r，lampinen jdifferential evolution-a practical approach to global optimization mberlin：springer-verlag press，2005【6】栾丽君、谭立静、牛奔。一种基于粒子群优化算法和差分进化算法的新型混合全局优化算法j。信息与控制,2007,36(6):708-714。【7】zhang m, luo w j, wang x f. differential evolution with dynamic stochastic selection for constrained optimizationj. in- formation sciences: an international journal, 2008, 178(15):3043-3074.附件：附表一：第一问罐容标定表h(mm)v(l)h(mm)v(l)h(mm)v(l)h(mm)v(l)01.6179559760.31608.90397160.621821.5834270.933082.5712950.013.4121632480.32643.14407120.631863.5029230.943119.7753930.026.0526473020.33677.8771720.641905.4569160.953156.6008010.039.6389360080.34713.08182230.651947.4336170.963193.0293550.0414.259541530.35748.73724690.661989.4210430.973229.042790.0519.994279420.36784.82363690.672031.40760.983264.6217810.0626.915234920.37821.32205290.682073.3812080.993299.7465190.0735.092077530.38858.21403880.692115.33007913334.3960360.0844.58693970.39895.48191150.72157.2424251.013368.5487830.0955.458765120.4933.10837440.712199.1062631.023402.1818590.167.763018830.41971.07680740.722240.9097081.033435.2713970.1181.551783870.421009.3709770.732282.6407771.043467.7918870.1296.874534290.431047.9752290.742324.2872961.053499.7166590.13113.77823180.441086.8740060.752365.8371871.063531.0171110.14132.30742230.451126.0524270.762407.2782731.073561.6626120.15152.49834170.461165.4958050.772448.5982821.083591.6202110.16174.18285560.471205.1897420.782489.7845581.093620.8541550.17197.12267690.481245.1200330.792530.824731.13649.3254070.18221.19015850.491285.2730550.82571.7060441.113676.990870.19246.292730.51325.6350850.812612.4156491.123703.8021320.2272.35627690.511366.1927880.822652.9404981.133729.7045020.21299.31876340.521406.9329280.832693.2674521.143754.634010.22327.12656040.531447.8426510.842733.3829811.153778.5142250.23355.7325130.541488.9089130.852773.2734621.163801.2477430.24385.09410430.551530.1192490.862812.9251721.173822.6970620.25415.17297260.561571.4610.872852.3238111.183842.6023460.26445.93375190.571612.9217020.882891.454981.193860.8477380.27477.34407120.581654.4891780.892930.304091.23877.4752690.28509.37378250.591696.1512550.92968.8558720.29541.99438010.61737.8955630.913007.0949630.3575.17948430.611779.7103130.923045.005515附表二：无变位进油时理论容积值与实际值的偏差x轴 hy轴 vv 理论值影响159.02312.00322.88093.370%176.14362.00374.63133.372%192.59412.00426.36323.369%208.50462.00478.133.374%223.93512.00529.85013.369%238.97562.00581.6043.371%253.66612.00633.35033.371%268.04662.00685.07933.369%282.16712.00736.8453.372%296.03762.00788.57593.370%309.69812.00840.32713.371%323.15862.00892.05443.369%336.44912.00943.80063.369%349.57962.00995.54053.369%362.561012.001047.2953.370%375.421062.001099.0533.371%388.161112.001150.8063.372%400.791162.001202.5523.372%413.321212.001254.2923.372%425.761262.001306.033.371%438.121312.001357.7733.371%450.401362.001409.493.369%462.621412.001461.2343.369%474.781462.001512.9773.369%486.891512.001564.7363.370%498.951562.001616.4843.371%510.971612.001668.2413.371%522.951662.001719.9823.371%534.901712.001771.7283.371%546.821762.001823.4573.370%558.721812.001875.1913.370%570.611862.001926.9533.371%582.481912.001978.6773.370%594.351962.002030.4323.370%606.222012.002082.1963.371%618.092062.002133.953.372%629.962112.002185.6723.371%641.852162.002237.4313.371%653.752212.002289.1613.371%665.672262.002340.8853.370%677.632312.002392.6693.371%678.542315.832396.6033.370%690.532365.832448.3693.371%690.822367.062449.6193.370%702.852417.062501.3943.371%714.912467.062553.1133.371%727.032517.062604.8823.371%739.192567.062656.5893.370%751.422617.062708.3363.370%763.702666.982760.0083.371%764.162668.832761.9383.371%776.532718.832813.6633.370%788.992768.832865.4183.371%801.542818.832917.1683.371%814.192868.832968.9173.371%826.952918.833020.6653.371%839.832968.833072.413.371%852.843018.833124.1433.371%866.003068.833175.893.371%879.323118.833227.6333.371%892.823168.833279.3823.371%892.843168.913279.4583.371%906.533218.913331.1783.370%920.453268.913382.9363.371%934.613318.913434.6733.370%949.053368.913486.4243.371%963.803418.913538.1653.371%978.913468.913589.9193.371%994.433518.913641.6683.371%1010.433568.913693.4173.371%1026.993618.913745.1363.370%1044.253668.913796.8893.371%1062.373718.913848.6483.371%1081.593768.913900.3843.371%1102.333818.913952.1373.371%1125.323868.914003.8613.371%1152.363918.914055.6123.371%1193.493968.914107.363.371%附表三：倾斜变位时标高对容积的影响x 轴hy轴 v（实验）校正v（理论）影响411.29962.86975.999 1.3462%423.451012.861022.655 0.9578%438.331062.861080.358 1.6197%450.541112.861128.176 1.3576%463.901162.861180.947 1.5316%477.741212.861236.076 1.8782%489.371262.861282.737 1.5496%502.561312.791336.000 1.7373%514.691362.791385.278 1.6234%526.841412.731434.898 1.5449%538.881462.731484.302 1.4534%551.961512.731538.212 1.6566%564.401562.731589.690 1.6959%576.561612.731640.179 1.6735%588.741662.731690.897 1.6658%599.561712.731736.057 1.3437%611.621762.731786.490 1.3300%623.441812.731835.999 1.2674%635.581862.731886.910 1.2814%646.281912.731931.816 0.9880%658.591962.731983.501 1.0472%670.222012.732032.331 0.9645%680.632062.732076.025 0.6404%693.032112.732128.034 0.7192%704.672162.732176.800 0.6463%716.452212.732226.077 0.5996%727.662262.732272.883 0.4467%739.392312.732321.750 0.3885%750.902362.732369.572 0.2887%761.552412.732413.691 0.0398%773.432462.732462.741 0.0004%785.392512.732511.924 -0.0321%796.042562.732555.537 -0.2815%808.272612.732605.386 -0.2819%820.802662.732656.174 -0.2468%832.802712.732704.522 -0.3035%844.472762.732751.243 -0.4175%856.292812.732798.243 -0.5177%867.602862.732842.892 -0.6978%880.062912.732891.689 -0.7276%892.922962.732941.593 -0.7186%904.343012.732985.491 -0.9124%917.343062.733034.954 -0.9152%929.903112.733082.197