发布时间 : 星期四 文章2014年美赛数学建模A题论文更新完毕开始阅读8747b27df78a6529657d5336
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)