求助：求瞬时频率ifndq问题

cboboc

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

lcx36

回复 1 # cboboc 的帖子

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

a1runner

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

ohaiou

刚看到，谢谢

yrj01

在求出正交函数“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

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

lzqfhz

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

lzqfhz

回复 2 # lcx36 的帖子

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