DemoInteractivePeakFitter.m

% Demonstration of Interactive Peak Fitter with built-in signal generator.  
% You can change the character of the signal in lines 25-31.
%  T. C. O'Haver http://www.wam.umd.edu/~toh/spectrum/InteractivePeakFitter.htm
% Version 2.2, May 2008, adds peak area to FitResults (in FitAndPlot.m)
clear
global x
global y
global xo
global dx
global NumPeaks
global Shape
global Scan
global Width
global c
global FitResults
global start
global extra

close
format compact
warning off all
warning('OFF', 'Exiting:')

% Generate simulated signal
n=2000;  % maximum x-value
increment=2; % Difference between adjacent x values
x=[1:increment:n];
amp=[1 2 3 1 4 1 3 2 1 2 1];  % Amplitudes of the peaks (one entry for each peak)
pos=[100 130 200 400 600 630 800 850 900 1200 1700];   % Positions of the peaks (one entry for each peak)
wid=[30 40 50 30 20 20 30 40 50 100 200];   % Widths of the peaks (one entry for each peak)
noise=.01;
% A = matrix containing one of the unit-amplidude functions in each of its rows
clear A
A=zeros(length(pos),length(x));
for k=1:length(pos),
  A(k,:)=gaussian(x,pos(k),wid(k));  % Can be any peak function
end
TrueY=amp*A;  % Multiplies each row by the corresponding amplitude and adds them up
% TrueY=TrueY+lorentzian(x,0,1000); % Adds background 
y = TrueY + noise .* randn(size(x));

% Arrange the peak parameters of the simulated signal into a matrix
ColumnLabels='           Peak#   Position     Height     Width';
NumPeaks=length(pos);
for m=1:NumPeaks,
  if m==1,  
     ActualParameters=[m pos(m) amp(m) wid(m)];, 
  else
     ActualParameters=[ActualParameters; [m pos(m) amp(m) wid(m)]];
  end
end

% Display the peak parameters of the simulated signal 
ColumnLabels
ActualParameters

% Execute the Interactive Peak Fitter script
InteractivePeakFitter;