function [ delay Err ] = delay_ergodic( high_order, delay_max, input, output )
%DELAY_ERGODIC 此处显示有关此函数的摘要
%   遍历delay参数的可行域，找到使loss最小的参数
%   首先实现MISO的形式
%   此处显示详细说明
Nu = size(input,2);
Ny = size(output,2);
N = size(output,1);
n=high_order;
na = n;
Err=[];
nb=n*ones(1,Nu);
delay=zeros(Ny,Nu);
range=delay_max*ones(1,Ny);

for i=1:Ny
    fprintf('output%d delay selection :',i);
    for j1=1:range(i)
        for j2=1:range(i)
            for j3=1:range(i)
            delay(i,j)=d % for test
            nk = delay(i,:) + 1;
            TH=arx([output(:,i),input],[na,nb,nk]);
            ysim=idsim(TH,input);
            Err(j,i,d)=cov(ysim-output(:,i))/cov(output(:,i));
            FOE(j,i,d)=(N+(2+2*Nu)*d)/(N-(2+2*Nu)*d)*cov(output(:,i)-ysim);
        end
    end
%     [err_FOE,delay(i)]=min(FOE(:,i,:));
%     fprintf(' output%d delay = %d \n',i,delay(i));
end

end