Visualizza wikitesto di Calcolatore serie di Fourier


  {{Risorsa
| tipo = lezione
| materia1 = Sistemi di calcolo automatico per il sonar
| avanzamento = 100%
}}

-Processo di approssimazione per la determinazione dello spettro dei segnali periodici-

==La serie di Fourier con il supporto delle figure==

Prima di descrivere il calcolatore della serie di Fourier, al quale è dedicato il paragrafo successivo, esaminiamo  l'algoritmo  omonimo mostrato in figura 1

:<math>
\begin{align}
F(x) &= a_0+\sum_{n=1}^\infin\left( a_n \cos nx + b_n \sin nx \right) \\
&\begin{cases}
a_0 = \frac{1}{2\pi}\int_{-\pi}^{+\pi}f(x)dx \\
a_n = \frac{1}{\pi}\int_{-\pi}^{+\pi}f(x)\cos nx\ dx \\
b_n = \frac{1}{\pi}\int_{-\pi}^{+\pi}f(x)\sin nx\ dx
\end{cases} \qquad\qquad\qquad\qquad\qquad 1)
\end{align}
</math>

questo consente, per via esclusivamente analitica, lo studio dei fenomeni periodici  funzioni  del tempo per il computo del loro spettro in frequenza, spettro espresso
in figura 1 con il simbolo <math> F(x)</math> ( dove <math> x = 2 \cdot \pi \cdot f \cdot t )</math>.

Generalmente la funzione <math> F(x)</math> è espressa con il simbolo <math> G(\omega)</math>.

Un esempio di segnale periodico <math> f(t)</math>, onda rettangolare, espresso matematicamente è mostrato in figura 2a

:<math>f(2\pi f_0 t) =
\begin{cases}
-E: & \text{per} \ (2k-1)\pi < 2 \pi f_0 t < 2k\pi \\
+E: & \text{per} \ 2k\pi < 2 \pi f_0 t < (2k+1)\pi
\end{cases} \qquad\qquad 2a)
</math>

il cui andamento è tracciato in figura 2b
 
[[File:trf2b.jpg|thumb|left|400px|figura 2b]]
{{clear}}

L'impiego della serie di Fourier porta  alla determinazione dello spettro di frequenza a righe <math> G( \omega)</math>, dell'onda di figura 2b,così come mostrato in figura 3.

[[File:trf3.jpg|thumb|left|400px|figura 3]]
{{clear}}
 
Se dell'onda periodica da analizzare non si conosce la legge <math> f(t)</math>, come ad esempio per la curva di figura 4: 

[[File:trf4.jpg|thumb|left|400px|figura 4]]
{{clear}}

la serie di Fourier non è applicabile; è sostituibile, con buona approssimazione,
dall'algoritmo di figura 5 che prevede la determinazione dello spettro mediante l'elaborazione dei coefficienti <math> An'</math>  e <math>  Bn' </math> da <math> n </math> campioni di <math> f(t) </math> , messi a calcolo, rilevati da grafici o da presentazioni su oscilloscopio:

:<math>C'_n = A'_0 + \sum_{n=1}^{n=k}A'_n \cos (n 2 \pi f_0 t) + B'_n \sin (n 2 \pi f_0 t) \qquad\qquad\qquad 5)</math>

Si osservi che:

*Con il simbolo <math> C'n </math> s'intende l'ampiezza della n esima riga dello spettro di frequenza.

* I coefficienti di figura 5 sono indicati con <math> A'n </math> e <math> B'n</math>  per distinguerli da quelli diversi di figura 1.

* <math> A'o,  A'n,  B'n </math> sono determinati con processo numerico iterativo implementando nel programma di calcolo le funzioni esposte in figura 5a dove, ad esempio, <math> 20 </math>  campioni del periodo in analisi sono riportati, in apposita memoria del P.C, con i simboli <math> a(1).....a(20) </math> .

