Stimatore del tempo di attesa per finire un download
Appendice
(calcoli octave)


obiezione respinta logo
Segnalazioni rivolte al mondo femminile ma non solo, delle farmacie, degli ospedali e dei consultori che non rispettano la legge e che non offrono la pillola anticoncezionale, che non fanno IVG, che non danno la RU486, che trattano malissimo la paziente, che non danno la pillola dei 5 giorni dopo magari con la scusa che l'hanno finita. Ma anche le segnalazioni positive di buona sanità, di competenza e attenzione.
Mappa delle segnalazioni.



Firefox: Riprenditi il web




Statistiche Europee d'uso dei browser 17/10/2017 :
Chrome...56%
Firefox..18%
IE.......8%
Safari...7%

Fonte: StatCounter

Borsa valori della moneta Bitcoin (18/10/2017):

prezzo di 1 bitcoin:
4785.69 €
LTC
DRK


Fonte: Bitcoincharts.com

Sostieni Wikileaks!
Se hai qualche bitcoin fai una donazione a wikileaks all'indirizzo:
1HB5XMLmzFVj8ALj6 mfBsbifRoD4miY36v

Queste sono le donazioni fatte fin'ora: 4025 bitcoin.




L'intero sito kensan.it ha avuto questi Share:

facebook: 6684
twitter: -
google+: 0
linkedIn: 0




IL TUO 5 PER MILLE PER GLI OSPEDALI DI EMERGENCY codice fiscale:
971 471 101 55







Calcoli Octave Eseguiti per la Parte I:

N=30; n=100; l=0.1; x=poisson_rnd(1/l,1,N);
Si tratta di avere N campioni di poisson e di stimare le rimanenti n-N

ls=N/(sum(x)); tN=sum(x); tn=(n-N)/ls+tN; var_tn=sqrt((n-N)/(ls^2));
l è la densità degli arrivi, ls la sua stima, tn la media dell'arrivo finale, seguita dalla varianza normalizzata.

t_fin=3000; t=[0:t_fin];
fn=prod((ones(n-N,1)*ls))/prod(1:n-N-1)*prod(ones(n-N-1,1)*t).*exp(-ls*t);
plot(t,fn)
Plottaggio della densità di probabilità di tn

x0=x-1/ls; var_xN=sqrt(1/(N-1)*sum(x0.*x0)-1/(n*(n-1))*(sum(x0))^2)
Varianza fatta con i campioni

Calcoli Octave Eseguiti per la Parte II:


d_t1=0.00001;t1=[d_t1:d_t1:d_t1*t_fin];
Vettore tempo per la v.a. ß

fsn=prod((ones(N,1)*ls))/prod(1:N-1)./prod(ones(N+1,1)*t1).*exp(-ls./t1);
plot(t1,fsn);
Densità di ß e plottaggio

printf("Integrale di fsn(t1): %e errore rispetto all'1: %e\n",sum(fsn*d_t1),1-sum(fsn*d_t1));
Verifica di correttezza sulla densità

semilogy(t1(90:3000),fsn(90:3000));
disegno di una parte della densità di ß

med_landa=sum(t1.*fsn*d_t1)*N; pot_landa=sum(t1.^2.*fsn*d_t1)*N^2;var_landa=pot_landa-med_landa^2; var_norm_landa=sqrt(var_landa); printf("valori di landa medio: %f, potenza: %f, varianza: %f e varianza normalmizzata: %f\n",med_landa,pot_landa,var_landa,var_norm_landa);
Calcolo dei momenti di ß

vet_ls=[]; for i=[1:N]; vet_ls(i)=i/sum(x(1:i)); end;
plot([1:N],vet_ls,'xg');hold;plot([1:N],(ls-var_norm_landa)*ones(N,1));plot([1:N],(ls+var_norm_landa)*ones(N,1)); hold; printf("landa vera: %f, landa stimata: %f, varianza normalizzata di landa: %f\n",l,ls,var_norm_landa);
Plottaggio del grafico mostrato nella Parte II (Figura 1)
Ciao Anonimo, commenta questo articolo!

NON ci sono ancora COMMENTI per questo articolo :

commenti abilitati per gli anonimi
Commenti sperimentali
by kensan & Mp
kensan logo Licenza Creative Commons 3.0
I miei testi sono sotto la Licenza "Creative Commons 3.0 Italia": se sei interessato a pubblicare i miei articoli leggi le note aggiuntive (Licenza di kensan.it) dove troverai anche le attribuzioni dei diritti per tutte le immagini pubblicate.
Questo sito memorizza sul tuo pc uno o più cookie di tipo tecnico, leggi l'informativa estesa.
Kensan geek site

e-mail
e-mail cifrata
Sandro kensan

kensan.it Geek&Hacker site Owner, 1 ott 2003