channel.m

这是一个vblast程序

function [rec_sig,channel_matrix,noise_var]=channel(tran_sig,Eb_NO_dB,R_Antenna,mod_num,code_rate)
Eb_NO=10^(Eb_NO_dB/10);  
[T_Antenna,sig_length]=size(tran_sig);
noise_var = 1/(log2(mod_num)*Eb_NO*code_rate);       
rec_sig=zeros(R_Antenna,sig_length);
channel_matrix=(randn(R_Antenna,T_Antenna,sig_length)+j*randn(R_Antenna,T_Antenna,sig_length))/sqrt(2);
for i1=1:sig_length
    rec_sig(:,i1)=channel_matrix(:,:,i1)*tran_sig(:,i1);
end
noise=sqrt(1/2*noise_var)*(randn(R_Antenna,sig_length)+j*randn(R_Antenna,sig_length));
rec_sig=rec_sig+noise;