get_noise.m

来自「一个UWB仿真程序包」· M 代码 · 共 34 行

function noise = get_noise(len,noise_var,samples_per_chip,TC,FS_CONT,RX_BW)


n_run = 5;

ALPHA_N = ceil(len/2/samples_per_chip/n_run); %make it integer multiple
T     = TC;
Nt    = T*FS_CONT;
Scale = (2*RX_BW)/FS_CONT;

%try to get close to a power of 2
N = ALPHA_N * Nt;
target_N = 2^ceil(log2(N));
ALPHA_N = ceil(target_N/Nt);

noise = zeros(1,ALPHA_N*2*samples_per_chip*n_run);

for i=1:n_run
    %noise = sqrt(noise_var)*(randn(1,length(rx_uoi)));
    [noise_r noise_i] = ...
        randn_bl(ALPHA_N,T,Nt,RX_BW,noise_var*Scale);
    
    start = (i-1)*2*ALPHA_N*samples_per_chip+1;
    noise(start:start+ALPHA_N*samples_per_chip-1)=noise_r;
    noise(start+ALPHA_N*samples_per_chip:start+2*ALPHA_N*samples_per_chip-1)=noise_i;
    
    %noise = [noise noise_r noise_i];
    clear noise_r;
    clear noise_i;
end

noise = noise(1:len);

return;