如何找出一组数据的过零点？？

lierr

比如采样到一组正弦波信号，由于采样率的关系，未必有真正意义上的零点
只是在接近零点处有值，如何找到其过零点
我自己想了个办法，就是[a,b]=min(abs(yy));
但是这样只能找到若干过零点的最小一个值，有没有现成的命令可以把所有零点都找出来？？

eight

到信号处理版块，寻找 EMD 工具箱的程序，然后打开 extr.m（或者 extr 函数）就可以看到

[ 本帖最后由 ChaChing 于 2009-12-7 21:37 编辑 ]

lierr

不好意思，以前没有用过emd，完全不懂，能否具体说说，谢谢了

[ 本帖最后由 ChaChing 于 2009-12-7 21:37 编辑 ]

eight

我不是建议你使用 EMD，你也没有必要理解它，我的意思是 EMD 的实现过程中用到这个方面的知识，你找到 extr 函数就能看到你的问题的解答了（extr 函数用于求一个信号的所有零点和极值点）

[ 本帖最后由 ChaChing 于 2009-12-7 21:38 编辑 ]

lierr

看了一些帖子，我现在要下emd的工具箱，可是你发的那个地址我打不开，也许我这里是教育网的缘故，好久没有能用的代理了
能不能给俺发一下lh900@163.com

[ 本帖最后由 ChaChing 于 2009-12-7 21:46 编辑 ]

eight

我贴一下吧：

1. 
   
2. m = length(x); x1=x(1:m-1); x2=x(2:m);
   
3. indzer = find(x1.*x2<0);
   
4. 
   
5. if any(x == 0)
   
6.    iz = find( x==0 ); indz = [];
   
7.    if any(diff(iz)==1)
   
8.      zer = x == 0;
   
9.      dz = diff([0 zer 0]);
   
10.      debz = find(dz == 1);
    
11.      finz = find(dz == -1)-1;
    
12.      indz = round((debz+finz)/2);
    
13.    else
    
14.      indz = iz;
    
15.    end
    
16.    indzer = sort([indzer indz]);
    
17. end
    
18. 
    

复制代码

[ 本帖最后由 ChaChing 于 2009-12-7 21:47 编辑 ]

mulan

原帖由 lierr 于 2007-6-5 22:02 发表

                               
登录/注册后可看大图

看了一些帖子，我现在要下emd的工具箱，可是你发的那个地址我打不开，也许我这里是教育网的缘故，好久没有能用的代理了
能不能给俺发一下lh900@163.com

给你推荐个加速器吧，亿目，和教育网直通车类似，现在是免费的

lierr

原帖由 eight 于 2007-6-5 22:06 发表

                               
登录/注册后可看大图

我贴一下吧：

m = length(x);
x1=x(1:m-1);
x2=x(2:m);
indzer = find(x1.*x2

这个程序有点问题，它找的是零点左侧的点，有时候右侧的点会更接近零点一些

我改了一下，不知道改的好不好，俩百分号之间是我加的
function pzero=fz(x)
m = length(x);
x1=x(1:m-1);
x2=x(2:m);
indzer = find(x1.*x2<0);
%%
n=length(indzer);
for i=1:n
    if abs(x(indzer(i)))>abs(x(indzer(i)+1))
        indzer(i)=indzer(i)+1;
    end
end
%%
if any(x == 0)
   iz = find( x==0 );
   indz = [];
   if any(diff(iz)==1)
     zer = x == 0;
     dz = diff([0 zer 0]);
     debz = find(dz == 1);
     finz = find(dz == -1)-1;
     indz = round((debz+finz)/2);
   else
     indz = iz;
   end
   indzer = sort([indzer indz]);
end
pzero=indzer;

lierr

原帖由 mulan 于 2007-6-5 22:24 发表

                               
登录/注册后可看大图

给你推荐个加速器吧，亿目，和教育网直通车类似，现在是免费的

好啊，万分感谢！！

leeking30

There may has something wrong in the code!
if the input is x=[1 0 -2 -3 2 -5 3 0 2 0     0   5   ]
the reture value should be pzero=[2 5 7 8 10 11]
but using the code above, the reture value
indzer =

     2     4     5     6     8    11
where is 10  ?? Obviously it is zero!

something wrong with my chinese input, Sorry!

I just make a little modification of the code above
here's my codes

function pickzero=fz(x)
m = length(x);
x1=x(1:m-1);
x2=x(2:m);
indzer = find(x1.*x2<=0);
n=length(indzer);
for i=1:n
    if abs(x(indzer(i)))>abs(x(indzer(i)+1))
        indzer(i)=indzer(i)+1;
    end
end
pickzero=indzer

[ 本帖最后由 leeking30 于 2008-10-29 06:48 编辑 ]