Problema con Matlab[Pippo Help me]
Inviato: lun mag 08, 2006 10:11 amHo 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
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