lyvv 发表于 2009-5-21 15:53

有人做过瞬时相位吗? 是不是和hilbert有关

本帖最后由 wdhd 于 2016-9-18 13:56 编辑


上图是别人做的瞬时相位的结果


但是没有算法,在网上到找搜索,好像是和hilbert有关,又找了一个hilbert的算法,但是结果和上图不一致


上图是在网上找的算法,差异非常大。



网上hilbert算法的C代码


//n和N我都是带入数组的长度
void hilbert(long n, float delta[], float kappa[],long N)
{   
float d1, d2, d3, d4;   
int i1, i2;   
for (i1 = 0; i1 < n; i1++)   
{   
kappa = 0.;   
for (i2 = 1; i2 < n; i2++)   
{   
   d1 = (i1+i2<N)? delta:0.;   
   d2 = (i1-i2>=0)? delta: 0.;   
   d3 = (i1+i2+1<N)? delta: 0.;   
   d4 = (i1-i2-1>=0)? delta: 0.;   
   kappa -= 0.5 * (d1-d2) / i2 + 0.5 * (d3 - d4) / (i2+1);   
}   
kappa /= PI;   
}   
}

  希望有高手能告诉怎么得到第一张图底部图的波形??不甚感激~!

qq 1960064

[ 本帖最后由 lyvv 于 2009-5-21 16:02 编辑 ]

lyvv 发表于 2009-5-21 15:58

一份数据文件
和一个可以查看波形的简单程序 (把后缀名改成zip)

程序截图


[ 本帖最后由 lyvv 于 2009-5-21 16:00 编辑 ]

lyvv 发表于 2009-5-22 09:20

有人知道么?
页: [1]
查看完整版本: 有人做过瞬时相位吗? 是不是和hilbert有关