Modellistica di Impianti e sistemi 2 2007/2008

lionel ha scritto:

Logue ha scritto:

lionel ha scritto:

rik ha scritto:

anche noi per la matrice ancora non abbiamo deciso niente.
I metodi iterativi non ci convincono, perchè ne abbiamo provati un po' e non convergono quasi mai.
Stiamo pensando di integrare matlab.
Stavamo anche provando a usare gauss, ma abbiamo un problema.
Dato che la matrice Q deve avere tutti i termini noti = 0, con gauss ci viene sempre la soluzione banale (tutti 0 ), e non è molto utile.
Per chi ha provato a usare gauss, ha avuto lo stesso problema? Come avete risolto?

In una delle righe della matrice (dopo che l'hai trasposta) devi sostituire la condizione di normalizzazione, cioè che la somma delle probabilità di stato sia 1.
Ovviamente devi cambiare anche il termine noto.
Noi abbiamo utilizzato una libreria java (mtj) e, con un pò di pazienza, siamo riusciti a risolvere sistemi con più di 300000 stati.

Una piccola curiosità, quanto impiegate a risolvere la matrice di queste dimensioni con mtj?

Su una notebook con 1 GB di ram ci mettiamo due-tre ore. Però nel tempo devi includere anche generazione stati e costruzione matrice di transizione.
Credo comunque che con Matlab la risoluzione si velocizza di moooolto!!

Ma la matrice che gestite è sparsa o piena?

La matrice è sparsa. Se conosci la libreria usiamo il tipo di dato FlexCompRowMatrix.

E' stato svolto a lezione il calcolo degli indici per l'infinite server?
Io negli appunti presi a lezioni, nè sui libri lo ritrovo, voi come li avete calcolati ?

ciao

eraclite ha scritto:

E' stato svolto a lezione il calcolo degli indici per l'infinite server?
Io negli appunti presi a lezioni, nè sui libri lo ritrovo, voi come li avete calcolati ?

ciao

quoto

eraclite ha scritto:

E' stato svolto a lezione il calcolo degli indici per l'infinite server?
Io negli appunti presi a lezioni, nè sui libri lo ritrovo, voi come li avete calcolati ?

ciao

Stesso identico problema, non sappiamo come gestirlo.

lionel ha scritto:

La matrice è sparsa. Se conosci la libreria usiamo il tipo di dato FlexCompRowMatrix.

Ma a voi mtj risolve tutte le matrici con la stessa combinazione di IterativeSolver-Preconditioner senza trovare valori negativi come soluzione ?


ciao