EoRImaging / eppsilon

eppsilon - error propagated power spectrum with interleaved observed noise
BSD 2-Clause "Simplified" License
5 stars 4 forks source link

dft with true z #105

Closed nicholebarry closed 4 years ago

nicholebarry commented 5 years ago

I think maybe the true dft isn't quite right (but I could totally be wrong).

Right now it's z_exp_ztrue = exp(1.*complex(0,1)*matrix_multiply(comov_dist_los-min(comov_dist_los), kz_mpc_orig_trim, /btranspose)) I'm suggesting z_exp_ztrue = exp(-1.*complex(0,1)*matrix_multiply(reverse(comov_dist_los-min(comov_dist_los)), kz_mpc_orig_trim, /btranspose))

They don't give the same transform, I checked.

nicholebarry commented 5 years ago

Here is the real part of the FT matrix for a regularized z, the original true z, and the proposed true z. Notice how the regularized z is more similar to the proposed true z (look for a left-right "shift"). The images are similar for the imaginary part. Screen Shot 2019-06-28 at 3 52 27 pm

Here is the PS difference between reg z and the original true z, for reference. reg_minus_true

Here is the PS difference between reg z and the proposed true z. Notice how the red and blue flipped in the window? That seems like a logical result to flipping the dft direction. However, the majority of the difference seems to mimic the one above, indicating that the overall de-regularization of the dft obscures this direction effect. reg_minus_prop

Here is the PS difference between the original true z and the proposed true z. MWA_limit2019_fullimg_ch9-126_averemove_swbh_dencorr__btl_noalltv_noocc4_minus_btl_noalltv_noocc4_2dkpower

The difference PS didn't really help me. But the matrix comparison (first image) gives me confidence that this is right.

bhazelton commented 5 years ago

Very interesting. I need to think about this a little more.

bhazelton commented 4 years ago

closed by #107