MSF ACE

deri

% Raspunsul unui integrator y(t) = b*udot(t)

% la semnalul de intare u(t) definit prin u(t1)=u1, ...., u(tn)=un.

% Intrari: vect_der=[t1 .... tn]

% vecu_der=[u1 .... un]

% b

% Iesiri: vecy_der=[y1 y2 ... yn-1]

vect_der=input('Vectorul momentelor de timp: ');

vecu_der=input('Vectorul valorilor intrarii: ');

b=input('Constanta derivatorului, b=');

vecy_der=(b*diff(vecu_der))./diff(vect_der);

subplot(2,1,1)

plot(vect_der(1:length(vect_der)),vecu_der(1:length(vecu_der)))

xlabel('timp');

ylabel('intrare');

subplot(2,1,2), plot(vect_der(1:length(vect_der)-1),vecy_der)

xlabel('timp');

ylabel('iesire');

hpri_rez

function hdot=hpri_rez(twork,hwork)

% Apelat de "ode23" in programul "rezervor".

global t q aria nrint qout;

for j=1:nrint

if (t(j)<=twork(1)) & (twork(length(twork))<=t(j+1))

debin=q(j);

end

end

hdot=(debin-qout)/aria;

inte_cst

% Raspunsul unui integrator a*ydot(t)=u(t), y(0)=y0,

% la o succesiune de trepte de amplitudini u1, u2,...,un

% aplicate pe intervalele de timp [0,t1), [t1,t2),...,[tn-1,tn]

% Apeleaza "ode23" cu functia "ypri_cst"

% Intrari: a, y0,

% vect_int=[t1 t2 .... tn]

% vecu_int=[u1 u2 .... un]

% Iesiri: timp, udetimp, ydetimp

clear

clc

vect_int=input('Vectorul momentelor de timp: ');

vecu_int=input('Vectorul valorilor intrarii: ');

yini=input('Conditia initiala y0=');

a=input('Constanta integratorului, a=');

global t u alfa nrint

t=vect_int;

%u=vecu_int;

alfa=a;

nrint=length(vecu_int);

%t=[0 t];

%y(1)=y0;

timp=[];

udetimp=[];

ydetimp=[];