I need to isolate certain values of a matrix based on whether one of the values fall in a certain range
3 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Wynand
el 19 de Oct. de 2023
Comentada: Dyuman Joshi
el 19 de Oct. de 2023
I am writing a word recognition program that takes an array from 'audioread' and applies a few methods to it to get a frequency spectrum. This spectrum is saved into a cell aray with each cell containing the frequency peaks for a quater of the word at a time.
So for each word it returns a 4x1 cell array. Each of these cells need to be further subdivided to frequency bands
An example of one such a cell follows
% ans = 2×2
% 24.9843526628176 32
% 5.44020124026728 62
The first column indicates amplitude and the second indicates the index. I need to split the matrix into divisions that are 30 wide each. So all 2nd column values from 0 to 30 must be together and all from 30 to 60 etc. All the way up to 300.
The following code exert is the function that I wrote for this application for context to the question.
It is probably flawed and unoptimal so if there are glaring issues in it feel free to include commentary in any answers.
Thank you
function FINAL_OUTPUT = FFT_Output(OG_Audio, Filter_count, Min_peak_distance)
% Setup
Sectional_Length = round(length(OG_Audio)/4);
BM_Filtered_Signal = (blackman(length(OG_Audio))).*OG_Audio;
First_Filter_Parameter = ((1/(Sectional_Length/100))*ones(1,(Sectional_Length/100)));
Word_section = cell(4,1);
FFT_Word = cell(4,1);
FINAL_OUTPUT = cell(4,1);
for i = 1:4
% Word sections
Word_section{i} = BM_Filtered_Signal(((Sectional_Length*(i-1))+1):(Sectional_Length*i));
% Moving average filter
j=0;
while(j<=Filter_count)
Word_section{i} = filter(First_Filter_Parameter, 1, Word_section{i});
j=j+1;
end
% FAST FOURIER TRANSFORM
FFT_Word{i} = abs(fft(Word_section{i},length(Word_section{i})));
FFT_Word{i} = FFT_Word{i}(1:(round(length(FFT_Word{i})/(2))));
% Peak detection
[M,I] = findpeaks(FFT_Word{i}, "MinPeakDistance",Min_peak_distance, "MinPeakHeight",(0.1*max(FFT_Word{i}))); % Identifying peak frequencies for comparisons
FINAL_OUTPUT{i} = [M,I];
end
end
0 comentarios
Respuesta aceptada
Dyuman Joshi
el 19 de Oct. de 2023
Editada: Dyuman Joshi
el 19 de Oct. de 2023
Here's a general approach - Discretize the data according to the groups and then split then into cells.
bins = 0:30:300;
idx = discretize(data, bins)
out = splitapply(@(x) {x}, data, idx)
2 comentarios
Más respuestas (0)
Ver también
Categorías
Más información sobre Blackman en Help Center y File Exchange.
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!