channel_est.m

OFDM系统仿真

function [ich10,qch10]=channel_est(ich9,qch9,pilot_symbol)
ich10 = zeros(size(ich9));
qch10 = zeros(size(qch9));
for m=1:4
    pilot_symbol_re=ich9(:,9*m-8);
    pilot_symbol_im=qch9(:,9*m-8);
    rx_pilot_symbol=pilot_symbol_re+i*pilot_symbol_im;
    rx_pilot_symbol=reshape(rx_pilot_symbol,1,768);
    h=rx_pilot_symbol./pilot_symbol;
    h=reshape(h,768,1);
%     coeff = [1/8, 1/4, 1/2, 1, 1/2, 1/4, 1/8];
    
    f = [0, 0.25, 0.5, 1];
    amp = [1 ,1 ,0 ,0];
    coeff= remez(16, f, amp);
    h1 = conv(h, coeff);
    h2 = h1((length(h1)-768)/2:(length(h1)-768)/2+767);
    for n=1:8
        ich10(:,9*m+n-8)= real((ich9(:,9*m+n-8)+j*qch9(:,9*m+n-8))./h2);
        qch10(:,9*m+n-8)= imag((ich9(:,9*m+n-8)+j*qch9(:,9*m+n-8))./h2);
    end
end