:<math>
\begin{align}
& A'_0 = \frac{1}{2\pi} \left[ 0,5a(1) + a(2) + a(3) + \dots + a(19) + 0,5a(20) \right] \frac{2\pi}{19} \qquad\qquad\qquad\qquad\qquad\qquad 5a) \\
& A'_n = \frac{1}{\pi} \left[ 0,5a(1)\cos(nwt_1) + a(2)\cos(nwt_2) + \dots + 0,5a(20)\cos(nwt_{20}) \right] \frac{2\pi}{19} \\
& B'_n = \frac{1}{\pi} \left[ 0,5a(1)\sin(nwt_1) + a(2)\sin(nwt_2) + \dots + 0,5a(20)\sin(nwt_{20}) \right] \frac{2\pi}{19} \\
& C'_n = \left({A'_n}^2 + {B'_n}^2\right)^\frac{1}{2}
\end{align}
</math>

* Il modulo della <math> n </math> esima  riga dello spettro è indicato da  <math>  C'n</math> ; nel calcolatore l valori <math>  C'n </math>  sono valutati senza l'addizione di <math>  A'o </math> che viene calcolato separatamente.  

Vedremo di seguito, in dettaglio, l'impiego degli algoritmi di figure 5 e 5a contestualmente all'esame del pannello operativo del calcolatore.

== Il pannello operativo del calcolatore ==
Il pannello operativo, mostrato in figura 6, è  generato da un  file eseguibile disponibile all'indirizzo:  [http://www.sonar-info.info/wikifourier/seriefourier.html  wikiFourier ] oppure
[https://github.com/2021dtc/prova20/archive/main.zip  seriefourier]
 
Il file esguibile consente il calcolo e la presentazione grafica  dello spettro dell'onda da analizzare.

[[File:tfr6.jpg|thumb|left|1000px|figura 6]]
{{clear}}

In figura 6, con una serie di numeri in rosso, sono distinte le 5 sezioni funzionali: 

'''Sezione 1 (Variabili di calcolo)''' 

La sezione contiene il pulsante ingresso variabili  e  4 textbox  nelle quali inserire i seguenti dati: 

*Durata del periodo dell'onda da analizzare espressa con <math>  T </math> in secondi.

*Il numero <math> m = n </math>  dei campioni da prelevare dal periodo dell'onda da analizzare.

*Il numero <math> k </math>  delle righe dello spettro che si vogliono calcolare.

*L'incremento di calcolo <math> dm </math>  dal quale dipendono sia la precisione che il tempo calcolo.
 
Al lancio del programma la sezione è impostata, in via preliminare, con i valori: 

<math> T = 0.01 \ S </math>; <math> \quad m = n =  20 </math>; <math>\quad  k = 10 </math>; <math> \quad dm = 0.01</math> .

'''Sezione 2: (Campioni nel periodo)'''

La sezione contiene:
*Un  label, a sinistra, che indica il numero progressivo dei campioni ricavati dal periodo da analizzare e un textbox, a destra, per l'inserzione delle ampiezze dei diversi campioni.

*Un pulsante "Ingresso campioni" per l'inserzione progressiva dei valori sopra citati;il pulsante presenta colore verde nella fase d'introduzione campioni, cambia colore in arancio una volta che i campioni dell'onda sono pari al valore impostato di <math> m = n </math> 

*Un pulsante " Azzeramento memorie" per la ripetizione, se necessario, dell'inserzione dati.  

'''Sezione 3: ( Presentazione spettro)'''

La sezione contiene:

*Il pulsante "Calcolo dello spettro" per l'avvio del programma di calcolo per la  presentazione grafica del modulo delle righe dello spettro tracciate come visione indicativa dello stesso.

*Un listato numerico delle frequenze e relative ampiezze normalizzate delle righe dello spettro. 
L'azione di questo comando provoca inoltre il tracciamento del periodo dell'onda inserito a campioni nella sezione 2 per un controllo della correttezza del suo profilo.

'''Sezione 4: (Listato valori calcolati)'''

Nel listato compaiono:

*Le <math>  n </math>  coppie calcolate nella forma, ad esempio per una generica <math> 3^a </math>  coppia, <math>  F = 1300 \ Hz</math>  <math> \ C'3 = 0.65432 </math> ; ciascuna coppia definisce una riga dello spettro sia in frequenza che in ampiezza normalizzata. 

*In alto, sopra il listato, il valore calcolato di <math>C'o</math> ( componente continua dell'onda ove sia presente )


'''Sezione 5: (Area di presentazione grafici)''' 

I grafici sono tracciati a solo scopo illustrativo dato che i valori che li generano sono già contenuti nei listati della sezione 4.

*In alto lo spettro dell'onda con le <math> m = n </math>  righe messe a calcolo; le ascisse del tracciato non sono calibrate in frequenza

*In basso il profilo del periodo dell'onda definito con tanti punti quanti sono gli <math> m = n </math>  campioni messi a calcolo; le ascisse del tracciato non sono calibrate in tempo.

==Primo  esempio di calcolo della serie di Fourier con controllo analitico==

Questo esempio, che interessa l'analisi di un'onda rettangolare, si avvale di un riscontro analitico in quanto lo spettro dell'onda in esame è definito dalla caratteristica serie di figura 7 che servirà come controllo dati elaborati dal calcolatore:

:<math>y=\frac{4}{\pi}E\left(\cos x - \frac{1}{3}\cos 3x + \frac{1}{5}\cos 5x - \frac{1}{7} \cos 7x + \dots \right) \qquad\qquad\qquad 7)</math>


