clear all

fordat2mat

for k=1:1:length(yval)
cas=(xval{k,1})';
amp=(yval{k,1})';9


figure('Name',strcat([num2str(k),'/',num2str(length(yval)),' ',num2str(signal{k,1})]),'Position', [5 180 600 400])
plot(cas,amp)

%% dialog %%%
  dlgtitle = 'Nastaveni pro vypocet';  
  prompt = {'Pocatek frekvencniho pasma pro vypocet RMS (Hz):',
  'Konec frekvencniho pasma pro vypocet RMS (Hz):',
  'Pocatecni cas vyhodnocovani dat (s):',
  'Konecny cas vyhodnocovani dat (s):'};
  dims = [1 1 1 1];
  definput = {'10','1000', min(cas), max(cas)};
  answer = inputdlg(prompt,dlgtitle,dims,definput);
 
 fminimum{k,1}=str2double(answer(1));
 fmaximum{k,1}=str2double(answer(2));
 tmin{k,1}=str2double(answer(3));
 tmax{k,1}=str2double(answer(4));
 
 close all
%%%%%%%%%%%%%%%%%%%%%%
end

wtb=waitbar(0,'Pocitam...');
for k=1:1:length(yval)


 matice=[xval{k,1}',yval{k,1}'];
 cas=matice(matice(:,1)>=tmin{k,1} & matice(:,1)<tmax{k,1},1);
 amp=matice(matice(:,1)>=tmin{k,1} & matice(:,1)<tmax{k,1},2);

ns=length(cas);
dt=(cas(end)-cas(1))/(ns-1); 
fs=1/dt;
 

[freqf,ampf]=furiertrans(cas,amp);  %% fft, vystup finalni spektrum
[ampfrange]=specrange(freqf,ampf,fminimum{k,1},fmaximum{k,1}); %% oriznuti spektra na pozadovanou hodnotu

%% PREVOD ZE ZRYCHLENI NA RYCHLOST - LAPLACE a/s=a/(2*pi*f)=v
ampf(2:end)=ampf(2:end)./(2*pi*freqf(2:end));
ampfrange(2:end)=ampfrange(2:end)./(2*pi*freqf(2:end));  %% od 2 hodnoty, protoze 1. je DC složka


mohutnost_vyrezu_ze_spektra=sqrt((sqrt(sum(ampfrange(2:end).^2))/sqrt(2))^2+(ampfrange(1)/2)^2)*1000; %% rychlost z vyrezu frekvenci
 

% fprintf('%d-%ds pro cele namerene frekvencni spektrum:\n',cas(1),cas(end)+dt);
%mohutnost_ze_signalu=sqrt(sum(amp.^2)/ns); %% zrychleni ze signalu
%mohutnost_ze_spektra=sqrt((sqrt(sum(ampf(2:end).^2))/sqrt(2))^2+(ampf(1)/2)^2)*1000; %% rychlost ze signalu


index(k,1)=k;
%RMS{1,1}(k,1)=mohutnost_vyrezu_ze_signalu;
RMS{1,2}(k,1)=mohutnost_vyrezu_ze_spektra;
%RMS{1,3}(k,1)=mohutnost_ze_signalu;
%RMS{1,4}(k,1)=mohutnost_ze_spektra;

waitbar(k/length(yval),wtb)




clear mohutnost_vyrezu_ze_signalu mohutnost_vyrezu_ze_spektra mohutnost_ze_signalu mohutnost_ze_spektra;
end
ukladani_mohutnost
close (wtb)
msgbox('Hotovo');


