gusucode.com > 混沌时间序列进行深度学习matlab编程 程序源码 > Main_Volterra.m
% 混沌时间序列的 volterra 预测(一步预测) -- 主函数 % 使用平台 - Matlab6.5 / Matlab7.0 % 作者:陆振波,海军工程大学 % 欢迎同行来信交流与合作,更多文章与程序下载请访问我的个人主页 % 电子邮件:luzhenbo@yahoo.com.cn % 个人主页:http://luzhenbo.88uu.com.cn clc clear close all %-------------------------------------------------- % 公共参数 k1 = 6000; % 前面的迭代点数 k2 = 2000; % 后面的迭代点数 (总样本数) %-------------------------------------------------- % 产生混沌序列 sigma = 10; % Lorenz 方程参数 a b = 8/3; % b r = 34; % c y = [-1,0,1]; % 起始点 (1 x 3 的行向量) h = 0.01; % 积分时间步长 z = LorenzData(y,h,k1+k2,sigma,r,b); X = z(k1+1:end,1); %-------------------------------------------------- X = normalize_a(X,1); % 信号归一化到均值为0,振幅为1 tau = 1; % 时延 m = 3; % 嵌入维数 p = 3; % Volterra阶数 %-------------------------------------------------- train_num = 500 % 训练样本数 test_num = 1500 % 测试样本数 %-------------------------------------------------- % 混沌序列的相空间重构 (phase space reconstruction) x_train = X(1:train_num); x_test = X(train_num+1:train_num+test_num); [xn_train,dn_train] = PhaSpaRecon(x_train,tau,m); [xn_test,dn_test] = PhaSpaRecon(x_test,tau,m); %------------------------------------------------------ % 训练与预测 % 训练算法求权矢量 [Wn,err_mse1] = volterra_train_lu(xn_train,dn_train,p); err_mse1_std = err_mse1/var(X) % 测试样本的一步预测 dn_pred = volterra_test(xn_test,p,Wn); err2 = dn_test - dn_pred; err_mse2 = sum(err2.^2)/length(err2); err_mse2_std = err_mse2/var(X) %---------------------------------------------------- % 结果作图 figure; subplot(211); plot(1:length(err2),dn_test,'r+:',1:length(err2),dn_pred,'bo-'); title('真实值(+)与预测值(o)') subplot(212); plot(err2,'k'); title('预测绝对误差')