要使k为常数，则α、β应使等式成立，即为该方程的根，又因我们为了计算时方便，近似地取等效面积代替，则需用α或β中的一个来表示偏差。 计算得??1.1350。再将α代入V对H的含参数的积分函数式，可得??3.50。
图f
由上图可知推出的函数非常近似于实际数据，在误差允许的范围内，即验证了α、β的值正确，??1.1350，??3.50。
由图f我们也可以看到，通过模型画出的曲线与实际检测数据模拟的曲线偏差相当小，这就验证了模型的正确性。
由于我们的模型建立是通过积分方法又经过matlab作数据处理，同时也参照了相关资料而综合得出的，对于小误差范围内，可以推广到主体为椭圆柱的储油罐，从而可以验证方法的可靠性。
罐体变位后油位高度间隔10cm的罐容标定值如下表： 表3 罐体变位后油位高度间隔10cm的罐容标定值 油位高相应油油位高相应油油位高相应油油位高相应油度间隔量/m^3 度间隔量/m^3 度间隔量/m^3 度间隔量/m^3 /m 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.0135 1.1288 2.4736 4.0316 5.7864 7.7216 9.8207 12.067 /m 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 14.445 16.937 19.528 22.201 24.939 27.726 30.546 33.382 /m 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 36.218 39.037 41.823 44.56 47.231 49.82 52.31 54.685 /m 2.4 2.5 2.6 2.7 2.8 2.9 3 56.929 59.024 60.956 62.707 64.26 65.6 66.71
16
参考文献
[1]胡庆波，储油罐计量系统误差分析及对策，油气田地面工程，第29卷第5期：59-60页，2010.05 [2]陈佰军 左振滨，编制卧式驻油罐容量表的计算机模型，黑龙江八一农垦大学学报，第7卷第2期：55-60，1993.12
[3]王郑耀，卧式加油罐剩余油料体积的计算，西安交通大学 理学院，2004.08 [4]姜起源 谢金星 叶俊，数学模型，北京：高等教育出版社，2003
17
附录
图a程序：
plot(Book12,Book11,'-.b') >>hold on
>>ezplot(v1,[0,1.2])
>> xlabel('油位高度/m');ylabel('油量/m^3');legend('原始数据','推导函数');
图b程序： >> hold on
>> plot(Book11,Book12,'-.b')
>>b1=polyfit(Book12,Book11,2);b01=polyval(b1,Book12);plot(Book12,b01,'r') >>
b1=polyfit(Book12,Book11,3);b01=polyval(b1,Book12);plot(Book12,b01,'-y') >>
b1=polyfit(Book12,Book11,4);b01=polyval(b1,Book12);plot(Book12,b01,':m') >> xlabel('x');ylabel('y');legend('原始数据','2阶多项式拟合','3阶多项式拟合','4阶多项式拟合'); 图c程序：
plot(Book32,Book31,'-.r') >> hold on
>> ezplot(v2,[0,1.2]) >>
b1=polyfit(Book32,Book31,3);b01=polyval(b1,Book32);plot(Book32,b01,':g')
>> title('倾斜变位进油');xlabel('油位高度/m');ylabel('油量/m^3');legend('原始数据','推导函数','3阶多项式拟合'); 图d程序：
>>g=[0，-0.05，-0.1，-0.15，-0.2]; >> ezplot(v2+g,[0,1.2]) 图e程序：
plot(booka,bookb,'-.b') >> hold on >>
b1=polyfit(booka,bookb,3);b01=polyval(b1,booka);plot(booka,b01,'r') 图f程序：
plot(booka,bookb,'-.b') >> hold on >>
b1=polyfit(booka,bookb,3);b01=polyval(b1,booka);plot(booka,b01,'r') >> ezplot(v,[0,3])
18