%讀取一列風速值，求出湍流積分尺度Lu，然后繪制出實測風速的功率譜，和理論計算的功率譜
clear
clc
?a=load（‘data_1.txt‘）;?
ave=mean（a）;a1=a-ave;[mn]=size（a）;???%?ave=平均風速，a1=脈動風速，m，n分布是行列
temp=zeros（2*mn）;a2=temp;R=zeros（mn）;???
a2（1:m）=a1;
for?i=1:m+1
????temp（i:i+m-1）=a1;
????R（i）=sum（temp.*a2）;
????temp=zeros（2*mn）;
end?
R=R/m;
R1=sum（a1.*a1）/m;
pt=R/R1;
Lu=2.6*36*sum（pt（1:3065））;

N=72648;F=1/36;U=2.6;u_10=2.17;
x=a;???????????????%文檔第一列是時間數據，第二列才是電壓數據
y=fft（x）;
mag=abs（y）;???%*2/72648;????%看例子里面都是這么寫的。。
mag=mag（1:N/2）;
P