karman谱拟合

%读取一列风速值，求出湍流积分尺度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