資源簡介
對信號進(jìn)行EMD分解并繪制希爾伯特時(shí)頻譜的程序

代碼片段和文件信息
%?hilbert譜分析
clear;
clc;
t=1:200;
fs=100;
t=t*1/fs;
%?x=sin(2*pi*20*t);
%?y=0.4*sin(2*pi*40*t+140);
%z=cos(2*pi*5*t)+0.5*cos(2*pi*10*t)+0.7*cos(2*pi*25*t);
%?z=exp(-0.3*t).*cos(pi*t.*t+2*pi*t);
%z=0.3*exp(j*8*t)+1*exp(j*16*t);
[fnamepname]=uigetfile(‘0a.wav‘‘Open?Wave?File‘);
file=[pnamefname];
[zfsbits]=wavread(file);
%N=length(y);
%z=?awgn(y50);
imf=emd(z);
?emd_visu(z1:length(z)imf)
[Aftt]=hhspectrum(imf);
[imtt]=toimage(Af);
disp_hhs(im[]fs);??%?繪制hilbert時(shí)頻譜
??????colormap(flipud(gray));%這條語句加上注釋改變圖形
title(‘hilbert時(shí)頻譜‘);
colorbar;
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????618??2009-05-08?16:12??hsa.m
-----------?---------??----------?-----??----
??????????????????618????????????????????1
評論
共有 條評論