%% PRIMER IZRISA ANALITIČNE FUNKCIJE v MATLABU / OCTAVE
% Želimo izrisati krivuljo toka kot funkcijo časa v obliki enačbe i(t)=I0*exp(-t/Tau)*sin(omega*t)
% kjer so I0, Tau in omega konstante, t pa čas
 
% Določitev konstant in opis enot
I0= 10 % tok v [A]
Tau=0.2  % časovna konstanta v [s]
omega=20 % kotna frekvenca v [1/s]
 
% Določitev časovnih točk, kjer bomo računali funkcijo
tzac=0 % začetni čas
tkon=1 % končni čas
Nkorakov=100 % korak
t=linspace(tzac,tkon,Nkorakov);  % generiram vektor časa v [s] 
% po potrebi glej help linspace za nastavitve
 
% Določitev funkcije toka v Matlab sintaksi. 
% Bodi pozoren na dodajanje pike "." pred simboli za množenje in deljenje vektorjev
i=I0*exp(-t/Tau).*sin(omega*t);
 
 
% Izris funkcije
plot(t,i,'Linewidth',1.5)
xlabel('Cas   /s')
ylabel('Tok   /A')
 
 
% Za izpis tabele i in t
[t;i]'  % kreiramo matriko z dvema kolonama
 
%% Dodatne možnosti izrisa, odkomentiraj in preiskusi
% plot(t,i,'-*') % s simbolom in polno črto
% plot(t,i,'g*') % z zelenimi simboli
% grid   % doda mrežo
% title('Tok v odvnisnosti od časa, prehodni pojav')
 
% break % komentiraj to vrstico, če želiš izvršiti program do konca
%% Dodatni izris dveh krivulj za različne parametre
figure;  % kreiram nov graf
 
i1=1*exp(-t/0.2).*sin(20*t);
i2=2*exp(-t/0.5).*sin(20*t);
i3=3*exp(-t/0.2).*sin(10*t);
plot(t,i1,t,i2,t,i3)
 
% Prikaz s črticami,  glej help plot za druge moznosti
plot(t,i1,'--',t,i2,'.',t,i3,'.-')
 
xlabel('Cas   /s')
ylabel('Tok   /A')
legend('i1','i2','i3')  % izris legende

%break

% 2.	Del Matlab kode za zamik funkcije v t osi

t2=t-0.2;
t3=t+0.3;
i2=1*exp(-t2/0.2).*sin(20*t2);
i3=1*exp(-t3/0.2).*sin(20*t3);
plot(t,i1,t,i2,t,i3)
legend('i1','i2','i3')  % izris legende