【数据分析】基于matlab实现齿轮箱振动数据分析

MATLAB软件广泛应用于自动控制,信号分析,数理统计等工程领域.对于信号处理分析最为广泛.本文重点介绍齿轮箱信号的分析方法.

function varargout = Mainchilun(varargin)

% MAINCHILUN MATLAB code for Mainchilun.fig

% MAINCHILUN, by itself, creates a new MAINCHILUN or raises the existing

% singleton*.

%

% H = MAINCHILUN returns the handle to a new MAINCHILUN or the handle to

% the existing singleton*.

%

% MAINCHILUN('CALLBACK',hObject,eventData,handles,...) calls the local

% function named CALLBACK in MAINCHILUN.M with the given input arguments.

%

% MAINCHILUN('Property','Value',...) creates a new MAINCHILUN or raises the

% existing singleton*. Starting from the left, property value pairs are

% applied to the GUI before Mainchilun_OpeningFcn gets called. An

% unrecognized property name or invalid value makes property application

% stop. All inputs are passed to Mainchilun_OpeningFcn via varargin.

%

% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one

% instance to run (singleton)".

%

% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help Mainchilun

% Last Modified by GUIDE v2.5 09-Sep-2022 23:12:37

% Begin initialization code - DO NOT EDIT

gui_Singleton = 1;

gui_State = struct('gui_Name', mfilename, ...

'gui_Singleton', gui_Singleton, ...

'gui_OpeningFcn', @Mainchilun_OpeningFcn, ...

'gui_OutputFcn', @Mainchilun_OutputFcn, ...

'gui_LayoutFcn', [] , ...

'gui_Callback', []);

if nargin && ischar(varargin{1})

gui_State.gui_Callback = str2func(varargin{1});

end

if nargout

[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});

else

gui_mainfcn(gui_State, varargin{:});

end

% End initialization code - DO NOT EDIT

% --- Executes just before Mainchilun is made visible.

function Mainchilun_OpeningFcn(hObject, eventdata, handles, varargin)

% This function has no output args, see OutputFcn.

% hObject handle to figure

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

% varargin command line arguments to Mainchilun (see VARARGIN)

% Choose default command line output for Mainchilun

handles.output = hObject;

% Update handles structure

guidata(hObject, handles);

% UIWAIT makes Mainchilun wait for user response (see UIRESUME)

% uiwait(handles.figure1);

% --- Outputs from this function are returned to the command line.

function varargout = Mainchilun_OutputFcn(hObject, eventdata, handles)

% varargout cell array for returning output args (see VARARGOUT);

% hObject handle to figure

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure

varargout{1} = handles.output;

% --- Executes on selection change in listbox1.

function listbox1_Callback(hObject, eventdata, handles)

% hObject handle to listbox1 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

% Hints: contents = cellstr(get(hObject,'String')) returns listbox1 contents as cell array

% contents{get(hObject,'Value')} returns selected item from listbox1

% --- Executes during object creation, after setting all properties.

function listbox1_CreateFcn(hObject, eventdata, handles)

% hObject handle to listbox1 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles empty - handles not created until after all CreateFcns called

% Hint: listbox controls usually have a white background on Windows.

% See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

set(hObject,'BackgroundColor','white');

end

% --- Executes on button press in pushbutton1.

function pushbutton1_Callback(hObject, eventdata, handles)

% hObject handle to pushbutton1 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

axes(handles.axes1) ;

cla reset

fs=2000*2.56; % 采样率

load data.mat

t=0:1/fs:(length(signal(1,:))-1)/fs; %以时间为横坐标

signal_type = get(handles.signal,'Value') ;

switch signal_type

case 1

signaldo=signal(1,:);

case 2

signaldo=signal(2,:);

case 3

signaldo=signal(3,:);

case 4

signaldo=signal(4,:);

end

signalnew=wden(signaldo,'minimaxi','s','mln',5,'sym5'); %小波分解去噪

axes(handles.axes1);

plot(t,signaldo,'-b');hold on;

plot(t,signalnew,'-r');

xlim([0,11]);ylabel('加速度');xlabel('时间（s）');

legend('去噪前','去噪后');

feature(1)=max(signalnew); %特征1 max

feature(2)=min(signalnew); %特征2 min

feature(3)=mean(signalnew); %特征3 mean

feature(4)=var(signalnew); %特征4 方差

feature(5)=std(signalnew); %特征5 均方根

set(handles.edit1,'String',num2str(feature(1)));

set(handles.edit2,'String',num2str(feature(2)));

set(handles.edit3,'String',num2str(feature(3)));

set(handles.edit4,'String',num2str(feature(4)));

set(handles.edit5,'String',num2str(feature(5)));

% --- Executes on button press in pushbutton2.

function pushbutton2_Callback(hObject, eventdata, handles)

% hObject handle to pushbutton2 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

% --- Executes on button press in pushbutton3.

function pushbutton3_Callback(hObject, eventdata, handles)

% hObject handle to pushbutton3 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

axes(handles.axes1) ;

cla reset

fs=2000*2.56; % 采样率

load data.mat

t=0:1/fs:(length(signal(1,:))-1)/fs; %以时间为横坐标

signal_type = get(handles.signal,'Value') ;

switch signal_type

case 1

signaldo=signal(1,:);

case 2

signaldo=signal(2,:);

case 3

signaldo=signal(3,:);

case 4

signaldo=signal(4,:);

end

signalnew=wden(signaldo,'minimaxi','s','mln',5,'sym5'); %小波分解去噪

signalnew=signalnew./max(abs(signalnew)); %归一化

axes(handles.axes1);

plot(t,signalnew,'-k');

xlim([0,11]);ylabel('加速度');xlabel('时间（s）');

feature(1)=max(signalnew); %特征1 max

feature(2)=min(signalnew); %特征2 min

feature(3)=mean(signalnew); %特征3 mean

feature(4)=var(signalnew); %特征4 方差

feature(5)=std(signalnew); %特征5 均方根

set(handles.edit1,'String',num2str(feature(1)));

set(handles.edit2,'String',num2str(feature(2)));

set(handles.edit3,'String',num2str(feature(3)));

set(handles.edit4,'String',num2str(feature(4)));

set(handles.edit5,'String',num2str(feature(5)));

function edit1_Callback(hObject, eventdata, handles)

% hObject handle to edit1 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of edit1 as text

% str2double(get(hObject,'String')) returns contents of edit1 as a double

% --- Executes during object creation, after setting all properties.

function edit1_CreateFcn(hObject, eventdata, handles)

% hObject handle to edit1 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.

% See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

set(hObject,'BackgroundColor','white');

end

function edit2_Callback(hObject, eventdata, handles)

% hObject handle to edit2 (see GCBO)

% eventdata reserved - to be defined in a future version of MATLAB

% handles structure with handles and user data (see GUIDATA)

[1]刘俊, 杨俊华, 杨梦丽,等. 基于MATLAB与VC的风电机组齿轮箱振动信号采集与分析系统[J]. 广东电力, 2013, 26(6):7.