mridemo_t1w_principle.m

function mridemo_t1w_principle(do_export_gif)

if (nargin < 1), do_export_gif = 0; end

% Define timline
T_sim = 0.605;

my_timeline = timeline(T_sim);


% Define pulse sequence

rfs = {...
    rf('x', 90,  50e-3, 5e-3), ...
    rf('x', 90,  150e-3, 5e-3), ...
    rf('x', 90,  250e-3, 5e-3), ...
    rf('x', 90,  350e-3, 5e-3), ...
    };

grads = {};

acqs = {...
    acq( 160e-3, 5e-3), ...
    acq( 260e-3, 5e-3), ...
    acq( 360e-3, 5e-3), ...
    };

my_pulse_seq = pulse_sequence(rfs, grads, acqs, my_timeline);


% Setup spin system and run the simulation!
t1_list = [360 50] * 1e-3;
t2_list = [30 30] * 1e-3;

m0_list = [1 1];

n_arrow = 1;

b0_fun = @(n) zeros(n, 1);


% Setup plot functions
l_str = {'Long T1 and short T2', 'Short T1 and short T2'};

my_plot_engine = mrisim_plot_engine(l_str);

my_plot_engine.plot_timeline.do_extend_rf_plot = 1;
my_plot_engine.plot_timeline.do_plot_acq = 0;

my_plot_engine.do_export_gif = do_export_gif;


for c = 1:2

    % Set up and init spin system
    my_spin_system = spin_system(...
        m0_list(c), ...
        t1_list(c), ...
        t2_list(c), ...
        b0_fun, ...
        n_arrow);

    % Set the identity of the system to control plotting
    my_spin_system.c_system = c;

    % Setup the simulator
    my_mri_sim = mrisim(...
        my_pulse_seq, ...
        my_spin_system);

    % Run simulation
    my_mri_sim.simulate(my_plot_engine);

end