Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- using Plots;
- using DSP;
- # read I/Q samples
- data = reinterpret(Complex{Float32}, read("/home/domagoj/gqrx_20190101_165640_97192000_1800000_fc.raw"));
- # create and 24th order Butterworth bandpass filter to sample array
- responsetype = Bandpass(50, 150000; fs=1800000);
- prototype = Butterworth(24);
- db = filt(digitalfilter(responsetype, prototype), data);
- # generate and plot 4096-FFT results
- #db = data;
- windowSize = 1024;
- n = floor(Int, size(db,1)/windowSize);
- for i in 1:n
- plot(db[(i-1)*windowSize+1:i*windowSize])
- ylims!(-1,1)
- xlims!(-1,1)
- gui();
- sleep(0.04)
- end
- # EOF
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement