I am trying to re-analyze the data that Koji took last night.
I think that my script is just pulling out the I and Q data for each port, and each degree of freedom, calculating the magnitude from sqrt( I**2 + Q**2 ) and the phase from atan2( I / Q ). No calibration.
If I print out the results, I get:
Sensing Matrix, units = cts/ct, phase in degrees
MICH Mag MICH Phase PRCL Mag PRCL Phase
AS55_I 1.627E-02 62.063 4.189E-03 68.344
AS55_Q 2.073E-02 -105.353 1.983E-02 66.361
REFL11_I 8.165E+02 -112.624 2.441E+00 77.911
REFL11_Q 2.712E+02 -112.650 7.065E-01 -127.093
REFL33_I 8.028E+00 -112.154 6.282E-02 70.990
REFL33_Q 5.490E-02 -165.912 9.908E-03 61.269
REFL55_I 8.347E+00 -112.085 2.146E-02 78.928
REFL55_Q 3.003E-01 -151.652 7.924E-02 87.153
If, however, I take the raw values that are stored in the data file, for one row (say, REFL33_Q) and calculate by hand (same formulas), I get different results:
MICH Mag MICH Phase PRCL Mag PRCL Phase
REFL33_Q 9.9E-03 28.89 5.46E-02 -103.8
Contrast that with Koji's uncalibrated transfer function result from elog 8611:
MICH Mag MICH Phase PRCL Mag PRCL Phase
REFL33Q 1.8665e-5 71.1204 1.6310e-4 -141.73
I am currently confused, and need to re-look at my script, as well as make sure I am actually measuring the things I think I am.
EDIT: This has been fixed, in that my 2 calculations agree with one another. I have crossed out the incorrect numbers, and put correct numbers below. I still don't agree with Koji, but at least I agree with myself.
The phase issue: I needed to calculate the phase with "ATAN2(I,Q)", which I did when I calculated by hand, but the script had "atan2(Q,I)". This has been fixed.
The magnitude issue: They match, but my "pretty print" script labels MICH as PRCL, and vice versa. Doh.
Corrected values:
Sensing Matrix, units = cts/ct, phase in degrees
PRCL Mag PRCL Phase MICH Mag MICH Phase
AS55_I 1.627E-02 27.937 4.189E-03 21.656
AS55_Q 2.073E-02 -164.647 1.983E-02 23.639
REFL11_I 8.165E+02 -157.376 2.441E+00 12.089
REFL11_Q 2.712E+02 -157.350 7.065E-01 -142.907
REFL33_I 8.028E+00 -157.846 6.282E-02 19.010
REFL33_Q 5.490E-02 -104.088 9.908E-03 28.731
REFL55_I 8.347E+00 -157.915 2.146E-02 11.072
REFL55_Q 3.003E-01 -118.348 7.924E-02 2.847
|