Verfasst am: 13.04.2016, 15:05
Titel: Maximale Amplitude der Schwingung
Hallo liebe Leute,
ich hoffe ihr könnt mir weiterhelfen. Es geht hierbei um Messdatenausertung für meine Bachelorarbeit und ich wäre euch über jede Hilfe dankbar.
Den Plot der Messdaten könnt ihr auf den Fotos im Anhang sehen. Ausgewertet wird der Bereich mit Stick-Slip-Effekt.
Gesucht ist in diesem Bereich die maximale Amplitude. Mein Betreuer hat vorgeschlagen eine umhüllende Funktion (besser eine Ausgleichskurve) zu erzeugen und so aus dem maximalen vertikalen Abstand die maximale Amplitude zu bestimmen.
Hintergrund: Die Aufgabe ist die auftretenden Reibschwingungen (Stick-Slip) bei der Schraubenmontage zu untersuchen. Die Schrauben werden mit unterschiedlicher Drehzahl maschinell angezogen. Für die Auswertung ist zunächst die maximale Amplitude in diesem Bereich ausschlaggebend.
Gibt es eine Möglichkeit eine Ausgleichskurve (oder Umhüllende) darzustellen und daraus den maximalen Abstand (Amplitude) zu bestimmen?
Code für Plot
Code:
%% Import data from text file. % Script for importing data from the following text file:
%
% H:\Testmessung 1.BIN
%
% To extend the code to different selected data or a different text file, % generate a function instead of a script.
%% Read columns of data as strings: % For more information, see the TEXTSCAN documentation.
formatSpec = '%s%s%s%s%s%s%s%[^\n\r]';
%% Open the text file.
fileID = fopen(filename,'r');
%% Read columns of data according to format string. % This call is based on the structure of the file used to generate this % code. If an error occurs for a different file, try regenerating the code % from the Import Tool.
dataArray = textscan(fileID, formatSpec, 'Delimiter', delimiter, 'ReturnOnError', false);
%% Convert the contents of columns containing numeric strings to numbers. % Replace non-numeric strings with NaN.
raw = repmat({''},length(dataArray{1}),length(dataArray)-1);
for col=1:length(dataArray)-1
raw(1:length(dataArray{col}),col) = dataArray{col};
end
numericData = NaN(size(dataArray{1},1),size(dataArray,2));
for col=[1,2,3,4,5,6,7] % Converts strings in the input cell array to numbers. Replaced non-numeric % strings with NaN.
rawData = dataArray{col};
for row=1:size(rawData, 1);
% Create a regular expression to detect and remove non-numeric prefixes and % suffixes.
regexstr = '(?<prefix>.*?)(?<numbers>([-]*(\d+[\.]*)+[\,]{0,1}\d*[eEdD]{0,1}[-+]*\d*[i]{0,1})|([-]*(\d+[\.]*)*[\,]{1,1}\d+[eEdD]{0,1}[-+]*\d*[i]{0,1}))(?<suffix>.*)';
try
result = regexp(rawData{row}, regexstr, 'names');
numbers = result.numbers;
% Detected commas in non-thousand locations.
invalidThousandsSeparator = false;
ifany(numbers=='.');
thousandsRegExp = '^\d+?(\.\d{3})*\,{0,1}\d*$';
ifisempty(regexp(thousandsRegExp, '.', 'once'));
numbers = NaN;
invalidThousandsSeparator = true;
end end % Convert numeric strings to numbers. if ~invalidThousandsSeparator;
numbers = strrep(numbers, '.', '');
numbers = strrep(numbers, ',', '.');
numbers = textscan(numbers, '%f');
numericData(row, col) = numbers{1};
raw{row, col} = numbers{1};
end catch me
end end end
%% Replace non-numeric cells with NaN
R = cellfun(@(x) ~isnumeric(x) && ~islogical(x),raw); % Find non-numeric cells
raw(R) = {NaN}; % Replace non-numeric cells
%% Clear temporary variables clearvars filename delimiter formatSpec fileID dataArray ans raw col numericData rawData row regexstr result numbers invalidThousandsSeparator thousandsRegExp me R;
%% Plot plot(Time,Mges,'r') hold on
grid on
plot(Time,Mkopf,'g')
%plot(Time,Mgew,'b') % plot(Time,Fax,'k')
wenn du Hilfe erwartest, solltest du uns auch mitteilen was du versucht hast und was nicht funktioniert hat. Sprich ausführbaren Code inkl. Beispieldaten und evtl. Fehlermeldungen bzw. warum das Ergebnis nicht den Erwartungen entspricht.
Gruß DSP
Einstellungen und Berechtigungen
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
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.
RSS