From 46405a45a5c708170e5ddb4bb139a66755ec8b85 Mon Sep 17 00:00:00 2001 From: Jens Ahrensfeld Date: Mon, 10 Aug 2020 17:49:00 +0000 Subject: [PATCH] - improved git-svn-id: http://moon:8086/svn/matlab/trunk@132 801c6759-fa7c-4059-a304-17956f83a07c --- mc_pd_eval3.m | 17 +++++++++-------- 1 file changed, 9 insertions(+), 8 deletions(-) diff --git a/mc_pd_eval3.m b/mc_pd_eval3.m index ff609b1..e243138 100644 --- a/mc_pd_eval3.m +++ b/mc_pd_eval3.m @@ -46,25 +46,26 @@ N = Nc*Nspc*slave_oversample; assert (N == length(x)); # Create slave -pc = 0.5 + perr; +phase_lo = 0.0 + perr; dfreq = freq/slave_oversample; -dpc = 0; +omega = 0; sample_count = slave_oversample; for n=1:N, - perr = phase_det(pc, x(n)); + phase = mod(phase_lo - 0.5, 1); + perr = phase_det(phase, x(n)); sample_count = sample_count - 1; if sample_count == 0, sample_count = slave_oversample; else perr = 0; end - _lo(n) = mod(pc - 0.5, 1); - dpc = dfreq - (klag*lag + klead*perr); - pc = mod(pc + dpc, 1); + _lo(n) = phase_lo; + omega = dfreq - (klag*lag + klead*perr); + phase_lo = mod(phase_lo + omega, 1); lag = lag + perr; _perr(n) = perr; - _dpc(n) = dpc; + _omega(n) = omega; end function perr = phase_det(lo, x) @@ -76,7 +77,7 @@ subplot(3, 1, 1) plot(1:N, _lo, 1:N, x); legend('Lo', 'x'); grid subplot(3, 1, 2) -plot(1:N, _dpc); legend('f'); grid +plot(1:N, _omega); legend('f'); grid subplot(3, 1, 3) plot(1:N, _perr); legend('p'); grid