From cdabfe0a798fe9b157124acbb384cad9632655b4 Mon Sep 17 00:00:00 2001
From: ghaddab <ghaddab@eurecom.fr>
Date: Wed, 26 Feb 2014 16:17:37 +0000
Subject: [PATCH] git-svn-id: http://svn.eurecom.fr/openair4G/trunk@5091
818b1a75-f10b-46b9-bf7c-635c3b92a50f
---
.../USERSPACE/OCTAVE/rx_spec_multicard_sync.m | 170 ++++++++++++++++++
.../USERSPACE/OCTAVE/txsig_multicard_sync.m | 121 +++++++++++++
2 files changed, 291 insertions(+)
create mode 100755 targets/ARCH/EXMIMO/USERSPACE/OCTAVE/rx_spec_multicard_sync.m
create mode 100644 targets/ARCH/EXMIMO/USERSPACE/OCTAVE/txsig_multicard_sync.m
diff --git a/targets/ARCH/EXMIMO/USERSPACE/OCTAVE/rx_spec_multicard_sync.m b/targets/ARCH/EXMIMO/USERSPACE/OCTAVE/rx_spec_multicard_sync.m
new file mode 100755
index 00000000000..18541edd269
--- /dev/null
+++ b/targets/ARCH/EXMIMO/USERSPACE/OCTAVE/rx_spec_multicard_sync.m
@@ -0,0 +1,170 @@
+dual_tx=0;
+card0=0;
+card1=1;
+limeparms;
+active_rf = zeros(1,4);
+active_rf = [1 1 1 1];
+autocal = [1 1 1 1];
+resampling_factor = [2 2 2 2];
+
+%rf_mode = (RXEN+TXEN+TXLPFNORM+TXLPFEN+TXLPF25+RXLPFNORM+RXLPFEN+RXLPF25+LNA1ON+LNAMax+RFBBNORM)*[1 1 1 1];
+rf_mode = (RXEN+TXLPFNORM+TXLPFEN+TXLPF25+RXLPFNORM+RXLPFEN+RXLPF25+LNA1ON+LNAMax+RFBBNORM)*active_rf;
+rf_mode = rf_mode+((DMAMODE_RX)*active_rf);
+%freq_rx = 2540000000*[1 1 1 1];
+freq_rx = 2600000000*active_rf;
+%freq_rx = 2680000000*active_rf; %+ 2540000000*[0 1 0 0]; % + 859500000*[0 0 1 0];
+%freq_rx = 1912600000*[1 1 1 1];
+%freq_rx = 859500000*[1 1 1 1];
+freq_tx = freq_rx;
+tx_gain = 0*active_rf; %[1 1 1 1];
+rx_gain = 0*active_rf; %1 1 1 1];
+rf_local= rf_local*active_rf; %[1 1 1 1];
+rf_rxdc = rf_rxdc*active_rf; %[1 1 1 1];
+%rf_vcocal=rf_vcocal_859*[1 1 1 1];
+rf_vcocal=rf_vcocal_19G*active_rf; %1 1 1 1];
+eNBflag = 0;
+tdd_config = DUPLEXMODE_FDD + TXRXSWITCH_TESTRX;
+%tdd_config = DUPLEXMODE_FDD + TXRXSWITCH_LSB;
+syncmode0 = SYNCMODE_MASTER;
+syncmode1 = SYNCMODE_SLAVE;
+rffe_rxg_low = 31*active_rf; %[1 1 1 1];
+rffe_rxg_final = 63*active_rf; %[1 1 1 1];
+rffe_band = B19G_TDD*active_rf; %[1 1 1 1];
+printf('Configuring Card 0\n');
+oarf_config_exmimo(card0,freq_rx,freq_tx,tdd_config,syncmode0,rx_gain,tx_gain,eNBflag,rf_mode,rf_rxdc,rf_local,rf_vcocal,rffe_rxg_low,rffe_rxg_final,rffe_band,autocal,resampling_factor)
+sleep(1);
+printf('Configuring Card 1\n');
+oarf_config_exmimo(card1,freq_rx,freq_tx,tdd_config,syncmode1,rx_gain,tx_gain,eNBflag,rf_mode,rf_rxdc,rf_local,rf_vcocal,rffe_rxg_low,rffe_rxg_final,rffe_band,autocal,resampling_factor)
+
+gpib_card=0; % first GPIB PCI card in the computer
+gpib_device=28; % this is configured in the signal generator Utilities->System
+cables_loss_dB = 6; % we need to account for the power loss between the signa
+power_dBm = -95;
+
+%gpib_send(gpib_card,gpib_device,['POW ' int2str(power_dBm+cables_loss_dB) 'dBm']);
+%gpib_send(gpib_card,gpib_device,'OUTP:STAT ON'); % activate output
+%vlen=76800;
+%v = 5000* [ zeros(1,1000) ones(1,500) zeros(1,500) ones(1,1) zeros(1,vlen-1000-1000-1)]; v=[v;v;v;v]'; oarf_send_frame(1,v,16);
+%sleep(1)
+
+s=oarf_get_frame(0);
+%s1=oarf_get_frame(1);
+
+%s = [s s1];
+
+if (resampling_factor(1)== 2)
+length1 = 76800;
+f0 = (7.68*(0:length(s(:,1))/4-1)/(length(s(:,1))/4))-3.84;
+axis([-3.84,3.84,40,200]);
+elseif (resampling_factor(1) == 1)
+length1 = 153600;
+f0 = (15.36*(0:length(s(:,1))/2-1)/(length(s(:,1))/2))-7.68;
+axis([-7.68,7.68,40,200]);
+else
+length1 = 307200;
+f0 = (30.72*(0:length(s(:,1))-1)/(length(s(:,1))))-15.36;
+axis([-15.36,15.36,40,200]);
+endif
+
+
+if (resampling_factor(2)== 2)
+length2 = 76800;
+f1 = (7.68*(0:length(s(:,2))/4-1)/(length(s(:,2))/4))-3.84;
+axis([-3.84,3.84,40,200]);
+elseif (resampling_factor(2) == 1)
+length2 = 153600;
+f1 = (15.36*(0:length(s(:,2))/2-1)/(length(s(:,2))/2))-7.68;
+axis([-7.68,7.68,40,200]);
+else
+length2 = 307200;
+f1 = (30.72*(0:length(s(:,2))-1)/(length(s(:,2))))-15.36;
+axis([-15.36,15.36,40,200]);
+endif
+
+if (resampling_factor(3)== 2)
+length3 = 76800;
+f2 = (7.68*(0:length(s(:,3))/4-1)/(length(s(:,3))/4))-3.84;
+axis([-3.84,3.84,40,200]);
+elseif (resampling_factor(3) == 1)
+length3 = 153600;
+f2 = (15.36*(0:length(s(:,3))/2-1)/(length(s(:,3))/2))-7.68;
+axis([-7.68,7.68,40,200]);
+else
+length3 = 307200;
+f2 = (30.72*(0:length(s(:,3))-1)/(length(s(:,3))))-15.36;
+axis([-15.36,15.36,40,200]);
+endif
+
+if (resampling_factor(4)== 2)
+length4 = 76800;
+f3 = (7.68*(0:length(s(:,4))/4-1)/(length(s(:,4))/4))-3.84;
+axis([-3.84,3.84,40,200]);
+elseif (resampling_factor(4) == 1)
+length4 = 153600;
+f3 = (15.36*(0:length(s(:,4))/2-1)/(length(s(:,4))/2))-7.68;
+axis([-7.68,7.68,40,200]);
+else
+length4 = 307200;
+f3 = (30.72*(0:length(s(:,4))-1)/(length(s(:,4))))-15.36;
+axis([-15.36,15.36,40,200]);
+endif
+
+length_slave = 76800;
+
+spec0 = 20*log10(abs(fftshift(fft(s(1:length1,1)))));
+spec1 = 20*log10(abs(fftshift(fft(s(1:length2,2)))));
+spec2 = 20*log10(abs(fftshift(fft(s(1:length3,3)))));
+spec3 = 20*log10(abs(fftshift(fft(s(1:length4,4)))));
+spec4 = 20*log10(abs(fftshift(fft(s(1:length_slave,5)))));
+spec5 = 20*log10(abs(fftshift(fft(s(1:length_slave,6)))));
+spec6 = 20*log10(abs(fftshift(fft(s(1:length_slave,7)))));
+spec7 = 20*log10(abs(fftshift(fft(s(1:length_slave,8)))));
+
+axis_time = [1:76800*4];
+
+
+clf
+if (active_rf(1) == 1)
+hold off;
+plot(f0',spec0,'r');
+figure
+plot(f0',spec4,'b');
+endif
+%hold on
+if (active_rf(2) == 1)
+hold off;
+figure;
+plot(f1',spec1,'b');
+endif
+%hold on
+if (active_rf(3) == 1)
+hold off;
+figure;
+plot(f2',spec2,'g');
+endif
+%hold on
+if (active_rf(4) == 1)
+hold off;
+figure;
+plot(f3',spec3,'m');
+endif
+%hold on
+hold off;
+figure;
+plot(real(s(1:length1,1)),'r');
+hold on
+plot(real(s(1:length_slave,5)),'b');
+hold off;
+figure;
+plot(imag(s(1:length1,1)),'r');
+hold on
+plot(imag(s(1:length_slave,5)),'b');
+
+
+%len=2500;off=0; chan=[1 5]; hold off ; figure; plot(real(s(off+1:off+len,chan)),'-o',"markersize",2) ; hold off ; figure; plot(imag(s(off+1:off+len,chan)),'r-o',"markersize",2);
+
+%hold off
+%axis([-3.84,3.84,40,160]);
+%gpib_send(gpib_card,gpib_device,'OUTP:STAT OFF'); % activate output
+%legend('Antenna Port 0','Antenna Port 1','Antenna Port 2','Antenna Port 3');
+grid
diff --git a/targets/ARCH/EXMIMO/USERSPACE/OCTAVE/txsig_multicard_sync.m b/targets/ARCH/EXMIMO/USERSPACE/OCTAVE/txsig_multicard_sync.m
new file mode 100644
index 00000000000..30a8ad0d4c0
--- /dev/null
+++ b/targets/ARCH/EXMIMO/USERSPACE/OCTAVE/txsig_multicard_sync.m
@@ -0,0 +1,121 @@
+%fc = 2660000000;
+fc = 1907600000;
+%fc = 859.5e6;
+
+rxgain=0;
+txgain=25;
+eNB_flag = 0;
+card = 0;
+active_rf = [1 1 0 0];
+autocal = [1 1 1 1];
+resampling_factor = [2 2 2 2];
+limeparms;
+rf_mode = (RXEN+TXEN+TXLPFNORM+TXLPFEN+TXLPF25+RXLPFNORM+RXLPFEN+RXLPF25+LNA1ON+LNAMax+RFBBNORM) * active_rf;
+rf_mode = rf_mode + (DMAMODE_RX + DMAMODE_TX)*active_rf;
+%rf_mode = RXEN+TXEN+TXLPFNORM+TXLPFEN+TXLPF25+RXLPFNORM+RXLPFEN+RXLPF25+LNA1ON+LNAByp+RFBBLNA1;
+%rf_local= [8253704 8253704 8257340 8257340]; %eNB2tx %850MHz
+%rf_local= [8255004 8253440 8257340 8257340]; % ex2 700 MHz
+rf_local = [8254744 8255063 8257340 8257340]; %eNB2tx 1.9GHz
+%rf_local = [8257292 8257300 8257340 8257340]; %ex2 850 MHz
+%rf_local = rf_local * chan_sel;
+rf_rxdc = rf_rxdc * active_rf;
+%rf_vcocal = rf_vcocal_859 * chan_sel;
+rf_vcocal = rf_vcocal_19G * active_rf;
+%rf_vcocal = rf_vcocal_26G_eNB * chan_sel;
+rxgain = rxgain*active_rf;
+txgain = txgain*active_rf;
+freq_tx = fc*active_rf;
+freq_rx = freq_tx;
+%freq_rx = freq_tx-120000000*chan_sel;
+%freq_tx = freq_rx+1920000;
+tdd_config = DUPLEXMODE_FDD + TXRXSWITCH_TESTTX;
+syncmode0 = SYNCMODE_MASTER;
+syncmode1 = SYNCMODE_SLAVE;
+rffe_rxg_low = 61*[1 1 1 1];
+rffe_rxg_final = 61*[1 1 1 1];
+rffe_band = B19G_TDD*[1 1 1 1];
+
+oarf_config_exmimo(0, freq_rx,freq_tx,tdd_config,syncmode0,rxgain,txgain,eNB_flag,rf_mode,rf_rxdc,rf_local,rf_vcocal,rffe_rxg_low,rffe_rxg_final,rffe_band,autocal,resampling_factor);
+sleep(1);
+oarf_config_exmimo(1, freq_rx,freq_tx,tdd_config,syncmode1,rxgain,txgain,eNB_flag,rf_mode,rf_rxdc,rf_local,rf_vcocal,rffe_rxg_low,rffe_rxg_final,rffe_band,autocal,resampling_factor);
+amp = pow2(14)-1;
+n_bit = 16;
+
+s = zeros(76800*4,4);
+
+select = 1;
+
+switch(select)
+
+case 1
+ s(:,1) = floor(amp * (exp(1i*2*pi*(0:((76800*4)-1))/7680)));
+ s(:,2) = floor(amp * (exp(1i*2*pi*(0:((76800*4)-1))/7680)));
+ s(:,3) = floor(amp * (exp(1i*2*pi*(0:((76800*4)-1))/7680)));
+ s(:,4) = floor(amp * (exp(1i*2*pi*(0:((76800*4)-1))/7680)));
+
+case 2
+ s(38400+128,1)= 80-1j*40;
+ s(38400+138,1)= -80+1j*40;
+
+case 3
+ for i=0:(12*5-1)
+ s((38401+640*i):(38400+640*(i+1)),1)=floor(linspace(-127,128,640)); %+1j*floor(linspace(128,-127,640));
+ end
+
+case 4
+ pss0_f0=[0,0,0,0,0,0,0,0,0,0,32767,0,-26120,-19785,11971,-30502,-24020,-22288,32117,6492,31311,9658,-16384,-28378,25100,-21063,-7292,-31946,20429,25618,14948,29158,11971,-30502,31311,9658,25100,-21063,-16384,28377,-24020,22287,32117,6492,-7292,31945,20429,25618,-26120,-19785,-16384,-28378,-16384,28377,-26120,-19785,-32402,4883,31311,-9659,32117,6492,-7292,-31946,32767,-1,25100,-21063,-24020,22287,-32402,4883,-32402,4883,-24020,22287,25100,-21063,32767,-1,-7292,-31946,32117,6492,31311,-9659,-32402,4883,-26120,-19785,-16384,28377,-16384,-28378,-26120,-19785,20429,25618,-7292,31945,32117,6492,-24020,22287,-16384,28377,25100,-21063,31311,9658,11971,-30502,14948,29158,20429,25618,-7292,-31946,25100,-21063,-16384,-28378,31311,9658,32117,6492,-24020,-22288,11971,-30502,-26120,-19785,32767,0,0,0,0,0,0,0,0,0,0,0];
+
+ pss0_f = pss0_f0(1:2:length(pss0_f0)) + sqrt(-1)*pss0_f0(2:2:length(pss0_f0));
+
+ pss0_f = [0 pss0_f(37:72) zeros(1,512-73) pss0_f(1:36)];
+ pss0_t = ifft(pss0_f);
+ pss0_t = [pss0_t((512-127):512) pss0_t];
+ pss0_t = [pss0_t zeros(1,512+36) pss0_t];
+ pss0_max = max(max(abs(real(pss0_t))),max(abs(imag(pss0_t))));
+
+ pss0_t_fp_re = floor(real(8192*pss0_t/pss0_max));
+ pss0_t_fp_im = floor(imag(8192*pss0_t/pss0_max));
+
+ s(38400+(1:length(pss0_t_fp_re)),1) = 2*floor(pss0_t_fp_re) + 2*sqrt(-1)*floor(pss0_t_fp_im);
+ s(38400+(1:length(pss0_t_fp_re)),2) = 2*floor(pss0_t_fp_re) + 2*sqrt(-1)*floor(pss0_t_fp_im);
+ s(38400+(1:length(pss0_t_fp_re)),3) = 2*floor(pss0_t_fp_re) + 2*sqrt(-1)*floor(pss0_t_fp_im);
+ s(38400+(1:length(pss0_t_fp_re)),4) = 2*floor(pss0_t_fp_re) + 2*sqrt(-1)*floor(pss0_t_fp_im);
+case 5
+ x=1:76800;
+
+ s(:,1) = 1i*8*(mod(x-1,4096)); % + 1i*(8*(mod(x-1,4096)));
+ s(:,2) = 8*(mod(x-1,4096)) + 1i*(8*(mod(x-1,4096)));
+ s(:,3) = 8*(mod(x-1,4096)) + 1i*(8*(mod(x-1,4096)));
+ s(:,4) = 8*(mod(x-1,4096)) + 1i*(8*(mod(x-1,4096)));
+% s(:,4) = 8*(rem(x-1,2048))-2**15 + 1i*(8*(rem(x-1,2048))-2**15);
+
+% s(:,4) = 8*(mod(x-1,4096)) + 1i*8*(mod(x-1,4096));
+
+ %s(:,1) = 8*floor(mod(x-1,76800)/75) + 1i*8*floor(mod(x-1,76800)/75);
+ %s(:,2) = 8*floor(mod(x-1,76800)/75) + 1i*8*floor(mod(x-1,76800)/75);
+ %s(:,3) = 8*floor(mod(x-1,76800)/75) + 1i*8*floor(mod(x-1,76800)/75);
+ %s(:,4) = 8*floor(mod(x-1,76800)/75) + 1i*8*floor(mod(x-1,76800)/75);
+
+otherwise
+ error('unknown case')
+endswitch
+
+%s = s*2 - 2**15 -1i*(2**15);
+s = s*2;
+
+
+%s(38400:end,1) = (1+1j);
+%s(38400:end,2) = (1+1j);
+
+sleep (1)
+%keyboard
+
+oarf_send_frame(0,s,n_bit);
+sleep(1);
+oarf_send_frame(1,s,n_bit);
+%r = oarf_get_frame(card);
+figure(1)
+hold off
+plot(real(s(:,1)),'r')
+%hold on
+%plot(imag(s(:,2)),'b')
--
GitLab