VMD确定分解层数的程序.zip_VMD的分解层数_naturally6y7_vmd 分解层数_vmd分解层数_vmd分解程序

%%机械故障诊断 vemd clear,clc,close all N = 1200; fs=12000; n = 0:N-1; t=n/fs; num = 50; %% VMD分解 % some sample parameters for VMD alpha = 1000; % moderate bandwidth constraint tau = 0; % noise-tolerance (no strict fidelity enforcement) K = 5; % 4 modes DC = 0; % no DC part imposed init = 0; % initialize omegas uniformly tol = 1e-8; %--------------- Run actual VMD code mmm = 30; % rrr=8; %%正常数据,HP=2,rmp=1750 load Normal2.mat; for c1 = 1:num normal_data(:,c1) = X099_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(normal_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); 以前用SE的程序 e_normal(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_normal(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_normal(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_normal(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_normal(5,c1) = gqiyi_shang(imf5,mmm); e_normalt(c1)=1; end %% 12k Drive End Bearing Fault Data IR0.007 HP=2,rmp=1750 load IR007_2.mat; for c1 = 1:num ir007_data(:,c1) = X107_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(ir007_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_ir007(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_ir007(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_ir007(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_ir007(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_ir007(5,c1) = gqiyi_shang(imf5,mmm); e_ir007t(c1)=2; end %% 12k Drive End Bearing Fault Data IR0.014 HP=2,rmp=1750 load IR014_2.mat; for c1 = 1:num ir014_data(:,c1) = X171_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(ir014_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_ir014(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_ir014(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_ir014(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_ir014(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_ir014(5,c1) = gqiyi_shang(imf5,mmm); e_ir014t(c1)=3; end %% 12k Drive End Bearing Fault Data IR0.021 HP=2,rmp=1750 load IR021_2.mat; for c1 = 1:num ir021_data(:,c1) = X211_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(ir021_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_ir021(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_ir021(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_ir021(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_ir021(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_ir021(5,c1) = gqiyi_shang(imf5,mmm); e_ir021t(c1)=4; end %% 12k Drive End Bearing Fault Data B0.007 HP=2,rmp=1750 load B007_2.mat; for c1 = 1:num b007_data(:,c1) = X120_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(b007_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_b007(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_b007(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_b007(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_b007(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_b007(5,c1) = gqiyi_shang(imf5,mmm); e_b007t(c1)=5; end %% 12k Drive End Bearing Fault Data B0.014 HP=2,rmp=1750 load B014_2.mat; for c1 = 1:num b014_data(:,c1) = X187_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(b014_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_b014(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_b014(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_b014(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_b014(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_b014(5,c1) = gqiyi_shang(imf5,mmm); e_b014t(c1)=6; end %% 12k Drive End Bearing Fault Data B0.021 HP=2,rmp=1750 load B021_2.mat; for c1 = 1:num b021_data(:,c1) = X224_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(b021_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_b021(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_b021(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_b021(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_b021(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_b021(5,c1) = gqiyi_shang(imf5,mmm); e_b021t(c1)=7; end %% 12k Drive End Bearing Fault Data OR0.007@6 HP=2,rmp=1750 load OR007@6_2.mat; for c1 = 1:num or007_data(:,c1) = X132_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(or007_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_or007(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_or007(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_or007(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_or007(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_or007(5,c1) = gqiyi_shang(imf5,mmm); e_or007t(c1)=8; end %% 12k Drive End Bearing Fault Data OR0.014@6 HP=2,rmp=1750 load OR014@6_2.mat; for c1= 1:num or014_data(:,c1) = X199_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(or014_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_or014(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_or014(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_or014(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_or014(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_or014(5,c1) = gqiyi_shang(imf5,mmm); e_or014t(c1)=9; end %% 12k Drive End Bearing Fault Data OR0.021@6 HP=2,rmp=1750 load OR021@6_2.mat; for c1 = 1:num or021_data(:,c1) = X236_DE_time(c1*N:c1*N+N-1); [u, u_hat, omega] = VMD(or021_data(:,c1), alpha, tau, K, DC, init, tol); imf= u'; imf1 = imf(:,1); % r = 0.2*std(imf1); e_or021(1,c1) = gqiyi_shang(imf1,mmm); imf2 = imf(:,2); % r = 0.2*std(imf2); e_or021(2,c1) = gqiyi_shang(imf2,mmm); imf3 = imf(:,3); % r = 0.2*std(imf3); e_or021(3,c1) = gqiyi_shang(imf3,mmm); imf4 = imf(:,4); % r = 0.2*std(imf4); e_or021(4,c1) = gqiyi_shang(imf4,mmm); imf5 = imf(:,5); % r = 0.2*std(imf5); e_or021(5,c1) = gqiyi_shang(imf5,mmm); e_or021t(c1)=10; end ee=[e_normal';e_ir007';e_ir014';e_ir021';e_b007';e_b014';e_b021';e_or007';e_or014';e_or021']; t