%?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時頻譜
??????colormap（flipud（gray））;%這條語句加上注釋改變圖形
title（‘hilbert時頻譜‘）;
colorbar;

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

?????文件????????618??2009-05-08?16:12??hsa.m

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

??????????????????618????????????????????1