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個文件信息