lierr 发表于 2007-6-5 21:09

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

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

eight 发表于 2007-6-5 21:14

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

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

lierr 发表于 2007-6-5 21:29

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

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

eight 发表于 2007-6-5 21:31

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

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

lierr 发表于 2007-6-5 22:02

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

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

eight 发表于 2007-6-5 22:06

我贴一下吧：
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();
   debz = find(dz == 1);
   finz = find(dz == -1)-1;
   indz = round((debz+finz)/2);
   else
   indz = iz;
   end
   indzer = sort();
end



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

mulan 发表于 2007-6-5 22:24

原帖由 lierr 于 2007-6-5 22:02 发表 http://www.chinavib.com/forum/images/common/back.gif



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

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

lierr 发表于 2007-6-5 23:06

原帖由 eight 于 2007-6-5 22:06 发表 http://www.chinavib.com/forum/images/common/back.gif


我贴一下吧：

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();
   debz = find(dz == 1);
   finz = find(dz == -1)-1;
   indz = round((debz+finz)/2);
   else
   indz = iz;
   end
   indzer = sort();
end
pzero=indzer;

lierr 发表于 2007-6-5 23:07

原帖由 mulan 于 2007-6-5 22:24 发表 http://www.chinavib.com/forum/images/common/back.gif


给你推荐个加速器吧，亿目，和教育网直通车类似，现在是免费的
好啊，万分感谢！！

leeking30 发表于 2008-10-29 05:24

There may has something wrong in the code!
if the input is x=
the reture value should be pzero=
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 编辑 ]
页: [1]
查看完整版本: 如何找出一组数据的过零点？？