91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 0.97M
    文件類型: .zip
    金幣: 1
    下載: 0 次
    發布日期: 2021-03-28
  • 語言: Matlab
  • 標簽: matlab??

資源簡介


代碼中的plot_hht函數實現了Hilbert-Huang變換(HHT)，HHT就是Hilbert-Huang Transform（希爾伯特黃變換），美國航天航空局黃鱷教授發明的，大概的過程是：先對信號進行經驗模態分解（Empirical mode decomposition －EMD），得出本征模態函數（IMF intrinsic mode function），再對本征模態函數進行希爾伯特變換，從而過進一步得該信號的希爾伯特譜、時頻能量譜等，以便對信號進行分析，據介紹對非線性及非平穩信號有較好的分析和處理效果。

資源截圖

代碼片段和文件信息 

function?imf?=?emd（x）
%?Empiricial?Mode?Decomposition?（Hilbert-Huang?Transform）
%?imf?=?emd（x）
%?Func?:?findpeaks

x???=?transpose（x（:））;
imf?=?[];
while?~ismonotonic（x）
???x1?=?x;
???sd?=?Inf;
???while?（sd?>?0.1）?|?~isimf（x1）
??????s1?=?getspline（x1）;
??????s2?=?-getspline（-x1）;
??????x2?=?x1-（s1+s2）/2;
??????
??????sd?=?sum（（x1-x2）.^2）/sum（x1.^2）;
??????x1?=?x2;
???end
???
???imf{end+1}?=?x1;
???x??????????=?x-x1;
end
imf{end+1}?=?x;

%?FUNCTIONS

function?u?=?ismonotonic（x）

u1?=?length（findpeaks（x））*length（findpeaks（-x））;
if?u1?>?0?u?=?0;
else??????u?=?1;?end

function?u?=?isimf（x）

N??=?length（x）;
u1?=?sum（x（1:N-1）.*x（2:N）?u2?=?length（findpeaks（x））+length（findpeaks（-x））;
if?abs（u1-u2）?>?1?u?=?0;
else??????????????u?=?1;?end

function?s?=?getspline（x）

N?=?length（x）;
p?=?findpeaks（x）;
s?=?spline（[0?p?N+1][0?x（p）?0]1:N）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????958087??2007-01-25?14:59??HHT.pdf
?????文件?????????965??2008-04-24?13:36??plot_hht.m
?????文件?????????893??2008-03-13?11:49??emd.m
?????文件?????????146??2007-08-26?00:42??findpeaks.m
?????文件??????102892??2007-04-28?15:57??Hum.wav

評論

共有 條評論