Changeset 9452 for ossiedev/branches/hvolos/components/MIMODec/MIMODec.cpp

Timestamp:

07/16/09 22:53:04 (4 years ago)

Author:

hvolos

Message:

Added alamouti STBC

Files:

• ossiedev/branches/hvolos/components/MIMODec/MIMODec.cpp (modified) (5 diffs)

ossiedev/branches/hvolos/components/MIMODec/MIMODec.cpp

@@ -10 +10 @@
 #include "MIMODec.h"
 
+using namespace arma;
+
 MIMODec_i::MIMODec_i(const char *uuid, omni_condition *condition) : 
     Resource_impl(uuid), component_running(condition) 
 {
-    dataOut_0 = new standardInterfaces_i::complexShort_u("Symbols_OUT");
-    dataIn_0 = new standardInterfaces_i::complexShort_p("MUX_MIMO_Symbols_IN");
+    dataOut = new standardInterfaces_i::complexShort_u("Symbols_OUT");
+    dataIn = new standardInterfaces_i::complexShort_p("MUX_MIMO_Symbols_IN");
 
     //Create the thread for the writer's processing function 
@@ -26 +28 @@
 MIMODec_i::~MIMODec_i(void)
 {   
-    delete dataOut_0;
-    delete dataIn_0;
+    delete dataOut;
+    delete dataIn;
 }
 
@@ -43 +45 @@
     CORBA::Object_var p;
 
-    p = dataOut_0->getPort(portName);
+    p = dataOut->getPort(portName);
 
     if (!CORBA::is_nil(p))
         return p._retn();
 
-    p = dataIn_0->getPort(portName);
+    p = dataIn->getPort(portName);
 
     if (!CORBA::is_nil(p))
@@ -99 +101 @@
             CORBA::Short simple_temp;
             props[i].value >>= simple_temp;
-            simple_0_value = simple_temp;
+            Nt = simple_temp;
         }
 
         if (strcmp(props[i].id, "DCE:74a38a7c-3b3e-11de-b7e0-0000e80014f9") == 0)
         {
+            const char * prop_str;
+            props[i].value >>= prop_str;
+            MIMOmethod=prop_str;
+        }
+        
+        if (strcmp(props[i].id, "DCE:a1026f5c-3b3e-11de-b7e0-0000e80014f9") == 0)
+        {
+            CORBA::Short simple_temp;
+            props[i].value >>= simple_temp;
+            Mr = simple_temp;
+        }
+        
+        if (strcmp(props[i].id, "DCE:17cabf40-3b3f-11de-b7e0-0000e80014f9") == 0)
+        {
+            const char * prop_str;
+            props[i].value >>= prop_str;
+            Modmethod=prop_str;
+        }
+      }
 }
 
@@ -110 +131 @@
     DEBUG(3, MIMODec, "ProcessData() invoked")
 
-    PortTypes::ShortSequence I_out_0, Q_out_0;
+    PortTypes::ShortSequence I_out, Q_out;
 
 
-    PortTypes::ShortSequence *I_in_0(NULL), *Q_in_0(NULL);
-    CORBA::UShort I_in_0_length, Q_in_0_length;
+    PortTypes::ShortSequence *I_in(NULL), *Q_in(NULL);
+    CORBA::UShort Data_length;
 
     while(1)
     {
-        dataIn_0->getData(I_in_0, Q_in_0);
+       short  int NH=4;
+        dataIn->getData(I_in, Q_in);
 
-        I_in_0_length = I_in_0->length();
-        Q_in_0_length = Q_in_0->length();
+        Data_length = I_in->length();
 
-        I_out_0.length(); //must define length of output
-        Q_out_0.length(); //must define length of output
+
+        I_out.length((Data_length-NH)/2); //must define length of output
+        Q_out.length((Data_length-NH)/2); //must define length of output
 
         /*insert code here to do work*/
+        
+       
+        cx_mat H(2,2);
+        for (int i=0;i<2;i++){
+                for (int j=0;j<2;j++){
+                H(i,j)=cx_double((*I_in)[i*2+j],(*Q_in)[i*2+j])/10000.0;
+                //std::cout<<"Iin:"<<(*I_in)[i*2+j];
+                }
+        }
+        //std::cout<<std::endl;
+        
+        /*Alamouti
+        for i=1:2:length(RX_symbols1),
+          r1=RX_symbols1(i);
+          r2=RX_symbols1(i+1);
+          RX_symbols(i)=r1*conj(H(1,1))+conj(r2)*H(1,2);
+          RX_symbols(i+1)=r1*conj(H(1,2))-conj(r2)*H(1,1);
+        end;
+        */
+
+        cx_double r1;
+        cx_double r2;
+        cx_mat rx_symbols(((Data_length-NH)/2),1);
+        short int j=0;
+        //std::cout<<"Decode"<<std::endl;
+        std::cout<<"Hdec="<<H;
+        for (short int i=NH; i<Data_length;i+=4){
+                r1=cx_double((*I_in)[i],(*Q_in)[i]);
+                r2=cx_double((*I_in)[i+2],(*Q_in)[i+2]);
+               // std::cout<<((Data_length-NH)/2)<<" i:"<<i<<" j"<<j<<std::endl;
+                rx_symbols(j,0)=r1*conj(H(0,0))+conj(r2)*H(0,1);
+                rx_symbols(j+1,0)=r1*conj(H(0,1))-conj(r2)*H(0,0);                                            
+                j+=2;
+        } 
+        
+        for (short int i=0; i<((Data_length-NH)/2);i++){
+        I_out[i]=(CORBA::Short)real(rx_symbols(i,0));
+        Q_out[i]=(CORBA::Short)imag(rx_symbols(i,0));
+        std::cout<<"Symbol"<<rx_symbols(i,0)<<std::endl;
+        std::cout<<"I:"<<I_out[i]<<"Q:"<<Q_out[i]<<std::endl;
+        }
 
 
-
-
-
-
-        dataIn_0->bufferEmptied();
-        dataOut_0->pushPacket(I_out_0, Q_out_0);
+        dataIn->bufferEmptied();
+        dataOut->pushPacket(I_out, Q_out);
     }
 }