Matlab隐式情节

Question

如何根据以下等式绘制 h / L 和 h / L0 之间的曲线：

  

[（h / L）{（（Dk ^ 4 / g）+ 1）tanh（kh）}] ^（1/2）=（h / L0）^（1/2）{1-（TU /小时）（H / L）}

我自己尝试了但是我无法获得正确的曲线。我的代码如下：

D=10^11;
g=9.8;
T=1;
U=5;
h=1;
a=T*U/h;
L= (-3:0.01:3);
L0=(-3:0.01:3);
y=h./L;
x=h./L0;
k=(2*pi/h).*y;
y= (x.*(1-a.*y).^2)./(((16*D*(pi^4)*(y.^4)/(g*h^4))+1).*tanh(2*pi.*y));
plot(x,y,'b');
xlabel('h/L0');
ylabel('h/L');

此处g是由重力引起的加速度，k是波数。 D是刚性因子，等于 D = EI / d 其中 d 是密度，而 EI 是刚性和 k = 2pi / L 。