maki/openairinterface5g
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

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

Browse Source
This commit is contained in:
Florian Kaltenberger
2013-11-12 14:19:24 +00:00
parent f60f1d0e70
commit 0f333f9db8
1 changed files with 213 additions and 194 deletions
Download Patch File Download Diff File Expand all files Collapse all files
targets/PROJECTS/TDDREC/runmeas_full_duplex.m
+213 -194
View File
@@ -1,198 +1,217 @@
# % Author: Mirsad Cirkic, Florian Kaltenberger % Author: Mirsad Cirkic, Florian Kaltenberger
# % Organisation: Eurecom (and Linkoping University) % Organisation: Eurecom (and Linkoping University)
# % E-mail: mirsad.cirkic@liu.se % E-mail: mirsad.cirkic@liu.se
if(paramsinitialized && ~LSBSWITCH_FLAG) if(paramsinitialized && ~LSBSWITCH_FLAG)
disp(["\n\n------------\nThis code is, so far, only written for single runs. Multiple " ... disp(['\n\n------------\nThis code is, so far, only written for single runs. Multiple ' ...
"runs will overwrite the previous measurement data, i.e., the " ... 'runs will overwrite the previous measurement data, i.e., the ' ...
"data structures are not defined for multiple runs. You will need to " ... 'data structures are not defined for multiple runs. You will need to ' ...
"add code in order to save the intermediate measurements and the " ... 'add code in order to save the intermediate measurements and the ' ...
"corresponding timestamps.\n------------"]) 'corresponding timestamps.\n------------'])
N=76800; N=76800;
M=4; M=4;
indA=find(active_rfA==1); indA=find(active_rfA==1);
indB=find(active_rfB==1); indB=find(active_rfB==1);
Nanta=length(indA); Nanta=length(indA);
Nantb=length(indB); Nantb=length(indB);
if(Nanta!=1) if(Nanta~=1)
error("Node A can only have one antenna active\n"); error('Node A can only have one antenna active\n');
endif end
Niter=1; Niter=1;
if(Niter!=1) if(Niter~=1)
error("We should only use one get_frame at each run.\n"); error('We should only use one get_frame at each run.\n');
endif end
Nmeas = 10; Nmeas = 10;
# %% ------- Prepare the signals for A2B ---------- %% %% ------- Prepare the signals for A2B ---------- %%
signalA2B=zeros(N,4); signalA2B=zeros(N,4);
signalB2A=zeros(N,4); signalB2A=zeros(N,4);
ia=1; ib=1; ia=1; ib=1;
Db2a_T=[]; Db2a_T=[];
for i=1:4 for i=1:4
if(indA(ia)==i) if(indA(ia)==i)
[Da2b_T, tmps]=genrandpskseq(N,M,amp); [Da2b_T, tmps]=genrandpskseq(N,M,amp);
signalA2B(:,i)=tmps*2; %make sure LSB is 0 (switch=tx) signalA2B(:,i)=tmps*2; %make sure LSB is 0 (switch=tx)
signalB2A(:,i)=repmat(1+1j,76800,1); %make sure LSB is 1 (switch=rx) signalB2A(:,i)=repmat(1+1j,76800,1); %make sure LSB is 1 (switch=rx)
if(length(indA)> ia) ia=ia+1; endif if(length(indA)> ia) ia=ia+1; end
endif end
% if(indB(ib)==i) % if(indB(ib)==i)
% % This part could be improved by creating fully orthogonal sequences % % This part could be improved by creating fully orthogonal sequences
% [tmpd, tmps]=genrandpskseq(N,M,amp); % [tmpd, tmps]=genrandpskseq(N,M,amp);
% signalB2A(:,i)=tmps*2; % signalB2A(:,i)=tmps*2;
% signalA2B(:,i)=repmat(1+1j,76800,1); % signalA2B(:,i)=repmat(1+1j,76800,1);
% Db2a_T=[Db2a_T tmpd]; % Db2a_T=[Db2a_T tmpd];
% if(length(indB)> ib) ib=ib+1; endif % if(length(indB)> ib) ib=ib+1; end
% endif % end
endfor end
#%%------------Prepare the signals for B2A---------------%% %%------------Prepare the signals for B2A---------------%%
for i=1:4 for i=1:4
if(indB(ib)==i) if(indB(ib)==i)
[tmpd, tmps]=genrandpskseq(N,M,amp); [tmpd, tmps]=genrandpskseq(N,M,amp);
signalB2A(:,i)=tmps*2; %make sure LSB is 0 (switch=tx) signalB2A(:,i)=tmps*2; %make sure LSB is 0 (switch=tx)
signalA2B(:,i)=repmat(1+1j,76800,1); %make sure LSB is 1 (switch=rx) signalA2B(:,i)=repmat(1+1j,76800,1); %make sure LSB is 1 (switch=rx)
Db2a_T=[Db2a_T tmpd]; Db2a_T=[Db2a_T tmpd];
if(length(indB)> ib) ib=ib+1; endif if(length(indB)> ib) ib=ib+1; end
endif end
endfor end
if (!chanest_full) if (~chanest_full)
signalB2A(1:38400,3)=0; signalB2A(1:38400,3)=0;
signalB2A(38401:end,2)=0; signalB2A(38401:end,2)=0;
Db2a_T(1:60,302:end) = 0; Db2a_T(1:60,302:end) = 0;
Db2a_T(61:end,1:301) = 0; Db2a_T(61:end,1:301) = 0;
end end
Da2b_R=zeros(Niter*120,Nantb*301,Nmeas); Da2b_R=zeros(Niter*120,Nantb*301,Nmeas);
Db2a_R=zeros(Niter*120,Nanta*301,Nmeas); Db2a_R=zeros(Niter*120,Nanta*301,Nmeas);
for meas=1:Nmeas for meas=1:Nmeas
# %% ------- Node A to B transmission ------- %% %% ------- Node A to B transmission ------- %%
oarf_send_frame(card,signalA2B,n_bit); oarf_send_frame(card,signalA2B,n_bit);
%keyboard %keyboard
sleep(0.01); sleep(0.01);
receivedA2B=oarf_get_frame(card); receivedA2B=oarf_get_frame(card);
%oarf_stop(card); %not good, since it does a reset %oarf_stop(card); %not good, since it does a reset
sleep(0.01); sleep(0.01);
#%%----------Node B to A transmission---------%% %%----------Node B to A transmission---------%%
oarf_send_frame(card,signalB2A,n_bit); oarf_send_frame(card,signalB2A,n_bit);
%keyboard %keyboard
sleep(0.01); sleep(0.01);
receivedB2A=oarf_get_frame(card); receivedB2A=oarf_get_frame(card);
%oarf_stop(card); %not good, since it does a reset %oarf_stop(card); %not good, since it does a reset
# %% ------- Do the A to B channel estimation ------- %% %% ------- Do the A to B channel estimation ------- %%
for i=0:119; for i=0:119;
ifblock=receivedA2B(i*640+[1:640],indB); ifblock=receivedA2B(i*640+[1:640],indB);
ifblock(1:128,:)=[]; ifblock(1:128,:)=[];
fblock=fft(ifblock); fblock=fft(ifblock);
fblock(1,:)=[]; fblock(1,:)=[];
fblock(151:360,:)=[]; fblock(151:360,:)=[];
Da2b_R((Niter-1)*120+i+1,:,meas)=vec(fblock); Da2b_R((Niter-1)*120+i+1,:,meas)=vec(fblock);
endfor end
HA2B=repmat(conj(Da2b_T),Niter,Nantb).*Da2b_R(:,:,meas); HA2B=repmat(conj(Da2b_T),Niter,Nantb).*Da2b_R(:,:,meas);
phasesA2B=unwrap(angle(HA2B)); phasesA2B=unwrap(angle(HA2B));
if(mean(var(phasesA2B))>0.5) if(mean(var(phasesA2B))>0.5)
disp("The phases of your estimates from A to B are a bit high (larger than 0.5 rad.), something is wrong."); disp('The phases of your estimates from A to B are a bit high (larger than 0.5 rad.), something is wrong.');
endif end
chanestsA2B(:,:,meas)=reshape(diag(repmat(Da2b_T,Niter,Nantb)'*Da2b_R(:,:,meas))/size(Da2b_T,1),301,Nantb); chanestsA2B(:,:,meas)=reshape(diag(repmat(Da2b_T,Niter,Nantb)'*Da2b_R(:,:,meas))/size(Da2b_T,1),301,Nantb);
#fchanestsA2B=zeros(512,Nantb); %fchanestsA2B=zeros(512,Nantb);
#for i=1:Nantb %for i=1:Nantb
# fchanestsA2B(:,i)=[0; chanestsA2B([1:150],i,meas); zeros(210,1); chanestsA2B(151:301,i,meas)]; % fchanestsA2B(:,i)=[0; chanestsA2B([1:150],i,meas); zeros(210,1); chanestsA2B(151:301,i,meas)];
#endfor %end
fchanestsA2B(:,:,meas)=[zeros(1,Nantb); chanestsA2B([1:150],:,meas); zeros(210,Nantb); chanestsA2B(151:301,:,meas)]; fchanestsA2B(:,:,meas)=[zeros(1,Nantb); chanestsA2B([1:150],:,meas); zeros(210,Nantb); chanestsA2B(151:301,:,meas)];
tchanestsA2B(:,:,meas)=ifft(fchanestsA2B(:,:,meas)); tchanestsA2B(:,:,meas)=ifft(fchanestsA2B(:,:,meas));
%% ------- Do the B to A channel estimation ------- %% %% ------- Do the B to A channel estimation ------- %%
for i=0:119; for i=0:119;
ifblock=receivedB2A(i*640+[1:640],indA); ifblock=receivedB2A(i*640+[1:640],indA);
ifblock(1:128,:)=[]; ifblock(1:128,:)=[];
fblock=fft(ifblock); fblock=fft(ifblock);
fblock(1,:)=[]; fblock(1,:)=[];
fblock(151:360,:)=[]; fblock(151:360,:)=[];
Db2a_R((Niter-1)*120+i+1,:,meas)=fblock.'; Db2a_R((Niter-1)*120+i+1,:,meas)=fblock.';
endfor end
HB2A=conj(repmat(Db2a_T,Niter,1)).*repmat(Db2a_R(:,:,meas),1,Nantb); HB2A=conj(repmat(Db2a_T,Niter,1)).*repmat(Db2a_R(:,:,meas),1,Nantb);
phasesB2A=unwrap(angle(HB2A)); phasesB2A=unwrap(angle(HB2A));
if(mean(var(phasesB2A))>0.5) if(mean(var(phasesB2A))>0.5)
disp("The phases of your estimates from B to A are a bit high (larger than 0.5 rad.), something is wrong."); disp('The phases of your estimates from B to A are a bit high (larger than 0.5 rad.), something is wrong.');
endif end
if (chanest_full) if (chanest_full)
chanestsB2A(:,:,meas)=zeros(301,Nantb); chanestsB2A(:,:,meas)=zeros(301,Nantb);
inds=repmat([1:Nantb]',1,301); inds=repmat([1:Nantb]',1,301);
for ci=1:301; for ci=1:301;
data=Db2a_T(:,ci+[0:Nantb-1]*301); data=Db2a_T(:,ci+[0:Nantb-1]*301);
rec=Db2a_R(:,ci,meas); rec=Db2a_R(:,ci,meas);
chanestsB2A(ci,:,meas)=(inv(data'*data)*data'*rec).'; chanestsB2A(ci,:,meas)=(inv(data'*data)*data'*rec).';
endfor end
else else
chanestsB2A(:,:,meas)=reshape(diag(repmat(Db2a_T,Niter,1)'*repmat(Db2a_R(:,:,meas),1,Nantb)/(Niter*60)),301,Nantb); chanestsB2A(:,:,meas)=reshape(diag(repmat(Db2a_T,Niter,1)'*repmat(Db2a_R(:,:,meas),1,Nantb)/(Niter*60)),301,Nantb);
end end
#fchanestsB2A=zeros(512,Nantb); %fchanestsB2A=zeros(512,Nantb);
#for i=1:Nantb %for i=1:Nantb
# fchanestsB2A(:,i)=[0; chanestsB2A([1:150],i); zeros(210,1); chanestsB2A(151:301,i)]; % fchanestsB2A(:,i)=[0; chanestsB2A([1:150],i); zeros(210,1); chanestsB2A(151:301,i)];
#endfor %end
fchanestsB2A(:,:,meas) = [zeros(1,Nantb); chanestsB2A([1:150],:,meas); zeros(210,Nantb); chanestsB2A(151:301,:,meas)]; fchanestsB2A(:,:,meas) = [zeros(1,Nantb); chanestsB2A([1:150],:,meas); zeros(210,Nantb); chanestsB2A(151:301,:,meas)];
tchanestsB2A(:,:,meas)=ifft(fchanestsB2A(:,:,meas)); tchanestsB2A(:,:,meas)=ifft(fchanestsB2A(:,:,meas));
end end
%% -- Some plotting code -- %% (you can uncomment what you see fit) %% -- Some plotting code -- %% (you can uncomment what you see fit)
received = receivedB2A; received = receivedB2A;
phases = phasesB2A; phases = phasesB2A;
tchanests = [tchanestsA2B(:,:,end), tchanestsB2A(:,:,end)]; tchanests = [tchanestsA2B(:,:,end), tchanestsB2A(:,:,end)];
fchanests = [fchanestsA2B(:,:,end), fchanestsB2A(:,:,end)]; fchanests = [fchanestsA2B(:,:,end), fchanestsB2A(:,:,end)];
clf clf
figure(1) figure(1)
for i=1:4 for i=1:4
subplot(220+i);plot(20*log10(abs(fftshift(fft(received(:,i)))))); subplot(220+i);plot(20*log10(abs(fftshift(fft(received(:,i))))));
endfor end
figure(2) figure(2)
t=[0:512-1]/512*1e-2; t=[0:512-1]/512*1e-2;
plot(t,20*log10(abs(tchanests))) plot(t,20*log10(abs(tchanests)))
xlabel('time')
ylabel('|h|')
legend('A->B1','A->B2','A->B3','B1->A','B2->A','B3->A');
%legend('A->B1','A->B2','B1->A','B2->A');
figure(4)
plot(20*log10(abs(fchanests)));
ylim([40 100])
xlabel('freq')
ylabel('|h|')
legend('A->B1','A->B2','A->B3','B1->A','B2->A','B3->A');
%legend('A->B1','A->B2','B1->A','B2->A');
if (0)
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') xlabel('time')
ylabel('phase') ylabel('|h|')
ylim([-pi pi]) legend('A->B1','A->B2','A->B3','B1->A','B2->A','B3->A');
drawnow; %legend('A->B1','A->B2','B1->A','B2->A');
pause(0.1)
endfor figure(4)
phavar=var(phases); plot(20*log10(abs(fchanests)));
plotphavar=[]; ylim([40 100])
for i=0:Nantb-1 xlabel('freq')
plotphavar=[plotphavar; phavar([1:301]+i*301)]; ylabel('|h|')
endfor legend('A->B1','A->B2','A->B3','B1->A','B2->A','B3->A');
plot([1:150 362:512],plotphavar,'o'); %legend('A->B1','A->B2','B1->A','B2->A');
%ylim([0 pi])
xlabel('subcarrier') if (0)
ylabel('phase variance') figure(3)
end 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)
end
phavar=var(phases);
plotphavar=[];
for i=0:Nantb-1
plotphavar=[plotphavar; phavar([1:301]+i*301)];
end
plot([1:150 362:512],plotphavar,'o');
%ylim([0 pi])
xlabel('subcarrier')
ylabel('phase variance')
end
%% estimate F matrix assuming it is diagonal for sanity checking
Fhatloc = zeros(Nmeas,301,Nantb);
for t=1:Nmeas
for s=1:301
ya=chanestsB2A(s,:,t);
yb=chanestsA2B(s,:,t);
Fhatloc(t,s,:)=(yb.*conj(ya))./(ya.*conj(ya));
end
end
figure(5)
plot_style={'rx','bo','gs'};
hold off
for n=1:Nantb
plot((squeeze(Fhatloc(:,:,n))),plot_style{n})
hold on
end
else else
if(LSBSWITCH_FLAG) error("You have to unset the LSB switch flag (LSBSWITCH_FLAG) in initparams.m.\n") if(LSBSWITCH_FLAG) error('You have to unset the LSB switch flag (LSBSWITCH_FLAG) in initparams.m.\n')
else error("You have to run init.params.m first!") else error('You have to run init.params.m first!')
endif end
endif end