边界谱问题求助！

li03152034625 · 发表于 2011-5-7 17:26

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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集
%%%%%%%%%%%%%求边际谱
[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('幅值');

我想把fs=1000; %fs为采样频率  ，改为fs=10000，请问函数其他地方要怎么改啊？谢谢。

chenlu1986 · 发表于 2011-5-14 19:18

回复 1 # li03152034625 的帖子

什么意思？？

li03152034625 · 发表于 2011-5-16 19:35

回复 2 # chenlu1986 的帖子

就是想把采样频率变大，不知道程序其他地方是否需要相应的修改。

Cena马 · 发表于 2011-5-17 10:48

其他地方不用改！

li03152034625 · 发表于 2011-5-20 09:02

回复 4 # Cena马的帖子

我知道啦十分感谢

账号		自动登录	找回密码
密码			我要加入

[HHT] 边界谱问题求助！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。