finding zero crossings in a signal

39 次查看（过去 30 天）
SSG_newbiecoder
SSG_newbiecoder 2018-1-6
编辑： Star Strider 2018-1-9
I need to find out the zero crossing points in a signal. How can I achieve this in matlab? I tried the code below but it is not working for my signal.Can anybody help me out?
t = [1:0.01:5]; % Time Vector
y = -2+2.*sin(2*pi*t); % Signal
zci = @(v) find(v(:).*circshift(v(:), [-1 0]) <= 0); % Returns Zero-Crossing Indices Of Argument Vector
zx = zci(y); % Approximate Zero-Crossing Indices
figure(1)
plot(t, y, '-r')
hold on
plot(t(zx), y(zx), 'bp')
hold off
grid
legend('Signal', 'Approximate Zero-Crossings')

请先登录，再回答此问题。

回答（2 个）

Star Strider
Star Strider 2018-1-6
编辑：Star Strider 2018-1-9
The code works correctly. The amplitude of your signal goes from -4 to 0, so the ‘zero-crossings’ in your signal are correctly detected at the highest value of each cycle.
If you want the indices where the signal crosses the -2 value, the ‘zci’ call should be:
zx = zci(y + 2); % Approximate Zero-Crossing Indices
EDIT #1 — (08 Jan 2018 at 13:00 UCT)
‘My issue is that my signal need not have an array index with zero as a value so I'm supposed to take the value as close to zero as possible and consider the index of that value as the zero crossing point.’
My code does exactly that!
EDIT #2 — (09 Jan 2018 at 00:36 UCT)
If you want to find the midpoint crosses of your waveform, consider using the midcross (link) function.
Birdman
Birdman 2018-1-6
To find their time values, type
t=t(y==0)
To find their values in y vector, type
y=y(y==0)
  2 个评论
SSG_newbiecoder
SSG_newbiecoder 2018-1-8
My issue is that my signal need not have an array index with zero as a value so I'm supposed to take the value as close to zero as possible and consider the index of that value as the zero crossing point.
Birdman
Birdman 2018-1-8
Ok, then you may change the condition as following
idx=find(y==0);
y(idx+1)
y(idx-1)
This will result in the closest value to zero in your y vector which is dependent on t vector. The result is(for both):
-0.0039 -0.0039 -0.0039 -0.0039
What if your t vector's increment change? For instance:
t=1:0.001:5;
y=-2+2.*sin(2*pi*t);
Then, when the same procedure applies,
idx=find(y==0);
y(idx+1)
y(idx-1)
The result is:
1.0e-04 *
-0.3948 -0.3948 -0.3948 -0.3948
which is different than before. Hope this helps.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Pulse and Transition Metrics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by