estimateChannelError = [];
    estimateCutChannelError = [];

    for j = 1 : 3
        [basebandChannel, channel_response] = baseband_channel(transmitter_output_data, channel);
        basebandChannel = awgn(basebandChannel, noise(i), 'measured');
        pilot_ach = get_pilot(basebandChannel(1:68160), 2);

        for p = 1 : 10
            estimateChannel_fdomain = pilot_ach(:, p) ./ pilot(:, p);
            estimateChannel_tdomain = ifft(estimateChannel_fdomain, N);
            estimateChannel_tdomain = [estimateChannel_tdomain(1 : CP); zeros([(N - CP) 1])];
            estimateCutChannel_tdomain = fft(estimateChannel_tdomain, N);
            channel_response_fft = fft(channel_response, 128).';
            estimateChannel_fdomain_error_part(p) = sum( abs(channel_response_fft - estimateChannel_fdomain).^2);
            estimateCutChannel_tdomain_error_part(p) = sum( abs(channel_response_fft - estimateCutChannel_tdomain).^2);
        end
        estimateChannelError = [estimateChannelError estimateChannel_fdomain_error_part];
        estimateCutChannelError = [estimateCutChannelError estimateCutChannel_tdomain_error_part];
    end

    estimateChannel_fdomain_error(i) = mean(estimateChannelError);
    estimateCutChannel_tdomain_error(i) = mean(estimateCutChannelError);