求助：求瞬时频率ifndq问题

cboboc · 发表于 2010-3-10 08:14

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

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

x

请问各位：
在国立中央大学的求瞬时频率的程序中，求相位函数为什么用这个式子dq(i)=sqrt(1-nvimf(i)*nvimf(i));？根据是什么呢？

lcx36 · 发表于 2010-11-19 09:39

回复 1 # cboboc 的帖子

我也有同样问题，希高手指点。
而且，求的dq其实在后面又没有用到，不知何故?
还有，omgcos(i)=abs(nvimf(i+1)-nvimf(i-1))*0.5/sqrt(1-nvimf(i)*nvimf(i));是根据什么公式或理论来的，也有点晕。

a1runner · 发表于 2010-11-23 21:44

可以访问他们的网站看，dq(i)=sqrt(1-nvimf(i)*nvimf(i))是一种更好的球瞬时频率的方法

ohaiou · 发表于 2011-3-29 10:47

刚看到，谢谢

yrj01 · 发表于 2011-5-13 15:35

在求出正交函数“dq(i)“之后的下述语句是怎么求瞬时频率的？哪里是通过求相位函数进而对其求导来计算频率的？
rangebot=-rangetop;
for i=2:(vlength-1),
if nvimf(i)>rangebot & nvimf(i) < rangetop
      omgcos(i)=abs(nvimf(i+1)-nvimf(i-1))*0.5/sqrt(1-nvimf(i)*nvimf(i));
else
      omgcos(i)=-9999;
end
end
omgcos(1)=-9999;
omgcos(vlength)=-9999;

jj=1;
for i=1:vlength,
if omgcos(i)>-1000
      ddd(jj)=i;
      temp(jj)=omgcos(i);
      jj=jj+1;
end
end
temp2=spline(ddd,temp,dd);
omgcos=temp2;

for i=1:vlength,
omega(i)=omgcos(i);
end

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

帮你们顶一下，偶没有这个方法所以不懂得，呵呵

lzqfhz · 发表于 2012-6-25 14:29

LZ 能否给分享下他的瞬时频率求解的程序。
lzq@bit.edu.cn

lzqfhz · 发表于 2012-6-26 11:15

回复 2 # lcx36 的帖子

你只要对cos(\theta)求导数就能看出他使用的算法了。最后得到的omgcos是\theta的导数。lzq@bit.edu.cn

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

[HHT] 求助：求瞬时频率ifndq问题

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