[数学建模]灵敏度分析记录

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/qq_37360631/article/details/86606924

1、灵敏度调试

function y=mfun(L)
g=9.8;      
m=75;
pi=3.1415926;
u=0.02;
%L=2;
a=18*pi/180;
k=0.5;
v0=60;
E0=1/2*m*v0*v0;

r=5;
E_edge=(1-2.*u.*L./r).*(E0.*exp(-2.*pi.*u)-3.*m.*g.*r.*u.*cos(a).*(exp(-2.*pi.*u)+1))./(4.*u.*u+1)-m.*g.*L.*cos(a);
h_max=(1./2.*k).*log(1+2.*k.*E_edge./m./g);
plot(L,h_max,'r');
hold on;
r=4;
E_edge=(1-2.*u.*L./r).*(E0.*exp(-2.*pi.*u)-3.*m.*g.*r.*u.*cos(a).*(exp(-2.*pi.*u)+1))./(4.*u.*u+1)-m.*g.*L.*cos(a);
h_max=(1./2.*k).*log(1+2.*k.*E_edge./m./g);
plot(L,h_max,'b');
hold on;
r=3;
E_edge=(1-2.*u.*L./r).*(E0.*exp(-2.*pi.*u)-3.*m.*g.*r.*u.*cos(a).*(exp(-2.*pi.*u)+1))./(4.*u.*u+1)-m.*g.*L.*cos(a);
h_max=(1./2.*k).*log(1+2.*k.*E_edge./m./g);
plot(L,h_max,'g');
hold on;

for L=0:1:10
    E_edge=(1-2.*u.*L./r).*(E0.*exp(-2.*pi.*u)-3.*m.*g.*r.*u.*cos(a).*(exp(-2.*pi.*u)+1))./(4.*u.*u+1)-m.*g.*L.*cos(a);
    h_max=(1./2.*k).*log(1+2.*k.*E_edge./m./g);
    plot(L,h_max,'b*');
end

grid on;  							%加表格，刻度均匀
legend('r=5','r=4','r=3');			%图例标注
xlabel('L');						%横纵坐标
ylabel('Maximum height');
set(gca,'XTick',0:5:60);			%修改x轴的坐标值，效果是0-0,5-1,10-2...
set(gca,'XTicklabel',{'0','1','2','3','4','5','6','7','8','9','10'});

end

2、

function y=mfun(v0)
pi=3.14159;
g=9.8;      
m=75;
u=0.02;
k=0.5;
%v0=60;
w1=0.5;
w2=0.5;
u_front=0.8;
u_back=0.01;
a=18;
L=2;
S=5;
r=2.8;
R=0.4;
qq=exp(-pi.*u);
qqq=m*g*r*cosd(a);
qq=exp(-pi.*u);
qqq=m*g*r*cosd(a);
h_min=1;            %最低是1m
h_max=8;
a_min=180;          %基础动作180    
a_max=1260;         %高难度1260

u=0.02;
E0=1./2.*m.*v0.*v0;
Eb=(E0*qq)-qqq*(u*qq-2*u*u-1)/(4*u*u+1);
E_edge=(1-2*u*L/r)*((Eb-u*m*g*S)*qq-((2*u*u+1)*qq+u)*qqq/(4*u*u+1))-m*g*L*cosd(a);
h=(1/2*k)*log(1+2*k*E_edge/m/g);
plot(v0,h,'r');
hold on;

u=0.12;
E0=1./2.*m.*v0.*v0;
Eb=(E0*qq)-qqq*(u*qq-2*u*u-1)/(4*u*u+1);
E_edge=(1-2*u*L/r)*((Eb-u*m*g*S)*qq-((2*u*u+1)*qq+u)*qqq/(4*u*u+1))-m*g*L*cosd(a);
h=(1/2*k)*log(1+2*k*E_edge/m/g);
plot(v0,h,'g');
hold on;

u=0.22;
E0=1./2.*m.*v0.*v0;
Eb=(E0*qq)-qqq*(u*qq-2*u*u-1)/(4*u*u+1);
E_edge=(1-2*u*L/r)*((Eb-u*m*g*S)*qq-((2*u*u+1)*qq+u)*qqq/(4*u*u+1))-m*g*L*cosd(a);
h=(1/2*k)*log(1+2*k*E_edge/m/g);
plot(v0,h,'b');
hold on;

legend('u=0.02','u=0.12','u=0.22');			%图例标注
end