http://www.gomatlab.de/index.php goMatlab.de - Deutschsprachiges Forum für die Benutzung der Software MATLAB / Simulink. Hilfe bei Bedienungsfragen oder bei der Umsetzung technischer Probleme. de http://www.gomatlab.de/index.php http://www.gomatlab.de/templates/vSpiralFull/images/_header_logo.gif 3 http://www.gomatlab.de/viewtopic.php?p=204690#204690 Tue, 13 Aug 2024 17:46:13 +0200 Ach dank dir Harald. Danke noch mla. Und sehr fleißig was du hier so weg schaffst 3 http://www.gomatlab.de/viewtopic.php?p=204650#204650 Mon, 05 Aug 2024 12:45:11 +0200 Moin, ah, das hatte ich übersehen. Was ist der Output von der Zeile bevor du die Triggereinstellungen vornimmst? triggerinfo(vid2) Ich habe die Vermutung, dass es an dieser Zeile liegt, da "DeviceSpecific" wahrscheinlich keine korrekte Konfiguration ist: triggerconfig(vid2, 'hardware', 'DeviceSpecific', 'DeviceSpecific'); Gruß Steve 0 http://www.gomatlab.de/viewtopic.php?p=204544#204544 Wed, 21 Feb 2024 09:23:51 +0100 Hallo, ich versuche aktuell ein paar Schritte einer Veröffentlichung nachzuvollziehen und komme nicht weiter. Es liegt eine Matrix vor (siehe Anhang). Diese soll nun mittels 2D FFT analysiert und in einem Bubble Plot dargestellt werden. Die Bubblegröße stellt dabei die Ampliduden der Ordnungen dar. Wenn ich die 2D FFT nun auf die Matrix anwende, komme ich zunächst erstmal nicht auf die korrekten Ordnungen (die habe ich zuvor stellenweise geprüft mittels 1D-FFT) und auch nicht auf die korrekten Amplitudenwerte. Könnt ihr mir da helfen? load("mat_signal_force.mat") signal = mat_signal_force; figure() surf(signal', 'EdgeColor','none') xlabel('x - Datenpunkte Ort') ylabel('y - Datenpunkte Zeit') L = size(signal); P2 = abs(fft2(signal))/(L(1)*L(2)); P1 = P2(1:fix(L(1)/2)+1,1:fix(L(2)/2)+1); P1(2:end-1) = 2*P1(2:end-1); [X,Y] = meshgrid(0:L(1)/2, 0:L(2)/2); figure() surf(X,Y,P1') xlim([0 50]) ylim([0 50]) 0 http://www.gomatlab.de/viewtopic.php?p=204419#204419 Thu, 24 Aug 2023 14:15:58 +0200 Hallo an alle, aktuell beschäftige ich mich mit der Auswertung einer abgehängten Decke die per Laservibrometer aufgrund Ihrer Resonanzfrequenzen untersucht wird. Hierbei werden Datensätze im ASCII Format importiert, Fouriertransformiert, energetisch gemittelt und zurück in den Zeitbereich transofrmiert, um hier über die Funktion Findpeaks später den mittleren Dämpfungsgrad zu ermitteln. Nun scheint aber bei der IFFT ein Fehler aufzutreten. Die Zeitachse scheint nach linksverschoben. Sämtliche Lösungsansätze bringen mich nicht weiter und eine mehrmalige Fehlersuche und analyse brachte mich nicht darauf wo der Fehler steckt. Vielleicht kann mir jemand von euch auf die Sprünge helfen. Danke schonmal im vorraus!! Fabi Hier ist der Code %%extrahieren der Geschwindigkeits- und Zeitwerte des ASCII Datensatzes data_G1= data_v1.data(1:end,2); data_T= data_v1.data (1:end,1); %% FFT Berechnung fa = 1/diff(data_T (1:2)); % Abtastfrequenz (1/T) Ts = 1/fa; % Abtastrate fn = fa/2.56 ; % Nyquistfrequenz %Ermitteln der Abtastpunkte für FFT N = length(data_G1); % gewünschte FFT-Länge df = fa/(N-1); % Frequenzauflösung f=(-fa/2:df:fa/2); %Frequenzachse um Nullpunkt zentriert for i = 1:numVectors % Name des Vektors vectorName = sprintf('data_G%d', i); %erzeugen eines Arrays % Lade den Vektor timeVector = eval(vectorName); % Führe die FFT durch frequencyVector = fft(timeVector); frequencyVector = abs(fftshift(frequencyVector)); %frequencyVector=frequencyVector(1:N/2); % Speichere den frequenzbasierten Vektor frequencyVectors(i, :) = frequencyVector; % % Zeige den frequenzbasierten Vektor an % disp(['Frequenzbasierten Vektor ', vectorName, ':']); % disp(frequencyVec ... 1 http://www.gomatlab.de/viewtopic.php?p=204141#204141 Wed, 29 Mar 2023 08:25:00 +0200 Im Anhang befindet sich noch die 2 Files 0 http://www.gomatlab.de/viewtopic.php?p=203620#203620 Mon, 31 Oct 2022 12:03:14 +0100 Hallo, ich habe die Messdaten eines Beschleunigungssensors und möchte daraus die 1. Eigenfrequenz ermitteln und das Signal im Bereich der ersten Eigenfrequenz anschließend wieder in den Zeitbereich transformieren um aus dem geglätteten Beschleuniguns-Zeitschrieb die Dämpfung zu bestimmen (es handelt sich um die Daten einer als Einfeldträger gelagerte Platte). Aktuell mache ich das, indem ich alle Amplituden außerhalb des Frequenzbereiches zu 0 setzt und anschließend eine IFFT durchführe. Allerdings scheint mir diese nicht Richtig, da die Dämpfungswerte zu gering sind und außerdem zu stark abhängig vom Frequenzbereich. Wie wäre das besser zu lösen? Danke! clear; %-------------------------------------------------------------------------- % E I N G A B E W E R T E %-------------------------------------------------------------------------- dateiname = 'Versuch_02_20220928_001.xlsb'; kanal_kraft = 1; kanal_beschleunigung = 2; minf = 3; % minimale Frequenz für Darstellung maxf = 50; % maximale Frequenz für Darstellung %-------------------------------------------------------------------------- % Split Impulshammerschläge t1 = (1:4700); t2 = (4701:8900); t3 = (8901:11800); t4 = (11801:17000); t5 = (17001:21200); t6 = (21201:30247); werte = xlsread('Versuch_02_20220928_001.xlsb'); zeit = werte(t4,1); kraft = werte(t4,1+kanal_kraft); beschleunigung = werte(t4,1+kanal_beschleunigung); fs = round(1/mean(diff(zeit))); % Abtastrate in Hz (Durchschnittliche Aufzeichnungen je Sekunde) n = length(kraft); % Länge verwendete Daten freq = linspace(0,(n-1)/n,n)*fs; % Frequenz F = fft(kraft)/n; A = fft(beschleunigung)/n; close all %Grafische Darstellung der Messdaten %Plot Impulshammerschläge figure plot( ... 1 http://www.gomatlab.de/viewtopic.php?p=202879#202879 Wed, 04 May 2022 14:21:03 +0200 Der UDP Stream wird ja ähnlich wie der Audiostream über ein erstelltes Objekt behandelt, welches abschließend released wird. laufen diese Objekte eigenständig für sich? DAnn wäre mein PRoblem ja fast gelöst, oder? -> Denn die UDP Datenverarbeitung könnte ich ja auch nach der Aufnahme machen, sofern der Buffer groß genug ist. 1 http://www.gomatlab.de/viewtopic.php?p=202577#202577 Tue, 22 Mar 2022 04:48:05 +0100 Wenn das Filter nicht online funktionnieren wird, konntest du ein FIR low pass mit filterDesign probieren (ich siehe etwa 10 Hz). Siehe Mathworks, wenn du das Signal Processing TB hast. 0 http://www.gomatlab.de/viewtopic.php?p=202336#202336 Sat, 05 Feb 2022 21:36:11 +0100 Wie kann man Positionsmessdaten mithilfe eines Kalman Filters filtern? Es handelt sich um eine statische Positionsmessung, das heißt das Objekt bewegt sich nicht, die Messung ist jedoch ungenau, sodass es Abweichungen bei den Messpunkten gibt. Welche Ansätze gibt es einen Kalman filter in Matlab zu implementieren, um diese Messdaten zu filtern? 0 http://www.gomatlab.de/viewtopic.php?p=201382#201382 Thu, 23 Sep 2021 11:50:12 +0200 Hallo zusammen, ich bin noch recht neu in Matlab und möchte folgende Aufgabe lösen: Ich habe eine Ãœbertragungsfunktion eines elektrischen Bauteils bereits gemessen und diese Messwerte liegen in der Form Frequenz--> Rea- und Imaginärteil vor. Ich habe also eine Tabelle, bei der die Frequenz von z.B. 1 kHz bis 10 kHz aufgelistet ist mit insgesamt 1601 Datenpunkten, welche linear zwischen 1 kHz und 10 kHz verteilt sind. Das habe ich mehrere Male, immer über Frequenzdekaden mit wechselnder Anzahl an Datenpunkten (201, 1601 und 2048). Jetzt habe ich drei Spannungssignale gleichzeitig über eine gewisse Zeitspanne mit einem Oszilloskop aufgenommen. Diese Aufnahmen bilden Bereiche von z.B. 50 ns/ div ab, also eine Betrachtungsdauer von insgesamt 500 ns und die Messung umfasst 1250 Werte. Ich habe also einen linearen Abstand zwischen den Messwerten von 400 Picosekunden. Die gemessenen Spannungssignale möchte ich im Frequenzbereich mit der Ãœbertragungsfunktion multiplizieren, um ein Ausgangssignal berechnen zu können. Mein Ziel ist es also, die FFT der Spannungssignale in Matlab so zu gestalten, dass ich einen Frequenzvektor erzeuge, der dem der Ãœbertragungsfunktion angepasst ist und ich das Ergebnis der FFT bei einer bestimmten Frequenz (komplexe Zahl) mit der komplexen Zahl der Ãœbertragungsfunktion bei der identischen Frequenz multiplizieren kann. Zur FFT: Es hakt an der Stelle, an der ich den Frequenzvektor erstellen will. Die FFt gibt mir für jeden Messwert den komplexen Fourier Koeffizienten aus. Nur weiß ich nicht, auf was sich die komplexen Fourier-Koeffizienten beziehen. Mir fehlt die konkrete Frequenz, bzw. ich benötige die exakt die Frequenzen der Ãœbertragungsfunktion für die FFT. Nach der FFT soll jeder Frequenz, die es auch in der Ãœbertragungsfunktion gibt, eine komplexe Zahl zugeordnet werden. Ich hoffe, dass das Problem verständlich erklärt wurde. Ich habe keine Probleme mit dem Einlesen der Daten aus ...