www.pudn.com > maldicode.zip > cPeakQuantification.m

 
clear all; clc 
load Windows 
load casecirr_3030_prepro 
 
m  = size(sample_tr,2); 
 
% TAKE THE MAXIMUM OF THE INTENSITIES IN THE IDENTIFIED WINDOW 
 
sample = [sample_tr]; 
max_win = []; 
for i = 1: length(win_tick) 
    if win_tick(i,1) == win_tick(i,2) 
        a=sample(win_tick(i,1),:); 
    else 
        [a in] = max(sample(win_tick(i,1):win_tick(i,2),: )); % sum/max 
        x1 = min(in); x2 = max(in); 
        wint(i,1) = win_tick(i,1)+x1-1;   
        wint(i,2) = win_tick(i,1)+x2-1; 
    end 
    max_win = [max_win;  a]; 
end 
wind = MZ(wint); 
sample_tr = max_win; 
 
sample = [sample_pr]; 
max_win2 = []; 
for i = 1: length(wint) 
    if wint(i,1) == wint(i,2) 
        a2=sample(wint(i,1),:); 
    else 
        a2 = max(sample(wint(i,1):wint(i,2),: )); % sum/max 
    end 
    max_win2 = [max_win2;  a2]; 
end 
 
sample_pr = max_win2; 
 
windows = wind; 
win_tick = wint; 
save casecirr_249wind sample_tr sample_pr label_tr label_pr windows win_tick;