实现计算波浪力学中波长的超越方程函数，matlab程序，可以直接调用.zip

function[L]=ldis(d,T) %利用弥散关系求波长 %L=波长 %d=水深 %T=波周期 g=9.8; %重力加速度m/s^2 sigma=2*pi/T; y=sigma^2*d/g; % d1=0.6666666666; d2=0.3555555555; d3=0.1608465608; d4=0.0632098765; d5=0.0217540484; d6=0.0065407983; % Tem=y/(1+d1*y+d2*y^2+d3*y^3+d4*y^4+d5*y^5+d6*y^6)+y^2; kd=sqrt(Tem); k=kd/d; L1=2*pi/k; tol=0.0001; %精度 Lo=g*T^2/(2*pi); L=Lo*tanh(2*pi*d/L1); while abs(L/L1-1)>tol; L1=(L+L1)/2; L=Lo*tanh(2*pi*d/L1); end