How to apply pca to audio files (pitch features)?

10 次查看（过去 30 天）

显示更早的评论

W.Khan 2019-11-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/489474-how-to-apply-pca-to-audio-files-pitch-features

评论： Walter Roberson 2019-11-20

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

Hello,

I applied pitch feature extraction to audio files. Now I want to apply PCA (for dimensionality reduction) to the pitch features. How can I do that? The problem is that I cannot save the features in csv or mat file. Also, the new audio feature replace the previous audio features thus giving me only feature vector i.e; comp_PCA1 gives me feature vector of 32600x1. But my files are 94 so the output should be 32600x94.

My code is given below:

clc
clear all
path_to_directory= dir ('**/*.wav'); % path to the directory 
for i = 1:length(path_to_directory)
    path= strcat(path_to_directory(i).folder, '\', path_to_directory(i).name);
    [signal, Fs]= audioread(path); % audio read here
     pitch = pwelch(signal);   % pitch feature extraction part
     [coeff,score] = pca(pitch); % pca for dimensionality reduction
     Comp_PCA1 = score(:,1);
% %      y (i,:)= coeff;   
end

7 个评论
显示 5更早的评论隐藏 5更早的评论

Ridwan Alam 2019-11-20

Great! Please mark this problem as solved or accept Walter's answer. Thanks!

Walter Roberson 2019-11-20

(I had only posted as a Comment before, so it could not be Accepted.)

请先登录，再进行评论。

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

采纳的回答

Walter Roberson 2019-11-20

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/489474-how-to-apply-pca-to-audio-files-pitch-features#answer_402329

在 MATLAB Online 中打开

path_to_directory= dir ('**/*.wav'); % path to the directory 
for i = 1:length(path_to_directory)
    filename = fullfile(path_to_directory(i).folder, path_to_directory(i).name);
    [signal, Fs]= audioread(path); % audio read here
     pitch = pwelch(signal);   % pitch feature extraction part
     [coeff,score] = pca(pitch); % pca for dimensionality reduction
     this_Comp = score(:,1);
     nR = size(this_Comp,1);
     if i == 1
         Comp_PCA1 = this_Comp;
     elseif nR <= size(Comp_PCA1,1)
         Comp_PCA1 = [Comp_PCA1(1:nR,:), this_Comp];
     else
         Comp_PCA1(:,i) = this_Comp(1:size(Comp_PCA1,1));
     end
end