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代码问题可私信交流。
部分理论引用网络文献,若有侵权联系博主删除。