|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
初涉EMD内容。在网上找到一些关于hht边际谱的内容,不大理解。用matlab如何求解边际谱?在论坛以前的帖子里找到一段小程序,有些地方不明白,想请教下大家:
clear;
fs=1000; %fs为采样频率;
N=1000; %采样点数(如何定?)
t=1/fs:1/fs:1;%
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=[y1;y2;zeros(size(y1))]; %IMF集(集合中为什么还要追加0行?)
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);%
[E,tt1]=toimage(A,fa,tt,length(tt));%
E=flipud(E);%(为什么要水平翻转?)
plot(tt1,E);
hold on;
for k=1:size(E,1)%
bjp(k)=sum(E(k,:))*1/fs;%(bjp是什么量?)
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
hold off; |
|