LTE link level simulation

维也纳大学LTE链路级仿真代码，适合LTE初学者了解全系统架构

function varargout = jCapacity（Hf varargin）
% 
% [CTxWF CNoTx PfTxWF P N] = jCapacity（Hf P N）
%
% Computes the capacity in bits/s/Hz （bps/Hz） for frequency-selective 
% wireless channels defined as a set of parallel channels in the frequency
% domain. The capacity is computed with channel knowledge at the
% transmitter side by using water filling power allocation over frequency
% （CTxWF） and without channel knowledge at the transmitter side （except 
% the knowledge of the SNR） thus the power is equally allocated over all
% frequencies （CNoTx）.
%
% Input parameters:
%   Hf : attenuation values （channel coefficients）. The dimension of Hf is
%     assumed like this:
%       size（hn） = [L nRA nTA]
%     where:
%       L is the number of frequency bins （number of parallel channels in the
%         frequency domain.
%       nRA = number of receive antennas.
%       nTA = number of transmit antennas.
%   P : sum power in linear units. The total transmit power. 
%     It is very important to note that in the case of capacity withouth
%     channel knowledge at the transmitter side all frequencies allocate
%     the same amount of power equal to: P / size（Hf 1）
%   N : noise power in linear units.
%
% Output parameters:
%   CTxWF : resulting capacity in bps/Hz with channel knowledge at the
%     transmitter side （water filling capacity）.
%   CNoTx : resulting capacity in bps/Hz without channel knowledge at the
%     transmitter side.
%   PfTxWF: transmit power per frequency bin i.e. the water filling
%     coefficients. Restriction: sum（PfTxWF） <= P. The transmit power
%     coefficients for the case when no channel knowledge at the
%     transmitter is available are always equal to P / size（Hf 1）.
%   PN can be recovered again.
%
% [CTxWF CNoTx PfTxWF P N] = jCapacity（Hf SNR）
%
% The same as before but specifying the SNR instead of the power of
% the noise and the total transmit power.
%
% [CTxWF CNoTx PfTxWF P N] = jCapacity（Hf）
%
% Assumes SNR = 1 which means that the SNR is somehow included in the
% attenuation values instead of being normalized.
%

% vargargout initialisation
for index=1:nargout
    varargout{index} = -1;
end

if （nargin == 1）
    % No input parameter for the SNR. It means SNR == 1.
    P = 1;
    N = 1;
elseif （nargin == 2）
    % The input parameter is the SNR value.
    P = varargin{1};
    N = 1;
elseif （nargin == 3）
    % The input parameters are （by this order） the total transmit power and
    % the noise power both in linear units.
    P = varargin{1};
    N = varargin{2};
else
    error（‘ Wrong input parameters. ‘）;
end

L = size（Hf 1）;
nRA = size（Hf 2）;
nTA = size（Hf 3）;

% MIMO => SVD => spatial water filling
gs = zeros（Lmin（nRAnTA））;

for frqIndex = 1:L
    G = squeeze（Hf（frqIndex : :））;
    singValues = svd（G）;
    gs（frqIndex:） = singValues;
end

Hf = reshape（gs1[]）;



% The following code is valid

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-10-04 17:57  LTE_link_Level_1.7_r1089\
     目录           0  2011-10-04 17:57  LTE_link_Level_1.7_r1089\+channels\
     文件        2941  2011-10-04 16:36  LTE_link_Level_1.7_r1089\+channels\channelMatrixTrace.m
     文件        1872  2011-10-04 16:36  LTE_link_Level_1.7_r1089\+channels\uplinkChannel.m
     目录           0  2011-10-04 17:57  LTE_link_Level_1.7_r1089\+channel_resampling\
     文件        2978  2011-10-04 16:36  LTE_link_Level_1.7_r1089\+channel_resampling\nearestNeighborInterpolator.m
     文件        6706  2011-10-04 16:36  LTE_link_Level_1.7_r1089\+channel_resampling\sincInterpolator.m
     目录           0  2011-10-04 17:57  LTE_link_Level_1.7_r1089\+network_elements\
     文件        7373  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\AdaptiveFeedbackScheduler.m
     文件        3893  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\bestCqiSchedulerMIMO.m
     文件        4525  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\bestCqiSchedulerMIMOHARQ.m
     文件         861  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\clock.m
     文件        3247  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\eesmAverager.m
     文件       11490  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\eNodeB.m
     文件        6964  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\fixedScheduler.m
     文件         985  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\harqProcess.m
     文件       26424  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\lteScheduler.m
     文件        4218  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\miesmAverager.m
     目录           0  2011-10-04 17:57  LTE_link_Level_1.7_r1089\+network_elements\MI_data\
     文件      157173  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\BICMC_16.mat
     文件      157105  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\BICMC_4.mat
     文件       81514  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\BICMC_64.mat
     文件         492  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\BICM_k_d.mat
     文件         478  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\BICM_k_d_MSE.mat
     文件         466  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\Efficiency_k_d_MSE.mat
     文件       16416  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\MI_curves.fig
     文件       79431  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\MI_data\MI_curves.png
     文件        9685  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\ProportionalFairScheduler.m
     文件        4857  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\roundRobinScheduler.m
     文件         410  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\sinrAverager.m
     文件        6910  2011-10-04 16:35  LTE_link_Level_1.7_r1089\+network_elements\UE.m
............此处省略309个文件信息