3pulse.m

%%  This simulation depicts a 3-pulse sequence with a typical array of spins

clear all

%% Creates a matrix of spins along the z axis

for i=1:250
    Spins.degrees(i,:) = [0,0,1];
end

Spins.offset = linspace(-16000000,16000000,250);
Spins.Tone = 0.0002;
Spins.Ttwo = 0.000002;


%% Sets the parameters and function for first pi/2 pulse

Sim.direction = 'y';
Sim.angle = 90;
Sim.framerate = 20;
Sim.nframes = 30;

Spins.degrees = pulse(Spins,Sim);

%% Sets the evolution parameters and function between pulses

Sim.length = 0.0000002;
Sim.nframes = 50;

Spins.degrees = evolution(Spins,Sim);

%% Sets the parameters and function for second pi/2 pulse

Sim.direction = 'y';
Sim.angle = 90;
Sim.nframes = 30;

Spins.degrees = pulse(Spins,Sim);

%% Sets the evolution parameters and function between pulses

Sim.length = 0.000001;
Sim.nframes = 100;

Spins.degrees = evolution(Spins,Sim);

%% Sets the parameters and function for first pi/2 pulse

Sim.direction = 'y';
Sim.angle = 90;
Sim.nframes = 30;

Spins.degrees = pulse(Spins,Sim);

%% Sets the evolution parameters and function after pulses

Sim.length = 0.000003;
Sim.nframes = 200;

Spins.degrees = evolution(Spins,Sim);