Commit cdabfe0a authored by ghaddab's avatar ghaddab

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

parent a03898b0
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
%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')
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment