Changeset 4009 in Sophya for trunk/AddOn

Timestamp:

Jun 23, 2011, 11:18:28 PM (14 years ago)

Author:

ansari

Message:

code combinaison I-Q (codage fonction CombineIQPairs) , Reza 23/06/2011

File:

• trunk/AddOn/TAcq/corrcrtadc.cc (modified) (1 diff)

trunk/AddOn/TAcq/corrcrtadc.cc

@@ -294 +294 @@
 }
 
+/* --Fonction-- */
 int CombineIQPairs(TMatrix< complex<r_4> >& cfour)
 {
-  // to be coded ...
+  // We assume the following channel numbering for the I-Q pairs
+  // Row pairs : (0,4) (1,5) (2,6) (3,7) (8,12) (9,13) (10,14) (11,15) ... are assumed to be 
+  // IQ pairs in the matrix cfour 
+  // We compute the I rows by the left side band, and the Q rows by the right side band 
+
+  // Left = a I + b Q
+  // Right = c I + d Q 
+  complex<r_4> za(1.,0.);
+  complex<r_4> zb(0.,1.);
+  complex<r_4> zc(1.,0.);
+  complex<r_4> zd(0.,-1.);
+
+  int nboards = cfour.NRows()/4;    
+  for(int jb=0; jb<nboards-1; jb++)  {
+    for(int c=0; c<4; c++) {
+      sa_size_t ii = jb*4+c;  // sa_size_t are long integers 
+      sa_size_t iq = ii+4;
+      //  Combine the two I,Q fourier coefficient vectors to form the left and right sid bands 
+      TVector< complex<r_4> > leftband  = za*cfour.Row(ii)+zb*cfour.Row(iq); 
+      TVector< complex<r_4> > rightband = zc*cfour.Row(ii)+zd*cfour.Row(iq); 
+      // replace the two I,Q vectors by left and right computed side bands 
+      cfour.Row(ii)=leftband;
+      cfour.Row(iq)=rightband;
+    }
+  }
   return 0;
 }