A题讲评数学建模解析.ppt

2021/2/7 * 2021/2/7 * α h Y Z X 0 V V α x y 0 z h L d h-r 油浮子坐标 (0, h-r, L/2-d) 过油浮子的油平面方程 坐标变换 2021/2/7 * x y 0 z h r s R V α 油位高度h的储油区域为D 油位高度h时的储油量V 定义示性函数 2021/2/7 * 利用MATLAB的三重积分程序计算储油量 triplequad(@(x,y,z)(-x*cos(a)*sin(b)+y*cos(a)* cos(b)+z*sin(a)=cos(a)*cos(b)*(h-r)+sin(a)* (L/2-d)).*(x.^2+ y.^2 =r^2).*(z=-s-sqrt(R^2-x.^2-y.^2)).*(z=s+sqrt(R^2-x.^2-y.^2)),-r,r,-r, r,-(s+R),s+R) triplequad (@function, x1,x2, y1,y2, z1,z2) 被积函数 function 在长方体中的三重积分: 用于function 是示性函数 I(x,y,z)的情况: 计算时间较长;计算精度受积分域大小的影响. 在上述坐标旋转下罐体纵向倾斜角α为负(β正负无关). 2021/2/7 * 1）罐体变位后储油量的计算（重积分）: 两端有油和一端有油区域的划分; 积分限的确定; 积分顺序的选择。 对学生论文缺陷的评述 2）用所谓“等效高度”:得出不变位时高度与体积的关系V(h);变位后根据(主)截面面积不变得到等效高度与原高度的关系h =f(h);将它代入原式V(h) 得到变位后高度与体积的关系V (h) 。 2021/2/7 * 3）辨识变位参数所用的数据和准则不对:用附件2中的 “显示油量”与模型在“显示油高”下的计算油量作拟合。 如果模型正确,用这种错误方法得到的变位参数应该是α=0, β =0 。 ——反映对该问题的理解根本错误 4）含混地叙述:“用数据和最小二乘得到 α= , β = ” “显示油量”是罐体未变位时的计量值. 2021/2/7 * 5）对问题（1）用实验数据修正模型投入的精力过多（解释误差、用多项式拟合误差、用比例系数修正模型等）. 6）对问题（1）的要求“罐体变位后对罐容表的影响”未作讨论. 7）对问题（2）的参数辨识不做检验和灵敏性分析 。 8）问题（1）,（2）未按要求给出罐容表. 9）从文献中引用工程算法、近似公式的问题. 10）用随机模拟(Monte-Carlo)方法计算重积分的问题. 2021/2/7 * 赛题评析 全国大学生数学建模竞赛 2010年A题 2021/2/7 * 题目 题目的背景与分析 问题（1）的通常做法 问题（2）的通常做法 一种建立储油量模型的新方法 对学生论文的评述 题目的启示 2021/2/7 * 【研究背景】通常加油站都有若干个储存燃油的地下储油罐,并且一般都有与之配套的“油位计量管理系统”,采用流量计和油位计来测量进/出油量与罐内油位高度等数据,通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算,以得到罐内油位高度和储油量的变化情况。 图1 是一种典型的储油罐尺寸及形状示意图,其主体为圆柱体,两端为球冠体 2010A题:储油罐的变位识别与罐容表标定 2021/2/7 * 油 油浮子 出油管 油位探测装置 注油口 检查口 地平线 2m 6m 1m 1m 3 m 油位高度 图1 储油罐正面示意图 油位探针 2021/2/7 * 许多储油罐在使用一段时间后,由于地基变形等原因,使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位）,从而导致罐容表发生改变。 按照有关规定,需要定期对罐容表进行重新标定。。 图2是其罐体纵向倾斜变位的示意图, 图3是罐体横向偏转变位的截面示意图。 2021/2/7 * 油位探针 α 地平线 图2 储油罐纵向倾斜变位后示意图 油 油浮子 出油管 油位探测装置 注油口 检查口 水平线 2021/2/7 * 图3 储油罐截面示意图 （b）横向偏转倾斜后正截面图 β 地平线垂直线 油位探针 （a）无偏转倾斜的正截面图 油位探针 油位探测装置 3m 2021/2/7 * (b) 小椭圆油罐截面示意图 α 油 油浮子 出油管 油位探针 注油口 水平线 2.05mcm 0.4m 1.2m 1.2m 1.78m (a) 小椭圆油罐正面示意图 图4 小椭圆型油罐形状及尺寸示意图 2021/2/7 * 【建模目的】用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 1 为了掌握罐体变位后对罐容表的影响,利用如图4的小椭圆