基于Lipschitz李式指数的随机信号特征识别和故障检测matlab仿真

1.程序功能描述

2.测试软件版本以及运行结果展示

3.核心程序

4.本算法原理

4.1 Lipschitz李式指数定义与性质

4.2 Lipschitz李式指数的估计

4.3 Lipschitz李式指数在信号特征识别与故障检测中的应用

5.完整程序

1.程序功能描述

基于Lipschitz李式指数的随机信号特征识别和故障检测.

2.测试软件版本以及运行结果展示

MATLAB2013B版本运行

3.核心程序

% --- 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)
global Signal
global Lipschitz

%    ߶ С   任
Ss      = 1;   
Ls      = 32;
wt      = cwt(Signal,Ss:Ls,'gaus2');
%ģ     
Max_Pos = func_WMMT(wt,10);
%  С   仯   ÿһ   Ķ Ӧ     ֵ    
[Max_List,Ptr,Lengths] = func_find_Max_nlevel(Max_Pos,10);
%    Lipschitzָ   
Lipschitz = func_cal_Lipschitz(Ptr,Max_List,wt,Lengths); 


axes(handles.axes2);
for k=1:length(Ptr)
	vec = Max_List(:,Ptr(k):Ptr(k) + Lengths(k)-1);
    plot(vec(2,:),log2(vec(1,:)));
    hold on;
end
title('ģ        ');
xlabel('u');
ylabel('log2');
hold off
axes(handles.axes3);
plot(Lipschitz,'r');
axis([0,length(Lipschitz),0,10]);
title('Lipschitzָ  ');
hold off
% --- 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)
%        ָ   ļ      Թ    źż  
global Signal
global Lipschitz

cnt = 0;
PP  = [];
for i = 1:length(Lipschitz)
    if abs(Lipschitz(i)) >= 4
       cnt = cnt + 1; 
       PP(cnt) = i;
    end
end
Len = length(Signal);
axes(handles.axes4);
if isempty(PP) == 0
    START = min(PP);
    ENDS  = max(PP);
    plot(Signal);title('   ϼ  ');
    hold on
    plot(14*START,-4:0.01:4,'r');
    hold on
    plot(Len,-4:0.01:4,'r');
else
    plot(Signal);title(' ޹   ');    
end
hold off
% --- Executes on button press in pushbutton4.
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton4 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
clc;
clear;
close all;


% --- Executes on button press in pushbutton6.
function pushbutton6_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton6 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global Signal
global Lipschitz

KHI = func_KHI(Lipschitz);
axes(handles.axes1);
plot(KHI,'b-o');
title('KHIָ  ');

FGP1 = func_FGP1(Lipschitz);
axes(handles.axes3);
plot(FGP1,'b-o');
title('FGP1ָ  ');


FGP2 = func_FGP2(Lipschitz);
axes(handles.axes2);
plot(FGP2,'b-o');
title('FGP2ָ  ');


FGP3 = func_FGP3(FGP2,0.5);
axes(handles.axes4);
plot(FGP3,'b-o');
title('FGP3ָ  ');
16_017m