Davenport风速谱

%表面阻力系数K=0.02，V10=28.4m/s,f——频率，Svf——Davenport风速谱
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)—频率的对数');
ylabel('log(S(n))—风速功率谱的对数');

基于Davenport风速谱单点时程模拟

分享离散频谱的能量重心校正法程序