Лаба 10 / part2

function part2()
global hMenu1 hA2 hA3 hA4 RR3 RR
 
r = get(hMenu1,'Value');
if r == 1
    axes(hA2)
    cla
    axes(hA3)
    cla    
    axes(hA4)
    cla
else
    RR = RR3(:, r-1);
    axes(hA2)
    cla
    hold on
    for i = 1:450
        x(1) = i;
        x(2) = i;
        R(1) = 0;
        R(2) = RR(i);
        plot(x, R);
        hold on
    end
    maxRR = max(RR)*1.2;
    set(hA2, 'Ylim',[0 maxRR])
    
    axes(hA3)
    dH = 0.05;
    X = 0:dH:maxRR;
    H = histc(RR,X)
    SH = sum(H);
    PH= H/SH*100;
    bar(X,PH, 'histc')
    maxh = max(PH)*1.2;
    set(hA3, 'Ylim',[0 maxh])
    
    axes(hA4)
    NRR = length(RR);
    plot( RR(1:NRR - 1), RR(2:NRR), 'o');
    maxr = max(RR)*1.2;
    set(hA4, 'Ylim',[0 maxRR],'Xlim',[0 maxRR])
    
    RRsr = mean(RR);
    RRmin = min(RR);
    RRmax = max(RR);
    dRR = RRmax - RRmin;
    [AMo, iMo] = max(PH);
    Mo = iMo*dH;
    axes(hA3)
    xt = maxRR/20;
    yt = maxh;
    dyt = maxh/15;
    text (xt, yt - dyt, ['RRsr = ' num2str(RRsr) ' c'])
    text (xt, yt - 2*dyt, ['RRmin = ' num2str(RRmin) ' c'])
    text (xt, yt - 3*dyt, ['RRmax = ' num2str(RRmax) ' c'])
    text (xt, yt - 4*dyt, ['dRR = ' num2str(dRR) ' c'])
    text (xt, yt - 5*dyt, ['Mo = ' num2str(Mo) ' c']) 
end
return
end