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

1 内容介绍

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

2 部分代码

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)

3 运行结果

4 参考文献

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

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机、雷达通信、无线传感器等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。