blind_equalization.m

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

function [avg_ber,snr_db]=blindequalization(num_of_frames,cnst,SNR_POINTS,tap_length)
    frame_length=cnst*tap_length;
    seed=15;
    snr_db=zeros(1,SNR_POINTS);
    avg_ber=zeros(1,SNR_POINTS);
    for iframe=1:num_of_frames
        iframe
        b=bit_gen(frame_length,seed);
        s=BPSK(b);
    
       for isnr=1:SNR_POINTS
           snr_db(isnr)=(isnr-1)*3+10;
           channel_output=channel(s,snr_db(isnr));
           equalized_output=equalizor(channel_output,tap_length);
           b_hat=BPSK_DE(equalized_output);
           ber = ber_compute(b, b_hat, frame_length);
           avg_ber(isnr) = ber/iframe + avg_ber(isnr)*(iframe-1)/iframe;
       end
    end