noise_estimation.m

mimo2x2天线选择系统的全系统matlab程序,先前的是dsp程序.

function sigma =  noise_estimation (  rA_data_I, rA_data_Q, rB_data_I, rB_data_Q,...
                                      tr1_I, tr1_Q , tr2_I, tr2_Q, H_hat, L );
                                  
                                  
                                    
Y = [ rA_data_I + j*rA_data_Q; rB_data_I + j*rB_data_Q ];

G_ZF = inv(H_hat'*H_hat)*H_hat';

Y = G_ZF*Y;

S = [ tr1_I+ j* tr1_Q ; tr2_I+j*tr2_Q ];

for k=1:2
Si(k,:) = interp(S(k,:),L);
end

% figure
% plot(real(Y(1,:)))
% hold on;
% plot(real(Si(1,:)),'g');

sigma = mean(std((Y-Si)'));