Algoritmo Transformada

Discreta de Fourier.

M. en C. y M. en E. Manuel Aguila Muñoz

Expo Instrumentos Digitales para la asignatura de Mediciones Eléctricas

IPN/ESIME Zacatenco/Ingeniería Eléctrica/Electrotecnia

Junio -2022
Algoritmo Transformada Discreta de Fourier.
Variantes: Algoritmo
• Ciclo completo
Algoritmo Transformada Discreta de Fourier.
Variantes: Algoritmo
• Medio ciclo
Aplicación TDF Ciclo completo
Señal de entrada: h1 h3, w1=2*pi*f1, f1=60 [Hz], nmpc=12, nc=2
k v(t)
1er Ventana = . ∗ ! " ∗ # ° + 8.5∗ ! #∗" ∗ °
1 9.226273229
2 19.31263991
3 10.00342317
4 1.824243019
5 7.37703559
6 7.142711675
7 -9.22627323
8 -19.3126399
9 -10.0034232
10 -1.82424302
11 -7.37703559
12 -7.14271167
13 9.226273229
14 19.31263991
15 10.00342317
16 1.824243019
17 7.37703559
18 7.142711675
19 -9.22627323
20 -19.3126399
21 -10.0034232
22 -1.82424302
23 -7.37703559
24 -7.14271167
Aplicación TDF Ciclo completo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1
1 19.3126399 0.5 0.8660254
2 10.0034232 0.8660254 0.5
3 1.82424302 1 2.833E-16
4 7.37703559 0.8660254 -0.5
5 7.14271167 0.5 -0.8660254
6 -9.2262732 5.666E-16 -1
7 -19.31264 -0.5 -0.8660254
8 -10.003423 -0.8660254 -0.5
9 -1.824243 -1 -1.838E-16
10 -7.3770356 -0.8660254 0.5
11 -7.1427117 -0.5 0.8660254
Aplicación TDF Ciclo completo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1 0 9.22627323
1 19.3126399 0.5 0.8660254 9.65631996 16.7252368
2 10.0034232 0.8660254 0.5 8.66321859 5.00171158
3 1.82424302 1 2.833E-16 1.82424302 5.1681E-16
4 7.37703559 0.8660254 -0.5 6.38870023 -3.6885178
5 7.14271167 0.5 -0.8660254 3.57135584 -6.1857698
6 -9.2262732 5.666E-16 -1 -5.228E-15 9.22627323
7 -19.31264 -0.5 -0.8660254 9.65631996 16.7252368
8 -10.003423 -0.8660254 -0.5 8.66321859 5.00171158
9 -1.824243 -1 -1.838E-16 1.82424302 3.3525E-16
10 -7.3770356 -0.8660254 0.5 6.38870023 -3.6885178
11 -7.1427117 -0.5 0.8660254 3.57135584 -6.1857698
Aplicación TDF Ciclo completo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1 0 9.22627323 0 1.5377122
1 19.3126399 0.5 0.8660254 9.65631996 16.7252368 1.60938666 4.32525167
2 10.0034232 0.8660254 0.5 8.66321859 5.00171158 3.05325642 5.15887027
3 1.82424302 1 2.833E-16 1.82424302 5.1681E-16 3.35729693 5.15887027
4 7.37703559 0.8660254 -0.5 6.38870023 -3.6885178 4.4220803 4.5441173
5 7.14271167 0.5 -0.8660254 3.57135584 -6.1857698 5.01730627 3.51315567
6 -9.2262732 5.666E-16 -1 -5.228E-15 9.22627323 5.01730627 5.05086788
7 -19.31264 -0.5 -0.8660254 9.65631996 16.7252368 6.62669293 7.83840734
8 -10.003423 -0.8660254 -0.5 8.66321859 5.00171158 8.0705627 8.67202594
9 -1.824243 -1 -1.838E-16 1.82424302 3.3525E-16 8.3746032 8.67202594
10 -7.3770356 -0.8660254 0.5 6.38870023 -3.6885178 9.43938657 8.05727297
11 -7.1427117 -0.5 0.8660254 3.57135584 -6.1857698 10.0346125 7.02631135
Aplicación TDF Ciclo completo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1 0 9.22627323 0 1.5377122 1.5377122
1 19.3126399 0.5 0.8660254 9.65631996 16.7252368 1.60938666 4.32525167 4.61496776
2 10.0034232 0.8660254 0.5 8.66321859 5.00171158 3.05325642 5.15887027 5.99469075
3 1.82424302 1 2.833E-16 1.82424302 5.1681E-16 3.35729693 5.15887027 6.15511048
4 7.37703559 0.8660254 -0.5 6.38870023 -3.6885178 4.4220803 4.5441173 6.34064635
5 7.14271167 0.5 -0.8660254 3.57135584 -6.1857698 5.01730627 3.51315567 6.125
6 -9.2262732 5.666E-16 -1 -5.228E-15 9.22627323 5.01730627 5.05086788 7.11931377
7 -19.31264 -0.5 -0.8660254 9.65631996 16.7252368 6.62669293 7.83840734 10.2641945
8 -10.003423 -0.8660254 -0.5 8.66321859 5.00171158 8.0705627 8.67202594 11.8464347
9 -1.824243 -1 -1.838E-16 1.82424302 3.3525E-16 8.3746032 8.67202594 12.0556216
10 -7.3770356 -0.8660254 0.5 6.38870023 -3.6885178 9.43938657 8.05727297 12.4105466
11 -7.1427117 -0.5 0.8660254 3.57135584 -6.1857698 10.0346125 7.02631135 12.25
Aplicación TDF Ciclo completo

