example.m

% The following code illustrates some use cases of MrVox2d

addpath(genpath(fullfile(pwd, 'mrvox')));

%% FID example
fprintf('\n\t>> Simulating a FID on simple voxel...\n\n');
[sa, sphi] = VoxelSim2D_do_one('config/voxpar_simple.txt', 'config/seqpar_fid.txt');
fprintf('\n\n[Press any key to continue]\n\n')
pause
close all

%% Spin echo example
fprintf('\n\t>> Simulating a spin echo on simple voxel...\n\n');
[sa, sphi] = VoxelSim2D_do_one('config/voxpar_simple.txt', 'config/seqpar_se.txt');
fprintf('\n\n[Press any key to continue]\n\n')
pause
close all

%% Spin echo example with strong susceptibility and reduced diffusion
fprintf('\n\t>> Simulating a spin echo on voxel with higher susceptibility and lower diffusion...\n\n');
[sa, sphi] = VoxelSim2D_do_one('config/voxpar_simple_highsusc_lowdiff.txt', 'config/seqpar_se.txt');
fprintf('\n\n[Press any key to continue]\n\n')
pause
close all

%% Stimulated echo example
fprintf('\n\t>> Simulating a spin echo on simple voxel...\n\n');
[sa, sphi] = VoxelSim2D_do_one('config/voxpar_simple_highsusc_lowdiff.txt', 'config/seqpar_stim.txt');
fprintf('\n\n[Press any key to continue]\n\n')
pause
close all

%% Narrow diffusion gradient example with hindered diffusion and cells
fprintf('\n\t>> Simulating a narrow diffusion gradient on simple voxel...\n\n');
[sa, sphi] = VoxelSim2D_do_one('config/voxpar_simple_hindered_cell.txt', 'config/seqpar_narrowpulse_singlerefoc.txt');
fprintf('\n\n[Press any key to continue]\n\n')
pause
close all

%% Dictionary generation example
fprintf('\n\t>> Simulating a dictionary on a simplistic voxel with varying vessel radius...\n\n');
dico = GenLookUp2D('config/voxpar_dico.txt', 'config/seqpar_GESFIDE.txt');

% Just to get input parameters for plotting legend
Model   = ReadModel('config/voxpar_dico.txt');
figure,plot(abs(dico')), xlabel('Time (ms)'), ylabel('MR signal magnitude')
title('Dictionary')
legend(cellstr(num2str(Model.geo.vasc.R'*1e6, 'R=%-dum')))