%---------------------------DEFINITIVNI VERZE
%skript vygenerovani prubehu tuhosti, pruznosti a hmotnosti pro mnoho stejnych poli
%pomoci funkci bpl_f (bez pridavnych lan) a spl_f (s pridavnymi lany)
%vystupnimi parametry skriptu jsou vektory:

%mc (prubeh hmotnosti vedeni kg),
%ec (prubeh pruznosti vedeni m/N), 
%kc (prubeh tuhosti vedeni N/m)
%lc (vektor vzdalenosti pro predchozi vektory m)



%---------------------------------------------------------vypocet mc, ec, kc a lc
n=length(x)+1;							%pocet elementu (vesaky+zavesy)
x1=[x,[0,0,0,0,0,0,0,0,0]];         	%rozsireni vektoru x doprava
for a=1:5;								%rozkopirovani vektoru do peti radku
   x1(a,:)=x1(1,:);						%docasny vektor odlehlosti
   n1(a)=n(1);                          %docasny vektor poctu prvku
   xpl1(a,:)=xpl;                       %docasny vektor pridavaku
end

l=0;
for a=1:n(1)							%vypocet delky rozpeti z vesakovych odlehlosti
   l=x1(1,a)+l;
end

np=l/delta;						%pocet prvku vystupnich matic na jedno rozpeti
%np=single(np);
%np=double(np);
if np~=round(np);
    np=round(np);
end

if xpl==[0,0]					%volame vypocet parametru docasnymi vektory pro pet stejnych poli
   [mc, ec, lc, a]=bpl_f(n1,x1,xpl1,mtd,mnl,mpl,Ftd,Fnl,Fpl,delta,5); %bez pridavaku - v matici pridavnych lan same nuly
else
   [mc, ec, lc, a]=spl_f(n1,x1,xpl1,mtd,mnl,mpl,Ftd,Fnl,Fpl,delta,5); %s pridavakem
end

clear a;
clear x1;                       %vymazani docasnych vektoru
clear n1;
clear xpl1;



%---------------------------------vyber prostredniho rozpeti, ktere se rozkopiruje
a=1;
while lc(a)<2*l						    %zjistime na kterem miste ziskanych vektoru
   a=a+1;								%zacina prostredni (treti) rozpeti
end

mc=mc(a:np+a-1);							%ze vseho vezmeme jenom prostredni rozpeti
ec=ec(a:np+a-1);
lc=[delta/2:delta:l];


%------------------------------------------------rozkopirovani na pozadovany pocet rozpeti
for a=0:pocet_poli-1	
    
       
   mc(a*np+1:(a+1)*np)=mc(1:np);			%vaha v kg
   ec(a*np+1:(a+1)*np)=ec(1:np);			%pruznost v m/N
   lc(a*np+1:(a+1)*np)=lc(1:np)+a*l;        %odlehlost pro predchozi vektory
end

%lc=lc./1000;                             %prevod odlehlosti na metry
kc=1./ec;                          %prevod na tuhost

clear a;