home/comp_snr_ts.md

@@ -62,5 +62,19 @@ plt.close()
 - PSD of the whitened time series in pycbc
 As expected, PSD ~ 1
 
+```
+
+d_pycbc = d_gwpy.to_pycbc()
+psd = d_pycbc.psd(1)
+
+w = d_pycbc.whiten(1, 1)
+psd = w.psd(1)
+
+plt.loglog(psd.sample_frequencies, psd.data, label='psd of the whitened timeseries')
+plt.xlim(xmin=10)
+plt.legend()
+plt.savefig('test.png')
+plt.close()
+```
 ![Screen_Shot_2021-03-11_at_4.25.13_PM](uploads/f238606e0dd7f52ff3bdccf3ed8dfe67/Screen_Shot_2021-03-11_at_4.25.13_PM.png)