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)+