PlotNormalizedProfiles.m

function PlotNormalizedProfiles(arrays, handle, method)
% PlotNormalizedProfiles plots a cell array of [n x m] matrices to the 
% provided axes handle. For each cell array, the elements (n) within the 
% first and last columns (m) are plotted, normalized based on method. If 
% n > 2, the first two columns will be multiplied together as the 
% horizontal axis values.
% 
% Method can be one of the following strings:
%   'none' for no normalization
%   'first' to normalize to the first value
%   'max' to normalize to the max value
%   'mean' to normalize to the average value
%
% Example usage:
%  profiles = {[1:10;rand(1,10)]', [1:10;rand(1,10)]'};
%  PlotNormalizedProfiles(profiles, figure, 'max');
%
% Author: Mark Geurts, mark.w.geurts@gmail.com
% Copyright (C) 2022 Aspirus, Inc.
%
% This program is free software: you can redistribute it and/or modify it 
% under the terms of the GNU General Public License as published by the  
% Free Software Foundation, either version 3 of the License, or (at your 
% option) any later version.
%
% This program is distributed in the hope that it will be useful, but 
% WITHOUT ANY WARRANTY; without even the implied warranty of 
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General 
% Public License for more details.
% 
% You should have received a copy of the GNU General Public License along 
% with this program. If not, see http://www.gnu.org/licenses/.

% Normalize values based on method input
s = 0;
for i = 1:length(arrays)
    if ~isempty(arrays{i})
        if s == 0 
            s = i;
        end
        % Available options are 'first', 'max', and 'mean'
        switch method
            case 'first'
                arrays{i}(:,end) = arrays{i}(:,end)/arrays{i}(1,end);
            case 'max'
                arrays{i}(:,end) = arrays{i}(:,end)/max(arrays{i}(:,end));
            case 'mean'
                arrays{i}(:,end) = arrays{i}(:,end)/mean(arrays{i}(:,end));
        end
    end
end

% Plot first non-empty array
axes(handle);
if size(arrays{s},2) > 2
    plot(handle, sqrt(arrays{s}(:,1).^2 + arrays{s}(:,2).^2), ...
        arrays{s}(:,end));
else
    plot(handle, arrays{s}(:,1), arrays{s}(:,end));
end
hold(handle, 'on');

% Loop through remaining arrays, plotting non-empty ones
for i = 2:length(arrays)
    if ~isempty(arrays{i})
        if size(arrays{i},2) > 2
            plot(handle, sqrt(arrays{i}(:,1).^2 + ...
                arrays{i}(:,2).^2), arrays{i}(:,end));
        else
            plot(handle, arrays{i}(:,1), arrays{i}(:,end));
        end
    else
        plot(handle, NaN,NaN);
    end
end

% Finish plot
hold(handle, 'off');

% Clean up 
clear s i;