附录
(1)
第一问
程序 1
function weibainwei
%未变位时罐体模型
y=xlsread( 数学建模第二期培训 第二题 附件 1.xls' ,1, 'c2:c79' );
x=xlsread( 数学建模第二期培训 第二题 附件 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( 数学建模第二期培训 第二题 附件 1.xls' ,3, 'd2:d54' );
y=xlsread( 数学建模第二期培训 第二题 附件 1.xls' ,3, 'c2:c54' );
x=x/1000;
y=y/1000;
for i=1:53
if x(i)>0&x(i)<(l-d)*tan(thta)
N1=(x(i)-b+d*tan(thta))/b;
M1=sqrt(1-N1^2);
y1(i)=(pi*a*b/2)*(l-b/tan(thta))-(a*b^2/tan(thta))*(M1^3/3-M1-N1*asin
(N1))-y(i);