aboutsummaryrefslogtreecommitdiff
path: root/octave/doppler_spread_ut.m
diff options
context:
space:
mode:
Diffstat (limited to 'octave/doppler_spread_ut.m')
-rw-r--r--octave/doppler_spread_ut.m51
1 files changed, 0 insertions, 51 deletions
diff --git a/octave/doppler_spread_ut.m b/octave/doppler_spread_ut.m
deleted file mode 100644
index f7d96b0..0000000
--- a/octave/doppler_spread_ut.m
+++ /dev/null
@@ -1,51 +0,0 @@
-% doppler_spread_ut.m
-% David Rowe Jan 2016
-%
-% Unit test script for doppler_spread
-
-f = 1;
-Fs = 8000;
-N = Fs*10;
-
-[spread states] = doppler_spread(f, Fs, N);
-
-% use spreading samples to modulate 1000Hz sine wave
-% You can listen to this with: sine1k_1Hz.raw
-
-% $ play -t raw -r 8000 -s -2
-s = cos(2*pi*(1:N)*1000/Fs);
-s = s .* spread;
-s = real(s)*5000;
-fs = fopen("sine1k_1Hz.raw","wb"); fwrite(fs,s,"short"); fclose(fs);
-
-% Some plots
-
-x = states.x; y = states.y; b = states.b;
-
-H = freqz(b,1,x);
-
-figure(1)
-clf
-subplot(211)
-plot(x,y,';target;')
-title('Gaussian Filter Freq Resp Lin');
-legend('boxoff');
-subplot(212)
-plot(x,20*log10(y),';target;')
-hold on;
-plot(x,20*log10(y),'g+;actual;')
-hold off;
-axis([0 f*10/2 -60 0])
-title('Gaussian Filter Freq Resp dB');
-xlabel('Freq (Hz)');
-legend('boxoff');
-
-figure(2);
-subplot(211)
-plot(abs(spread))
-title('Spreading Function Magnitude');
-subplot(212)
-plot(s)
-title('1000Hz Sine Wave');
-xlabel('Time (samples)')
-