Problema con Matlab[Pippo Help me]

[litie]

Ho un problema per un'esame.

Si tratta di un problema di farmacocinetica,dosi ripetute in un modello monocompartimentale.....

I parametri farmacocinetici sono questi

V(Volume di distribuzione) = 5litri
Dose=500mg
K01=1.2 ore-1

Rappresento il sistema spazio degli stati

dq/dt=-K01q(t) +u(t),
y(t)=q(t)/V;
q(0)=0;


in pratica u(t) deve essere un segnale periodico impulsivo.Le somministrazioni avvengono ad intervalli regolari TAU.

Come genero questo segnale con Matlab?

Ho provato cos?

A=[-K01];B=[500];C=[1/V];D=[0];
sys=(A,B,C,D);
[u,t]=gensig('pulse',6,72,0.01);
[y,t]=lsim(sys,u,t);

dove, 6=TAU ? l'interrvallo tra dosi,72 ? la durata della cura e 0.01(ipotizzato da me) ? il passo di quatizzazione.
Dovrebbe essere corretto ma col fischio che funziona,y(t) non ? sballata.....la concentrazione media nel plasma nel caso di dosi ripetute la si calcola cos?:

Yss=DOSE/K01*TAU*V


magari sbaglio ad utilizzare la lsim....se guardo u(t) generato con gensig mi sembra OK ,cio? ho un impulso ad intervallo regolari....per? poi y(t) calcolata con lsim ? sballata.....

help

[Pippo]

Litie, ti rispondo con MP