From ff6a4fa2f8c5bd8a87a517efb5bbaa8ace2d82aa Mon Sep 17 00:00:00 2001 From: drowe67 Date: Thu, 19 Apr 2018 23:50:28 +0000 Subject: [PATCH] updated fading sample generation git-svn-id: https://svn.code.sf.net/p/freetel/code@3509 01035d8c-6547-0410-b346-abe4f91aad63 --- codec2-dev/octave/tcohpsk.m | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/codec2-dev/octave/tcohpsk.m b/codec2-dev/octave/tcohpsk.m index 59ddae01..076a9589 100644 --- a/codec2-dev/octave/tcohpsk.m +++ b/codec2-dev/octave/tcohpsk.m @@ -64,8 +64,8 @@ randn('state',1); % select which test ---------------------------------------------------------- -test = 'compare to c'; -%test = 'awgn'; +%test = 'compare to c'; +test = 'awgn'; %test = 'fading'; % some parameters that can be over ridden, e.g. to disable parts of modem @@ -734,9 +734,11 @@ endfunction % function to write float fading samples for use by C programs -function write_noise_file(raw_file_name, Fs, len_samples) - [spread spread_2ms hf_gain] = init_hf_model(Fs, len_samples); - hf_gain +function write_noise_file(raw_file_name, Fs, dopplerSpreadHz, len_samples) + spread = doppler_spread(dopplerSpreadHz, Fs, len_samples); + spread_2ms = doppler_spread(dopplerSpreadHz, Fs, len_samples); + hf_gain = 1.0/sqrt(var(spread)+var(spread_2ms)); + % interleave real imag samples inter = zeros(1,len_samples*4); -- 2.25.1