第三次培训 matlab代码

cumcm 2010A题 储油罐的变位识别与罐容表标定

通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。

请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。

（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。

（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度 ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。

附件1：小椭圆储油罐的实验数据
附件2：实际储油罐的检测数据

未倾斜

v2.m

function y = v2(h)
%计算未倾斜的v-h方程
l = 240;a = 89; b = 60;
t = (h - b)/b;
s = ( pi/2 + t.*sqrt(1 - t.^2) + asin(t) )*a*b;
y = s.*l;
y = double(y); 

v2_apply.m

clear, clc
%读取数据h
filename = '问题A附件1：实验采集数据表.xls';
sheet = 1;
xlRange1 = 'G2:G79';
h= xlsread(filename,sheet,xlRange1);
h = h / 10; %单位转化为厘米
V = v2(h)./1000; %单位转化为升

sheet = 1;
xlRange2 = 'E2:E79';
xlswrite(filename,V,sheet,xlRange2)

第一问 纵向变位

volumn.m

function