%%%離線驗證程序
%%%10和15修正、選擇、arma、SVM_10_15和ARMA平均、SVM_c和ARMA平均
%%%by?7/24

%???clc;clear?all;close?all;
%?%load?model_power;
%?%?load?model_power;
Q?=?20;%訓練天數(shù)
Q2?=?100;%arma的訓練數(shù)據
sp_in?=?2.25;
sp_out?=?25;
power_max?=?1504.15087;
power_min?=?10.20950758;
string=‘1700270A‘;
excel_name?=?‘tianjindagang.xlsx‘;
fengji_num?=?33;%---要預測的表格的風機數(shù)
day_excel?=?54;%---表格里面的天數(shù)
qsrq=1;
qsyue=3;
ri?=20;%%要預測的日期
yue=4;
sample_lag=5;
flag_nwp1?=?ones（311）;%3.25（18）3.26（12），4.7（18）4.224.234.24
flag_nwp2?=?ones（651）;
flag_nwp?=?[flag_nwp1;flag_nwp2];
Hz?=?[];

%?[N1TX]?=?xlsread（‘環(huán)港00風速1.xlsx‘）;
%?[N2TX]?=?xlsread（‘環(huán)港00風速2.xlsx‘）;
%?[N3TX]?=?xlsread（‘環(huán)港00風速3.xlsx‘）;
%?speed_nwp_1=（N1（293:98））‘;
%?speed_nwp_2=（N2（293:98））‘;
%?speed_nwp_3=（N3（293:98））‘;
%獲取數(shù)據
[rangerange_nwp]?=?excel_range（1fengji_num6day_excel）;
[NTX]?=?xlsread（excel_namerange）;%讀取平均風速、平均功率、平均溫度、風機數(shù)
[N_kongTX]?=?xlsread（excel_name[‘A1:A‘?num2str（day_excel*288-2）]）;
DATE?=?TX（:1）;
count_N?=?size（N）;
[rangerange_nwp]?=?excel_range（1+8fengji_num6day_excel）;%讀取18點的nwp
[range_zsrange_nwp_zs]?=?excel_range（1+9fengji_num6day_excel）;%讀取18點的nwp

[Nwp_288TXwp_288]?=?xlsread（excel_namerange_nwp）;
[Nwp_288_zsTXwp_288_zs]?=?xlsread（excel_namerange_nwp_zs）;
Nwp_288_zs=Nwp_288_zs/fengji_num;
?%要將nwp從288點變成96點
???count?=?size（Nwp_288）;
???Nwp?=?Nwp_288（1:3:count（1））;
???count_zs?=?size（Nwp_288_zs）;
???Nwp_zs?=?Nwp_288_zs（1:3:count（1））;

N_cc?=?[N（:count_N（2））N（:count_N（2）-3）N（:count_N（2）-2）N（:count_N（2）-1）];
[Xname1]=powercurvechazhi（string）;
v?=?1;
for?i=-47:240
if?yue==3||yue==4
read_datetime?=?char（DATE（（（yue-qsyue）*31+（ri-qsrq））*288+i1））;
else?
read_datetime?=?char（DATE（（31+30+（ri-qsrq））*288+i1））;
end
?[time_arrayfilename_excelfilename_txttxt_time]?=?dateparting_ceshi（read_datetime）;

%%判斷時刻
run_flag?=?is_15minutes（time_array）;

if?run_flag?==?1
?????if?yue==3
????????time=ri-（qsrq-1）-1;
????elseif?yue==4
??????time=31-（qsrq-1）+ri-1;
????elseif?yue==5
??????time=30+31-（qsrq-1）+ri-1;
????end
????if?time_array（3）==ri
????????n=time;
????else
????????
????????%%提取風機運行數(shù)量
????????n=time-1;
????end
????m=12*time_array（4）+time_array（5）/5+1+1;%預測開始的時間
????run_num?=?N（n*288+m-1count_N（2））;%%現(xiàn)在時刻的運行臺數(shù)?
????sp_5?=?N（n*288+m-1-2count_N（2）-3）;
????sp_10?=?N（n*288+m-1-1-47:n*288+m-1-1count_N（2）-3）;
????sp_15?=?N（n*288+m-1-47:n*288+m-1count_N（2）-3）;
????
????Wspc?=?N_cc（（m+288*n）:（m+288*n+47）2）;%真實風速
????sp_nwp_48=Nwp_288（（m+288*n）:（m+288*n+47）1）;
????Power_nwp_48_zs=Nwp_288_zs（（m+288*n）:（m+288*n+47）1）;
????F_Power?=?N_cc（（m+288*n）:（m+288*n+47）3）;%真實發(fā)電量
????D_pre=N_cc（（m+288*n-48）:（m+288*n-1）2）;?
????D_pre_qushi?=??D_pre?-?detrend（D_pre）;
????D_pre_gui?=?maxmin（D_pre）;
????z_b2?=?N_cc（m-9+288*n:m-1+288*n2）;%前9點
?%%預測?
%?????[b_2sp_15_averF_PowerD1]?=?predictspeed_ceshi_zuichu（QN_cctime_arrayqsyueyueqsrqrisample_lag）;
[b_2D1bn16_1_D_prebn16_1_D_pre_gui]?=?net_shenjing（

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

?????文件??????12725??2012-08-13?22:26??awayline_SVM_ARMA.m

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

????????????????12725????????????????????1