Frequency spectrum of a sound signal

51 visualizaciones (últimos 30 días)
Lok Hua Lee
Lok Hua Lee el 20 de Oct. de 2019
Comentada: Keira Duffy el 20 de Sept. de 2021
Hi guys, I would like to know some hints on how to plot frequency spectrum of magnitude and phase spectra of an audio signal in both omega and frequency as x-axis parameter (plot separately). Thanks. My code is as below and i'm not sure what's going on.
[y,fs] = audioread('test.wave');
N = length(y);
t = (0:N-1/fs);
n = (0:N-1);
y = y(:,1);
% spectral analysis
w = hanning(N, 'periodic');
[X, f] = periodogram(y, w, N, fs, 'power');
X = 20*log10(sqrt(X)*sqrt(2));
% plot the signal spectrum
figure(1);
subplot(2,1,1);
semilogx(f, X, 'r');
xlim([0 max(f)]);
grid on;
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14);
title('Amplitude spectrum of the signal');
xlabel('Frequency, Cycles/Second');
ylabel('Magnitude, dB');
subplot(2,1,2);
semilogx((2*pi*f), X, 'r');
xlim([0 max(2*pi*f)]);
grid on;
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14);
title('Amplitude spectrum of the signal');
xlabel('Angular Frequency, Radian/Sample');
ylabel('Magnitude, dB');

Respuestas (1)

Kaashyap Pappu
Kaashyap Pappu el 23 de Oct. de 2019
Using a test sound file of my own, I was able to generate the plots attached. The function “audioread” works if the suffix of the sound file name is ‘.wav’ not ‘.wave’.
A similar question has been addressed here.
Hope this helps!
  2 comentarios
Lok Hua Lee
Lok Hua Lee el 23 de Oct. de 2019
Thanks for the reply, the following code is what i think should be working fine. But im not sure if the graph should be the same for both frequency and angular frequency?
clear, clc, close all;
[y,fs] = audioread('test.wav');
N = length(y); % Length of vector y, number of samples
Y = fft(y,N); % Fourier transform of y
F = ((0:1/N:1-1/N)*fs); % Frequency vector
w = 2*pi*F; % Angular frequency vector
magnitudeY = abs(Y); % Magnitude of the FFT
phaseY = unwrap(angle(Y)); % Phase of the FFT
figure (1);
subplot(2,1,1);
plot(F, magnitudeY);
grid on;
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14);
xlabel('Frequency, Hz');
ylabel('Magnitude, dB');
title('Magnitude spectrum of sound wave in frequency');
subplot(2,1,2);
plot(w, magnitudeY);
grid on;
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14);
xlabel('Angular Frequency, rad/sample');
ylabel('Magnitude, dB');
title('Magnitude spectrum of sound wave in angular frequency');
figure (2);
subplot(2,1,1);
plot (F, phaseY);
grid on;
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14);
xlabel('Frequency, Hz');
ylabel('Phase angle, radian');
title('Phase spectrum of sound wave in frequency');
subplot(2,1,2);
plot (w, phaseY);
grid on;
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14);
xlabel('Angular Frequency, rad/sample');
ylabel('Phase angle, radian');
title('Phase spectrum of sound wave in angular frequency');
Keira Duffy
Keira Duffy el 20 de Sept. de 2021
When you do magnitudeY = abs(Y); you automatically assumed that it is in decibels but it is not. Dont you need to take the log10 to convert it?

Iniciar sesión para comentar.

Categorías

Más información sobre Vibration Analysis en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by