2010年全国研究生数学建模竞赛优秀论文选-附录B 程序清单.doc

附录B程序清单

functionB1B2

%2010年全国研究生数学建模B题

%1小题2小题主要程序。用来计算模型及识别其中的系数。

%程序分为

%a程序运行参数选择；

%b符号计算解出模型的特解；

%c选择截面积、速度、投放高度及试验数据；

%d基于包映射的二分法搜索

%等四部分。

%吉林大组。2010年9月。

clc;

%程序运行参数选择

flgSharp=1;%重物形状选择。1大实心方体2大空心方体3小空心方体4小空心蜂巢5大空心蜂巢6大实心蜂巢7小实心蜂巢8大三角9小三角

flgV=1;%重物速度选择。1代表v=0.55，2代表v=0.47，3代表v=0.4。

flgH=1;%重物高度选择。1代表H=0，2代表H=0.05，3代表H=0.12。

flgPose=1;%重物姿态选择。1代表平放，2代表立放，3代表竖放。

%选择说明：

%选择1大实心方体时，flgV、flgH可任意选择。

%选择2大空心方体3小空心方体4小空心蜂巢5大空心蜂巢6大实心蜂巢7小实心蜂巢8大三角9小三角时，

%flgV、flgH只能选择[1,1]或[2,2]，其余数据没有填入。

%另外，选择8大三角9小三角时，flgPose只能选择1平放。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%符号计算解出模型的特解

sy=dsolve(m*D2sy=m*g-pw*g*Tj-A1*pw*c1*Dsy^2-A2*c2*Dsy,Dsy(0)=g*sqrt(2*H/g),sy(0.04)=0);

sx=dsolve(m*D2sx=A3*pw*c3*(v-Dsx)^2+A4*c4*(v-Dsx),Dsx(0)=0,sx(0.04)=0);

%常数部分赋值

g=9.8;%#ok*NASGU%重力加速度

pw=1000;%#okNASGU%水密度

ps=2300;%重物密度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%选择截面积、速度、投放高度及试验数据

ifflgSharp==1

strSharp=大实心方体%#okNASGU,NOPRT%在commandwindows中输出形状

Tj=0.08*0.08*0.04;%重物体积

ifflgV==1

v=0.55%#okNASGU,NOPRT%在commandwindows中输出速度

ifflgH==1

H=0%#okNASGU,NOPRT%在commandwindows中输出投放高度

ifflgPose==1

strPose=平放%#okNASGU,NOPRT%在commandwindows中输出姿态

%各面截面积

A1=(0.08*0.08);

A2=(0.08*0.04*4);

A3=(0.08*0.04);

A4=(0.08*0.08*2+0.08*0.04*2);

%试验数据

y55=[27,2525.25,2220,1815.8,13107.166666667,2.733333333,0.166666667,-0.416666667]/100;

x55=[-1.833333333,-1.6,-1.25,-0.666666667,0.15,0.85,2.716666667,4,5.633333333,8.666666667,101213.1]/100;

end

ifflgPose==2

strPose=立放%#okNASGU,NOPRT%在commandwindows中输出姿态

%各面截面积

A1=(0.08*0.04);

A2=(0.08*0.08*2+0.08*0.04*2);

A3=(0.08*0.08);