Goertzel Algorithmus (Frequenz + Phase)

[TiGü]

Zitat von stoxx:

och man .. ich mach das nicht jeden Tag

Was auch immer das für eine FFT ist.

Kann ich dich hauen?

Tut mir leid, aber kann es sein das dir ein nicht unerheblicher Teil der mathematischen Grundlagen fehlen?

Zitat von stoxx:

hast du, mein letzter Wunsch für heute ..

nochmal 10 und 20 Hz?
mit etwa gleich hoher Amplitude?
auch logarithmisch. sieht das bild auch so gut aus?

Machts dir selbst (ja, Scilab ist Freeware, installiere es endlich):

Delphi-Quellcode:

clf()
clc

sample_rate = 20000;
t = 0 : 1 / sample_rate : 1;

a1 = 5;
f1 = 10;
w1 = 2 * %pi * f1;
y1 = a1 * sin(w1 * t);

a2 = 500;
f2 = 20;
w2 = 2 * %pi * f2;
y2 = a2 * sin(w2 * t);

y3= y1 + y2;

N = size(t,'*')
f = sample_rate * (0 : N / 2) / N;
n = size(f, '*');

spekt = fft(y3);

playsnd(y3, sample_rate)

subplot(211)
plot(t, y3)

subplot(212)
plot2d("nl",f , (abs(2 / N * spekt(1 : n))))