仅提供参考!
Ee=@(t)(2.94+1)/2+(2.94-1)/2*(1+12/t).^(-0.5)
func=@(t)(t<=1).*(60*(log(8./t+t/4))./sqrt(Ee(t)))+(t>=1).*(120*pi./(t+1.393+0.667*log(t+1.444))./sqrt(Ee(t)))
fplot(func,[0,8])
xlabel('t');ylabel('Z0')


Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!