%從txt中統(tǒng)計(jì)數(shù)據(jù)，每個(gè)月單獨(dú)統(tǒng)計(jì)，用的前5個(gè)月份的txt
%最終統(tǒng)計(jì)結(jié)果放在data下，data的第n列表示第n型號(hào)設(shè)備每天的數(shù)量
clear;
clc;
data1=readtxt（1‘data_2015_1.txt‘）;
data2=readtxt（2‘data_2015_2.txt‘）;
data3=readtxt（3‘data_2015_3.txt‘）;
data4=readtxt（4‘data_2015_4.txt‘）;
data5=readtxt（5‘data_2015_5.txt‘）;
data=[data1;data2;data3;data4;data5];

%下邊處理都是用的8型號(hào)的數(shù)據(jù)
%剔除壞值，保證數(shù)據(jù)在均值上下3個(gè)標(biāo)準(zhǔn)差范圍內(nèi)
data8=data（:8）.‘;
len=length（data8）;
result=removedata（data8）;
%濾波
?Wc=0.9;??????????%數(shù)字帶寬????????????????????????????????
?[ba]=butter（4Wc）;
?result_filter=filter（baresult）;
?figure（3）
?plot（1:length（result_filter）result_filter‘.-r‘）
?hold?on
?plot（1:length（result_filter）result‘.-g‘）
?hold?off
?
?%arma模型定階，其中的ij表示的是pq
?s?=?1;?%周期是1
x?=?result_filter;%初始數(shù)據(jù)的錄入，把濾波之后的數(shù)據(jù)作為輸入
n?=?14;?%預(yù)報(bào)的個(gè)數(shù)
m1?=?length（x）;?%原始的數(shù)據(jù)的個(gè)數(shù)
k=0;
w?=?x;?
m2?=?1;
for?i?=?1:9
????for?j?=?1:2
????????ToEstMd?=?arima（i0j）;?%指定模型的結(jié)構(gòu)
????????k?=?k?+?1;
????????R（k）?=?i;
????????M（k）?=?j;
????????[EstMdEstParamCovLogLinfo]?=?estimate（ToEstMdw‘）;%模型擬合
????????numParams?=?sum（any（EstParamCov））;%計(jì)算擬合參數(shù)的個(gè)數(shù)
????????[aic（k）bic（k）]?=?aicbic（LogLnumParamsm2）;%計(jì)算aic，bic
????end
end
fprintf（‘RMAICBIC的對(duì)應(yīng)值如下\n%f‘）;%顯示計(jì)算結(jié)果
check??=?[R‘M‘a(chǎn)ic‘bic‘]

%根據(jù)上面結(jié)果，p，q為93時(shí)aic，bic較小，用arma（9，3）做預(yù)測(cè)
www=result_filter（1:40）;
toEstMd?=?arima（903）;%構(gòu)建arma（93）模型
[EstMdEstParamCovLogLinfo]?=?estimate（toEstMdwww‘）;%將上面模型用在www上
f=?forecast（EstMd14‘Y0‘www‘）;%對(duì)www進(jìn)行預(yù)測(cè)，預(yù)測(cè)14個(gè)點(diǎn)
figure（4）
plot（1:54result_filter（1:54）‘.-‘）
hold?on
plot（41:54f‘.-r‘）
hold?off
grid?on

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????1758??2018-03-11?12:41??test\codecraft.m

?????文件????????703??2018-03-10?17:20??test\data.mat

?????文件???????8459??2018-01-10?16:42??test\data_2015_1.txt

?????文件??????49174??2018-01-10?16:42??test\data_2015_12.txt

?????文件??????12074??2018-02-27?20:47??test\data_2015_2.txt

?????文件???????8903??2018-02-10?21:16??test\data_2015_3.txt

?????文件??????13885??2018-02-10?21:16??test\data_2015_4.txt

?????文件??????19945??2018-02-10?21:16??test\data_2015_5.txt

?????文件??????36591??2018-01-10?16:42??test\data_2016_1.txt

?????文件????????448??2018-03-09?21:22??test\readtxt.m

?????文件????????292??2018-03-10?22:19??test\removedata.m

?????目錄??????????0??2018-03-11?12:42??test

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

???????????????152232????????????????????12