Dämpfung von Signal - Einhüllende berechnen - Mein MATLAB Forum

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Dämpfung von Signal - Einhüllende berechnen

Gast09

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 03.06.2012, 13:29 Titel: Dämpfung von Signal - Einhüllende berechnen

Hallo, ich habe folgendes Problem: Ich habe ein Signal, wie in meinem Anhang zu sehen (als Bild und Datenpunkte). Ich möchte nun die Einhüllende bilden, um das Dämpfungsverhalten über das logarithmische Dekrement zu bestimmen. (Habe mehrere Signale, welche ich vergleichen muss.) Leider habe ich bisher weder von MATLAB noch von Signalverarbeitung allgemein besonders viel Ahnung, brauche diese Auswertung aber für eine Arbeit, die ich schreiben muss. Kann mir jemand sagen, wie ich da vorgehen muss? Danke!

106a.xls

Beschreibung:

Download

Dateiname:

106a.xls

Dateigröße:

38.5 KB

Heruntergeladen:

1164 mal

DSP

Forum-Meister


	Beiträge: 2.117

	Anmeldedatum: 28.02.11

	Wohnort: ---

	Version: R2014b

Verfasst am: 04.06.2012, 11:24 Titel:

Damit kannst du dir den Dämpfungsverlauf deines Signals ansehen. Ich habe noch deine Datenvektoren t und y in einen .mat file gepackt um sie in matlab verwenden zu können. Handelt es sich hier um die Impulsantwort eines Systems? Bist du dir auch sicher, dass das Abtasttheorem eingehalten worden ist? Fa = 500 Hz scheint mir etwas wenig.

Code:

% Zeitbereich
% ----------------------------------
load data_y_t;

fa = 1/diff(t(1:2)); % Abtastfrequenz
fn = fa/2; % Nyquistfrequenz
N = length(y); % gewünschte FFT-Länge (N=2^x, sonst wird der
DFT-Algorithmus verwendet!)
% Momentan sind es 276 Messwerte und somit wird die DFT verwendet.
% Du kannst N auch auf 512 Werte erhöhen, wodurch die Freq.-Auflösung besser wird. Die FFT() hängt dann automatisch Nullen an y vor der Transformation an um auf 512 Messwerte zu kommen. Das ändert aber nichts am Ergebnis.
df = fa/N; % Frequenzauflösung

% Graphische Darstellung
% ----------------------
% max. Amplitude zur Skalierung der graphischen Darstellung feststellen:
max_y = max(abs(y))*1.1;
fig = figure(1);
plot(t,y,'.-')
%axis([0 N -max_y max_y])
title('Datensatz')
ylabel('Amplitude')
xlabel('Time')
grid

% Frequenzbereich
% ----------------------------------

% Berechnung der FFT
% ------------------
H = fft(y, N);
% Berechnung des Amplitudengangs aus dem komplexen Frequenzvektor H:
amplH = abs(H);
% Graphische Darstellung
% ----------------------
% Frequenzvektoren (werden bei der graphischen Darstellung benötigt):
x_fn = 0 : df : fn-df;

% Darstellung des interessierenden Frequenzbereichs des
% Amplitudengangs (0...fn) und
% daran angepasste Amplitudenskalierung (Normierung auf N/2):
amplitudengang = [amplH(1)/N; amplH(2:N/2)/(N/2)]; % DC-Bin auf N normieren!
fig = figure(fig+1);
plot(x_fn, amplitudengang, 'b.-')
%axis([0 fn 0 a*1.1])
title('Amplitudengang')
ylabel('Amplitude')
xlabel(['Auflösung: ',num2str(df),' Hz Frequenz in Hz'])
grid

% Ausgabe in dB
% ------------------
fig = figure(fig+1);
plot(x_fn, 20*log10(amplitudengang))
%axis([0 fn -100 20*log10(a)+3])
title('Amplitudengang')
ylabel('Amplitude in dB')
xlabel(['Auflösung: ',num2str(df),' Hz Frequenz in Hz'])
grid

Funktion ohne Link?

Willst du auch den Phasenverlauf sehen?

data_y_t.zip

Beschreibung:

In Matlab work Verzeichnis laden

Download

Dateiname:

data_y_t.zip

Dateigröße:

2.57 KB

Heruntergeladen:

820 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.