clc
clear
close all
a = dir( [ 'Bearing1_3/' '*.csv']);
val = xlsread(['Bearing1_3\' a(2200).name]);
figure
plot(val(:,5),'k-')
hold on
plot(val(:,6),'r-')
hold on
legend('水平信号','垂直信号')
ylabel('振幅')
signal=(val (:,5));
%----------------- Initialization
maxAlpha=20000; %compactness of mode
tau=0;%time-step of the dual ascent
tol=1e-6; %tolerance of convergence criterion;
stopc=4;%the type of stopping criteria
%-------------- SVMD function
[u,uhat]=svmd(signal,maxAlpha,tau,tol,stopc);
figure
for ii =1:size(u,1)
if ii ==1
title('svmd分解')
end
subplot(size(u,1),1,ii)
plot(u(ii,:),'b-')
end
figure
plot(signal)
hold on
plot(sum(u))
legend('原始信号','重构信号')
figure
plot(abs(uhat))
m= 2;
t = 2;
Scale = 10;
x = sum(u);
MPE = MPerm(x',m,t,Scale)