11.07.2015 Views

Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...

Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...

Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

10.STIMA DEI PARAMETRI DELLA TURBOLENZA ATMOSFERICA—————————————————————⎯⎯⎯⎯⎯⎯⎯⎯————————APPENDICE 10.1Data un vettore x di N elementi, la routine FOURIER realizza la DFT del segnale x e produce i vettorixr e xi che sono rispettivamente la parte reale e la parte immaginaria della DFT di x. Se ier non èuguale a zero, N non è una potenza di 2.! =================================================================Subroutine FOURIER(N,x,xr,xi,ier)! =================================================================! Calcolo della DFT del vettore x!! INPUT:! N = numero di elemnti del vettore (potenza di 2)! x = vettore di cui si vuole la DFT! OUTPUT:! xr = parte reale della DFT! xi = parte immaginaria della DFT! ier = 0 (nessun problema)! 1 (N non e' potenza di 2)! Da: Proakis e Manolakis (1992)! ----------------------------------------------------------------REAL, ALLOCATABLE :: wr(:)REAL, ALLOCATABLE :: wi(:)REAL x(N),xi(N),xr(N)!ier = 1DO i=1,15m = in2 = 2**iIF(N.EQ.n2) Thenier = 0EXITEndifEndDOIF(ier.EQ.1) RETURNALLOCATE(wr(N),wi(N))!DO i=1,Nxr(i) = x(i)xi(i) = 0.EndDO!Call WTABLE(N,wr,wi)Call FFT(N,xr,xi,wr,wi)!DEALLOCATE(wr,wi)RETURNEND!! =================================================================Subroutine FFT(N,xr,xi,wr,wi)! =================================================================REAL xr(N),xi(N),wr(N),wi(N)n2 = NDO i=1,m————————————————————————————————————————- 500 -

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!