¿y la segunda ventana?
Aplicación TDF Ciclo completo
Señal de entrada: h1 h3, w1=2*pi*f1, f1=60 [Hz], nmpc=12, nc=2
k v(t)
2da Ventana = . ∗ ! " ∗ # ° + 8.5∗ ! #∗" ∗ °
1 9.226273229
2 19.31263991
3 10.00342317
4 1.824243019
5 7.37703559
6 7.142711675
7 -9.22627323
8 -19.3126399
9 -10.0034232
10 -1.82424302
11 -7.37703559
12 -7.14271167
13 9.226273229
14 19.31263991
15 10.00342317
16 1.824243019
17 7.37703559
18 7.142711675
19 -9.22627323
20 -19.3126399
21 -10.0034232
22 -1.82424302
23 -7.37703559
24 -7.14271167
Aplicación TDF Ciclo completo
2da Ventana
k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp
0 19.3126399 0 1 0 19.3126399 0 3.21877332 3.21877332
1 10.0034232 0.5 0.8660254 5.00171158 8.66321859 0.8336186 4.66264308 4.73657688
2 1.82424302 0.8660254 0.5 1.5798408 0.91212151 1.0969254 4.81466334 4.93803892
3 7.37703559 1 2.833E-16 7.37703559 2.0899E-15 2.32643133 4.81466334 5.34726713
4 7.14271167 0.8660254 -0.5 6.18576976 -3.5713558 3.35739296 4.21943736 5.39219242
5 -9.2262732 0.5 -0.8660254 -4.6131366 7.990187 2.58853685 5.5511352 6.125
6 -19.31264 5.666E-16 -1 -1.094E-14 19.3126399 2.58853685 8.76990851 9.14394982
7 -10.003423 -0.5 -0.8660254 5.00171158 8.66321859 3.42215545 10.2137783 10.7718343
8 -1.824243 -0.8660254 -0.5 1.5798408 0.91212151 3.68546225 10.3657985 11.0014731
9 -7.3770356 -1 -1.838E-16 7.37703559 1.3557E-15 4.91496818 10.3657985 11.471996
10 -7.1427117 -0.8660254 0.5 6.18576976 -3.5713558 5.94592981 9.77057256 11.4375771
11 9.22627323 -0.5 0.8660254 -4.6131366 7.990187 5.17707371 11.1022704 12.25
Aplicación TDF Ciclo completo

