# downsampling / reducing fft data for log plot

17 views (last 30 days)

Show older comments

Hello,

I want to plot the data from an fft, however as the source data is very large consequently my fft data consists of one million samples. I want to plot the data in a log-log plot, but of course not all of the samples are necessary.

How can I downsample / reduce the number of samples/frequencies of the fft result while keeping the quality of the plot? If I remove for example every 4th sample it won't have an effect on the high frequencies, but I will lose precious samples on the low frequency range.

Thanks for your help!

To make it a bit more clear: The fft results have a sampling rate of 100Hz: 100Hz, 200Hz, 300Hz,...1e8Hz When I plot the spectrum logarithmically the range between 100 and 1000Hz has 9 samples (ok) between 1000 and 10000Hz has 90 samples (ok) between 10000 and 100000Hz has 900 samples (still ok, starts getting too much) between 1e7 and 1e8Hz has 900000 samples (way too much)

So I want to reduce the number of samples, especially in the higher frequencies, while keeping the low frequencies more intact, as a strong reduction of samples will change the plot significantly.

##### 2 Comments

Walter Roberson
on 30 Oct 2012

### Answers (4)

Wayne King
on 31 Oct 2012

Edited: Wayne King
on 31 Oct 2012

If you have the Signal Processing Toolbox, you can use goertzel.m to compute the DFT at specified frequencies. For example:

Fs = 1000;

t = 0:1/Fs:1-1/Fs;

x = cos(2*pi*100*t)+1/2*cos(2*pi*200*t)+1/4*cos(2*pi*300*t);

The indices in the DFT for 100,200, and 300 Hz are 101,201,and 301

idx = [101 201 301];

xdft = goertzel(x,idx);

A better example would be: Resolve the 1.24 kHz and 1.26 kHz components in the following noisy cosine which also has a 10 kHz component. The data is sampled at 32 kHz.

Fs = 32e3; t = 0:1/Fs:2.96;

x = cos(2*pi*t*10e3)+cos(2*pi*t*1.24e3)+cos(2*pi*t*1.26e3)...

+ randn(size(t));

N = (length(x)+1)/2;

f = (Fs/2)/N*(0:N-1); % Generate frequency vector

indxs = find(f>1.2e3 & f<1.3e3); % Find frequencies of interest

X = goertzel(x,indxs);

plot(f(indxs)/1e3,20*log10(abs(X)/length(X)));

title('Mean Squared Spectrum');

xlabel('Frequency (kHz)');

ylabel('Power (dB)');

grid on;

set(gca,'XLim',[f(indxs(1)) f(indxs(end))]/1e3);

where you select a range of frequencies.

You can also just plot a range of frequencies in the fft() output. You don't have to plot the entire output.

Wayne King
on 31 Oct 2012

To make it a bit more clear: The fft results have a sampling rate of 100Hz: 100Hz, 200Hz, 300Hz,...1e8Hz When I plot the spectrum logarithmically the range between 100 and 1000Hz has 9 samples (ok) between 1000 and 10000Hz has 90 samples (ok) between 10000 and 100000Hz has 900 samples (still ok, starts getting too much) between 1e7 and 1e8Hz has 900000 samples (way too much)

You do not mean sampling rate. You mean frequency resolution. Specifically, you mean the spacing between DFT bins is 100 Hz.

Have you seen my suggestion of goertzel.m?

##### 0 Comments

Dr. Seis
on 31 Oct 2012

A possible workaround since you want to compute and retain ALL frequencies, but only plot a subset. You can try plotting a subset of your data so that you keep more low frequency points and ignore more high frequency points by:

N = 1000; % Number of samples you want to plot

f_begin = 1e2; f_end = 1e8;

log10_df = (log10(f_end)-log10(f_begin))/N;

new_f = 10 .^ (log10(f_begin) : log10_df : log10(f_end));

Then plot:

plot(new_f,interp1(OLD_Frequencies_1e2_to_1e8, OLD_Amplitudes_1e2_to_1e8, new_f)

If you don't want to interpolate, then you can either try to find the index associated with the closest matching frequency OR possibly try to take an average of the amplitude data around the matching frequency.

##### 0 Comments

### See Also

### Categories

### Products

### Community Treasure Hunt

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

Start Hunting!