function [anoise] = VelocityNoise(N1, delta)

% Spectrum not available online. Assume similar spectrum as that of Acceleration
% noise..but assumed to be orders of magnitude larger than accel. noise

anoise = randn(N1, 1); 
anoise = fft (anoise);

df = (1/(N1*delta)); 

for k1 = 2:(N1/2) + 1
    fk = (k1-1)*df; 
    
    tf = (10^(-10)) * sqrt(2.5/6)/(fk);
    
    anoise(k1) = tf * anoise(k1);
end

anoise(1) = anoise(1)*(10^(-19)); 

for k1 = (N1/2)+2: N1
    anoise(k1) = conj(anoise( N1 - k1 +2)); 
end

%imaginary part in the order of 10^-34, 
%real part 10^-19. imaginary part : computational artifact
anoise = real(ifft(anoise));