PmWikiFisica | NMR / OnlineAnalysis browse

< Quick introduction to gtknmr | Index | Echo with two equal pulses >

start a new terminal or tab in the terminal

cd to the data directory

invoke matlab

  matlab7 (go back to matlab6, as of 10/2011)

in the matlab command window issue the command

  setup_nmr (not needed any more, as of 10/2011)

to initialize the nmr analysis routines and then

  parmaecho filename.pna -T -F

to fast Fourier transform (FFT) and produce amplitudes from FFT vs. zone number.

A GUI appears with signals and buttons

+Z -Z moves forward and backward through zones,

the cursors select:

zero before this cursor put outside the window to the left
center of filter window, when full,
point where baseline is evaluated for suntraction, place to the right, outside the window

apod selects the filter width and shape, (see below). The signal will be multiplied by the selected filter window, black line, prior to FFT,

done calculates FFT

The button apod opens the filter GUI.

Select:

shape: lorenzian, gaussian, hypergaussian, rectangular;
mode: half (only positive times) or full,
width in Hz (narrower width implies longer decay time of the black envelope)
done to finish

When done in the time selections the FFT amplitude (square root of square real plus square imaginary parts of FFT) appears in a new GUI

zoom in

cursors set the interval where the FFT amplitude integral is calculated and stored in the output amplitude variable (NMRamp), one value for each zone

center of interval
one extreme of the interval

done performs the calculation, saves analysis parameters and exits

The results are stored in a number of matlab variables, see them with

  whos

Among these NMRamp contains the calculated integrals in an array, one element for each zone

  parmaecho filename.pna -p -s7 -ATXAT

repeats the last analysis with stored parameters, only on zones from 7 to the last, outputs NMRamp and one acquisition parameter value, e.g. TxAt, for each analysed zone in the array NMRmisc

Therefore the plot of amplitudes e.g. vs. transmitter attenuation is obtained by

  plot(NMRmisc,NMRamp,'o')

< Quick introduction to gtknmr | Index | Echo with two equal pulses >