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

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

 

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

 

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

 

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

 

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

 

我贴一下吧：

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

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

 

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

我改了一下，不知道改的好不好，俩百分号之间是我加的
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;

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 编辑 ]