From: drowe67 <drowe67@01035d8c-6547-0410-b346-abe4f91aad63>
Date: Thu, 22 Dec 2016 06:27:29 +0000 (+0000)
Subject: used parabolic interpolator for phase synth
X-Git-Url: http://git.whiteaudio.com/gitweb/?a=commitdiff_plain;h=a8653d7ce4441a213f9410745c0ffcbf0307f5f6;p=freetel-svn-tracking.git

used parabolic interpolator for phase synth

git-svn-id: https://svn.code.sf.net/p/freetel/code@2944 01035d8c-6547-0410-b346-abe4f91aad63
---

diff --git a/codec2-dev/octave/newamp.m b/codec2-dev/octave/newamp.m
index 3609b9eb..3a151203 100644
--- a/codec2-dev/octave/newamp.m
+++ b/codec2-dev/octave/newamp.m
@@ -128,8 +128,13 @@ function [phase Gdbfk s Aw] = determine_phase(model, f, ak)
   L       = min([model(f,2) max_amp-1]);
   Wo      = model(f,1);
 
-  mask_sample_freqs_kHz = (Fs/1000)*[0:Nfft/2]/Nfft;           % fft frequency grid (nonneg freqs)
-  Gdbfk = resample_mask(model, f, mask_sample_freqs_kHz);
+  sample_freqs_kHz = (Fs/1000)*[0:Nfft/2]/Nfft;           % fft frequency grid (nonneg freqs)
+  Am = model(f,3:(L+2));
+  AmdB = 20*log10(Am);
+  rate_L_sample_freqs_kHz = (1:L)*Wo*4/pi;
+  Gdbfk = interp_para(rate_L_sample_freqs_kHz, AmdB, sample_freqs_kHz);    
+
+  % Gdbfk = resample_mask(model, f, mask_sample_freqs_kHz);
 
   % optional input of aks for testing