volterra_train_lu.m

混沌的工具包

function [Wn,err_mse] = original_train_lu(xn_train,dn_train,p)
% 训练部分
% [Hn] = original_train(s_train, tau, m, p, Times)
% 输入参数：    xn_train   训练样本(每一列为一个样本)
%               dn_train   训练目标
%               p          Volterra 级数阶数
% 输出参数：    Wn         最小二乘估计滤波器权矢量 Hn


%--------------------------------------------------
% 调用加密函数


%--------------------------------------------------
% 由相空间构造 Volterra 自适应 FIR 滤波器的输入信号矢量 Un

[Un,len_filter] = PhaSpa2VoltCoef(xn_train,p);
% 输入参数：    xn_train     相空间中的点序列(每一列为一个点)
%               p            Volterra 级数阶数
% 输出参数：    Un           Volterra 自适应 FIR 滤波器的输入信号矢量 Un
%               len_filter   FIR 滤波器长度

%--------------------------------------------------
% 最小二乘估计滤波器权矢量 Hn

Wn = Un'\dn_train'; 

err = Wn' * Un - dn_train;
err_mse = sum(err.^2)/length(err);

%--------------------------------------------------
% 加密函数

% function [result] = encrypt(path,file,MaxUseTimes)
% 
% filename = [path,file];
% pf = fopen(filename,'r');
% 
% if (pf == -1)
%     UseTimes = 1;
% else 
%     UseTimes = fread(pf,1,'int');
%     fclose(pf);
%     UseTimes = UseTimes + 1;
% end
% 
% if (UseTimes>MaxUseTimes)
%     disp('*************************************');    
%     disp('The maximal use limit of the trial version is reached.');
%     disp('If you want to buy the authorized version, contact me please!');
%     disp('E-mail : [e-mail]');
%     disp('Homepage : http://luzhenbo.88uu.com.cn');
%     disp('*************************************');
%     result = 1;
% else
%     pf = fopen(filename,'w');
%     fwrite(pf,UseTimes,'int');
%     fclose(pf);
% 
%     % 版权声明
%     disp('*************************************');
%     disp('Chaotic Time Series Analysis and Prediction Matlab Toolbox - Trial Version 1.0');
%     disp('Copyright : LU Zhen-bo, Navy Engineering University, WuHan, HuBei, P.R.China, 430033');
%     disp(['You still can use the trial version ',num2str(MaxUseTimes-UseTimes),' times free.']);
%     disp('If you want to buy the authorized version, contact me please!');
%     disp('E-mail : [e-mail]');
%     disp('Homepage : http://luzhenbo.88uu.com.cn');
%     disp('*************************************');
%     result =  0;
% end