clear;
clc;
tic

year=2
%模擬的年限

for?l=1:year;
????
????%%%%%%%%%%%%%%%%%%%%%%%%%??第一步：數據導入與預處理??%%%%%%%%%%%%%%%%%%%
????
????SW0=load（‘windspeed.txt‘）;????%載入原始風速數據
????SW0=SW0‘;?????????????????????%將風速數據從橫列轉置成豎列
????SW0=SW0/10*36;????????????????%原始數據的風速單位為m/s，這里轉化為km/h
????N=size（SW02）;????????????????%把SW0的列數賦值給N
????mu=mean（SW0）;?????????????????%求SW0的平均值
????sigma=var（SW0）;
????sigma=sigma^0.5;??????????????%求SW0的標準差
????y=（SW0-mu）./sigma;????????????%數據預處理
????
????N1=8736;??????????????????????%N1為模擬的小時數
????WTGnum=25;????????????????????%WTGnum為風力發電機數量
????
????%%%%%%%%%%%%%%??第二步：根據AIC準則確定ARMA模型的階數??%%%%%%%%%%%%%%%%%%
????
????for?n=2:7;
????????m=armax（y‘[nn-1]）;
????????fai=-m.a;
????????theta=m.c;????????????????%把armax函數得到的參數，取出來
????????
????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????251??2017-04-09?20:59??windspeed.txt

?????文件??????69317??2013-01-19?19:43??load.txt

?????文件??????11176??2017-04-17?17:03??wind.m

-----------?---------??----------?-----??----

????????????????80744????????????????????3