function p= CWTobradaMNF(ctrlPrikaza,signal,kut,X,Fs,kutLMA,L,R,LMA,frekvencije,waveime)
% %% 
% 
% %%
% % uèitavanje signala i njihovih parametara
% 


%% Osnovna obrada

signal = signal(:)'; % row vector
len = length(signal);

MNF= [];
tMNF= [];

for i=1:R:len-L
    %x = signal(i:i+L-1); %.* window; ########## (i:i+L-1) ili (i+i+L-1)/2
    %X = cwt(x, skale, waveime);
    X_psd =X(:,i:i+L-1).*conj(X(:,i:i+L-1)); % skalogram 
    MNF{end+1}= meanfreqCWT(X_psd,frekvencije);
    tMNF{end+1}= (i+i+L-1)/2; %vremenski trenuci u kojima je izraèunata frekvencija
    
end


%% Izraèunavanje REGRESIJSKOG pravca

% izglaðivnaje izraèunatih frekvencija sa dvostrukim filtriranjem
smMNF= smoothMF(cell2mat(MNF),LMA);

% odreðivanje iznosa i vremenskih trenutaka lokalnih maksimuma
[lokMAX, tMAX]= lokMaximum(smMNF,cell2mat(tMNF),kutLMA,kut);

% linearna interpolacija u smislu najmanje kvadratne pogreške -LMSE
% pREG je regresijski pravac izraèunat za vremensku skalu u [ms]
% za skalu u [sekundama] treba množiti p(1) sa 1000
% pREG= p(1)*x + p(2)
p= polyfit(tMAX,lokMAX,1);
pREG= p(1)*[1:length(signal)]+ p(2);


%% PRIKAZ dobivenih velièina

if ctrlPrikaza==1
    
    figure,
    subplot(4,1,1); plot([1:length(signal)]/1000,signal/1000);
        xlabel('vrijeme [s]'); ylabel('[mV]'); 
        title(['EMG signal: L= ',num2str(L),', R= ',num2str(R)]);
    subplot(4,1,2); plot(cell2mat(tMNF)/1000,cell2mat(MNF)); 
        xlabel('vrijeme [s]'); ylabel('[Hz]'); 
        title('Srednja frekvencija - MNF');
    subplot(4,1,3); plot(cell2mat(tMNF)/1000,smMNF); hold on;
        plot(tMAX/1000,lokMAX,'rx','MarkerSize',10); 
        xlabel('vrijeme [s]'); ylabel('[Hz]'); 
        title('Izglaðena srednja frekvencija - MNF'); 
        hold off;
    subplot(4,1,4); plot([1:length(signal)]/1000,pREG,'k'); hold on; 
        plot(tMAX/1000,lokMAX,'rx','MarkerSize',10); 
        xlabel('vrijeme [s]'); ylabel('[Hz]'); 
        title('Regresijski pravac i lokalni maksimumi');
        str1= strcat('y= ',num2str(60*1000*p(1)),'[Hz/min]*x + ',num2str(p(2)));
        str2= 'lok. maksimumi';
        legend(str1,str2);
        hold off;

end


