1/f amplitude spectrum and slope

8 次查看（过去 30 天）

Mercede Erfanian 2021-2-23

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/753654-1-f-amplitude-spectrum-and-slope

评论： ANURAG SINGH 2021-2-24

Hi,

I would like to compute 1/f amplitude spectrum and loglog plot Or log10 .wav files and here is the code which does not work for many reasons. Any help is appreciated.

clear all;

close all;

clc;

[x,Fs] = audioread('*.wav');

Pxx = zeros(129,size(x,2));

for nn = 1:size(x,2)

[Pxx,F] = pwelch(x(:,nn),hamming(128),[],[],Fs,'psd');

end

PSD_x = 1./F(1:end).^2;

figure;

plot(log10(F),10*log10(PSD_x),'r','linewidth',4);

hold on;

plot(log10(F),10*log10(Pxx));

xlabel('loglog(Hz)'); ylabel('loglog(dB)');

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

采纳的回答

Mathieu NOE 2021-2-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/753654-1-f-amplitude-spectrum-and-slope#answer_631274

在 MATLAB Online 中打开

hello Mercede

Slightly modified code, simply give it a valid wav file name - it works as can be seen below - or where is your issue ?

clear all;
close all;
clc;
 
[x,Fs] = audioread('unfiltered_sound.wav'); % my test file name
 
% Pxx = zeros(129,size(x,2));
for nn = 1:size(x,2)
[Pxx,F] = pwelch(x(:,nn),hamming(128),[],[],Fs,'psd');
end
 
PSD_x = 1./F(1:end).^2;
figure;
% plot(log10(F),10*log10(PSD_x),'r','linewidth',4);
semilogx(F,10*log10(PSD_x),'r','linewidth',4);
hold on;
semilogx(F,10*log10(Pxx));
% plot(log10(F),10*log10(Pxx));
% xlabel('loglog(Hz)'); ylabel('loglog(dB)');
xlabel('(Hz)'); ylabel('(dB)');