Come si vede la serie analitica è formata soltanto da addendi dispari i cui coefficienti, espressi con frazioni,  sono gli analoghi dei termini <math> Cn </math> elaborati dal calcolatore.


Utilizziamo sempre i valori preliminari assegnati al calcolatore al momento dell'avvio,
come   mostra la figura 8 della sezione 1.
 
[[File:tfr8.jpg|thumb|left|300px|figura 8]]
{{clear}}


Quindi pigiamo il pulsante "Ingresso variabili". 
 
Impostiamo nella sezione 2, premendo " Ingresso campioni" (pulsante arancione),  per lo spettro di un'onda rettangolare quale quella  riportata in figura 2b, la seguente sequenza di campioni:

<math>1;1;1;1;1;1;1;1;1;1; -1;-1;-1;-1;-1;-1;-1;-1;-1;-1 </math>

relativa ad un periodo  dell'onda sopra citata fino a quando il pulsante stesso diventa arancione.

Dopo la pressione del pulsante " Calcolo dello spettro ", si ha sia il listato dei valori così come riportato in  figura 9 che i grafici di figura 10

[[File:tfr9.jpg|thumb|left|600px|figura 9]]
{{clear}}


La lista indica righe multiple di <math> 1000 \ Hz </math> dato che, essendo <math>T = 0.001 \ S </math>, la frequenza è <math>f = 1 / 0.001</math> = <math>1000 \ Hz </math>.  

[[File:tfr10.jpg|thumb|left|300px|figura 10]]
{{clear}}

 
Procediamo ora al controllo della lista di figura 9 impostando una tabella che riporta, tanto le ampiezze calcolate con l'espressione dello spettro teorico espresso in figura  7, che quelle del listato di figura 9:

:{| class="wikitable"
|-
! Righe>
! C'<sub>1</sub>
! C'<sub>2</sub>
! C'<sub>3</sub>
! C'<sub>4</sub>
! C'<sub>5</sub>
! C'<sub>6</sub>
|-
! Con il calcolatore
| 1,0000
| 0,0012
| 0,3333
| 0,0012
| 0,2000
| 0,0012
|-
! Per via analitica
| 1,0000
| 0
| 0,3333
| 0
| 0,2000
| 0
|}

Come si vede in tabella di tavolo i valori di <math>C'n </math> per <math> n </math> dispari, praticamente coincidono, per <math> n </math> pari il calcolatore commette un errore presentando ampiezze pari a <math>0.002 </math> invece che <math> 0 </math> come da processo analitico.

Gli errori si riducono ripetendo il computo con un passo di calcolo inferiore a <math>dm = 0.01 </math>, ad esempio <math>dm = 0.001</math>; naturalmente il processo di calcolo richiede un tempo superiore di sviluppo.

==Secondo esempio: calcolo della serie di Fourier con controllo analitico==

Questo esempio, che interessa l'analisi di un'onda a dente di sega, si avvale di un riscontro analitico in quanto lo spettro dell'onda in esame è definito dalla caratteristica serie di figura 11 che servirà come controllo dei
dati elaborati dal calcolatore:

:<math>y=\frac{2}{\pi}E\left(\sin x - \frac{1}{2}\sin 2x + \frac{1}{3}\sin 3x - \frac{1}{4} \sin 4x + \dots \right) \qquad\qquad\qquad 11)</math>

