小波预测风电功率

clear
vm=xlsread（‘C:\Users\Administrator\Desktop\123.xls‘）;    %取原始风速数据200个 
%用db3小波进行3层分解
[cl]=wavedec（vm3‘db4‘）;
%重构第1-5层逼近系数.
a3 = wrcoef（‘a‘cl‘db4‘3）;
%重构第1-3层细节系数
d3 = wrcoef（‘d‘cl‘db4‘3）;
d2 = wrcoef（‘d‘cl‘db4‘2）;
d1 = wrcoef（‘d‘cl‘db4‘1）;
d=[a3d3d2d1];

totalnumber=200;
totalnumber2=200;
%totalnumber2=8560;

v=d（1:totalnumber1:4）;
v1=d（totalnumber+1:totalnumber+totalnumber21:4）;


for m=1:4
%1标准化处理
average_v=mean（d（1:endm））; %平均值
var_v=var（d（1:endm））;      %方差
var_v_s=sqrt（var_v）;        %标准差
for i=1:totalnumber
     DX（i1）=（v（im）-average_v）/var_v_s;
end
vs=DX;
%2标准化处理
for i=1:totalnumber2
     DX1（i1）=（v1（im）-average_v）/var_v_s;
end
vs1=DX1（1:totalnumber21）;

z=iddata（vs）;                       %格式转换
z1=iddata（vs1）;
for p=1:10
    q=p-1;
    n=armax（z[p q]）;
    AIC=aic（n）;
end
for k=1:10
    if AIC（k）==min（AIC）
        p_final=k;
        q_final=p_final-1;
        break
    end
end

m_final=armax（z[p_final q_final]）;   %确定最终模型
%预测过程
p=