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@4440 818b1a75-f10b-46b9-bf7c-635c3b92a50f

Browse Source
This commit is contained in:
Florian Kaltenberger
2013-11-19 16:23:43 +00:00
parent 4f528b3805
commit ca8d901752
2 changed files with 12 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
targets/PROJECTS/TDDREC/runmeas_full_duplex.m
+9 -3
View File
@@ -65,6 +65,12 @@ if(paramsinitialized)
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);
chanestsA2B=zeros(301,Nantb,Nmeas);
tchanestsA2B=zeros(512,Nantb,Nmeas);
fchanestsA2B=zeros(512,Nantb,Nmeas);
chanestsB2A=zeros(301,Nantb,Nmeas);
tchanestsB2A=zeros(512,Nantb,Nmeas);
fchanestsB2A=zeros(512,Nantb,Nmeas);
for meas=1:Nmeas for meas=1:Nmeas
%% ------- Node A to B transmission ------- %% %% ------- Node A to B transmission ------- %%
@@ -115,9 +121,9 @@ if(paramsinitialized)
end 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.');
#end %end
if (chanest_full) if (chanest_full)
chanestsB2A(:,:,meas)=zeros(301,Nantb); chanestsB2A(:,:,meas)=zeros(301,Nantb);
targets/PROJECTS/TDDREC/runmeas_long_chanest.m
+3 -3
View File
@@ -120,9 +120,9 @@ if(paramsinitialized)
%% ------- Do the B to A channel estimation ------- %% %% ------- Do the B to A channel estimation ------- %%
HB2A=conj(Db2a_T.*repmat(Db2a_R,1,Nantb)); HB2A=conj(Db2a_T.*repmat(Db2a_R,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.');
#end %end
if (chanest_full) if (chanest_full)
chanestsB2A=zeros(301,Nantb); chanestsB2A=zeros(301,Nantb);