%從txt中統計數據，每個月單獨統計，用的前5個月份的txt
%最終統計結果放在data下，data的第n列表示第n型號設備每天的數量
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型號的數據
%剔除壞值，保證數據在均值上下3個標準差范圍內
data8=data（:8）.‘;
len=length（data8）;
result=removedata（data8）;
%濾波
?Wc=0.9;??????????%數字帶寬????????????????????????????????
?[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;%初始數據的錄入，把濾波之后的數據作為輸入
n?=?14;?%預報的個數
m1?=?length（x）;?%原始的數據的個數
k=0;
w?=?x;?
m2?=?1;
for?i?=?1:9
????for?j?=?1:2
????????ToEstMd?=?arima（i0j）;?%指定模型的結構
????????k?=?k?+?1;
????????R（k）?=?i;
????????M（k）?=?j;
????????[EstMdEstParamCovLogLinfo]?=?estimate（ToEstMdw‘）;%模型擬合
????????numParams?=?sum（any（EstParamCov））;%計算擬合參數的個數
????????[aic（k）bic（k）]?=?aicbic（LogLnumParamsm2）;%計算aic，bic
????end
end
fprintf（‘RMAICBIC的對應值如下\n%f‘）;%顯示計算結果
check??=?[R‘M‘aic‘bic‘]

%根據上面結果，p，q為93時aic，bic較小，用arma（9，3）做預測
www=result_filter（1:40）;
toEstMd?=?arima（903）;%構建arma（93）模型
[EstMdEstParamCovLogLinfo]?=?estimate（toEstMdwww‘）;%將上面模型用在www上
f=?forecast（EstMd14‘Y0‘www‘）;%對www進行預測，預測14個點
figure（4）
plot（1:54result_filter（1:54）‘.-‘）
hold?on
plot（41:54f‘.-r‘）
hold?off
grid?on

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

?????文件???????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