how to save frames of data into one struct?

조회 수: 10(최근 30일)
Jun W
Jun W 2018년 5월 13일
댓글: Walter Roberson 2018년 5월 25일
I have the following code to read one frame data--> n points-->find 10 peaks(from n points). say I want to read 10000 frames of data, then save all the (1000*10)peaks into a struct. how do I modify the code?
a=arduino
n=20;
for ii=1:n
v(ii)=readVoltage(a,'A0');
end
[pks,loc]=findpeaks(real(fft(v)),'npeaks',10);

채택된 답변

Walter Roberson
Walter Roberson 2018년 5월 13일
repeat_count = 10000;
my_data(repeat_count) = struct('pks', [], 'loc', []); %initialize for efficiency
n=20;
for count = 1 : repeat_count
v = zeros(1, n);
for ii=1:n
v(ii) =readVoltage(a,'A0');
end
[my_data(count).pks, my_data(count).loc] = findpeaks(real(fft(v)),'npeaks',10);
end
Note: if you want 10 peaks out of 20 points. then you will need to select index 2 through 11 from the fft result.
Selecting peaks from the entire fft result is not the right thing to do because fft of real data returns complex conjugate data symmetric around the middle of the output.
Because peak selection requires that the point before and the point after be less than or equal to the current point, the maximum number of peaks that can be located is half of the input size -- and if you are going for that many you might as well just take every second point.
Effectively the maximum number of peaks you can find from real(fft(v)) is length(v)/4, and you should probably be aiming for less than that.
  댓글 수: 6
Walter Roberson
Walter Roberson 2018년 5월 25일
? length of the signal would be length(v) like you are using now, so I am not sure what you are asking about?
I suggest, by the way,
P1 = P2(1:ceil(end/2));

댓글을 달려면 로그인하십시오.

추가 답변(0개)

Community Treasure Hunt

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

Start Hunting!

Translated by