研发埠

标题: matlab hht工具箱已经安装，看看我这个哪里的错误？ [打印本页]

作者: 欧阳丹 时间: 2013-5-28 13:25
标题: matlab hht工具箱已经安装，看看我这个哪里的错误？
>> clear;N=1024;n=1:N;fs=1024;t=n/fs;x=0.3*sin(2*pi*25 *t);y=sin(2*pi*50 *t);a=0.2* sin(2*pi*100*t);b= 0.01*sin(2*pi*150 *t); z=x+y+b+a;data=z;imf=emd(data); %对输入信号进行EMD分解 [A,f,t]=hhspectrum(imf); %对IMF分量求取瞬时频率与振幅：A：是每个IMF的振幅向量,f:每个IMF对应的瞬时频率，t：时间序列号[E,t,Cenf]=toimage(A,f); %将每个IMF信号合成求取Hilbert谱，E：对应的振幅值，Cenf：每个网格对应的中心频率这里横轴为时间，纵轴为频率 %即时频图(用颜色表示第三维值的大小)和三维图（三维坐标系：时间，中心频率，振幅） cemd_visu(data,1:length(data),imf); %显示每个IMF分量及残余信号--------------------------------------------disp_hhs(E); %希尔伯特谱----------------------------------------------------------%画出边际谱%N=length(Cenf);%设置频率点数 for k=1:size(E,1) bjp(k)=sum(E(k,

)*1/fs;endfigure(3);plot(Cenf(1,

*fs,bjp); % 作边际谱图进行求取Hilbert谱时频率已经被抽样成具有一定窗长的离散频率，所以此时的频率轴已经是中心频率xlabel(&#39;频率 / Hz&#39;);ylabel(&#39;幅值&#39;);??? Undefined command/function &#39;instfreq&#39;.Error in ==> hhspectrum at 45 f(i,

=instfreq(an(i,

&#39;,tt,l)&#39;;>>

欢迎光临研发埠 (http://bbs.yanfabu.com/)