wav file dft without fft
14 views (last 30 days)
Show older comments
I want wav audio file to dft conversion without fft
3 Comments
Rik
on 10 Jun 2020
Question title:
wav file dft without fft
Question body:
I want wav audio file to dft conversion without fft
Deleted comment (posted on 9 Jun 2020 at 3:31):
yes I want to use only DFT function don't use fft function
Accepted Answer
Ron Fredericks
on 10 Jun 2020
Edited: Ron Fredericks
on 11 Jun 2020
Hi Seung Hoon LEE:
I produced this MATLAB script to demonstrate calculation of a DFT starting from an mp3 file. Audio files typically have a large number of data points and my DFT function is very slow. So I use a loop reduction factor to speed up the DFT calculation. This parameter and frequency range can be adjusted near top of code. The audio file I use is an mp3 bird call that can be downloaded from xeno-canto.org.
%%
% PURPOSE: Compare fft with dft loop for an audio file
% AUDIO FILE: Use a public domain birdcall from xeno-canto.org
% RESULT: Produce 3 plots: 1) Time domain, 2) FFT, and 3) DFT
% DEVELOPER: http://www.biophysicslab.com/portfolio/matlab-projects/
% DATE: June 6, 2020
%
% REFERENCES:
% "Signal Processing Problems" hosted on Udemy.com
% by Dr. Mike X Cohen, sigprocMXC_SpectBirdcall.m
% "Master the Fourier transform and its applications" hosted on Udemy.com
% by Dr. Mike X Cohen, Fourier_DTFT.m
%
%%
% Load in birdcall (source: https://www.xeno-canto.org/403881).
[bc,fs] = audioread('XC403881.mp3');
% Configuration params for Fourier transforms.
freq_range = [0 8000]; % Hz
dft_loop_reduction = 100; % set to 1 for no reduction (Hint: takes a long time)
% let's hear it!
% soundsc(bc,fs)
n = length(bc);
hz = linspace(0,fs/2,floor(n/2)+1);
% Smooth birdcall audio audio for fft & dft
signal = detrend(bc(:,1))'; % transpose here for DFT loop later
% Create a time vector based on the data sampling rate.
timevec = (0:n-1)/fs;
% Plot the data from the two audio file channels.
figure(1), clf
subplot(311)
% Include a small offset for left and right audio channels.
plot(timevec,bsxfun(@plus,bc,[.2 0]))
xlabel('Time (sec.)')
title('Time domain')
set(gca,'ytick',[],'xlim',timevec([1 end]))
%% Compute & plot the power spectrum using MATLAB fft function.
bcpow_fft = abs(fft( signal )/n).^2;
subplot(312)
plot(hz,bcpow_fft(1:length(hz)),'linew',2)
xlabel('Frequency (Hz)')
ylabel('Power');
title(['Frequency domain using FFT with ' num2str(length(bcpow_fft)) ' points']);
set(gca,'xlim',freq_range)
% Make fft and dft y-limits be the same for comparison.
ylim_fft = get(gca,'ylim');
%% Compute & plot the power spectrum using DFT loop.
fourTime = (0:n-1)/n;
fCoefs = zeros(size(signal));
h = waitbar(0,'Please wait for DFT loop...');
% DFT loop is very inefficient, it is used here for demonstration only.
for fi=1:dft_loop_reduction:n
% Create complex sine wave.
csw = exp( -1i*2*pi*(fi-1)*fourTime );
% Compute dot product between sine wave and signal (Fourier coefficients).
fCoefs(fi) = sum( signal.*csw );
% GUI to show progress for long calculation times.
waitbar(fi/n)
end
close(h)
bcpow_dft = abs(fCoefs / n).^2;
subplot(313)
plot(hz,bcpow_dft(1:length(hz)),'linew',2)
xlabel('Frequency (Hz)')
ylabel('Power');
title(['Frequency domain using DFT loop with ' ...
num2str(floor(length(bcpow_dft)/dft_loop_reduction)) ' points']);
set(gca,'xlim',freq_range)
set(gca,'ylim',ylim_fft)
%% done.
Results...
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/313517/image.jpeg)
0 Comments
More Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!