s0 = I1z + I2z
-----p/2(I1x + I2x)----->
s1 = -I1y - I2y
-----W1t1I1z
+ W2t1I2z----->
s2 = -I1y cos(W1t1)
- I2y cos(W2t1)
+ I1x sin(W1t1)
+ I2x sin(W2t1)
-----pJ12t1(2I1zI2z)----->
s3 = [-I1y cos(W1t1)
- I2y cos(W2t1)
+ I1x sin(W1t1)
+ I2x sin(W2t1)]
cos(pJ12t1)
+ [I1xI2z cos(W1t1)
+ I1zI2x cos(W2t1)
+ I1yI2z sin(W1t1)
+ I1zI2y sin(W2t1)]
sin(pJ12t1)
-----p/2(I1x + I2x)----->
s4 = [-I1z cos(W1t1)
- I2z cos(W2t1)
+ I1x sin(W1t1)
+ I2x sin(W2t1)]
cos(pJ12t1)
+ [-I1xI2y cos(W1t1)
- I1yI2x cos(W2t1)
- I1zI2y sin(W1t1)
- I1yI2z sin(W2t1)]
sin(pJ12t1)
-----during acquisition----->
s5 = [-I1z cos(W1t1)
- I2z cos(W2t1)]
cos(pJ12t1)
+ [I1x sin(W1t1)
cos(W1t2) + I2x
sin(W2t1)
cos(W2t2)
+ I1y sin(W1t1)
sin(W1t2) + I2y
sin(W2t1)
sin(W2t2)] cos(pJ12t1)
cos(pJ12t2)
+ [I1yI2z sin(W1t1)
cos(W1t2) + I1zI2y
sin(W2t1)
cos(W2t2)
- I1xI2z sin(W1t1)
sin(W1t2) - I1zI2x
sin(W2t1)
sin(W2t2)] cos(pJ12t1)
sin(pJ12t2)
+ {expansion of [-I1xI2y cos(W1t1)
- I1yI2x cos(W2t1)]
sin(pJ12t1) leads
to 32 unobservable terms}
+ [-I1zI2y sin(W1t1)
cos(W2t2) - I1yI2z
sin(W2t1)
cos(W1t2)
+I1zI2x sin(W1t1)
sin(W2t2) + I1xI2z
sin(W2t1)
sin(W1t2)] sin(pJ12t1)
cos(pJ12t2)
+ [I2x sin(W1t1)
cos(W2t2) + I1x
sin(W2t1)
cos(W1t2)
+ I2y sin(W1t1)
sin(W2t2) + I2y
sin(W2t1)
sin(W1t2)] sin(pJ12t1)
sin(pJ12t2)