¿y la tercer ventana?
Aplicación TDF Ciclo completo
Señal de entrada: h1 h3, w1=2*pi*f1, f1=60 [Hz], nmpc=12, nc=2
k v(t)
3er Ventana = . ∗ ! " ∗ # ° + 8.5∗ ! #∗" ∗ °
1 9.226273229
2 19.31263991
3 10.00342317
4 1.824243019
5 7.37703559
6 7.142711675
7 -9.22627323
8 -19.3126399
9 -10.0034232
10 -1.82424302
11 -7.37703559
12 -7.14271167
13 9.226273229
14 19.31263991
15 10.00342317
16 1.824243019
17 7.37703559
18 7.142711675
19 -9.22627323
20 -19.3126399
21 -10.0034232
22 -1.82424302
23 -7.37703559
24 -7.14271167
Aplicación TDF Ciclo completo
3er Ventana
k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp
0 10.0034232 0 1 0 10.0034232 0 1.66723719 1.66723719
1 1.82424302 0.5 0.8660254 0.91212151 1.5798408 0.15202025 1.93054399 1.93652014
2 7.37703559 0.8660254 0.5 6.38870023 3.68851779 1.21680362 2.54529696 2.82119614
3 7.14271167 1 2.833E-16 7.14271167 2.0235E-15 2.40725557 2.54529696 3.50334354
4 -9.2262732 0.8660254 -0.5 -7.990187 4.61313661 1.07555774 3.31415306 3.4843127
5 -19.31264 0.5 -0.8660254 -9.65632 16.7252368 -0.5338289 6.10169253 6.125
6 -10.003423 5.666E-16 -1 -5.668E-15 10.0034232 -0.5338289 7.76892972 7.7872487
7 -1.824243 -0.5 -0.8660254 0.91212151 1.5798408 -0.3818087 8.03223652 8.04130595
8 -7.3770356 -0.8660254 -0.5 6.38870023 3.68851779 0.6829747 8.64698949 8.67391962
9 -7.1427117 -1 -1.838E-16 7.14271167 1.3126E-15 1.87342664 8.64698949 8.84760728
10 9.22627323 -0.8660254 0.5 -7.990187 4.61313661 0.54172881 9.41584559 9.43141656
11 19.3126399 -0.5 0.8660254 -9.65632 16.7252368 -1.0676578 12.2033851 12.25
Aplicación TDF Ciclo completo

¿y para el algoritmo de medio ciclo?

Aplicación TDF Medio ciclo
Señal de entrada: h1 h3, w1=2*pi*f1, f1=60 [Hz], nmpc=12, nc=2
k v(t) = . ∗ ! " ∗ # ° + 8.5∗ ! #∗" ∗ °
1 9.226273229
1er Ventana
2 19.31263991
3 10.00342317
4 1.824243019
5 7.37703559
6 7.142711675
7 -9.22627323
8 -19.3126399
9 -10.0034232
10 -1.82424302
11 -7.37703559
12 -7.14271167
13 9.226273229
14 19.31263991
15 10.00342317
16 1.824243019
17 7.37703559
18 7.142711675
19 -9.22627323
20 -19.3126399
21 -10.0034232
22 -1.82424302
23 -7.37703559
24 -7.14271167
Aplicación TDF Medio Ciclo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1
1 19.3126399 0.5 0.8660254
2 10.0034232 0.8660254 0.5
3 1.82424302 1 2.833E-16
4 7.37703559 0.8660254 -0.5
5 7.14271167 0.5 -0.8660254
6 -9.2262732 5.666E-16 -1
Aplicación TDF Medio Ciclo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1 0 9.22627323
1 19.3126399 0.5 0.8660254 9.65631996 16.7252368
2 10.0034232 0.8660254 0.5 8.66321859 5.00171158
3 1.82424302 1 2.833E-16 1.82424302 5.1681E-16
4 7.37703559 0.8660254 -0.5 6.38870023 -3.6885178
5 7.14271167 0.5 -0.8660254 3.57135584 -6.1857698
Aplicación TDF Medio Ciclo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1 0 9.22627323 0 3.07542441
1 19.3126399 0.5 0.8660254 9.65631996 16.7252368 3.21877332 8.65050334
2 10.0034232 0.8660254 0.5 8.66321859 5.00171158 6.10651285 10.3177405
3 1.82424302 1 2.833E-16 1.82424302 5.1681E-16 6.71459385 10.3177405
4 7.37703559 0.8660254 -0.5 6.38870023 -3.6885178 8.8441606 9.0882346
5 7.14271167 0.5 -0.8660254 3.57135584 -6.1857698 10.0346125 7.02631135
Aplicación TDF Medio Ciclo

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1 0 9.22627323 0 3.07542441 3.07542441
1 19.3126399 0.5 0.8660254 9.65631996 16.7252368 3.21877332 8.65050334 9.22993552
2 10.0034232 0.8660254 0.5 8.66321859 5.00171158 6.10651285 10.3177405 11.9893815
3 1.82424302 1 2.833E-16 1.82424302 5.1681E-16 6.71459385 10.3177405 12.310221
4 7.37703559 0.8660254 -0.5 6.38870023 -3.6885178 8.8441606 9.0882346 12.6812927
5 7.14271167 0.5 -0.8660254 3.57135584 -6.1857698 10.0346125 7.02631135 12.25
Aplicación TDF Medio ciclo
2da Ventana
k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp
0 19.3126399 0 1 0 19.3126399 0 6.43754664 6.43754664
1 10.0034232 0.5 0.8660254 5.00171158 8.66321859 1.66723719 9.32528617 9.47315375
2 1.82424302 0.8660254 0.5 1.5798408 0.91212151 2.19385079 9.62932667 9.87607784
3 7.37703559 1 2.833E-16 7.37703559 2.0899E-15 4.65286266 9.62932667 10.6945343
4 7.14271167 0.8660254 -0.5 6.18576976 -3.5713558 6.71478591 8.43887473 10.7843848
5 -9.2262732 0.5 -0.8660254 -4.6131366 7.990187 5.17707371 11.1022704 12.25

