[csidc:2004] - Re: Memorija by "Stojlljkovic Predrag" <lordps@tesla.rcub.bg.ac.yu>

Koliko sam ja skontao, algoritam minvar racuna inverznu matricu matrici R koja je dimenzija nxn gde je n broj mikrofona, pa to ispadne 16x16. Ali mislim da sam nasao nacin da resim problem memorije a da vreme izvrsavanja algoritma ostane prihvatljivo (negde oko n^4 instrukcija, pa nekoliko puta po toliko).
Ne razumem sta smatras pod numerickom stabilnoscu. Za izracunavanje determinante sam mislio da koristim algoritam svodjenja na trougaonu determinantu, posto je on najbolji u smislu iskoriscenja memorije i vremena izvrsavanja. E sad, on cesto deli elemente determinante, pa u tom smislu mozda nije bas najprecizniji, ali cu proveriti kad ga budem napravio koliko je odstupanje (videcu sta ce mi dati kao resenje kad izmnozim invertovanu i originalnu matricu).
Nadao sam se da cemo koristiti DSP sa 64K (mislim da je Gvozden pominjao TI5410) jer bih sa toliko memorije mogao da napravim visestruko brzi i precizniji algoritam, ali sad sta je tu je. Ako neko ima neku genijalnu ideju za izracunavanje determinante a da je algoritam precizniji od svodjenja na trougaonu (a u isto vreme i dovoljno brz), neka javi, a do tad cu koristiti ovaj.
            Pedja
  ----- Original Message ----- 
  From: Dusan Grujic 
  To: csidc@titan.etf.bg.ac.yu 
  Sent: Saturday, February 21, 2004 4:04 PM
  Subject: Re: [csidc] Memorija


  Imamo UKUPNO (za program i podatke :) 16k x 16 bit (instrukcije su 16-bitne). Cak i da ubacimo eksterni RAM, toliki broj pristupa eksternom ramu bi drasticno degradirao performanse. Jesi li siguran da je potreban toliki broj operacija za inverziju matrice 16 x 16?

  Inace, sa 16 mikrofona matrica kovarjanse nije 16 x 16, zbog spatial smoothing-a (K u simulaciji). Dimenzije matrice su zapravo (16-K+1) x (16-K+1), sto u nasem slucaju dodje 12 x 12.

  Posebno obrati paznju na numericku stabilnost algoritma. Ako se koristi najosnovniji algoritam, tj. izracunavanje determinante po definiciji, cak i za matrice 3 x 3 resenje moze biti netacno.

  Procitao sam u radu (onom sa SIMRAD-ovim sonarom) da bi matrica kovarjanse idealno trebalo da bude tzv. Toeplic - ova matrica (konjugovana simetrija po obe dijagonale). Trebalo bi da za ovaj tip matrice postoji algoritam koji je mnogo efikasniji od opsteg algoritma za inverziju matrice.

  Pozdrav,

  Gruja 



  -- 
  ___________________________________________________________
  Sign-up for Ads Free at Mail.com
  http://www.mail.com/?sr=signup

  ----------------------------------------------------------------- Belgrade CSIDC Team Home Page: http://titan.etf.bg.ac.yu/csidc/ ----------------------------------------------------------------- unsubscribe: minimalist@titan.etf.bg.ac.yu?subject=unsubscribe%20csidc -----------------------------------------------------------------