clear all;clc;
EbNo_Vec = [0:5:30];
BER_Vec=[];
for n = 1:length(EbNo_Vec);
    snr = EbNo_Vec(n)+10*log10(2);
    t=0:4e-5:1e2;
    
    sigma=1/(10^(snr/10));
    sim('lhw_VBLAST_MMSE');
    BER_Vec(n,:) = ErrorVec;
    format short e
    EbNo_Vec(n)
    ErrorVec
end;
semilogy (EbNo_Vec,BER_Vec(:,1), 'DisplayName', 'BER_Vec'); figure(gcf)
axis([0 30 1e-5 1]);
grid on