Nyquistgrenzfrequenz bestimmen

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

dabauer82

Forum-Century


	Beiträge: 184

	Anmeldedatum: 21.03.08

	Wohnort: ---

	Version: R2008a, R2010b

Verfasst am: 09.02.2011, 16:35 Titel: Nyquistgrenzfrequenz bestimmen

Hallo zusammen,

ich möchte einen Spagat wagen und mein Wissen aus der Regelungstechnik mit der fertigungstechnik verbinden Smile

Und zwar liegt mir ein optisch vermessener Rauhigkeitschrieb einer Oberfläche vor. Wenn ich diesen exportiere, habe ich die Möglichkeit, jeden Pixel oder nur jeden x-ten Pixel zu exportieren. Jeder Pixel bedeutet aber 2 mio Punkte, so dass ich an Rechnerkapazitäten stoßen. Meine Idee war daher, mittels dem Nyquist-Grenzfrequenzkriterium (http://de.wikipedia.org/wiki/Nyquist-Frequenz) jene Abtastfrequenz zu finden, für die gearde noch eine ausreichend genaue Abtastung und Rekonstruktion des Ausgangssignals möglich ist.

Diese Herangehensweise soll erst mal unabhängig von der Tatsache sein, dass diese Nyquist-Frequenz kein hinreichendes Kriterium für die Gültigkeit des reproduzierten Signals ist.

Mein Ansatz basiert stark auf einem Matlab-Beispiel - ich hab einfach noch nicht genug Ahnung Smile

Mittels FFT soll man angeblich aus dem Plot die Frequenz ablesen können. Allerdings bekomme ich kein so schönes Ergebnis wie in der Matlab-Hilfe gezeigt. Ich könnte nun vermuten, dass es aufgrund eines Rauschens oder fehlender Dämpfung zu wirr ist. Daher habe ich versucht den Befehl resample zu nutzen, aber dessen Realisierung liegt mir noch fern.

Bin ich überhaupt auf dem richtigen Weg oder ist eine andere Alternative sinnvoller?

Code:

%% Einlesen der Dateien
a = load('test_5x5.asc');
xa = a(:,1);
ya = a(:,2);
za = a(:,3);

la = length(a);

% Bestimmung der Frequenz eines Siganls mit FFT zwecks Nyquist
Fs = 10000;
NFFT = 2^nextpow2(la); % Next power of 2 from length of y
Za = fft(za,NFFT)/lb;
f = Fs/2*linspace(0,1,NFFT/2+1);
plot(f,2*abs(Za(1:NFFT/2+1))) ;
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

Funktion ohne Link?

test_5x5.zip

Beschreibung:

Download

Dateiname:

test_5x5.zip

Dateigröße:

231.74 KB

Heruntergeladen:

705 mal

Abtastfrequenz.png

Beschreibung:

Download

Dateiname:

Abtastfrequenz.png

Dateigröße:

3.44 KB

Heruntergeladen:

647 mal

Rauigkeitsprofil.png

Beschreibung:

Download

Dateiname:

Rauigkeitsprofil.png

Dateigröße:

2.28 KB

Heruntergeladen:

616 mal


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.