cpeakquantification.m

利用蚁群算法和PSO算法实现MALDI_TOF数据分析

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;