m=pi/40;
k=0.0075;
g=10;
t=2;
h=7;
a=exp(k*h/m);
for i=1:10
t=sqrt(m/(k*g))*log10(a+sqrt(a^2-1))/log10(exp(1));
end
for L=1:0.1:10
for t=0.5:0.05:2
v=10;
for i=1:10
v=exp(k/m*(L-m/k*log10(k/m*t+1/v)/log10(exp(1)))); %迭代法计算速度
hold on
mesh(L,t,v)
end
end
end