diff --git a/targets/PROJECTS/TDDREC/beamforming.m b/targets/PROJECTS/TDDREC/beamforming.m
index 1e8c68d65b..aeafa27396 100644
--- a/targets/PROJECTS/TDDREC/beamforming.m
+++ b/targets/PROJECTS/TDDREC/beamforming.m
@@ -3,6 +3,8 @@ signalA2B=zeros(N,4);
 signalB2A=zeros(N,4);
 signalB2A_prec_full=zeros(N,4);
 signalB2A_prec_diag=zeros(N,4);
+Da2b_R=zeros(Niter*120,Nantb*301);
+Db2a_R=zeros(Niter*120,Nanta*301);
 
 %%----------Node A to B transmission---------%%
 for i=1:4
diff --git a/targets/PROJECTS/TDDREC/calibration.m b/targets/PROJECTS/TDDREC/calibration.m
index abadc4b5ed..0ab5473827 100644
--- a/targets/PROJECTS/TDDREC/calibration.m
+++ b/targets/PROJECTS/TDDREC/calibration.m
@@ -18,16 +18,16 @@ end
 %% calculate full F matrix
 Fs = zeros(Nantb,Nantb,301);
 for s=1:301
-	[F, HA]=alterproj(CHB2A,CHA2B, s, Nantb, N_loc, Nmeas); 
-    Fs(:,:,s)=F;
-	%HAs{s}=HA;
+  [F, HA]=alterproj(CHB2A,CHA2B, s, Nantb, N_loc, Nmeas); 
+  Fs(:,:,s)=F;
+  %HAs{s}=HA;
 end
 
 %% calculate the diag F matrix based on last position
 Fds = zeros(Nantb,Nantb,301);
 for s=1:301
     CHA = squeeze(chanestsA2B(s,:,:));
-    CHB = aqueeze(chanestsB2A(s,:,:));
+    CHB = squeeze(chanestsB2A(s,:,:));
     Fds(:,:,s)=diag(diag(CHB*CHA')./diag(CHA*CHA'));
 end