git-svn-id: http://svn.eurecom.fr/openair4G/trunk@5204 818b1a75-f10b-46b9-bf7c-635c3b92a50f

targets/PROJECTS/CORRIDOR/emos_read_ca.m

 close all
 clear all
-nant1 = 1;
+
+enable_plots=2; %enables figures
+
+%% preload and init data
+addpath('../../../openair1/PHY/LTE_REFSIG');
+primary_synch; %loads the primary sync signal
+pss1_t = upsample(primary_synch0_time,4*4);
+pss2_t = upsample(primary_synch0_time,4*2);
+
+%load('E:\EMOS\corridor\ofdm_pilots_sync_2048_v7.mat');
+load('ofdm_pilots_sync_30MHz.mat');
+
+filename = 'E:\EMOS\corridor\eNB_data_20140319_133327.EMOS';
+
+nb_rb1 = 100; %this can be 25, 50, or 100
+num_carriers1 = 2048/100*nb_rb1;
+num_zeros1 = num_carriers1-(12*nb_rb1+1);
+prefix_length1 = num_carriers1/4; %this is extended CP
+ofdm_symbol_length1 = num_carriers1 + prefix_length1;
+frame_length1 = ofdm_symbol_length1*num_symbols_frame;
+useful_carriers1 = num_carriers1-num_zeros1-1;
+
+nb_rb2 = 50; %this can be 25, 50, or 100
+num_carriers2 = 2048/100*nb_rb2;
+num_zeros2 = num_carriers2-(12*nb_rb2+1);
+prefix_length2 = num_carriers2/4; %this is extended CP
+ofdm_symbol_length2 = num_carriers2 + prefix_length2;
+frame_length2 = ofdm_symbol_length2*num_symbols_frame;
+useful_carriers2 = num_carriers2-num_zeros2-1;
+
+nant1 = 2;
 nant2 = 2;
 samples_slot1 = 7680*2;
 samples_slot2 = 7680;
 samples_slot_agg = nant1*samples_slot1 + nant2*samples_slot2;
-nframes = 10;
+nframes = 50;
 slots_per_frame = 20;
 
-%% read data
-fid = fopen('E:\EMOS\corridor\eNB_data_20140311_173205.EMOS','r');
+d = dir(filename);
+nblocks = floor(d.bytes/(samples_slot_agg*slots_per_frame*nframes*4));
+PDP1_total = zeros(nblocks*nframes,useful_carriers1);
+PDP2_total = zeros(nblocks*nframes,useful_carriers2);
+
+%% main loop
+fid = fopen(filename,'r');
+block = 1;
+while ~feof(fid)
+    fprintf(1,'Processing block %d of %d',block,nblocks);
+    
+    %%
     [v,c]=fread(fid, 2*samples_slot_agg*slots_per_frame*nframes, 'int16',0,'ieee-le');
-fclose(fid);
+    if (c==0)
+        break
+    end
     v0 = double(v(1:2:end,:))+1j*double(v(2:2:end,:));
 
-%%
     v1 = zeros(samples_slot1*slots_per_frame*nframes,nant1);
-v2 = zeros(samples_slot1*slots_per_frame*nframes,nant2);
+    v2 = zeros(samples_slot2*slots_per_frame*nframes,nant2);
     for slot=1:slots_per_frame*nframes
-    v1((slot-1)*samples_slot1+1:slot*samples_slot1,1) = v0((slot-1)*samples_slot_agg+1:slot*samples_slot_agg-samples_slot2*nant2,1);
-    v2((slot-1)*samples_slot2+1:slot*samples_slot2,1) = v0((slot-1)*samples_slot_agg+samples_slot1+1:slot*samples_slot_agg-samples_slot2,1);
-    v2((slot-1)*samples_slot2+1:slot*samples_slot2,2) = v0((slot-1)*samples_slot_agg+samples_slot1+samples_slot2+1:slot*samples_slot_agg,1);
+        for a1=1:nant1
+            v1((slot-1)*samples_slot1+1:slot*samples_slot1,a1) = ...
+                v0((slot-1)*samples_slot_agg+(a1-1)*samples_slot1+1:...
+                   (slot-1)*samples_slot_agg+ a1   *samples_slot1,1);
+        end
+        for a2=1:nant2
+            v2((slot-1)*samples_slot2+1:slot*samples_slot2,a2) = ...
+                v0((slot-1)*samples_slot_agg+nant1*samples_slot1+(a2-1)*samples_slot2+1:...
+                   (slot-1)*samples_slot_agg+nant1*samples_slot1+ a2   *samples_slot2,1);
+        end
     end
     
