Frequenzmessung periodisches Signal?

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

zapman

Forum-Newbie


	Beiträge: 2

	Anmeldedatum: 18.04.08

	Wohnort: London

	Version: ---

Verfasst am: 21.05.2008, 11:54 Titel: Frequenzmessung periodisches Signal?

Hey!

Ich bin mit Simulink noch nicht all zu vertraut und habe folgendes Problem:

Wie bekomme ich eine Frequenzmessung oder Frequenzschätzung eines periodischen Signals (z.B. leicht verrauschter Sinus) hin?

Denke da an ZERO CROSSING DETECTION, um die Nulldurchgänge des Signals zu zählen und eventuell mit einem unverrauschten Referenzsignal zu vergleichen.

Kenne aber die benötigten Blöcke nicht genau bzw. deren Kombination.

Das Signal kann durch Tiefpass zweiter Ordnung gefiltert werden, andere Lösungen nehme ich aber auch gerne an! Wink

Besten Dank schonmal im Vorraus!

nschlange

Ehrenmitglied


	Beiträge: 1.320

	Anmeldedatum: 06.09.07

	Wohnort: NRW

	Version: R2007b

Verfasst am: 21.05.2008, 14:35 Titel:

Hi,

das klassische Werkzeug dafür ist eigentlich die Fourier-Transformation.
Dazu gibt es sicherlich einen Simulink-Block.
_________________

Viele Grüße
nschlange

"Chuck Norris ejakuliert fluessigen Stahl!"

Janvi

Forum-Anfänger


	Beiträge: 28

	Anmeldedatum: 20.05.08

	Wohnort: Rom

	Version: R2007a

Verfasst am: 04.06.2008, 09:55 Titel:

das klassische Werkzeug zur Frequenzmessung bei unbekannten Kurven ist die Autokorrelation. Ab 64 Punkten ist der Rechenaufwand über die fft nach Cooley/Tukey schneller weshalb i.d.R. nicht im Zeit sondern im Frequenzbereich gerechnet wird. D.h. du must zuerst vorwärts dann wieder zurücktransformieren. Die Multiplikation im Frequenzbereich entspricht dann der Faltung im Zeitbereich. Als Skript etwa so:

Code:

akf = abs(ifft(abs(fft(ripple)).^2)); % Fouriertransformation berechnen
akf = akf/1.0e+4; % zum Vergleich mit Eingangssignal skalieren
akf = minus(akf,min(akf)); % und Offset entfernen
ddf = diff(sign(diff(akf))); % Differential zeigt alle Extrempunkte
ddf(ddf==2)=0; % Tiefpunkte sind uninteressant
akf(ddf==0)=0; % alle anderen Werte auch
akf(n/2:n)=0; % zweite Haelfte sowieso. Moegliche Perioden
plot(akf,'r'); % mit verbleibenden Hochpunkten anzeigen
[val,i]=max(akf); % hoechster Wert ist Periode
plot(i,val,'bs','Markersize',20);
barh(val,i);

Funktion ohne Link?

Die Rechnung steht in der ersten Zeile, sodann kommt nur noch die Auswertung. Bestechender Vorteil der Methode ist, daß es unabhängig von der Kurvenform arbeitet und der Sinus auch um den Nullpunkt herum rauschen darf.


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.