From f2d9455f945808059e510b33c54b29f36ef257cc Mon Sep 17 00:00:00 2001
From: Florian Kaltenberger <florian.kaltenberger@eurecom.fr>
Date: Thu, 9 Jan 2014 11:09:57 +0000
Subject: [PATCH] adding a project directory for CORRIDOR

git-svn-id: http://svn.eurecom.fr/openair4G/trunk@4838 818b1a75-f10b-46b9-bf7c-635c3b92a50f
---
 targets/PROJECTS/CORRIDOR/emos_read.m | 43 +++++++++++++++++++++++++++
 1 file changed, 43 insertions(+)
 create mode 100644 targets/PROJECTS/CORRIDOR/emos_read.m

diff --git a/targets/PROJECTS/CORRIDOR/emos_read.m b/targets/PROJECTS/CORRIDOR/emos_read.m
new file mode 100644
index 0000000000..98b5eb00d9
--- /dev/null
+++ b/targets/PROJECTS/CORRIDOR/emos_read.m
@@ -0,0 +1,43 @@
+%% read data
+fid = fopen('E:\EMOS\corridor\eNB_data_20140108_185919.EMOS','r');
+[v,c]=fread(fid, 76800*4*100, 'int16',0,'ieee-le');
+v2 = double(v(1:2:end))+1j*double(v(2:2:end));
+fclose(fid);
+
+figure(1)
+plot(abs(fft(v2)))
+
+%% sync
+addpath('../../../openair1/PHY/LTE_REFSIG');
+primary_synch; %loads the primary sync signal
+
+pss_t = upsample(primary_synch0_time,4*4);
+[corr,lag] = xcorr(v2,pss_t);
+figure(2);
+plot(lag,abs(corr));
+
+%% frame start detection
+% TODO: do a proper peak detection. for now just take the max
+[m,i]=max(abs(corr));
+frame_start = lag(i);
+
+%% ofdm receiver
+load('../LOCALIZATION/ofdm_pilots_sync_2048.mat');
+ofdm_symbol_length = num_carriers + prefix_length;
+frame_length = ofdm_symbol_length*num_symbols_frame;
+useful_carriers = num_carriers-num_zeros-1;
+
+received = v2(frame_start:frame_start+frame_length);
+received_f = zeros(num_symbols_frame,useful_carriers);
+for i=0:num_symbols_frame-1;
+    ifblock=received(i*ofdm_symbol_length+[1:ofdm_symbol_length]);
+    ifblock(1:prefix_length)=[];
+    fblock=fft(ifblock);
+    received_f(i+1,:) = [fblock(end-useful_carriers/2+1:end); fblock(2:useful_carriers/2+1)];
+end
+
+%% channel estimation
+H=conj(f).*received_f;
+
+
+