clear?all;
close?all;
clc;
%?給定參數
fs=1000000;
N=2048;
n=0:N-1;
t=n/fs;
number=fix（log2（N））-1;
x=xlsread（‘870.xlsx‘）;
data=x（:1）;
%繪制實測信號時域波形
figure（1）;
plot（tdata）;xlabel（‘時間?t/s‘）;ylabel（‘幅值/mV‘）;
set（gca‘xlim‘[0?2.047*10^-3]）;?
%?some?sample?parameters?for?VMD
alpha?=?2000;???????%?moderate?bandwidth?constraint
tau?=?0;????????????%?noise-tolerance?（no?strict?fidelity?enforcement）
K?=?8;??????????????%?5?modes
DC?=?0;?????????????%?no?DC?part?imposed
init?=?1;???????????%?initialize?omegas?uniformly
tol?=?1e-7;
%---------------?Run?actual?VMD?code
figure（2）;
[u?u_hat?omega]?=?VMD（dataalpha?tau?K?DC?init?tol）;
subplot（size（u1）+121）;
plot（u（1:）‘k‘）;grid?on;
title（‘VMD分解‘）;
ylabel（{‘imf1‘;‘幅值/mV‘}）
set（gca‘xlim‘[0?2048]）;
subplot（size（u1）+122）;
plot（abs（fft（u（1:）））‘k‘）;
title（‘對應頻譜‘）;
for?i?=?2:size（u1）+1
????subplot（size（u1）+12i*2-1）;
????plot（u（i-1:）‘k‘）;grid?on;
????ylabel（{‘imfi‘;‘幅值/mV‘}）
set（gca‘xlim‘[0?2