channel_est.m

用于模拟时变非平稳的ARMA过程

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm
alpha= 1/2;
load updown_h_rusk.mat;
N= 2048;
M0= 28;
%figure(1);mesh(abs(h_I2))

MCH= 3;
Mmax= 1;
Lmax= 1;

hnm= h_I2(113:end, M0:M0+MCH);
hnm= real(hnm/max(max(abs(hnm))));

%figure(2);mesh(abs(hnm))
AMBI= {};
for tau= 0:MCH
   for taup= 0:MCH
	 [tau taup]
         Rttp= corr_est(hnm(:, tau+1), hnm(:, taup+1), Mmax, alpha);
         Attp= fft(Rttp);
         Attp= [Attp(N/2+1:N, :); Attp(1:N/2, :)];
%         mesh(abs(Attp(N/2+1-3*Lmax:N/2+1+3*Lmax, :)));drawnow
         AMBI{tau+1, taup+1}= Attp(N/2+1-3*Lmax:N/2+1+3*Lmax, :);
   end;
end;
filename= sprintf('rhannelAMBI%d-%d-%d', MCH, Mmax, Lmax)
save(filename, 'AMBI', 'M0', 'hnm', 'N');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm
alpha= 1/2;
load updown_h_rusk.mat;

MCH= 3;
M= 15;
L= 10;
filename= sprintf('channelAMBI%d-%d-%d', MCH, M, L)
load(filename);

A= zeros((MCH+1)*M*(2*L+1), (MCH+1)*M*(2*L+1));
for tau= 0:MCH
   for taup= 0:MCH
      [tau taup]
      for l= -L:L
         for lp= -L:L
            for m= 1:M
               for mp= 1:M
%                  [tau taup m mp l lp]
%                   [(m-1)*(2*L+1)*(MCH+1)+(L+l)*(MCH+1)+tau+1  (mp-1)*(2*L+1)*(MCH+1)+(L+lp)*(MCH+1)+taup+1]
                  A((m-1)*(2*L+1)*(MCH+1)+(L+l)*(MCH+1)+tau+1, (mp-1)*(2*L+1)*(MCH+1)+(L+lp)*(MCH+1)+taup+1)...
		      = AMBI{tau+1, taup+1}(3*L+1+l-lp, M+1+m-mp);
               end;
            end;
         end;
      end;
   end;
end;

a= zeros((MCH+1)*M*(2*L+1), (MCH+1));
for tau= 0:MCH
   for taup= 0:MCH
      for m= 1:M
         for l= -L:L
            a((m-1)*(2*L+1)*(MCH+1)+(L+l)*(MCH+1)+tau+1, taup+1)= AMBI{tau+1, taup+1}(3*L+1+l, M+1+m);
	 end;
      end;
   end;
end;

theta= -inv(A)*a;

CML= {};
for tau= 0:MCH
   for taup= 0:MCH
      [tau, taup]
      Aml= [[zeros(L, 1); 1; zeros(L, 1)] zeros(2*L+1, M)];
      for m= 1:M
         for l= -L:L
            Aml(L+1+l, m+1)= theta((m-1)*(2*L+1)*(MCH+1)+(L+l)*(MCH+1)+tau+1, taup+1);
            Aml= param_stabilize(Aml, N, .98, 1);
	 end;
      end;
      CML{tau+1, taup+1}= {Aml 1};
   end;
end;

filename= sprintf('channelPARA%d-%d-%d', MCH, M, L)
save(filename, 'CML', 'M0', 'hnm', 'N');