\documentclass[11pt]{article}
\usepackage{palatino,graphicx}
%\pagestyle{empty}
\pagestyle{headings}
\input hdr_jmp

\topmargin=-0.5in
\textheight=9in
\evensidemargin 0.0in
\oddsidemargin 0.0in
\setlength{\textwidth}{6.5in}

\begin{document}
\ 
\hdr{1}

\begin{center}
{\Large\em EE469B:  Assignment 5} \\[0.15in]
{\em Due Thursday Oct. 31}
\end{center}


\paragraph{1. Spinor Magnitude Constraint}  Starting from the expressions for $\alpha$ and $\beta$ as a function of a rotation angle $\theta$ about a unit vector ${\bf n} = (n_x, n_y, n_z)$, show that
\[
\alpha\alpha^* + \beta\beta^* = 1
\]

\paragraph{2. Composite Pulses} Often combinations of pulses are used to perform specific tasks more accurately.  One example is the sequence consisting  of a rotation of $\theta$ about the $x$ axis, followed by a rotation by the same $\theta$ about the $y$ axis.  Our goal in this case is to start from equilbirum ${\bf M} = (0,0,M_0)$ and to accurately make $M_z$ zero.  The magnetization is left somewhere in the transverse plane, but we don't care where.

\subparagraph{a)} Solve for $\alpha$ and $\beta$ for this pulse sequence.

\subparagraph{b)} Find a simple expression for the longitudinal magnetization. Solve for the $\theta$ that produces an $M_z=0$. Plot the $M_z$ as a function of $\theta$ for $0<\theta<\pi$.  
What range of $\theta$ will result in $M_z$ being less than $0.05M_0$?

\subparagraph{c)} Find a simple expression for the transverse magnetization. Plot $ M_{xy}$ as a function of $\theta$ for  $0<\theta<\pi$. 

\subparagraph{d)} 
Plot $\angle M_{xy}$ as a function of $\theta$ for  $0<\theta<\pi$.  Note that this function is approximately linear about $\theta = \pi/2$.  Hence, what this composite pulse does is trade off accuracy in $M_z$ for phase error in $M_{xy}$. One use for this pulse sequence is the measurement of the RF field strength, since the actual flip angle $\theta_a$ is approximately the negative of magnetization phase $-\angle M_{xy}$.  How far can $\theta_a$ vary from $\pi/2$ before the error in this approximation is 5$^\circ$?.

\vspace{0.125in}
\noindent
{\em Hint: } For these problems the half and double angle formulas are useful
\begin{eqnarray*}
  \sin \theta & = & 2 \sin (\theta/2) \cos (\theta/2) \\
  \cos \theta & = &   \cos^2(\theta/2) - \sin^2(\theta/2)
\end{eqnarray*}
Your answers shouldn't have any half angles.

\paragraph{3. Inverse SLR Transform}  On the web site are several m-files that implement the inverse SLR transform.  These are
\begin{verbatim}
b2a.m, mag2mp.m, ab2rf.m
\end{verbatim}
The first is \verb+b2a.m+. This takes a beta polynomial and returns the consistent, minimum phase, alpha polynomial.  It uses \verb+mag2mp.m+ to go from the magnitude profile of alpha to the minimum phase alpha.  Then \verb+ab2rf.m+ takes the alpha and beta polynomials, and returns the corresponding RF pulse.  
There is also an m-file, 
\begin{verbatim}
ab2inv.m
\end{verbatim}
that takes the output of the simulator \verb+abrm.m+, and returns the $M_z$ profile that is produced.

Given a suitably scaled beta polynomial (i.e. the passband amplitude is $\sin(\theta/2)$), then the RF pulse can be computed as
\begin{verbatim}
>> a = b2a(b);
>> RF = ab2rf(a,b);
\end{verbatim}
One way to think about the inverse SLR transform is that it corrects for the nonlinearity of the Bloch Equation.  We can design the beta polynomial using Fourier arguments, just as a small tip angle pulse.  Then scale the beta polynomial to $\sin(\theta/2)$, and apply the inverse SLR transform to produce the RF pulse with that beta profile.

\subparagraph{a)}  Design a TBW=8 windowed sinc RF pulse, and scale it to a $\pi$ radian rotation.  Plot the RF pulse scaled to Gauss.  Assume a pulse length of 8 ms, and a gradient strength of 0.425 G/cm.  Simulate and plot the inversion profile.  Choose an interesting range of spatial locations, and compute the inversion profile with
\begin{verbatim}
>> mz = ab2inv(abrm(rf,x));
\end{verbatim}

\subparagraph{b)}  Now, use the windowed sinc waveform as the beta polynomial.  First, scale it to the proper value for an inversion.  Then find the corresponding minimum phase (and minimum power) alpha polynomial using \verb+b2a.m+.  Next compute the RF pulse using \verb+ab2rf.m+.  Plot the SLR inversion pulse and the windowed sinc inversion pulse from part (a), both scaled to Gauss.  By what factor has the peak amplitude increased?

\subparagraph{c)} Simulate the inversion profile of the SLR inversion pulse from part (b), and plot it along with the inversion profile of the windowed sinc pulse of part (a).  

\end{document}