回复 #15 rc-hw-0002 的帖子
能知道你的qq吗？[ 本帖最后由 zhlong 于 2007-8-22 15:14 编辑 ] 本帖最后由 VibInfo 于 2016-10-17 15:12 编辑
原帖由 rc-hw-0002 于 2007-4-21 23:05 发表
画边际谱
function ms=HHTms(x,N)
% Input-
% x - 2-D matrix x(k,n) of the HHT spectrum
% N - 分辨率
% Output-
% ms - vector ms(k) that specifies the margi ...
你这里的w=1:length(ms);w=w/N;是何用意呀？ 画边际谱
function ms=HHTms(x,N)
% Input-
% x - 2-D matrix x(k,n) of the HHT spectrum
% N - 分辨率
% Output-
% ms - vector ms(k) that specifies the marginal spectrum surf(h) shading interp
if(nargin<2)
N=1;
end
n=size(x);
k=n(1);
n=n(2);
ms=sum(x')'/n;
w=1:length(ms);
w=w/N;
%画边际谱：
figure
plot(w,ms);
xlabel('\fontsize{18}频率 / Hz');
ylabel('\fontsize{18}幅值');
legend('\fontsize{18}希尔伯特边际谱');
楼主的程序里N是分辨率，那么1/N就是采样周期了，单位是时间t/s,而w=w/N就是周期乘以一个序列，单位还是s,这和你的图对不上吧？不知道我这么理解对不对？
所以我想知道楼主的w=1:length(ms);w=w/N;是想实现什么的？
谢谢
回复 楼主 rc-hw-0002 的帖子
好贴呀！！！[ 本帖最后由 大鹏之举 于 2009-4-29 22:17 编辑 ] 本帖最后由 VibInfo 于 2016-10-17 15:12 编辑
原帖由 zhangnan3509 于 2007-4-21 23:20 发表
%----- Construct the ploting matrix
clear p;
h1=zeros(npt-1,nyy+1);
p=round(nyy*(omg-wmn)/dw)+1;
for j1=1:npt-1
for i1=1:knb
ii1=p(j1,i1);
h1(j1,ii1)=h1(j1,ii1)+a(j1,i1);
end
end
请问这里的p求出来的是什么频率呢？为什要这样求呢 本帖最后由 VibInfo 于 2016-10-17 15:12 编辑
原帖由 zhangnan3509 于 2007-4-21 23:20 发表
%----- Define the results
%w=linspace(wmn,wmx,ny-1)';
%xs=linspace(t0,t1,nx)';
h1=flipud(rot90(h1));
h1=h1(1:ny-1,:);
还有这里，%w=linspace(wmn,wmx,ny-1)';不是求频率轴上的坐标吗？是不是应该为nyy-1呢？ 还有，请问各位高手，什么是三点光滑？ 回复 13 # zhangnan3509 的帖子
请问程序里
for i=1:knb
a(:,i)=filtfilt(filtr,1,a(:,i));
omg(:,i)=filtfilt(filtr,1,omg(:,i));
end
为什么要循环两次啊？