Come si vede la serie analitica è formata  da addendi pari e dispari i cui coefficienti, espressi con frazioni,  sono gli analoghi dei termini <math>Cn</math> elaborati dal calcolatore.


Utilizziamo sempre i valori di base assegnati al calcolatore al momento dell'avvio,
come   mostra la figura 12 della sezione 1.

[[File:tfr8.jpg|thumb|left|600px|figura 12]]
{{clear}}


quindi pigiamo il pulsante "Ingresso variabili". 

Impostiamo nella sezione 2,  premendo " Ingresso campioni" ( pulsante giallo ),  i campioni d'ampiezza di  un'onda a dente di sega  quale quella riportata in figura 13 

[[File:tfr13.jpg|thumb|left|300px|figura 13]]
{{clear}}



la seguente sequenza di campioni:

<math> -1;-0.9;-0.8;-0.7;-0.6;-0.5;-0.4;-0.3;-0.2; 0.1;0;0.1;0.2;0.3;0.4;0.5;0.6;0.7:0.8:0.9  </math>

relativa ad un periodo <math>T</math> dell'onda sopra citata fino a quando
il pulsante stesso diventa arancione.

Dopo la pressione del pulsante " Calcolo dello spettro ", si ha sia il listato dei valori così come riportato in  figura 14 che i grafici di figura 15:
 
[[File:tfr14.jpg|thumb|left|600px|figura 14]]
{{clear}}


La lista indica righe multiple di <math>1000 Hz </math> dato che, essendo <math>T = 0.001 \ S </math>, <math>f = 1/ 0.001 </math>= <math>1000 \ Hz</math>. 

[[File:tfr15.jpg|thumb|left|300px|figura 15]]
{{clear}}


Procediamo ora al controllo della lista di figura 14 impostando una tabella che riporta, tanto le ampiezze dello spettro teorico espresso in figura  11 che quelle della lista nella figura citata:

:{| class="wikitable"
|-
! Righe>
! C'<sub>1</sub>
! C'<sub>2</sub>
! C'<sub>3</sub>
! C'<sub>4</sub>
! C'<sub>5</sub>
! C'<sub>6</sub>
|-
! Con il calcolatore
| 1,0000
| 0,5001
| 0,3335
| 0,2501
| 0,2001
| 0,1667
|-
! Per via analitica
| 1,0000
| 0,5000
| 0,3333
| 0,2500
| 0,2000
| 0,1666
|}

Come si vede in tabella di tavolo i valori di <math> C' </math> da calcolatore sono quasi coincidenti con dagli analoghi analitici; l'errore massimo è dell'ordine del  <math>2 </math> per mille.

== Terzo esempio: calcolo della serie di Fourier di un'onda generica ==

Questo esempio, che interessa l'analisi di un'onda generica, non si avvale alcun riscontro analitico.

Il tracciato dell'onda sia quello riportato in figura 16:
 
[[File:trf16.jpg|thumb|left|300px|figura 16]]
{{clear}}


Si esegue la campionatura di un periodo con 21 campioni così come mostrato in figura 17:

[[File:trf17.jpg|thumb|left|300px|figura 17]]
{{clear}}


Si compila, in tavola, la lista delle ampiezze dei  campioni ricavandola da rilievi geometrici su figura 17:

:{| class="wikitable"
|-
! periodo
! ampiezza
|-
| t<sub>0</sub>
| 0,00
|-
| t<sub>1</sub>
| 0,40
|-
| t<sub>2</sub>
| 0,64
|-
| t<sub>3</sub>
| 0,75
|-
| t<sub>4</sub>
| 0,85
|-
| t<sub>5</sub>
| 0,99
|-
| t<sub>6</sub>
| 1,04
|-
| t<sub>7</sub>
| 0,86
|-
| t<sub>8</sub>
| 0,53
|-
| t<sub>9</sub>
| 0,21
|-
| t<sub>10</sub>
| 0,00
|-
| t<sub>11</sub>
| -0,02
|-
| t<sub>12</sub>
| -0,17
|-
| t<sub>13</sub>
| -0,38
|-
| t<sub>14</sub>
| -0,47
|-
| t<sub>15</sub>
| -0,40
|-
| t<sub>16</sub>
| -0,28
|-
| t<sub>17</sub>
| -0,26
|-
| t<sub>18</sub>
| -0,29
|-
| t<sub>19</sub>
| -0,22
|-
| t<sub>20</sub>
| -0,00
|}