3er Ventana
k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp
0 10.0034232 0 1 0 10.0034232 0 3.33447439 3.33447439
1 1.82424302 0.5 0.8660254 0.91212151 1.5798408 0.3040405 3.86108799 3.87304029
2 7.37703559 0.8660254 0.5 6.38870023 3.68851779 2.43360724 5.09059392 5.64239228
3 7.14271167 1 2.833E-16 7.14271167 2.0235E-15 4.81451114 5.09059392 7.00668709
4 -9.2262732 0.8660254 -0.5 -7.990187 4.61313661 2.15111547 6.62830612 6.96862539
5 -19.31264 0.5 -0.8660254 -9.65632 16.7252368 -1.0676578 12.2033851 12.25
Aplicación TDF Ciclo completo.
¿y sí cambia la frecuencia de las señales de referencia?, ej. n=3 ?

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 9.22627323 0 1
1 19.3126399 1 6.1257E-17
2 10.0034232 1.2251E-16 -1
3 1.82424302 -1 -1.838E-16
4 7.37703559 -2.45E-16 1
5 7.14271167 1 3.0629E-16
6 -9.2262732 3.6754E-16 -1
7 -19.31264 -1 -4.288E-16
8 -10.003423 -4.901E-16 1
9 -1.824243 1 -1.225E-15
10 -7.3770356 6.1257E-16 -1
11 -7.1427117 -1 1.1025E-15
Aplicación TDF Ciclo completo.
¿y sí cambia la frecuencia de las señales de referencia?, ej. n=3
1er Ventana
k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp
0 9.22627323 0 1 0 9.22627323 0 1.5377122 1.5377122
1 19.3126399 1 6.1257E-17 19.3126399 1.183E-15 3.21877332 1.5377122 3.56722028
2 10.0034232 1.2251E-16 -1 1.2256E-15 -10.003423 3.21877332 -0.129525 3.22137834
3 1.82424302 -1 -1.838E-16 -1.824243 -3.352E-16 2.91473282 -0.129525 2.91760931
4 7.37703559 -2.45E-16 1 -1.808E-15 7.37703559 2.91473282 1.09998094 3.11538528
5 7.14271167 1 3.0629E-16 7.14271167 2.1877E-15 4.10518476 1.09998094 4.25
6 -9.2262732 3.6754E-16 -1 -3.391E-15 9.22627323 4.10518476 2.63769315 4.87954578
7 -19.31264 -1 -4.288E-16 19.3126399 8.2813E-15 7.32395808 2.63769315 7.78445805
8 -10.003423 -4.901E-16 1 4.9023E-15 -10.003423 7.32395808 0.97045595 7.38797311
9 -1.824243 1 -1.225E-15 -1.824243 2.2348E-15 7.01991758 0.97045595 7.08667959
10 -7.3770356 6.1257E-16 -1 -4.519E-15 7.37703559 7.01991758 2.19996188 7.3565668
11 -7.1427117 -1 1.1025E-15 7.14271167 -7.875E-15 8.21036952 2.19996188 8.5
Aplicación TDF Ciclo completo.
¿y sí cambia la frecuencia de las señales de referencia?, ej. n=3
2da Ventana
k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp
0 19.3126399 0 1 0 19.3126399 0 3.21877332 3.21877332
1 10.0034232 1 6.1257E-17 10.0034232 6.1278E-16 1.66723719 3.21877332 3.62493883
2 1.82424302 1.2251E-16 -1 2.235E-16 -1.824243 1.66723719 2.91473282 3.35787839
3 7.37703559 -1 -1.838E-16 -7.3770356 -1.356E-15 0.43773126 2.91473282 2.94741854
4 7.14271167 -2.45E-16 1 -1.75E-15 7.14271167 0.43773126 4.10518476 4.1284562
5 -9.2262732 1 3.0629E-16 -9.2262732 -2.826E-15 -1.0999809 4.10518476 4.25
6 -19.31264 3.6754E-16 -1 -7.098E-15 19.3126399 -1.0999809 7.32395808 7.40610019
7 -10.003423 -1 -4.288E-16 10.0034232 4.2895E-15 0.56725625 7.32395808 7.34589284
8 -1.824243 -4.901E-16 1 8.9399E-16 -1.824243 0.56725625 7.01991758 7.04279933
9 -7.3770356 1 -1.225E-15 -7.3770356 9.0372E-15 -0.6622497 7.01991758 7.05108626
10 -7.1427117 6.1257E-16 -1 -4.375E-15 7.14271167 -0.6622497 8.21036952 8.2370348
11 9.22627323 -1 1.1025E-15 -9.2262732 1.0172E-14 -2.1999619 8.21036952 8.5
Aplicación TDF Ciclo completo.
¿y que pasa con: CD + h1 + h3?
=% . . ∗ ! " ∗ # ° + 8.5∗ ! #∗" ∗ °

