基于ARMR和白噪声特性模型及风速威布尔分布研究（Matlab代码实现）

df=1/T;%频率分辨率
f=(0:n-1)*df;
length_f=length(f);
for i=1:length_f
    S(i)=4*(sigma_u)^2*L/average_v/(1+70.8*(f(i)*L/average_v)^2)^(5/6);%Von Karman谱功率密度函数
end
R=real(ifft(S));%自相关与功率谱密度是傅立叶变换关系，求自相关
%求解AR模型参数
a(1)=R(2)/R(1);
a(2)=R(3)/R(1)-R(2)*R(2)/R(1)/R(1);
%求解MA模型参数
Ry=(-a(1)+a(2)*a(1))*R(1)+(1+a(1)^2+a(2)^2)*R(2)+(-a(1)+a(1)*a(2))*R(3)+(-a(2))*R(4);
b=Ry/sigma_noise;%beta=Ry/白噪声的方差
%脉动风速
v(1)=u(1);
v(2)=a(1)*v(1)+u(2)+b*u(1);
for k=3:n
    v(k)=a(1)*v(k-1)+a(2)*v(k-2)+u(k)+b*u(k-1);%ARMA模型
end
%风速
for k=1:n
    x(k)=average_v+v(k);
end