MATLAB: How can i restore complex numbers from abs(fft)

complexfftfrequencyifftrestoresound

I made this function

 {
function [y,AudioFft,fs,TimeDomain,FreqDomain] = ReadWavFile(filename)
[y,fs] = audioread(filename); %Load File
TimeDomain = 0:1/fs:(length(y)-1)/fs; % Time Domain
N = length(y); %wav file length
AudioFft = fft(y); % fft
AudioFft = AudioFft(1:N/2+1); %Deleting mirror image that was made by abs 
FreqInterval = fs/N; % Frequency Interval
FreqDomain = 0:FreqInterval:fs/2; %Freq domain (deleted mirror image so fs/2)
end
}

I am making equalizer by matlab

I made the graph(x domain is Freq and y domain is amplitude(that was made by abs(AudioFft), exactly plot(FreqDomain,abs(AudioFft)) that changed by moving of sliderbar of gui

I want to play sound that changed by moving of slider bar

So I think I have to restore abs and return to Time domain by ifft to use sound function(sound(y,fs)).

How can i do?

In short, How can i restore complex number?

Best Answer

"In short, How can i restore complex number?"

In short, "you can't" once you've thrown away the separate magnitude and phase, it's lost irretrievably. "Don't do that!"

Related Solutions

MATLAB: Can the amplitude be plotted in dB

Try this:

semilogx(Fv, 20*log10(abs(FT_af(Iv,2))*2)), xlabel('Frequency, Hz'),ylabel('magnitude');

That will convert the amplitudes to dB. The y-axis will reflect those changes, however the tick marks will not be in integer powers-of-10 such as with your x-axis and the semilogx call. You will have to format the y-axis ticks yourself.

MATLAB: Wav file dft without fft

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...

Best Answer

Related Solutions

MATLAB: Can the amplitude be plotted in dB

MATLAB: Wav file dft without fft

Related Question