2014年美赛数学建模A题论文

Team # 31725 Page 21 of 25

title('Evaluation Of Capability') figure %2?1??ò3?3μí????t

scatter(DifferentValueThree,PropertyAbnormalTwo,'.') xlabel('DifferentValueThree') ylabel('PropertyAbnormalTwo') title('Evaluation Of Capability') figure

for VariableFour=1:LengthConditionOne

FinalProperty(VariableFour)=0.3*PropertyNormal(VariableFour)+0.6*PropertyAbnormalOne(VariableFour)+0.1*PropertyAbnormalTwo(VariableFour); %?óè¨1y3ì

FinalDifferentValue(VariableFour)=(DifferentValueOne(VariableFour)+DifferentValueTwo(VariableFour)+DifferentValueThree(VariableFour))/3; end

scatter(FinalDifferentValue,FinalProperty,'.') xlabel('FinalDifferentValue') ylabel('FinalProperty')

title('Evaluation Of Capability')

Program Two: %??3μ?ó?ù?èμ??à1?±í′?ê? n=4;

Velocity=30:10:120;

Acceleration=[0.287,0.345,0.475,0.442,0.467,0.490,0.551,0.584,0.614,0.644]; p=polyfit(Velocity,Acceleration,n) double_x=0:0.01:1;

double_y=polyval(p,double_x);

plot(double_x,double_y,'-b',Velocity,Acceleration,'r','Markersize',20)

Program Three: %?￡?aé?·??ú??·?3ì originalvelocity=60:90;

velocity=originalvelocity./3.6; originaldifferentvalue=0:30;

differentvalue=originaldifferentvalue./3.6; h=1.2;%í?é? B=0.4; %

C=-0.015; o=1; %

SafeDistance=50; for i=1:31

Team # 31725 Page 22 of 25

Acceleration=velocity(i)^4*(-3.115e-08)+velocity(i)^3*(1.007e-05)+velocity(i)^2*-0.001165+velocity(i)*0.05992-0.7174;

?celeration=velocity(VariableOne)^4*(-3.115e-08)+velocity(VariableOne)^3*(1.007e-05)+velocity(VariableOne)^2*-0.001165+velocity(VariableOne)*0.05992-0.7174; %n=4;

%Velocity=30:10:120;

% Acceleration=[0.287,0.345,0.475,0.442,0.467,0.490,0.551,0.584,0.614,0.644]; % p=polyfit(Velocity,Acceleration,n) for j=1:31

chasetime=((sqrt(differentvalue(j)*differentvalue(j)+2*Acceleration*SafeDistance)-differentvalue(j))/Acceleration);

chasedistance=velocity(i)*chasetime+Acceleration*chasetime*chasetime/4;

ROC=h*h*sqrt(1+chasedistance*chasedistance)*sqrt(1+chasedistance*chasedistance)*sqrt(1+chasedistance*chasedistance)/chasedistance; %ROC is radius of curvature

PropertyNormal(o)=sqrt(127*ROC*(B+C))/chasetime; o=o+1;

%PropertyNormal measures the property of this scheme. end end

originalarray=reshape(PropertyNormal,31,31); middlearray=originalarray'; finalarray=max(middlearray); maxposition_x=differentvalue; maxposition_y=finalarray; n=2;

% n=intput('input:n=');

p=polyfit(maxposition_x,maxposition_y,n)

Program Four:

originalvelocity=60:90;

velocity=originalvelocity./3.6; originaldifferentvalue=0:30;

differentvalue=originaldifferentvalue./3.6; h=1.2;%í?é? B=0.4; %

C=-0.015; o=1; %

SafeDistance=50; for n=1:31

for m=1:31

Acceleration=velocity(m)^4*(-3.115e-08)+velocity(m)^3*(1.007e-05)+velocity(m)^2*-0.001

Team # 31725 Page 23 of 25

165+velocity(m)*0.05992-0.7174;

?celeration=velocity(VariableOne)^4*(-3.115e-08)+velocity(VariableOne)^3*(1.007e-05)+velocity(VariableOne)^2*-0.001165+velocity(VariableOne)*0.05992-0.7174; %n=4;

%Velocity=30:10:120;

% Acceleration=[0.287,0.345,0.475,0.442,0.467,0.490,0.551,0.584,0.614,0.644]; % p=polyfit(Velocity,Acceleration,n)

chasetime=((sqrt(differentvalue(n)*differentvalue(n)+2*Acceleration*SafeDistance)-differentvalue(n))/Acceleration);

chasedistance=velocity(m)*chasetime+Acceleration*chasetime*chasetime/4;

ROC=h*h*sqrt(1+chasedistance*chasedistance)*sqrt(1+chasedistance*chasedistance)*sqrt(1+chasedistance*chasedistance)/chasedistance; %ROC is radius of curvature

PropertyNormal(o)=sqrt(127*ROC*(B+C))/chasetime; o=o+1;

%PropertyNormal measures the property of this scheme. end end

originalarray=reshape(PropertyNormal,31,31); middlearray=originalarray'; finalarray=min(middlearray); maxposition_x=differentvalue; maxposition_y=finalarray; n=2;

% n=intput('input:n=');

p=polyfit(maxposition_x,maxposition_y,n)

Program Five:

x=[0.5184 0.9793 1.7857 2.8917 4.5507 5.4954 6.4171 6.8779 7.3387 7.6152 8.0069 8.2604 8.6290 9.1129];

Team # 31725 Page 24 of 25

y=[242.8655 244.2690 245.9064 246.3743 245.6725 243.0994 238.4211 231.1696 226.0234 220.8772 216.4327 210.1170 203.8012 197.4854]; maxposition_x=x; maxposition_y=y; n=4;

%n=intput('input:n=');

p=polyfit(maxposition_x,maxposition_y,n) x1=[0.5184 0.4724 0.5876 1.0484 1.6014 2.4078 3.2834 4.3664 5.5184 6.3710 7.5461 8.3295 9.0668]; y1=[243.0994 226.4912 215.4971 200.2924 193.2749 187.6608 181.3450 179.7076 181.1111 182.5146 187.8947 192.1053

Team # 31725 Page 25 of 25

197.7193]; maxposition_x=x1; maxposition_y=y1; n=6;

%n=intput('input:n=');

p=polyfit(maxposition_x,maxposition_y,n)