clc
clear?all
close?all
tic
load?130.mat;%讀取外圈故障數(shù)據(jù)，這個(gè)數(shù)據(jù)在VMD文件夾給你了
toc
tic
x=X130_DE_time（1:4096）;%從讀取的數(shù)據(jù)中取X130_DE_time數(shù)據(jù)作為要分解的信號(hào)
fs=12000;%采樣頻率，這個(gè)是信號(hào)處理里面要設(shè)定信號(hào)采集時(shí)候的頻率，是由采集信號(hào)傳感器決定，你不用管
N=4096;%采樣點(diǎn)數(shù)，這是采集點(diǎn)數(shù)，也就是信號(hào)的長(zhǎng)短
n=0:N-1;
f=n*fs/N;
T=1000;
n=0:N-1;
t=n/fs;

%?some?sample?parameters?for?VMD
%這個(gè)是我要優(yōu)化的VMD信號(hào)處理的算法，包含6個(gè)參數(shù)，其中alpha和K是要優(yōu)化的，其他都不用管，取默認(rèn)值
alpha?=?4000;????????%這叫懲罰因子，可以取100~2000以?xún)?nèi)
tau?=?0;????????????%?noise-tolerance?（no?strict?fidelity?enforcement）
K?=?5;??????????????%?這個(gè)是信號(hào)分解的個(gè)數(shù)，設(shè)定多少就分解成多少個(gè)
DC?=?0;?????????????%?no?DC?part?imposed
init?=?1;???????????%?initialize?omegas?uniformly
tol?=?1e-5;
%上面的參數(shù)設(shè)定好后再進(jìn)行下面的信號(hào)分解，設(shè)定K=8就是把信號(hào)分解成8個(gè)，如果alpha設(shè)定不同，也會(huì)影響分解的8個(gè)分量的頻率不一樣，所以這兩個(gè)參數(shù)都會(huì)影響分解結(jié)果
%我們的思路是當(dāng)alpha取100~2000，K取1~10，為了減少計(jì)算量，把a(bǔ)lpha取100200300，整百間隔的跳，也就是20個(gè)alpha值和10個(gè)K值組合

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件????2928192??2010-03-08?15:43??VMDtest\130.mat

?????文件???????2968??2018-03-06?12:53??VMDtest\main.m

?????文件???????4645??2017-03-06?21:32??VMDtest\VMD.m

?????文件????????852??2017-12-22?15:11??VMDtest\yyshang.m

?????目錄??????????0??2018-03-07?19:54??VMDtest

-----------?---------??----------?-----??----

??????????????2936657????????????????????5