Commit eea02488 authored by Ling Sun's avatar Ling Sun
Browse files

fix para names

parent 14064319
...@@ -18,22 +18,22 @@ format long g ...@@ -18,22 +18,22 @@ format long g
startFreq = str2num(freq) + str2num(idx) startFreq = str2num(freq) + str2num(idx)
% Set directories and file name % Set directories and file name
dirName = '/home/ling.sun/O1/Viterbi/fstat/'; fstatdir = '/home/ling.sun/O1/Viterbi/fstat/';
outdir = '/home/ling.sun/O1/Viterbi/output/'; outdir = '/home/ling.sun/O1/Viterbi/output/';
outfile = strcat(outdir,'Viterbi-', num2str(startFreq),'.dat') outfile = strcat(outdir,'Viterbi-', num2str(startFreq),'.dat')
% Set parameters % Set parameters
a0median = 1.805 a0median = 1.805 % median asini being searched
P = 68023.70496 g = [-0.8:0.05:0.8]; % asini grid (+/-80% uncertainty)
steps = 13 P = 68023.70496 % orbit period
steps = 13 % total steps
% Read Fstats % Read Fstats
tic tic
[f,X] = readfstats(dirName, startFreq, steps); [f,X] = readfstats(fstatdir, startFreq, steps);
toc toc
fmean = f(1)+0.5; % Take the mean value of each 1-Hz sub-band fmean = f(1)+0.5; % Take the mean value of each 1-Hz sub-band
g = [-0.8:0.05:0.8]; % asini grid (+/-80% uncertainty)
a0 = (1+g)*a0median; % asini values being searched a0 = (1+g)*a0median; % asini values being searched
for m=1:length(g) for m=1:length(g)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment