equalizer.m

来自「ofdm_64QAM_simulation,只是进行的初步的实验,大家仅仅供参考」· M 代码 · 共 26 行

% Equalizer Part
function equal(j*SubCarrier+1:(j+1)*SubCarrier)=kh_k
mapping=[ 7*d + i*7*d;
   -7*d + i*7*d;
   7*d - i*7*d;
   -7*d - i*7*d];
for sub=1:SubCarrier,
   if j==0,
      estimated_c(sub)=1;                                     % initial estimate of ICI
   end
   %k=sub;
   y_k=FFTData(j*SubCarrier+sub);
   z_k(sub)=estimated_c(sub)*y_k;                      % equalizer out put
   for b=1:4,
      metrics(sub,b)=(real(z_k(sub))-real(mapping(b,1)) ).^2 + (imag(z_k(sub))-imag(mapping(b,1))).^2;
   end;
   % khryu=[];
   [min_metrics decs]=min(metrics(sub,:));
   khryu=mapping(decs);
   e_k(sub)=khryu-z_k(sub);
   estimated_c(sub)=estimated_c(sub)+delta*e_k(sub)*conj(y_k);   % coefficients update
   mse(sub)=e_k(sub)^2;
end
kh_k=z_k.';
lms=[lms mse(dis)];