1er Ventana

k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp

0 3.72627323 0 1 0 3.72627323 0 0.62104554 0.62104554
1 13.8126399 0.5 0.8660254 6.90631996 11.9620971 1.15105333 2.61472838 2.85687386
2 4.50342317 0.8660254 0.5 3.90007887 2.25171158 1.80106647 2.99001365 3.49056185
3 -3.675757 1 2.833E-16 -3.675757 -1.041E-15 1.18844031 2.99001365 3.21754129
4 1.87703559 0.8660254 -0.5 1.6255605 -0.9385178 1.45936706 2.83359401 3.18731976
5 1.64271167 0.5 -0.8660254 0.82135584 -1.42263 1.5962597 2.59648901 3.04791735
6 -14.726273 5.666E-16 -1 -8.344E-15 14.7262732 1.5962597 5.05086788 5.29710405
7 -24.81264 -0.5 -0.8660254 12.40632 21.4883765 3.66397969 8.63226396 9.37767179
8 -15.503423 -0.8660254 -0.5 13.4263583 7.75171158 5.90170608 9.92421589 11.5464365
9 -7.324243 -1 -1.838E-16 7.32424302 1.346E-15 7.12241325 9.92421589 12.215516
10 -12.877036 -0.8660254 0.5 11.1518399 -6.4385178 8.98105324 8.85112959 12.6095921
11 -12.642712 -0.5 0.8660254 6.32135584 -10.948909 10.0346125 7.02631135 12.25
Aplicación TDF Ciclo completo.
¿y que pasa con: CD + h1 + h3?
=% . . ∗ ! " ∗ # ° + 8.5∗ ! #∗" ∗ °

2da Ventana
k Ent RefSen RefCos Ent*RefSen Ent*RefCos Preal Pimag Vp
0 13.8126399 0 1 0 13.8126399 0 2.30210665 2.30210665
1 4.50342317 0.5 0.8660254 2.25171158 3.90007887 0.37528526 2.9521198 2.97587808
2 -3.675757 0.8660254 0.5 -3.1832989 -1.8378785 -0.1552646 2.64580672 2.65035851
3 1.87703559 1 2.833E-16 1.87703559 5.3177E-16 0.15757471 2.64580672 2.65049485
4 1.64271167 0.8660254 -0.5 1.42263004 -0.8213558 0.39467972 2.50891408 2.53976808
5 -14.726273 0.5 -0.8660254 -7.3631366 12.7533267 -0.8325097 4.63446853 4.70864853
6 -24.81264 5.666E-16 -1 -1.406E-14 24.8126399 -0.8325097 8.76990851 8.80933413
7 -15.503423 -0.5 -0.8660254 7.75171158 13.4263583 0.45944221 11.0076349 11.0172189
8 -7.324243 -0.8660254 -0.5 6.34298052 3.66212151 1.51660563 11.6179885 11.7165588
9 -12.877036 -1 -1.838E-16 12.8770356 2.3664E-15 3.66277823 11.6179885 12.1816912
10 -12.642712 -0.8660254 0.5 10.9489095 -6.3213558 5.48759648 10.5644292 11.9046579
11 3.72627323 -0.5 0.8660254 -1.8631366 3.22704728 5.17707371 11.1022704 12.25
Algoritmo TDF. Ciclo completo
• Sí se consideran los caso de la tabla 1, ¿el
? ? ? algoritmo funciona correctamente?
• ¿en qué condiciones no funciona?
Medio Ciclo
• Sí se consideran los caso de la tabla 1, ¿el
? algoritmo funciona correctamente?
Tabla 1. Casos de prueba. • ¿en qué condiciones no funciona?
Proponga magnitudes:
+CD + h1 Para ambos algoritmos:
h1 + h2
• ¿Se puede determinar correctamente el valor
-CD + h1 +h2
+CD + h1 + h2 +h5
pico de h2?
**Al menos tres ventanas • ¿Se puede determinar correctamente el valor
pico de h5?
Pruebe con fm=1440 [Hz]