no_gaus.m

脉动风速功率谱估计

% 峰度
for i=1:1:144
   Kur_u_1(i)=kurtosis(u_1(i,:));
   Kur_v_1(i)=kurtosis(v_1(i,:));
end

for i=1:1:144
   Kur_u_2(i)=kurtosis(u_2(i,:));
   Kur_v_2(i)=kurtosis(v_2(i,:));
end

for i=1:1:144
   Kur_u_3(i)=kurtosis(u_3(i,:));
   Kur_v_3(i)=kurtosis(v_3(i,:));
end

for i=1:1:144
   Kur_u_4(i)=kurtosis(u_4(i,:));
   Kur_v_4(i)=kurtosis(v_4(i,:));
end

for i=1:1:144
   Kur_u_5(i)=kurtosis(u_5(i,:));
   Kur_v_5(i)=kurtosis(v_5(i,:));
end

for i=1:1:144
   Kur_u_6(i)=kurtosis(u_6(i,:));
   Kur_v_6(i)=kurtosis(v_6(i,:));
end

for i=1:1:144
   Kur_u_7(i)=kurtosis(u_7(i,:));
   Kur_v_7(i)=kurtosis(v_7(i,:));
end

for i=1:1:144
   Kur_u_8(i)=kurtosis(u_8(i,:));
   Kur_v_8(i)=kurtosis(v_8(i,:));
end

for i=1:1:144
   Kur_u_9(i)=kurtosis(u_9(i,:));
   Kur_v_9(i)=kurtosis(v_9(i,:));
end

for i=1:1:144
   Kur_u_10(i)=kurtosis(u_10(i,:));
   Kur_v_10(i)=kurtosis(v_10(i,:));
end

%斜度

for i=1:1:144
   ske_u_1(i)=skewness(u_1(i,:));
   ske_v_1(i)=skewness(v_1(i,:));
end

for i=1:1:144
   ske_u_2(i)=skewness(u_2(i,:));
   ske_v_2(i)=skewness(v_2(i,:));
end

for i=1:1:144
   ske_u_3(i)=skewness(u_3(i,:));
   ske_v_3(i)=skewness(v_3(i,:));
end

for i=1:1:144
   ske_u_4(i)=skewness(u_4(i,:));
   ske_v_4(i)=skewness(v_4(i,:));
end

for i=1:1:144
   ske_u_5(i)=skewness(u_5(i,:));
   ske_v_5(i)=skewness(v_5(i,:));
end

for i=1:1:144
   ske_u_6(i)=skewness(u_6(i,:));
   ske_v_6(i)=skewness(v_6(i,:));
end

for i=1:1:144
   ske_u_7(i)=skewness(u_7(i,:));
   ske_v_7(i)=skewness(v_7(i,:));
end

for i=1:1:144
   ske_u_8(i)=skewness(u_8(i,:));
   ske_v_8(i)=skewness(v_8(i,:));
end

for i=1:1:144
   ske_u_9(i)=skewness(u_9(i,:));
   ske_v_9(i)=skewness(v_9(i,:));
end

for i=1:1:144
   ske_u_10(i)=skewness(u_10(i,:));
   ske_v_10(i)=skewness(v_10(i,:));
end