% main_dsr_harxmr_chk

randn('seed',0); rand('seed',0)   % Setting seeds in prbs1 function.
R=prbs1(N,75,175);                % The reference signal

randn('seed',0)

L=2: J=2; g=0;
M=10;
for j=1:M
    
    [Y,U]=fun_clop_data(R);
    
    [a1,b1,d1,e1,c1,f1]=dsr(Y,U,L,g,J,1,2);
    hd=ddcgain(a1,b1,d1,e1)
end