Design and simulate a set of gSlider excitation pulses¶

%matplotlib notebook
import numpy as np
import sigpy as sp
import sigpy.mri as mr
import sigpy.plot as pl
import sigpy.mri.rf as rf
import scipy.signal as signal
import matplotlib.pyplot as pyplot

N = 512
G = 5
exFlip = 90*np.pi/180
tb = 12
d1 = 0.01
d2 = 0.01
phi = np.pi
dt = 1

pulses = rf.slr.dz_gslider_rf(N, G, exFlip, phi, tb, d1, d2, cancel_alpha_phs = True)
Mxy = np.zeros((np.size(np.arange(-2*tb, 2*tb, 0.01)), G), dtype = complex)

for Gind in range(1, G+1):
    
    [a, b] = rf.sim.abrm(pulses[:, Gind-1], np.arange(-2*tb, 2*tb, 0.01), True)
    Mxy[:, Gind-1] = 2*np.multiply(np.conj(a),b)
    

pl.LinePlot(Mxy.T)

pl.LinePlot(pulses.T)