add_gauss_noise.m

来自「一些关于空时编码的编码解码仿真程序」· M 代码 · 共 17 行

function [output] = add_gauss_noise(rec,noise)

[K,D,N,G,N_t,N_r,L,C]= header;
[s,g,v] = size(rec);
% nn = randn(s,g,v)*attn/2 + j*randn(s,g,v)*attn/2;



% rc = reshape(rec,1,N_r*N*K);
% rec_pwr = (rc*rc')/N_r/N/K;
% noise = rec_pwr/2*10^(-snr/10);
% 
n_p_db = 10*log10(noise);
ne = wgn(s,g*v, n_p_db,'complex');
 nn = reshape(ne,s,g,v);
% 
output = rec + nn;