function [k_space, t] = reconstructionNMR (pd, t1, t2, tr, te, excit)
% =========================================================================
% [k_space, t] = resonstructionNMR (pd, t1, t2, tr, te, excit)
% 
% funkce rekonstruujici obrazek NMR, jejimi vystupem bude virtualni fantom
% (matice rho, T_1 a T_2) a matice K (k-prostor), kde radky budou odpovidat
% jednotlivym excitacim 
% 
% PARAMETRY:
% * pd - Hustota spinu
% * t1 - Mr?zkova relaxacn? konstanta
% * t2 - Spinova relaxacn? konstanta
% * tr - Opakovac? perioda
% * te - Cas do mereni
% * excit - pocet excitaci
% 
% VYSTUP:
% * k_space - matice K, kde radky budou odpovidat jednotlivym excitacim
% * t - cas (odpovidajici radkum v matici K)
% 
% PRIKLAD:
%     load('fantom.mat');
%     tr = 5000; te = 1;
%     [k_space, t] = reconstructionNMR (pd, t1, t2, tr, te, 4);
%     figure; imagesc( log(abs(k_space)) ); colormap gray;
% 
% -------------------------------------------------------------------------

if nargin<5,  error('too few parameters');  end;
