Finding number of pulses in a data set?

14 次查看（过去 30 天）
Joey Parker
Joey Parker 2020-5-6
评论： Michael Soskind 2020-5-6
I have a bit_stream.mat data file where it contains a single row of 1's and 0's in a random alternating sequence (1x664 double in workspace). How do I go about finding the number of pulses in that data set?

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

回答（1 个）

Michael Soskind
Michael Soskind 2020-5-6
Hi Joey,
Looks like findpeaks should be your friend here.
load('bit_stream.mat');
plot(bit_stream)
[ones,loc_ones,w,p] = findpeaks(bit_stream);
numel(ones) % number of peaks with respect to
[pulses,loc_pulses] = findpeaks(bit_stream, 'MinPeakWidth', 2);
numel(pulses) % number of peaks with a value of 1 2 or more subsequent times
Best,
Michael
  4 个评论
显示 2更早的评论
Joey Parker
Joey Parker 2020-5-6
Thank you that makes perfect sense. I have an additionally question if it involves the same function but for finding the pulse period of that data set?
Michael Soskind
Michael Soskind 2020-5-6
So it depends, usually, mathematically that would come from the fft of the signal. Hoewever, given that the pulses do not have a particular period, as they are pretty random in the example data you provided, the sidebands you would want are quite small. However, you could still use the following code:
% Using FFT to find the most prevalent period
figure(); % defining a figure to plot the fft in
plot(fftshift(abs(fft(bit_stream)))); % pefrforming an fft of the data (shifting to have 'zero mean', but this )
% from the figure above, looks to be a period of 4, as there is a shift of four between each of the peaks
% More trivial method of taking the average of the distance between the different elements
figure(); % defining a figure to plot the differences in
plot(diff(loc_ones)); hold on; % plotting the differences of loc_ones using the diff function
plot(1:numel(loc_ones), repmat(mean(diff(loc_ones)),1,numel(loc_ones))); % plotting the average, calculated as mean(diff(loc_ones))
The second set of code plots that the average of the difference between the different functions is in fact close to 4, confirming the effectiveness of the fft. However, taking the average of the differences might be easier.
The mean function and diff function are all you really need, just thought it might be of interest to see the similarities of two different methods.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Fourier Analysis and Filtering 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by