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))—风速功率谱的对数');
