power_m1.m

脉动风速功率谱估计

     WPower_u_1=zeros(144,301);
     WPower_v_1=zeros(144,301);
     WPower_u_2=zeros(144,301);
     WPower_v_2=zeros(144,301);
     WPower_u_3=zeros(144,301);
     WPower_v_3=zeros(144,301);
     WPower_u_4=zeros(144,301);
     WPower_v_4=zeros(144,301);
     WPower_u_5=zeros(144,301);
     WPower_v_5=zeros(144,301);
     WPower_u_6=zeros(144,301);
     WPower_v_6=zeros(144,301);
     WPower_u_7=zeros(144,301);
     WPower_v_7=zeros(144,301);
     WPower_u_8=zeros(144,301);
     WPower_v_8=zeros(144,301);
     WPower_u_9=zeros(144,301);
     WPower_v_9=zeros(144,301);
     WPower_u_10=zeros(144,301);
     WPower_v_10=zeros(144,301);

for i=1:1:144
[x1,f] = pwelch(u_1(i,:),hamming(512),10,600,1);
x2 = pwelch(v_1(i,:),hamming(512),10,600,1);
y1 = pwelch(u_2(i,:),hamming(512),10,600,1);
y2 = pwelch(v_2(i,:),hamming(512),10,600,1);
z1 = pwelch(u_3(i,:),hamming(512),10,600,1);
z2 = pwelch(v_3(i,:),hamming(512),10,600,1);
s1 = pwelch(u_4(i,:),hamming(512),10,600,1);
s2 = pwelch(v_4(i,:),hamming(512),10,600,1);
h1 = pwelch(u_5(i,:),hamming(512),10,600,1);
h2 = pwelch(v_5(i,:),hamming(512),10,600,1);
w1 = pwelch(u_6(i,:),hamming(512),10,600,1);
w2 = pwelch(v_6(i,:),hamming(512),10,600,1);
k1 = pwelch(u_7(i,:),hamming(512),10,600,1);
k2 = pwelch(v_7(i,:),hamming(512),10,600,1);
d1 = pwelch(u_8(i,:),hamming(512),10,600,1);
d2 = pwelch(v_8(i,:),hamming(512),10,600,1);     
e1 = pwelch(u_9(i,:),hamming(512),10,600,1);
e2 = pwelch(v_9(i,:),hamming(512),10,600,1);
r1 = pwelch(u_10(i,:),hamming(512),10,600,1);
r2 = pwelch(v_10(i,:),hamming(512),10,600,1);     
  for j=1:1:301
    WPower_u_1(i,j)=x1(j);
    WPower_v_1(i,j)=x2(j);
    WPower_u_2(i,j)=y1(j);
    WPower_v_2(i,j)=y2(j);
    WPower_u_3(i,j)=z1(j);
    WPower_v_3(i,j)=z2(j);
    WPower_u_4(i,j)=s1(j);
    WPower_v_4(i,j)=s2(j);
    WPower_u_5(i,j)=h1(j);
    WPower_v_5(i,j)=h2(j);
    WPower_u_6(i,j)=w1(j);
    WPower_v_6(i,j)=w2(j);
    WPower_u_7(i,j)=k1(j);
    WPower_v_7(i,j)=k2(j);
    WPower_u_8(i,j)=d1(j);
    WPower_v_8(i,j)=d2(j);
    WPower_u_9(i,j)=e1(j);
    WPower_v_9(i,j)=e2(j);
    WPower_u_10(i,j)=r1(j);
    WPower_v_10(i,j)=r2(j);
  
  end
end
plot(f, WPower_u_1(1,:),f,WPower_v_1(1,:))