clc;
clear?all;
%?load（‘VMD_basic.mat‘）
%%取出第一行與信號相減
load（‘simulation.mat‘）
x_trunc（1:1:2048）?=?x（1:1:2048）;
%?self_sequence（1:3:512）?=?1;
%?self_sequence（2:3:512）?=?2;
%?self_sequence（3:3:512）?=?3;
[SFNFT]?=?spectrogram（x_trunc?hamming（128）1272561/139.975e6‘centered‘）;
figure
imagesc（TFabs（SFN））;

%?some?sample?parameters?for?VMD
alpha?=?1000;????????%?moderate?bandwidth?constraint
tau?=?0;????????????%?noise-tolerance?（no?strict?fidelity?enforcement）噪聲耐受，一般取0就行
K?=?5;??????????????%?3?modes
DC?=?0;?????????????%?no?DC?part?imposed
init?=?1;???????????%?initialize?omegas?uniformly
tol?=?1e-9;?????????%總信號擬合誤差允許程度
[u_after_d?u_hat?omega]?=?VMD（x_trunc?alpha?tau?K?DC?init?tol）;
%---------------?Visualization

for?k=1:K
[SFNFT]?=?spectrogram（u_after_d（k:）?hamming（128）1272561/139.975e6‘centered‘）;
figure
mesh（TFabs（SFN））;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????923??2018-12-09?15:28??test_difference.m

?????文件?????126072??2018-12-06?09:59??simulation.mat

?????文件???????4851??2018-12-09?17:58??vmd_fun.m

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

???????????????131846????????????????????3