Esercizio Matlab
In quanto non conosco nulla di Matlab ma necessito quest'esercizio per la tesi , vi chiedo gentilmente di svolgere questo esercizio,ovvero bisogna cambiare i dati dei grafici con i seguenti:
1° grafico invariato.
2° grafico con assunzione del farmaco random di giorno in giorno (es.1° giorno 3 dosi,2° giorno 1 dose,3° giorno 5 dosi(max 7 dosi in un giorno))
3° grafico con dose crescente (1° giorno 1 dose,2° 2 dosi ecc..)
4° grafico con dose decrescente (1°giorno 7 dosi...)
Grazie mille in anticipo,qui l'esercizio con dati diversi:
1° grafico invariato.
2° grafico con assunzione del farmaco random di giorno in giorno (es.1° giorno 3 dosi,2° giorno 1 dose,3° giorno 5 dosi(max 7 dosi in un giorno))
3° grafico con dose crescente (1° giorno 1 dose,2° 2 dosi ecc..)
4° grafico con dose decrescente (1°giorno 7 dosi...)
Grazie mille in anticipo,qui l'esercizio con dati diversi:
sis=tf([.7 1],conv([.1 1],conv([.5 1],[36 1]))); t=0:.01:24*7; u=t*0; T1=1:2400:24*4.5*100; u(T1)=1; y=lsim(sis,u,t); t=t/24; T1=T1/24/100; close subplot(2,1,1) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con una dose al giorno') xlabel('giorni') axis([0,max(t),0,1e-3]) t=0:.01:24*7; u=t*0; T1=1:1200:24*4.5*100; u(T1)=1; y=lsim(sis,u,t); t=t/24; T1=T1/24/100; subplot(2,1,2) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con due dosi al giorno') xlabel('giorni') axis([0,max(t),0,1e-3]) pause t=0:.01:24*7; u=t*0; T1=1:4800:24*4.5*100; u(T1)=2; y=lsim(sis,u,t); t=t/24; T1=T1/24/100; close subplot(2,1,2) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con una dose doppia al giorno') xlabel('giorni') axis([0,max(t),0,1.2e-3]) t=0:.01:24*7; u=t*0; T1=1:1200:24*4.5*100; u(T1)=1; y=lsim(sis,u,t); t=t/24; T1=T1/24/100; subplot(2,1,1) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con due dosi al giorno') xlabel('giorni') axis([0,max(t),0,1.2e-3]) pause
Risposte
Finito,funziona tutto.. Grazie infinite!!
Risolto.
sis=tf([.7 1],conv([.1 1],conv([.5 1],[36 1]))); close t=0:.01:24*7; u=t*0; T1=1:2400:24*7*100; u(T1) =1; y=lsim(sis,u,t); t=t/24; T1=T1/24/100; subplot(2,1,1) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con una dose al giorno') xlabel('giorni') axis([0,max(t),0,4e-3]) t=0:.01:24*7; u=t*0; T1=1:2400:24*7*100; u(T1)=randi(7,1,7); y=lsim(sis,u,t); t=t/24; T1=T1/24/100; subplot(2,1,2) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con n dose al giorno') xlabel('giorni') axis([0,max(t),0,4e-3]); t=0:.01:24*7; u=t*0; T1=1:2400:24*7*100; u(T1)=[1,2,3,4,5,6,7]; y=lsim(sis,u,t); t=t/24; T1=T1/24/100; figure() subplot(2,1,1) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con dose crescente') xlabel('giorni') axis([0,max(t),0,4e-3]) t=0:.01:24*7; u=t*0; T1=1:2400:24*7*100; u(T1)=[7,6,5,4,3,2,1]; y=lsim(sis,u,t); t=t/24; T1=T1/24/100; subplot(2,1,2) plot(t,y,T1,T1*0,'.') grid on title('Concentrazione plasmatica con dose decrescente') xlabel('giorni') axis([0,max(t),0,4e-3]);
Risolto.
Ottimo!