Davenport風(fēng)速譜

%表面阻力系數(shù)K=0.02，V10=28.4m/s,f——頻率，Svf——Davenport風(fēng)速譜
K=0.02;V10=28.4;
n=linspace(-3,2,101);
f=zeros(1,101);
for i=1:101
    f(i)=10^n(i);
end 
Svf=zeros(1,101);
for i=1:101
    m=(1200*f(i)/V10)^2;
    Svf(i)=4*K*V10^2*m/(f(i)*(1+m)^(4/3));
end

 for i=1:101
n(i)=n(i)+n(i)*1/4.5*abs(n(i)-2);
 end

plot(n,log(Svf));
xlabel('log(f)—頻率的對(duì)數(shù)');
ylabel('log(S(n))—風(fēng)速功率譜的對(duì)數(shù)');

基于Davenport風(fēng)速譜單點(diǎn)時(shí)程模擬

分享離散頻譜的能量重心校正法程序