马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
前几天再论坛上看到一个有关边际谱的程序,有的语句不太懂,请高手指点一下
clear;
fs=1000; %fs为采样频率; tspan=2;
t=1/fs:1/fs:tspan; N=length(t); %采样点数 y1=5*sin(2*pi*241*t);
y2=3*sin(2*pi*73*t);
y=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E) for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值'); [E,tt1]=toimage(A,fa,tt,length(tt));这句是什么意思,其中的E是不是指toimage.m中的im, E(k,:)又是什么,为什么 边际谱就是bjp(k)=sum(E(k,:))*1/fs |