-%%
+    if enable_plots>=2
         figure(1)
-plot(20*log10(abs(fftshift(fft(v1)))))
-figure(2)
-plot(20*log10(abs(fftshift(fft(v2)))))
+        plot(abs(fftshift(fft(v1))))
         
+        figure(2)
+        plot(abs(fftshift(fft(v2))))
+    end
     
-%% sync
-addpath('../../../openair1/PHY/LTE_REFSIG');
-primary_synch; %loads the primary sync signal
+    %% frame start detection
+    [corr,lag] = xcorr(v2(:,1),pss2_t);
+    %[m,idx]=max(abs(corr));
+    [m,idx]=peaksfinder(corr,frame_length2);
 
-pss_t = upsample(primary_synch0_time,4*2);
-[corr,lag] = xcorr(v2(:,1),pss_t);
-figure(3);
+    if (enable_plots>=2)
+        figure(20);
+        hold off
         plot(lag,abs(corr));
-
-%% frame start detection
-load('E:\EMOS\corridor\ofdm_pilots_sync_30MHz.mat');
-ofdm_symbol_length = num_carriers + prefix_length;
-frame_length = ofdm_symbol_length*num_symbols_frame;
-useful_carriers = num_carriers-num_zeros-1;
-
-% TODO: do a proper peak detection. for now just take the max
-[m,i]=max(abs(corr));
-frame_start = lag(i) - prefix_length;
-
-%% ofdm receiver
-received = v2(frame_start:frame_start+frame_length,1);
-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(2:useful_carriers/2+1); fblock(end-useful_carriers/2+1:end)];
+        hold on
+        plot(lag(idx),m,'ro')
     end
     
-%% channel estimation
-H=conj(f2).*received_f;
-Ht = ifft(H,[],2);
-PDP = mean(abs(Ht(2:end,:).^2),1);
+    %%
+    for i=1:size(idx,2)-1; % the last frame is not complite
+        fprintf(1,'.');
+        %frame_start2 = lag(i) - prefix_length2;
+        frame_start2 = lag(idx(i))-prefix_length2;
+        frame_start1 = frame_start2*2;
+        
+        % ofdm receiver
+        received_f1 = OFDM_RX(v1(frame_start1:frame_start1+frame_length1,:),num_carriers1,useful_carriers1,prefix_length1,num_symbols_frame);
+        received_f2 = OFDM_RX(v2(frame_start2:frame_start2+frame_length2,:),num_carriers2,useful_carriers2,prefix_length2,num_symbols_frame);
+        
+        % channel estimation (SISO)
+        H1=conj(squeeze(f1(1,:,:))).*received_f1(:,:,1);
+        H2=conj(squeeze(f2(1,:,:))).*received_f2(:,:,1);
+        H1t = ifft(H1,[],2);
+        H2t = ifft(H2,[],2);
+        PDP1 = mean(abs(H1t(2:end,:).^2),1);
+        PDP2 = mean(abs(H2t(2:end,:).^2),1);
+        PDP1_total((block-1)*nframes+i+1,:) = PDP1;
+        PDP2_total((block-1)*nframes+i+1,:) = PDP2;
+        
+        if enable_plots>=1
             figure(3)
-surf(20*log10(abs(Ht)))
+            surf((abs(H1t)))
             xlabel('time [OFDM symbol]')
             ylabel('delay time [samples]')
             zlabel('power [dB]')
+            title('H1t')
             shading interp
             figure(4)
-plot(10*log10(PDP))
+            plot(10*log10(PDP1))
             xlabel('delay time [samples]')
             ylabel('power [dB]')
+            title('PDP1')
+            figure(30)
+            surf((abs(H2t)))
+            xlabel('time [OFDM symbol]')
+            ylabel('delay time [samples]')
+            zlabel('power [dB]')
+            title('H2t')
+            shading interp
+            figure(40)
+            plot(10*log10(PDP2))
+            xlabel('delay time [samples]')
+            ylabel('power [dB]')
+            title('PDP2')
+        end
         
+    end
+    fprintf(1,'\n');
+    block = block+1;
+end
 
+fclose(fid);

targets/PROJECTS/CORRIDOR/peaksfinder.m

@@ -2,7 +2,7 @@ function [m,ind]=peaksfinder(corr,frame_length)
 
 
 
-threshold=45000;
+threshold=max(abs(corr))*0.75;
 consecutivePos=[];
 highCorrVal=find(abs(corr)>threshold);
 num_peak=0;