I dati di tavola s'inseriscono progressivamente in sezione 2 pigiando il pulsante "Ingresso campioni". 

Dopo la pressione del pulsante " Calcolo dello spettro ", si ha sia il listato dei valori così come riportato in figura 19 che i grafici di figura 20:

[[File:trf19.jpg|thumb|left|300px|figura 19]]
{{clear}}


La lista indica righe multiple di <math>1000 \ Hz</math> dato che, essendo
<math> T = 0.001 \ S </math>, <math>\quad f = 1/ 0.001 </math> = <math> 1000 \ Hz. </math>

[[File:trf20.jpg|thumb|left|300px|figura 20]]
{{clear}}

== Quarto esempio: calcolo della serie di Fourier di un'onda precostruita==

Per un controllo ulteriore del calcolatore analizziamo un'onda precostruita che è caratterizzata dalla presenza della sola <math> 4^a </math> armonica di ampiezza <math>C'4 = 0.3333.</math> 

L'onda in oggetto, generata da un'apposita funzione matematica, è riportata in figura 21:

[[File:trf21.jpg|thumb|left|300px|figura 21]]
{{clear}}


L'onda ha un  periodo <math> T = 0.001 \ S </math> e quindi una frequenza di <math> 1000 \ Hz</math>; il periodo viene campionato con le stesse modalità relative al periodo dell'esercizio precedente; il risultato della campionatura è riportato in tavola:

:{| class="wikitable"
|-
! periodo
! ampiezza
|-
| t<sub>0</sub>
| 0,00
|-
| t<sub>1</sub>
| 0,62
|-
| t<sub>2</sub>
| 0,78
|-
| t<sub>3</sub>
| 0,61
|-
| t<sub>4</sub>
| 0,63
|-
| t<sub>5</sub>
| 0,99
|-
| t<sub>6</sub>
| 1,26
|-
| t<sub>7</sub>
| 1,00
|-
| t<sub>8</sub>
| 0,39
|-
| t<sub>9</sub>
| -0,00
|-
| t<sub>10</sub>
| -0,00
|-
| t<sub>11</sub>
| 0,00
|-
| t<sub>12</sub>
| -0,38
|-
| t<sub>13</sub>
| -1,00
|-
| t<sub>14</sub>
| -1,26
|-
| t<sub>15</sub>
| -1,00
|-
| t<sub>16</sub>
| -0,63
|-
| t<sub>17</sub>
| -0,61
|-
| t<sub>18</sub>
| -0,78
|-
| t<sub>19</sub>
| -0,63
|-
| t<sub>20</sub>
| -0,00
|}

Seguendo la procedura d'inserzione dati e campioni analoga a quella del paragrafo precedente, dopo la pressione del pulsante " Calcolo dello spettro ", si ha sia il listato dei valori così come riportato in figura 23 che i grafici di figura 24:

[[File:trf23.jpg|thumb|left|600px|figura 23]]
{{clear}}


[[File:trf24.jpg|thumb|left|300px|figura 24]]
{{clear}}


L'esame dei dati, sia numerici ( riga <math>C'4</math> ) che grafici evidenziano, come era nelle aspettative, la presenza di una riga a <math>4000 \ Hz </math> di ampiezza <math>C'4 = 0.3133</math>. 

Le righe a livello molto basso presenti in figura 23, non visibili in figura 24, sono dovute all'approssimazione  dell’algoritmo di calcolo.

==Bibliografia==
*G. Moretti, ''Analisi matematica-vol.II -parte II'',Hoepli, Milano 1953.
*A. Papoulis, ''The Fourier integral and its applications''Mc Graw-Hill, New York 1962.
*Edit: Milton Abramowitz ..''Handbook of mathematical functios'', USA 1970.
* C. Del Turco, ''La matematica con il personal computer'', Editrice MODERNA La Spezia 1998.

[[Categoria:Lezioni di Sistemi di calcolo automatico per il sonar]]

[[Categoria:Calcolatori automatici per il sonar]]