equalizor.asv

来自「THIS IS A BLIND QUALIZATION ALGORITHM. 」· ASV 代码 · 共 26 行

function [equalized_output]=equalizor(channel_output,tap_length)
    a=0.0001;     %a is the step size
    vec_length=length(channel_output);
    C_cfts=zeros(tap_length,1);
    indx1=-(tap_length-1)/2:1:(tap_length-1)/2;
    
%     for ii=1:vec_length
%         
%         %indx2=indx1+(ii-1)*tap_length;
%         if(ii<=(tap_length-1)/2)
%             
%           equalized_output(ii)=channel_output(ii);
%         elseif(ii<=vec_length-(tap_length-1)/2)
%           indx2=indx1+ii;
%           equalized_output(ii)=sum(C_cfts.*(channel_output(indx2)'));
%           C_cfts=C_cfts-a*(equalized_output(ii)-sign(equalized_output(ii)))*conj(channel_output(indx2)');
%           %C_cfts is a column vector;
%         else
%             equalized_output(ii)=channel_output(ii);
%         end
%     end
    for ii=(tap_length-1)/2+1:vec_length-(tap_length-1)/2
        indx2=indx1+ii;
          equalized_output(ii)=sum(C_cfts.*(channel_output(indx2)'));
          C_cfts=C_cfts-a*(equalized_output(ii)-sign(equalized_output(ii)))*conj(channel_output(indx2)');