da_olt3.m

来自「用VC++编写的遗传算法源程序」· M 代码 · 共 33 行

function [ypred,x]=da_olt2(x,xold,coefs,filt,poly);
%
% Non-linear dynamic model
%

%
% Ensure the coefficients are in a column
%
[D L]=size(coefs);
if L > D
	coefs=coefs';
end

%
% Polynomialise the data
%

%
% Filter the data:
%
[D L]=size(x);
for i = 1:L
	x(1,i)=(xold(1,i)*filt(i))+(x(1,i)*(1-filt(i)));
end


%
% Do model
%
ypred=(x.^poly)*coefs;