matlab实现BP预测数据

基于matlab，读取txt文本数据，对文本数据进行处理，提取特征，然后基于BP神经网络进行数据预测。

%  将输入的数据进行转化为统计数据
close all;clear;
%% 加载数据 
path_train = ‘TrainData_2015.1.1_2015.2.19.txt‘;
path_test = ‘TestData_2015.2.20_2015.2.27.txt‘;
[traind_flavortrain_time] = readtxt（path_train）;
[testd_flavortest_time] = readtxt（path_test）;

%% 对数据进行预测
% % 当所有时刻中虚拟机为0的规格，预测的结果是这些数据也都是为0
sum_flavor1 = sum（train_flavor）;
% sum_flavor2 = sum（test_flavor）;
% ind_zero = find（sum_flavor1==0）;
% ind_pre（1:ind_zero） = 0;


% 对不为0的虚拟机规格进行数量预测
data_nuum = length（traind_flavor（:1））;
traind = traind_flavor（1:data_nuum - 21:）;
trainl = traind_flavor（22:data_nuum:）;

% data_nuum2 = length（testd_flavor（:1））;
% testd = testd_flavor（1:data_nuum2 - 6:）;
% testl = testd_flavor（7:data_nuum2:）;

imp=emd（sum_flavor1）;
[Aft] = hhspectrum（imp）;
% 求瞬时频率
%     infreq = fs * f;
    
    %求信号的边界谱bjp，E：对应的振幅值
    [Ett1]=toimage（Aftlength（t））;
    enery = E;
    E=flipud（E）;
    band = 0;   % 对应于边际谱值大于1的带宽
    num = 0;    % 用于寻找以0.3位阈值的截止位置
    for k=1:size（E1）
        bjp（k）=sum（E（k:））;
        if bjp（k）>1
            band = band+1;
        end
        if bjp（k）>0.5
            num = num+1;
        end
        if num == 20
            stop = k-20;
        end
    end
    
    area = polyarea（1:size（E1）bjp）;       %求频谱面积
    
    
    
%% 测试
test_out=sim（nettraind_flavor（（data_nuum - 20）:data_nuum:）‘）;
test_out = test_out‘;
test_out（test_out<10）=0;
test_out = round（test_out）;

%% 数据评价

N = 15;
for i =1:7
    s1（i） = sqrt（sum（（sum（test_out（i:）） - testd_flavor（i:））.^2）/N）;
    s2（i） = sqrt（sum（test_out（i:）.^2））;
    s3（i） = sqrt（sum（testd_flavor（i:）.^2））;
    score（i） = 1- （s1/（s2+s3））;
end

% BP的效果很差

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       8459  2018-01-10 16:42  matlab_BP\data\data_2015_1.txt

     文件      49174  2018-01-10 16:42  matlab_BP\data\data_2015_12.txt

     文件      12074  2018-02-27 20:47  matlab_BP\data\data_2015_2.txt

     文件       8903  2018-02-10 21:16  matlab_BP\data\data_2015_3.txt

     文件      13885  2018-02-10 21:16  matlab_BP\data\data_2015_4.txt

     文件      19945  2018-02-10 21:16  matlab_BP\data\data_2015_5.txt

     文件      36591  2018-01-10 16:42  matlab_BP\data\data_2016_1.txt

     文件      33224  2018-03-13 22:04  matlab_BP\data.emf

     文件      18778  2018-03-13 22:04  matlab_BP\data.fig

     文件      35997  2018-03-13 22:07  matlab_BP\data.png

     文件       1815  2018-03-13 22:12  matlab_BP\prediction.m

     文件       1430  2018-03-13 21:27  matlab_BP\prediction_BP.m

     文件       1671  2018-03-13 19:38  matlab_BP\readtxt.m

     文件       2643  2018-01-16 14:53  matlab_BP\TestData_2015.2.20_2015.2.27.txt

     文件      16125  2018-01-16 14:53  matlab_BP\TrainData_2015.1.1_2015.2.19.txt

     目录          0  2018-03-13 19:52  matlab_BP\data

     目录          0  2018-03-16 11:13  matlab_BP

----------- ---------  ---------- -----  ----

               260714                    17