Cedric Roux
56c66c09aa
Pre-ci test passed git-svn-id: http://svn.eurecom.fr/openair4G/trunk@4000 818b1a75-f10b-46b9-bf7c-635c3b92a50f
120 lines
4.0 KiB
Matlab
120 lines
4.0 KiB
Matlab
load_channel = 1;
|
|
|
|
SNR=10^(SNRdB/10);
|
|
% MMSE regularization factor
|
|
mu = NTx/SNR;
|
|
% noise variance
|
|
n_p=1/SNR;
|
|
|
|
% open all files
|
|
for k=1:NUser
|
|
fid(k) = fopen(filename_emos{sched_users(k),idx},'r');
|
|
if (fid(k)==-1)
|
|
error('[eval_data_MU] Error opening file %s', filename_emos{sched_users(k),idx});
|
|
end
|
|
end
|
|
|
|
|
|
% read the data in chunks of 1000 Frames
|
|
NFrames_max = 1000;
|
|
NFrames_read = zeros(1,NUser);
|
|
NFrames_Tx = max(framestamp_max(:,idx))-min(framestamp_min(:,idx))+1;
|
|
NFreq = 40;
|
|
|
|
chunks = 0;
|
|
NFrames_Tx_read = 0; %Number of Tx frames read from file
|
|
NFrames_Tx_proc = 0; %Number of Tx frames processed (with all crc's positive)
|
|
|
|
align_matrix = false(NUser,NFrames_Tx);
|
|
for u = 1:NUser
|
|
align_matrix(u,est(u,idx).framestamp_tx - min(framestamp_min(:,idx)) + 1) = true;
|
|
end
|
|
|
|
crc_matrix = false(NUser,NFrames_Tx);
|
|
for u = 1:NUser
|
|
crc_matrix(u,est(u,idx).framestamp_tx - min(framestamp_min(:,idx)) + 1) = ~est(u,idx).err_ind;
|
|
end
|
|
|
|
NFrames_pos = sum(all(crc_matrix));
|
|
|
|
DProf = zeros(NUser, NTx, 64);
|
|
|
|
bit_vec = [4 8 12 16];
|
|
|
|
CAP_4U_ZF_CVQ = zeros(NFreq,NFrames_pos,length(bit_vec));
|
|
CAP_4U_MMSE_CVQ = zeros(NFreq,NFrames_pos,length(bit_vec));
|
|
|
|
%while any(NFrames_read < NFrames(:,idx).')
|
|
%fk 30.1.2008: changed the following two line since it produced an endless loop
|
|
while any(NFrames_read < sum(align_matrix,2).')
|
|
|
|
%Tx_frames = (NFrames_Tx_read+1) : min(NFrames_Tx_read + NFrames_max, NFrames_Tx);
|
|
Tx_frames = (NFrames_Tx_read+1) : min(NFrames_Tx_read + NFrames_max,size(align_matrix,2));
|
|
[H, H_fq, dummy_est, dummy_gps, NFrames_read] = load_estimates_MU_lm_new(fid, NTx, length(Tx_frames), NFrames_read, align_matrix(:,Tx_frames), IFversion);
|
|
if isempty(H_fq)
|
|
warning('[eval_data_MU] No data read! Terminating loop.');
|
|
break
|
|
end
|
|
NFrames_Tx_read = NFrames_Tx_read + size(H_fq,4);
|
|
|
|
%use only measurements where all users have positive CRC
|
|
H = H(:,:,:,all(crc_matrix(:,Tx_frames)),:);
|
|
H_fq = H_fq(:,:,:,all(crc_matrix(:,Tx_frames)),:);
|
|
|
|
chunks = chunks + 1;
|
|
disp(['[eval_data_MU] Chunk = ' num2str(chunks) ', NFrames_read = ' num2str(NFrames_read) ', NFrames_Tx_read = ' num2str(NFrames_Tx_read) ]);
|
|
|
|
% normalize H per user
|
|
% the mean squared Frobenius norm has been precomputed in est(k,idx).Hnorm
|
|
% in order to get E(|h_{i,j,k}|^2) = 1, we need
|
|
% size(H_fq) = [NRx,NTx,NFreq,NFrames,Nuser]
|
|
for k=1:NUser
|
|
H(:,:,:,:,k) = H(:,:,:,:,k)*sqrt(NTx*NRx)/sqrt(mean(est(sched_users(k),idx).Hnorm));
|
|
H_fq(:,:,:,:,k) = H_fq(:,:,:,:,k)*sqrt(NTx*NRx)/sqrt(mean(est(sched_users(k),idx).Hnorm));
|
|
end
|
|
|
|
disp('[eval_data_MU] Calculating SDMA capacities');
|
|
for k2=1:size(H_fq,4)
|
|
%[HSf,antenna_index] = antenna_selection(H_fq(:,:,:,k2,:));
|
|
for k1=1:size(H_fq,3)
|
|
HS = squeeze(H_fq(1,:,k1,k2,:)).';
|
|
for b = 1:length(bit_vec)
|
|
[CAP_4U_MMSE_CVQ(k1,NFrames_Tx_proc + k2,b), Hq] = MISO_CVQ(HS,n_p,mu,bit_vec(b));
|
|
[CAP_4U_ZF_CVQ(k1,NFrames_Tx_proc + k2,b), Hq] = MISO_CVQ(HS,n_p,0,bit_vec(b));
|
|
end
|
|
end
|
|
end
|
|
|
|
%% calculate and plot the PDP and frequency response over time
|
|
% The function Dprofile requires that size(H)=[N_t,N_tau,N_Tx,N_Rx]
|
|
% but load estimates returns size(H1) = [NRx,NTx,NTau,NFrames]
|
|
disp('[eval_data_MU] PDP Calculation');
|
|
for k=1:NUser
|
|
for l=1:NTx
|
|
DProf(k,l,:) = DProf(k,l,:) + reshape(DProfile(squeeze(permute(H(:,l,:,:,k),[4 3 2 1 5]))),1,1,[]);
|
|
end
|
|
end
|
|
|
|
|
|
NFrames_Tx_proc = NFrames_Tx_proc + size(H_fq,4);
|
|
disp(['[eval_data_MU] Chunk = ' num2str(chunks), ', NFrames_Tx_proc = ' num2str(NFrames_Tx_proc) ]);
|
|
|
|
end
|
|
|
|
% close all files
|
|
for k=1:NUser
|
|
fclose(fid(k));
|
|
end
|
|
|
|
DProf = DProf ./ chunks;
|
|
|
|
%% save calcualted data
|
|
filename_mat = fullfile(filepath,sprintf('results_CVQ_idx_%d.mat',idx));
|
|
if strfind(version,'7.3')
|
|
save(filename_mat,'-V7','CAP*','DProf');
|
|
else
|
|
save(filename_mat,'DProf','CAP*');
|
|
end
|
|
|
|
|