数学建模2010A题 储油罐的变位识别与罐容表标定 程序.doc

数学建模2010A题 储油罐的变位识别与罐容表标定 程序 附录 （1）第一问 程序1 function weibainwei %未变位时罐体模型 y=xlsread(‘F:\数学建模第二期培训\第二题\A\附件1.xls’,1,’c2:c79’); x=xlsread(‘F:\数学建模第二期培训\第二题\A\附件1.xls’,1,’d2:d79’); x=x/1000; y=y/1000; a=0.89; b=0.6; l=2.45; for i=1:78 m=(x(i)-b)/b; y1(i)=a*b*l*(pi/2+m*sqrt(1-m*m)+asin(m))-y(i); end y1=y1’; plot(x,y1) a=polyfit(x,y1,1); yy=polyval(a,x); plot(x,y1,x,yy) 程序2 function bianwei %变位后罐体误差拟合a=0.89; b=0.6; l=2.45; d=0.4; thta=pi*4.1/360; x=xlsread(‘F:\数学建模第二期培训\第二题\A\附件1.xls’,3,’d2:d54’); y=xlsread(‘F:\数学建模第二期培训\第二题\A\附件1.xls’,3,’c2:c54’); x=x/1000; y=y/1000; for i=1:53 if x(i)gt;0amp;x(i)lt;(l-d)*tan(thta) N1=(x(i)-b+d*tan(thta))/b; M1=sqrt(1-N1); y1(i)=(pi*a*b/2)*(l-b/tan(thta))-(a*b/tan(thta))*(M1/3-M1-N1*asin(N1))-y(i); else if x(i)gt;(l-d)*tan(thta)amp;x(i)lt;(2*b-l*tan(thta)) N1=(x(i)-b+d*tan(thta))/b; N2=(x(i)-b+(d-l)*tan(thta))/b; M1=sqrt(1-N1); M2=sqrt(1-N2); y1(i)=pi*a*b*l/2-(a*b/tan(thta))*(M1/3-M2/3+N2*asin(N2)-N1*asin(N1)+M2-M1)-y(i); else y1(i)=0; end end end y1=y1’; plot(x,y1) a=polyfit(x,y1,2); yy=polyval(a,x); plot(x,y1,’bo’,x,yy,’b’) 程序3 function xiuzhenghou %修正后对比与误差分析 a=0.89; b=0.6; l=2.45; d=0.4; thta=pi*4.1/360; x=xlsread(‘F:\数学建模第二期培训\第二题\A\附件1.xls’,3,’d2:d54’); y=xlsread(‘F:\数学建模第二期培训\第二题\A\附件1.xls’,3,’c2:c54’); x=x/1000; y=y/1000; for i=1:53 if x(i)gt;0amp;x(i)lt;(l-d)*tan(thta) N1=(x(i)-b+d*tan(thta))/b; M1=sqrt(1-N1); y1(i)=(pi*a*b/2)*(l-b/tan(thta))-(a*b/tan(thta))*(M1/3-M1-N1*asin(N1))-(-0.5926*x(i)+0.8472*x(i)+0.1308); y2(i)=y1(i)-y(i); else if x(i)gt;(l-d)*tan(thta)amp;x(i)lt;(2*b-l*tan(thta)) N1=(x(i)-b+d*tan(thta))/b; N2=(x(i)-b+(d-l)*tan(thta))/b; M1=sqrt(1-N1); M2=sqrt(1-N2); y1(i)=pi*a*b*l/2-(a*b/tan(thta))*(M1/3-M2/3+N2*asin(N2)-N1*asin(N1)+M2-M1)-(-0.5926*x(i)+0.8472*x(i)+0.1308); y2(i)=y1(i)-y(i); else y1(i)=0; end end end y1=y1’; y2=y2’; figure(1) plot(x,y,’r’,x,y1,’b’) figure(2) bar(x,y2) 程序4 function shujuchansheng %产生数据表格 a=0.89; b=0.6; l=2.45; d=0.4; thta=pi*4.1/360; x=0.01:0.01:1.2 for i=