Science/psdlag-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/psdlag-agn.git synced 2024-11-24 01:45:05 +00:00
Code Issues Packages Projects Releases Wiki Activity

new frequency bin paradigm, should be more robust# Please enter the commit message for your changes. Lines starting

Browse Source
This commit is contained in:
caes 2017-01-30 03:09:59 -05:00
parent 11216e4401
commit 9fe84e07b3
1 changed files with 3 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
scripts/analyze_lightcurve.py
View File

@ -18,7 +18,7 @@ except getopt.GetoptError:
# Time resolution determined from inspection and testing. This script # Time resolution determined from inspection and testing. This script
# does not expect evenly spaced data in time. # does not expect evenly spaced data in time.
dt = 0.01 dt = 0.1
#### Get the psd for the first light curve #### #### Get the psd for the first light curve ####
@ -54,12 +54,9 @@ fqL = np.array([0.0049999999, 0.018619375, 0.044733049, 0.069336227, 0.10747115,
# fqL = np.concatenate(([0.5/seg_length], fqL)) # fqL = np.concatenate(([0.5/seg_length], fqL))
# #
f1 = 1/270. f1 = 1/175.
f2 = 0.5/dt f2 = 0.5/dt
fqL = np.hstack((np.array(0.5*f1),np.logspace(np.log10(0.9*f1),np.log10(0.3*f2),9),np.array(2*f2)))
fqL = np.hstack((np.array(0.5*f1),np.logspace(np.log10(0.9*f1),np.log10(.5*f2),7),np.array(2*f2)))
nfq = len(fqL) - 1 nfq = len(fqL) - 1
fqd = 10**(np.log10( (fqL[:-1]*fqL[1:]) )/2.) fqd = 10**(np.log10( (fqL[:-1]*fqL[1:]) )/2.)