ADCFREE

 compute the free diffusion coefficients and the free signal
 
 Input:
     1. bvalues
     2. DIFF_cmpts
     3. VOL_cmpts
     4. IC_cmpts
 
 Output:
     1. Sig_free
     2. ADC_free_allcmpts

0001 function [Sig_free,ADC_free_allcmpts] = ADCFREE(bvalues,DIFF_cmpts,VOL_cmpts,IC_cmpts)
0002 
0003 % compute the free diffusion coefficients and the free signal
0004 %
0005 % Input:
0006 %     1. bvalues
0007 %     2. DIFF_cmpts
0008 %     3. VOL_cmpts
0009 %     4. IC_cmpts
0010 %
0011 % Output:
0012 %     1. Sig_free
0013 %     2. ADC_free_allcmpts
0014     
0015 nexperi = length(bvalues);
0016 Ncmpt = length(DIFF_cmpts);
0017 
0018 VOL_allcmpts = 0;
0019 
0020 for icmpt = 1:Ncmpt
0021     VOL_allcmpts  = VOL_allcmpts + VOL_cmpts(icmpt);
0022 end
0023 
0024 for icmpt = 1:Ncmpt
0025     VF_cmpts(icmpt) = VOL_cmpts(icmpt)/VOL_allcmpts;
0026 end
0027 
0028 Sig_free = zeros(size(bvalues(:)));
0029 
0030 for icmpt = 1:Ncmpt
0031     Sig_free = Sig_free+IC_cmpts(1,icmpt)*VOL_cmpts(icmpt)*exp(-DIFF_cmpts(icmpt)*bvalues(:));
0032 end
0033 
0034 ADC_free_allcmpts = sum((IC_cmpts.*VF_cmpts)'.*DIFF_cmpts')./sum((IC_cmpts.*VF_cmpts)');
0035