Mercurial > hg > fxanalyse

comparison FXAnalyse.c @ 117:96ab53ab9496

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix computaion of expected frep change
author Daniele Nicolodi <daniele.nicolodi@obspm.fr>
date Fri, 05 Jul 2013 16:23:43 +0200
parents 4f56dbdad85d
children 43b35f4aae78
comparison
equal deleted inserted replaced
116:4f56dbdad85d 117:96ab53ab9496
1157 f_rep_minus = f_rep_minus / n_3; 1157 f_rep_minus = f_rep_minus / n_3;
1158 f_beat_Sr_minus = f_beat_Sr_plus / n_3; 1158 f_beat_Sr_minus = f_beat_Sr_plus / n_3;
1159 1159
1160 // check delta frep 1160 // check delta frep
1161 double delta_f_rep = f_rep_plus - f_rep_minus; 1161 double delta_f_rep = f_rep_plus - f_rep_minus;
1162 double expected = N1 / Ndiv * 2.0 * DeltakHz_3 * 1000.0; 1162 double expected = Ndiv * 2.0 * DeltakHz_3 * 1000.0 / N1;
1163 logmsg("delta frep=%g Hz expected=%g Hz", delta_f_rep, expected); 1163 logmsg("delta frep=%g Hz expected=%g Hz", delta_f_rep, expected);
1164 1164
1165 // compute N3 1165 // compute N3
1166 double delta_f_beat_Sr = f_beat_Sr_plus - f_beat_Sr_minus + 2.0 * N3/N1 * Ndiv * DeltakHz_3 * 1000; 1166 double delta_f_beat_Sr = f_beat_Sr_plus - f_beat_Sr_minus + 2.0 * N3/N1 * Ndiv * DeltakHz_3 * 1000;
1167 logmsg("delta fbeat=%g", delta_f_beat_Sr); 1167 logmsg("delta fbeat=%g", delta_f_beat_Sr);