Function

[dastew at sympatico.ca]

Date: Sat, 4 Jul 2009 20:22:09 -0300
From: mjtrangoni at gmail.com
To: help-octave at octave.org
Subject: Function








Hi, Im Mario from Argentina.



I want to make a graphic of a continuous function that change itselft
at certain place. Maybe someone know how to manage this case.



It's a current's graphic of a chage/discharge of a resistant-inductance
serial circuit.



The script file that I was making is that: 





##Expresion de i_a en el tiempo

# Define eje x - de 0 a 5 

t = linspace(0, 2, 500);

tiempo = 2;

figure, clf

  

#Funcion a graficar

#t(1:length(t)/2)

x1 = (24.-(24.*(e.^(-1.*(t(1:length(t)./2))))));

  

#t(length(t)/2:end)

#x2 =
(24.-(24.*(e.^(-1.*(t(250)))))).*(e.^(-1.*((t(length(t)./2+1:end)))));

  

#x2 = (a.*(e.^(-1.*((t(length(t)/2+1:end).-(t(1:length(t)/2)))))));

x2 = ((24.-(24.*(e.^(-1.)))).*(e.^(-1.*((t(length(t)/2+1:end))))));

  

  

#x2 =
((24.-(24.*(e.^(-1.)))).*(e.^(-1.*((t(length(t)/2+1:end).-(t(1:length(t)/2)))))));

#Ploteo

  

plot(t, [x1 x2], "b", "linewidth", 0.5)

  

size(t)

size(x1)

size(x2)

  

#plot(t, y(t), "b", "linewidth", 0.5);

#plot(t, ((24.-1)./1)*(1.-e.^(-t)), "b", "linewidth", 0.5);

  

#Grilla

grid on, hold on

  

#Titulo

title('Corriente de Armadura i_a', 'FontSize',20);

  

#Ejes

ylabel('i_a (A)','FontSize',16)

xlabel('t (seg)','FontSize',16)

  

## Imprime solido en eps

print "-S1680,1050" -solid -deps ia.eps;

  

  
Looking forward for an answer





Mario Trangoni

UTN FRC

Argentina






Is this what you were thinking?

t1=0:.2:10;
t2=0:.2:10;
t3=0:.2:20.2;
x1=24*(1-e.^(-t1));
x2=24*e.^(-t2);

x=[x1 x2];

plot(t3,x)



Doug


-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://www-old.cae.wisc.edu/pipermail/help-octave/attachments/20090705/e68159d4/attachment-0001.html