From 4380042cf804db6e9442d9011eb13fdd62e0e98b Mon Sep 17 00:00:00 2001
From: Cirkic <cirkic@eurecom.fr>
Date: Thu, 18 Jul 2013 13:40:28 +0000
Subject: [PATCH] Added the tdd rec. script files
git-svn-id: http://svn.eurecom.fr/openair4G/trunk@4037 818b1a75-f10b-46b9-bf7c-635c3b92a50f
---
targets/PROJECTS/TDDREC/genorthpskseq.m | 32 ++++++
targets/PROJECTS/TDDREC/genorthqpskseq.m | 40 +++++++
targets/PROJECTS/TDDREC/genrandpskseq.m | 24 +++++
targets/PROJECTS/TDDREC/initparams.m | 42 ++++++++
targets/PROJECTS/TDDREC/runmeas.m | 130 +++++++++++++++++++++++
5 files changed, 268 insertions(+)
create mode 100644 targets/PROJECTS/TDDREC/genorthpskseq.m
create mode 100644 targets/PROJECTS/TDDREC/genorthqpskseq.m
create mode 100644 targets/PROJECTS/TDDREC/genrandpskseq.m
create mode 100644 targets/PROJECTS/TDDREC/initparams.m
create mode 100644 targets/PROJECTS/TDDREC/runmeas.m
diff --git a/targets/PROJECTS/TDDREC/genorthpskseq.m b/targets/PROJECTS/TDDREC/genorthpskseq.m
new file mode 100644
index 0000000000..184485e871
--- /dev/null
+++ b/targets/PROJECTS/TDDREC/genorthpskseq.m
@@ -0,0 +1,32 @@
+function [carrierdata, s]=genorthqpskseq(Ns,N,amp)
+
+if(N!=512*150)
+ error('The sequence length must be 76800.');
+endif
+
+s = zeros(N,Ns);
+H=1; for k=1:log2(128) H=[H H; H -H]; end; H=H(:,1:120);
+i=1; while(i<size(H,1)) h=H(i,:); inds=find(h*H'!=0); H(inds,:)=[]; H=[h; H]; i=i+1; end
+Hc=H+sqrt(-1)*H;
+if(sum(Hc*Hc'-240*eye(8))>0) error("The code is not orhtogonal\n"); endif
+
+carrierdata=zeros(120,Ns*301);
+orthinds=1:8;
+for i=1:301
+ for k=1:Ns
+ randind=ceil((9-k)*rand());
+ carrierdata(:,i+(k-1)*301)=Hc(orthinds(randind),:)';
+ orthinds(randind)=[];
+ endfor
+endfor
+
+for k=1:Ns
+ for i=0:119
+ fblock=[0 carrierdata(i,1:150) zeros(1,210) carrierdata(i,151:301)];
+ ifblock=ifft(fblock,512);
+ block = [ifblock(end-127:end) ifblock]; # Cycl. prefix
+ s([1:640]+i*640,k)=floor(amp*block);
+ endfor
+endfor
+
+endfunction
\ No newline at end of file
diff --git a/targets/PROJECTS/TDDREC/genorthqpskseq.m b/targets/PROJECTS/TDDREC/genorthqpskseq.m
new file mode 100644
index 0000000000..b42d449d7f
--- /dev/null
+++ b/targets/PROJECTS/TDDREC/genorthqpskseq.m
@@ -0,0 +1,40 @@
+function [carrierdata, s]=genorthqpskseq(Ns,N,amp)
+
+if(N!=512*150)
+ error('The sequence length must be 76800.');
+endif
+
+s = zeros(N,Ns);
+H=1; for k=1:log2(128) H=[H H; H -H]; end; H=H(:,1:120);
+i=1; while(i<size(H,1)) h=H(i,:); inds=find(h*H'!=0); H(inds,:)=[]; H=[h; H]; i=i+1; end
+Hc=H+sqrt(-1)*H;
+if(sum(Hc*Hc'-240*eye(8))>0) error("The code is not orhtogonal\n"); endif
+
+carrierdata=zeros(120,Ns*301);
+
+inds=1:8;
+ind=8;
+for i=1:301
+ for k=1:Ns
+ carrierdata(:,i+(k-1)*301)=Hc(inds(ind),:)';
+ inds(ind)=[];
+ ind=ind-1;
+ if(ind==0)
+ inds=1:8;
+ ind=8;
+ endif
+ endfor
+endfor
+
+for k=1:Ns
+ frstgroup=(k-1)*301+[1:150]; # The first group of OFDM carriers
+ sndgroup=(k-1)*301+[151:301]; # The second group of OFDM carriers
+ for i=0:119
+ fblock=[0 carrierdata(i+1,frstgroup) zeros(1,210) carrierdata(i+1,sndgroup)];
+ ifblock=ifft(fblock,512);
+ block = [ifblock(end-127:end) ifblock]; # Cycl. prefix
+ s([1:640]+i*640,k)=floor(amp*block);
+ endfor
+endfor
+
+endfunction
\ No newline at end of file
diff --git a/targets/PROJECTS/TDDREC/genrandpskseq.m b/targets/PROJECTS/TDDREC/genrandpskseq.m
new file mode 100644
index 0000000000..8d0030c13c
--- /dev/null
+++ b/targets/PROJECTS/TDDREC/genrandpskseq.m
@@ -0,0 +1,24 @@
+function [carrierdata, s]=genrandpskseq(N,M,amp)
+
+if(mod(N,640)~=0)
+ error('The sequence length must be divisible with 640.');
+endif
+s = zeros(N,1);
+MPSK=exp(sqrt(-1)*([1:M]*2*pi/M+pi/M));
+
+% OFDM sequence with 512 FFT using randomly
+% generated 4-QAM and 128 cyclic prefix
+carrierdata=zeros(120,301);
+for i=0:(N/640-1)
+ datablock=MPSK(ceil(rand(301,1)*M));
+ for j=1:301
+ carrierdata(i+1,j)=datablock(j);
+ endfor
+ fblock=[0 datablock(1:150) zeros(1,210) datablock(151:301)];
+ ifblock=ifft(fblock,512);
+ % Adding cycl. prefix making the block of 640 elements
+ block = [ifblock(end-127:end) ifblock];
+ s([1:640]+i*640)=floor(amp*block);
+endfor
+
+endfunction
\ No newline at end of file
diff --git a/targets/PROJECTS/TDDREC/initparams.m b/targets/PROJECTS/TDDREC/initparams.m
new file mode 100644
index 0000000000..222683a9a2
--- /dev/null
+++ b/targets/PROJECTS/TDDREC/initparams.m
@@ -0,0 +1,42 @@
+clear
+paramsinitialized=false;
+limeparms;
+rxgain=30;
+txgain=25;
+eNB_flag = 0;
+card = 0;
+Ntrx=4;
+dual_tx=0;
+active_rfA=[1 0 0 0];
+active_rfB=[0 1 1 0];
+active_rf=active_rfA+active_rfB;
+
+%fc = 2660000000;
+fc = 1912600000; %1907600000;
+%fc = 859.5e6;
+
+
+%rf_mode=(RXEN+TXEN+TXLPFNORM+TXLPFEN+TXLPF25+RXLPFNORM+RXLPFEN+RXLPF25+LNA1ON+LNAByp+RFBBLNA1) * active_rf;
+autocal_mode=active_rf;
+rf_mode=(TXLPFNORM+TXLPFEN+TXLPF25+RXLPFNORM+RXLPFEN+RXLPF25+LNA1ON+LNAMax+RFBBNORM) * active_rf;
+tdd_config = DUPLEXMODE_FDD + TXRXSWITCH_TESTTX;
+%tdd_config = DUPLEXMODE_FDD + TXRXSWITCH_LSB;
+syncmode = SYNCMODE_FREE;
+rf_local = [8254744 8255063 8257340 8257340]; %eNB2tx 1.9GHz
+rf_vcocal=rf_vcocal_19G*active_rf;
+
+rffe_rxg_low = 61*active_rf;
+rffe_rxg_final = 61*active_rf;
+rffe_band = B19G_TDD*active_rf;
+
+rf_rxdc = rf_rxdc*active_rf;
+
+freq_rx = fc*active_rf;
+freq_tx = freq_rx; %+1.92e6;
+tx_gain = txgain*active_rf;
+rx_gain = rxgain*active_rf;
+oarf_config_exmimo(card, freq_rx,freq_tx,tdd_config,syncmode,rx_gain,tx_gain,eNB_flag,rf_mode,rf_rxdc,rf_local,rf_vcocal,rffe_rxg_low,rffe_rxg_final,rffe_band,autocal_mode);
+autocal_mode=0*active_rf; % Autocalibration is only needed the first time we conf. exmimo
+amp = pow2(14)-1;
+n_bit = 16;
+paramsinitialized=true;
diff --git a/targets/PROJECTS/TDDREC/runmeas.m b/targets/PROJECTS/TDDREC/runmeas.m
new file mode 100644
index 0000000000..56ddd9d9c7
--- /dev/null
+++ b/targets/PROJECTS/TDDREC/runmeas.m
@@ -0,0 +1,130 @@
+if(paramsinitialized)
+ N=76800;
+ M=4;
+ signalB2A=zeros(N,1);
+ signalA2B=zeros(N,4);
+ indA=find(active_rfA==1);
+ indB=find(active_rfB==1);
+ Nanta=length(indA);
+ Nantb=length(indB);
+ Niter=1;
+ # %% ------- Node A to B channel measurements ------- %%
+ # rf_mode_current = rf_mode + (DMAMODE_TX+TXEN)*active_rfA +(DMAMODE_RX+RXEN)*active_rfB;
+ # oarf_config_exmimo(card, freq_rx,freq_tx,tdd_config,syncmode,rx_gain,tx_gain,eNB_flag,rf_mode_current,rf_rxdc,rf_local,rf_vcocal,rffe_rxg_low,rffe_rxg_final,rffe_band,autocal_mode);
+ # [Da2b_T, signalA2B]=genrandpskseq(N,M,amp);
+ # signalA2B=repmat(signalA2B,1,4);
+ # oarf_send_frame(card,signalA2B,n_bit);
+ # fprintf('A to B: Frame sent to Tx buffer.\n')
+ # fprintf('A to B: Acquiring chan. measurements....\n')
+ # loop=true;
+ # chanests=[];
+ # Da2b_R=zeros(Niter*120,Nantb*301);
+ # Da2b_T=repmat(Da2b_T,Niter,Nantb);
+ # for iter=0:(Niter-1) % Do the chan. est until SNR is good enough. That is, do the "oarf_get_fram(card)
+ # receivedA2B=oarf_get_frame(card);
+ # fprintf('A to B: Frame grabbed from Rx buffer.\n')
+ # for i=0:119;
+ # ifblock=receivedA2B(i*640+[1:640],indB);
+ # ifblock(1:128,:)=[];
+ # fblock=fft(ifblock);
+ # fblock(1,:)=[];
+ # fblock(151:360,:)=[];
+ # Da2b_R(iter*120+i+1,:)=vec(fblock);
+ # endfor
+ # endfor
+ # H=conj(Da2b_T).*Da2b_R;
+ # phases=unwrap(angle(H));
+ # chanests=diag(Da2b_T'*Da2b_R)/size(Da2b_T,1);
+ # for i=0:(Nantb-1)
+ # fchanests(:,i+1)=[0; chanests(301*i+[1:150]); zeros(210,1); chanests(301*i+[151:301])];
+ # endfor
+ # tchanests=ifft(fchanests);
+
+ # %% -- Do some plotting
+ # clf
+ # figure(1)
+ # for i=1:4
+ # subplot(220+i);plot(20*log10(abs(fftshift(fft(receivedA2B(:,i))))));
+ # endfor
+
+ # figure(2)
+ # t=[0:512-1]/512*1e-2;
+ # plot(t,abs(tchanests))
+ # xlabel('time')
+ # ylabel('|h|')
+
+ # figure(3)
+ # % wndw = 50;
+ # % for i=1:5:Nantb*301 %# sliding window size
+ # % phamean = filter(ones(wndw,1)/wndw, 1, phases(:,i)); %# moving average
+ # % plot(phamean(wndw:end),'LineWidth',2);
+ # % title(['subcarrier ' num2str(i)]);
+ # % xlabel('time')
+ # % ylabel('phase')
+ # % ylim([-pi pi])
+ # % drawnow;
+ # % pause(0.1)
+ # % endfor
+ # phavar=var(phases);
+ # plotphavar=[];
+ # for i=0:Nantb-1
+ # plotphavar=[plotphavar; phavar([1:301]+i*301)];
+ # endfor
+ # plot([1:150 362:512],plotphavar,'o');
+ # %ylim([0 pi])
+ # xlabel('subcarrier')
+ # ylabel('phase variance')
+
+
+ # figure(4)
+ # plot(20*log10(abs(fchanests))), ylim([40 100])
+
+ # %end
+ # oarf_stop(card);
+ # fprintf(' done\n')
+
+ %% ------- Node B to A channel measurements ------- %%
+ rf_mode_current = rf_mode + (DMAMODE_TX+TXEN)*active_rfB +(DMAMODE_RX+RXEN)*active_rfA;
+ oarf_config_exmimo(card, freq_rx,freq_tx,tdd_config,syncmode,rx_gain,tx_gain,eNB_flag,rf_mode_current,rf_rxdc,rf_local,rf_vcocal,rffe_rxg_low,rffe_rxg_final,rffe_band,autocal_mode);
+ signalB2A=zeros(N,4);
+ Db2a_R=zeros(Niter*120,Nanta*301);
+ [Db2a_T,tmps]=genorthqpskseq(Nantb,N,amp);
+ ind=1;
+ for i=1:4
+ if(indB(ind)==i)
+ signalB2A(:,i)=tmps(:,ind);
+ if(length(indB)> ind) ind=ind+1; endif
+ endif
+ endfor
+ %receivedB2A=sum(signalB2A*diag([1 2 3 4]),2)+100*randn(N,1);
+
+ oarf_send_frame(card,signalB2A,n_bit);
+ if(length(indA)!=1) error("Node A can only have one antenna active\n"); endif
+ for iter=0:(Niter-1) % Do the chan. est until SNR is good enough. That is, do the "oarf_get_fram(card)
+ fprintf('B to A: Frame grabbed from Rx buffer.\n')
+ receivedB2A=oarf_get_frame(card);
+
+ for i=0:119;
+ ifblock=receivedB2A(i*640+[1:640],indA);
+ ifblock(1:128,:)=[];
+ fblock=fft(ifblock);
+ fblock(1,:)=[];
+ fblock(151:360,:)=[];
+ Db2a_R(iter*120+i+1,:)=fblock.';
+ endfor
+ endfor
+ oarf_stop(card);
+
+ # # -- Do the channel estimation
+ chanests=reshape(diag(repmat(Db2a_T,Niter,1)'*repmat(Db2a_R,1,Nantb)/(Niter*120)),301,Nantb)
+ H=conj(repmat(Db2a_T,Niter,1)).*repmat(Db2a_R,1,Nantb);
+ phases=unwrap(angle(H));
+ # for i=0:(Nantb-1)
+ # fchanests(:,i+1)=[0; chanests(301*i+[1:150]); zeros(210,1); chanests(301*i+[151:301])];
+ # endfor
+ # tchanests=ifft(fchanests);
+
+
+else
+ error('You have to run init.params.m first!')
+end
--
GitLab