Frage zu DSP. Welche Methode ist die beste?

Wie wär´s mit Görtzel oder ich würde es mal mit einigen Biquads 
hintereinander versuchen. Eine Güte von ca. 100 sollte doch mit 2-3 
Filtern 2. Ordnung zuverlässig zu machen sein. Die Frage wäre, wie 
schnell die Erkennung sein soll. Die Einschwingzeit ist eben auch von 
der Bandbreite abhängig. Etwas mehr als 16 Bit sollten es aber schon 
sein, sonst fällt das Signal irgendwann hinten runter.

Goertzel sagt aber wenig über die Amplitude aus.
Besser ist, was man analogtechnisch als Heruntermischen auf 
Zwischenfrequenz (und eventuell wieder hochmischen) ausführt.
Digital entspricht das dem Multiplizieren mit einem digitalen "local 
oscillator", der nahe neben der gesuchten Frequenz liegt. Mit I/Q- 
Mischung kann man auch gleich die Spiegelfrequenz unterdrücken.

Moin,

Christoph K. schrieb:
> Goertzel sagt aber wenig über die Amplitude aus.

Mit bisschen nachrechnen sagt der ziemlich viel ueber die Leistung aus - 
damit auch ueber die Amplitude. Guggstu hier:

https://en.wikipedia.org/wiki/Goertzel_algorithm#Power-spectrum_terms

Wenn die Trennschaerfe/"Nebenzipfeldaempfung" vom Goertzel reicht, ist 
das das Verfahren der Wahl. Erst wenn der Analysefilterfrequenzgang 
genauer (z.B. rechteckiger) sein muss, dann wuerd' ich mit 
'runtermischen und tiefpassfiltern anfangen.

Gruss
WK

Christoph K. schrieb:
> Goertzel sagt aber wenig über die Amplitude aus.
> Besser ist, was man analogtechnisch als Heruntermischen auf
> Zwischenfrequenz (und eventuell wieder hochmischen) ausführt.
> Digital entspricht das dem Multiplizieren mit einem digitalen "local
> oscillator", der nahe neben der gesuchten Frequenz liegt. Mit I/Q-
> Mischung kann man auch gleich die Spiegelfrequenz unterdrücken.

Aus reiner Neugier: wie würdest du denn den digitalen LO realisieren? 
mit DDS? wenn man eine einigermassen gute Signalqualität mit DDS 
erzielen will, muss das ganze ja ordentlich flott laufen.

Wie schnell muss denn die Erkennung sein ?

Görzel halte Ich hier nicht für zielführend. Ein FIR definierter Tiefe 
sollte die Anforderung erfüllen. Die Bandbreite lässt sich ja darüber 
leicht festlegen. Beim Beispiel oben mit schon gegebener Abtastrate 
nimmt man einfach eine entsprechende Anzahl an Perioden. Mit dem Fenster 
muss man ein wenig aufpassen.

Hi,

mit nem IIR Filter als elliptischen Bandpass geht das ganz gut. Anbei 
der Matlab Entwurf für die gewünschten Bandbreiten, 8.Ordnung, 3dB 
ripple und 60dB Dämpfung.Das Bild zeigt den Frequenzgang und die 
Impulsantwort. Das Filter schwingt 200ms oder so nach, schnell ist das 
nicht. Ausserdem sind die Beträge der Pole so nahe an eins, dass man 
eventuell mit 32Bit float Rechengenauigkeit nicht mehr hinkommt, sodass 
man auf 64 Bit upgraden muss. Das gilt auch für eine integer 
Implementation.

Cheers
Detlef


clear
fs=30000;
fsh=fs/2;
wp=[(5500-50)/fsh (5500+50)/fsh];
n=4;
rp= 3;
rs=60;
f=5000:6000;
[fb,fa]=ellip(n,rp,rs,wp);
[H f]=freqz(fb,fa,f,fs);
subplot(2,1,1)
plot(f,20*log10(abs(H)));
grid
subplot(2,1,2)
plot(filter(fb,fa,[1 zeros(1,10000)]))
return