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=