Esercizio distribuzione di Poisson

MrBlonde93
Buongiorno a tutti,
scrivo il testo di un esercizio di statistica su una poissoiniana di cui non riesco a venire a capo. Ecco il testo:
Per portare a termine una lavorazione di un pezzo è necessario utilizzare due macchine M1 e M2 in successione. La macchina Mi guasta mediamente Li volte in un anno e in presenza del guasto di una delle macchine la lavorazione si arresta.
1- Calcolare la probabilità che la lavorazione dei pezzi non si interrompa per almeno un tempo P0.
2- Tempo medio della prima interruzione
3- Supponendo che l2 sia minore di l1 e che si decida di affiancare a M2 una macchina N2 con le stesse caratteristiche, da far intervenire in caso di guasto. Calcolare la probabilità che la lavorazione dei pezzi non si interrompa per almeno tu tempo P0

Risposte
Lo_zio_Tom
L'esercizio in questione è [strike]piuttosto standard[/strike] interessante [strike]e non[/strike] però è possibile che tu non riesca a fare nulla.....posta ciò che hai fatto e vediamo di proseguire....

MrBlonde93
Messaggio modificato e testo corretto.
1- Applico banalmente Poisson con $Mi$ e li $P(a=P0) = ...$
2- Quella media sarà data dalla media di $( 1 - P(a=P0) )$
3- Le probabilità saranno vincolate tra loro, quindi il fatto che N2 entri in azione sarà legato al guasto di M2 (Bayes?), ma non ho altre idee in merito.

Lo_zio_Tom
1. non avendo inserico formule non ho capito come intendi risolverlo ma, dato che dici che è banale, pur dubitando della corretta soluzione, lo do per scontato.

2. non è così banale.

Se le due macchine hanno una distribuzione di poisson $Po(L_(i))$ espressa in guasti/anno il tempo medio di durata senza guasto è un interarrivo di poisson -> si distribuisce come un'esponenziale negativa.

Ad esempio supponiamo che la distribuzione di poisson delle due macchine sia questa $Po(12)$, identica per le due macchine e le due macchine siano indipendenti; ciò significa che abbiamo le due macchine che presentano "mediamente" 12 guasti all'anno (un guasto al mese)....come si distribuisce il tempo di funzionamento di ogni macchina? Ovviamente è un'esponenziale negativa di media $1/12$ di anno (un mese).

In pratica così: $f(x)=12e^(-12x)$ se il tempo è espresso in anni oppure così $f(x)=e^(-x)$ se il tempo è espresso in mesi.
Il tuo esercizio è generico perché indica $L_(i)$ (tra l'altro diverso fra una macchina e l'altra) ma il concetto non cambia.

Ora dato che abbiamo DUE macchine e il testo ci chiede il tempo medio del PRIMO guasto è come se ci chiedesse il tempo medio della variabile $min(X,Y)$ . Quindi prima devi calcolare la distribuzione del minimo (che dovresti saper fare, con l'ipotesi di indipendenza) e poi calcolarne la media.... :oops:

PS: l'ipotesi di indipendenza delle due macchine è necessaria, anche se non specificata dal testo.

Se la spiegazione ti è chiara attendo sviluppi, altrimenti fammi sapere.

3. più articolato ma sullo stesso genere... per il momento concentriamoci sul 2. :wink:

MrBlonde93
Grazie! Diciamo che non ho mai trovato esercizi del genere, il che mi ha un po' spiazzato. Sinceramente non ho mai avuto a che fare con $min (X,Y)$. Il ragionamento sull'esponenziale negativa l'ho capito ed era qualcosa di noto.

Lo_zio_Tom
Vediamo come risolvere l'esercizio:

Dunque indichiamo con $X$ la prima macchina e $Y_(1),Y_(2)$ le due macchine successive. Indichiamo con $Z=min(X,Y_(i))$. Osserviamo subito che, $min(X,Y)>a hArr X,Y>a$

$f(x)=theta_(1)e^(-theta_(1)x)$

$f(y_(i))=theta_(2)e^(-theta_(2)y)$

Con semplici calcoli si ottiene $F_(Z)(z)=1-e^(-z(theta_(1)+theta_(2))$

quindi

1. il circuito della lavorazione è così schematizzabile




Per cui la probabilità che la lavorazione duri per almeno $p_(0)$ tempo è come dire che il minimo tempo prima dell'interruzione sia maggiore di $p_(0)$

$1-F_(Z)(p_(0))=e^(-p_(0)(theta_(1)+theta_(2))$

2. derivando la $F_(Z)$ otteniamo

$f(z)=(theta_(1)+theta_(2))e^(-z(theta_(1)+theta_(2)))~ EXP(theta_(1)+theta_(2))$

Riconosciamo subito una distribuzione esponenziale di parametro $(theta_(1)+theta_(2))$ e quindi la media è $1/((theta_(1)+theta_(2))$

3. Inserendo una seconda macchina $Y_(2)$ sussidiaria in caso di rottura di $Y_(1)$ il circuito è della lavorazione è così schematizzabile



quindi può essere risolto così:

$Y=Y_(1)+Y_(2)$

$Y~ Gamma(2;theta_(2))$

e di conseguenza la probabilità che la lavorazione non si fermi per almeno un tempo $p_(0)$ sarà (standardizzando la Gamma)

$e^(-p_(0)theta_(1))P{chi_((4))^2>=2theta_(2)p_(0)}$

dove il valore della probabilità della $chi_((4))^2$ lo leggi sulle tavole




ciao

MrBlonde93
Grazie, non trovando l'argomento nelle dispense proverò a cercare online.

Rispondi
Per rispondere a questa discussione devi prima effettuare il login.