hashtags/lmms
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

New Spectrum Analyzer (#4950)

Browse Source 
Replace old spectrum analyzer by new one with higher resolution and
many new features.

Resolves #2847.
This commit is contained in:
Martin Pavelek
2019-07-17 22:45:26 +02:00
committed by Johannes Lorenz
parent 73c2c70d96
commit c3b4d5188a
42 changed files with 5329 additions and 751 deletions
Download Patch File Download Diff File Expand all files Collapse all files

277
src/core/fft_helpers.cpp
View File

@@ -2,6 +2,7 @@
* fft_helpers.cpp - some functions around FFT analysis
*
* Copyright (c) 2008-2012 Tobias Doerffel <tobydox/at/users.sourceforge.net>
* Copyright (c) 2019 Martin Pavelek <he29.HS/at/gmail.com>
*
* This file is part of LMMS - https://lmms.io
*
@@ -28,127 +29,194 @@
#include <cmath>
#include "lmms_constants.h"
/* returns biggest value from abs_spectrum[spec_size] array
returns -1 on error
*/
float maximum(float *abs_spectrum, unsigned int spec_size)
/* Returns biggest value from abs_spectrum[spec_size] array.
*
* return -1 on error, otherwise the maximum value
*/
float maximum(const float *abs_spectrum, unsigned int spec_size)
{
float maxi=0;
float maxi = 0;
unsigned int i;
if ( abs_spectrum==NULL )
return -1;
if (abs_spectrum == NULL) {return -1;}
if (spec_size <= 0) {return -1;}
if (spec_size<=0)
return -1;
for ( i=0; i<spec_size; i++ )
for (i = 0; i < spec_size; i++)
{
if ( abs_spectrum[i]>maxi )
maxi=abs_spectrum[i];
if (abs_spectrum[i] > maxi) {maxi = abs_spectrum[i];}
}
return maxi;
}
float maximum(const std::vector<float> &abs_spectrum)
{
return maximum(abs_spectrum.data(), abs_spectrum.size());
}
/* apply hanning or hamming window to channel
returns -1 on error */
int hanming(float *timebuffer, int length, WINDOWS type)
/* Normalize the array of absolute magnitudes to a 0..1 range.
* Block size refers to FFT block size before any zero padding.
*
* return -1 on error, 0 on success
*/
int normalize(const float *abs_spectrum, float *norm_spectrum, unsigned int bin_count, unsigned int block_size)
{
int i;
float alpha;
if ( (timebuffer==NULL)||(length<=0) )
return -1;
if (abs_spectrum == NULL || norm_spectrum == NULL) {return -1;}
if (bin_count == 0 || block_size == 0) {return -1;}
for (i = 0; i < bin_count; i++)
{
norm_spectrum[i] = abs_spectrum[i] / block_size;
}
return 0;
}
int normalize(const std::vector<float> &abs_spectrum, std::vector<float> &norm_spectrum, unsigned int block_size)
{
if (abs_spectrum.size() != norm_spectrum.size()) {return -1;}
return normalize(abs_spectrum.data(), norm_spectrum.data(), abs_spectrum.size(), block_size);
}
/* Check if the spectrum contains any non-zero value.
*
* return 1 if spectrum contains any non-zero value
* return 0 otherwise
*/
int notEmpty(const std::vector<float> &spectrum)
{
for (int i = 0; i < spectrum.size(); i++)
{
if (spectrum[i] != 0) {return 1;}
}
return 0;
}
/* Precompute an FFT window function for later real-time use.
*
* return -1 on error
*/
int precomputeWindow(float *window, unsigned int length, FFT_WINDOWS type, bool normalized)
{
unsigned int i;
float gain = 0;
float a0;
float a1;
float a2;
float a3;
if (window == NULL) {return -1;}
// constants taken from
// https://en.wikipedia.org/wiki/Window_function#AList_of_window_functions
switch (type)
{
case HAMMING: alpha=0.54; break;
default:
case RECTANGULAR:
for (i = 0; i < length; i++) {window[i] = 1.0;}
gain = 1;
return 0;
case BLACKMAN_HARRIS:
a0 = 0.35875;
a1 = 0.48829;
a2 = 0.14128;
a3 = 0.01168;
break;
case HAMMING:
a0 = 0.54;
a1 = 1.0 - a0;
a2 = 0;
a3 = 0;
break;
case HANNING:
default: alpha=0.5; break;
a0 = 0.5;
a1 = 1.0 - a0;
a2 = 0;
a3 = 0;
break;
}
for ( i=0; i<length; i++ )
// common computation for cosine-sum based windows
for (i = 0; i < length; i++)
{
timebuffer[i]=timebuffer[i]*(alpha+(1-alpha)*cos(2*F_PI*i/((float)length-1.0)));
window[i] = (a0 - a1 * cos(2 * F_PI * i / ((float)length - 1.0))
+ a2 * cos(4 * F_PI * i / ((float)length - 1.0))
- a3 * cos(6 * F_PI * i / ((float)length - 1.0)));
gain += window[i];
}
// apply amplitude correction
gain /= (float) length;
for (i = 0; i < length; i++) {window[i] /= gain;}
return 0;
}
/* compute absolute values of complex_buffer, save to absspec_buffer
take care that - compl_len is not bigger than complex_buffer!
- absspec buffer is big enough!
returns 0 on success, else -1 */
int absspec(fftwf_complex *complex_buffer, float *absspec_buffer, int compl_length)
/* Compute absolute values of complex_buffer, save to absspec_buffer.
* Take care that - compl_len is not bigger than complex_buffer!
* - absspec buffer is big enough!
*
* return 0 on success, else -1
*/
int absspec(const fftwf_complex *complex_buffer, float *absspec_buffer, unsigned int compl_length)
{
int i;
if ( (complex_buffer==NULL)||(absspec_buffer==NULL) )
return -1;
if ( compl_length<=0 )
return -1;
if (complex_buffer == NULL || absspec_buffer == NULL) {return -1;}
if (compl_length <= 0) {return -1;}
for (i=0; i<compl_length; i++)
for (i = 0; i < compl_length; i++)
{
absspec_buffer[i]=(float )sqrt(complex_buffer[i][0]*complex_buffer[i][0] + complex_buffer[i][1]*complex_buffer[i][1]);
absspec_buffer[i] = (float)sqrt(complex_buffer[i][0] * complex_buffer[i][0]
+ complex_buffer[i][1] * complex_buffer[i][1]);
}
return 0;
}
/* build fewer subbands from many absolute spectrum values
take care that - compressedbands[] array num_new elements long
- num_old > num_new
returns 0 on success, else -1 */
int compressbands(float *absspec_buffer, float *compressedband, int num_old, int num_new, int bottom, int top)
/* Build fewer subbands from many absolute spectrum values.
* Take care that - compressedbands[] array num_new elements long
* - num_old > num_new
*
* return 0 on success, else -1
*/
int compressbands(const float *absspec_buffer, float *compressedband, int num_old, int num_new, int bottom, int top)
{
float ratio;
int i, usefromold;
float j;
float j_min, j_max;
if ( (absspec_buffer==NULL)||(compressedband==NULL) )
return -1;
if (absspec_buffer == NULL || compressedband == NULL) {return -1;}
if (num_old < num_new) {return -1;}
if (num_old <= 0 || num_new <= 0) {return -1;}
if (bottom < 0) {bottom = 0;}
if (top >= num_old) {top = num_old - 1;}
if ( num_old<num_new )
return -1;
usefromold = num_old - (num_old - top) - bottom;
if ( (num_old<=0)||(num_new<=0) )
return -1;
if ( bottom<0 )
bottom=0;
if ( top>=num_old )
top=num_old-1;
usefromold=num_old-(num_old-top)-bottom;
ratio=(float)usefromold/(float)num_new;
ratio = (float)usefromold / (float)num_new;
// for each new subband
for ( i=0; i<num_new; i++ )
for (i = 0; i < num_new; i++)
{
compressedband[i]=0;
compressedband[i] = 0;
j_min=(i*ratio)+bottom;
j_min = (i * ratio) + bottom;
if ( j_min<0 )
j_min=bottom;
if (j_min < 0) {j_min = bottom;}
j_max=j_min+ratio;
j_max = j_min + ratio;
for ( j=(int)j_min; j<=j_max; j++ )
for (j = (int)j_min; j <= j_max; j++)
{
compressedband[i]+=absspec_buffer[(int)j];
compressedband[i] += absspec_buffer[(int)j];
}
}
@@ -158,84 +226,73 @@ int compressbands(float *absspec_buffer, float *compressedband, int num_old, int
int calc13octaveband31(float *absspec_buffer, float *subbands, int num_spec, float max_frequency)
{
static const int onethirdoctavecenterfr[] = {20, 25, 31, 40, 50, 63, 80, 100, 125, 160, 200, 250, 315, 400, 500, 630, 800, 1000, 1250, 1600, 2000, 2500, 3150, 4000, 5000, 6300, 8000, 10000, 12500, 16000, 20000};
static const int onethirdoctavecenterfr[] = {20, 25, 31, 40, 50, 63, 80, 100, 125, 160, 200, 250, 315, 400, 500, 630, 800, 1000, 1250, 1600, 2000, 2500, 3150, 4000, 5000, 6300, 8000, 10000, 12500, 16000, 20000};
int i, j;
float f_min, f_max, frequency, bandwidth;
int j_min, j_max=0;
int j_min, j_max = 0;
float fpower;
if ( (absspec_buffer==NULL)||(subbands==NULL) )
return -1;
if ( num_spec<31 )
return -1;
if ( max_frequency<=0 )
return -1;
if (absspec_buffer == NULL || subbands == NULL) {return -1;}
if (num_spec < 31) {return -1;}
if (max_frequency <= 0) {return -1;}
/*** energy ***/
fpower=0;
for ( i=0; i<num_spec; i++ )
fpower = 0;
for (i = 0; i < num_spec; i++)
{
absspec_buffer[i]=(absspec_buffer[i]*absspec_buffer[i])/FFT_BUFFER_SIZE;
fpower=fpower+(2*absspec_buffer[i]);
absspec_buffer[i] = (absspec_buffer[i] * absspec_buffer[i]) / FFT_BUFFER_SIZE;
fpower = fpower + (2 * absspec_buffer[i]);
}
fpower=fpower-(absspec_buffer[0]); //dc not mirrored
fpower = fpower - (absspec_buffer[0]); //dc not mirrored
/*** for each subband: sum up power ***/
for ( i=0; i<31; i++ )
for (i = 0; i < 31; i++)
{
subbands[i]=0;
subbands[i] = 0;
// calculate bandwidth for subband
frequency=onethirdoctavecenterfr[i];
frequency = onethirdoctavecenterfr[i];
bandwidth=(pow(2, 1.0/3.0)-1)*frequency;
bandwidth = (pow(2, 1.0/3.0)-1) * frequency;
f_min=frequency-bandwidth/2.0;
f_max=frequency+bandwidth/2.0;
f_min = frequency - bandwidth / 2.0;
f_max = frequency + bandwidth / 2.0;
j_min=(int)(f_min/max_frequency*(float)num_spec);
j_max=(int)(f_max/max_frequency*(float)num_spec);
j_min = (int)(f_min / max_frequency * (float)num_spec);
j_max = (int)(f_max / max_frequency * (float)num_spec);
if ( (j_min<0)||(j_max<0) )
if (j_min < 0 || j_max < 0)
{
fprintf(stderr, "Error: calc13octaveband31() in fft_helpers.cpp line %d failed.\n", __LINE__);
return -1;
}
for ( j=j_min; j<=j_max; j++ )
for (j = j_min; j <= j_max; j++)
{
if( j_max<num_spec )
subbands[i]+=absspec_buffer[j];
if (j_max < num_spec) {subbands[i] += absspec_buffer[j];}
}
} //for
return 0;
}
/* compute power of finite time sequence
take care num_values is length of timesignal[]
returns power on success, else -1 */
float signalpower(float *timesignal, int num_values)
/* Compute power of finite time sequence.
* Take care num_values is length of timesignal[]
*
* return power on success, else -1
*/
float signalpower(const float *timesignal, int num_values)
{
if ( num_values<=0 )
return -1;
if (num_values <= 0) {return -1;}
if( timesignal==NULL )
return -1;
if (timesignal == NULL) {return -1;}
float power=0;
for ( int i=0; i<num_values; i++ )
float power = 0;
for (int i = 0; i < num_values; i++)
{
power+=timesignal[i]*timesignal[i];
power += timesignal[i] * timesignal[i];
}
return power;

1
src/gui/MainWindow.cpp
View File

@@ -625,6 +625,7 @@ SubWindow* MainWindow::addWindowedWidget(QWidget *w, Qt::WindowFlags windowFlags
SubWindow *win = new SubWindow(m_workspace->viewport(), windowFlags);
win->setAttribute(Qt::WA_DeleteOnClose);
win->setWidget(w);
if (w->sizeHint().isValid()) {win->resize(w->sizeHint());}
m_workspace->addSubWindow(win);
return win;
}

1
src/gui/SubWindow.cpp
View File

@@ -30,6 +30,7 @@
#include <QMdiArea>
#include <QMoveEvent>
#include <QPainter>
#include <QScrollBar>
#include "embed.h"

11
src/gui/widgets/EffectView.cpp
View File

@@ -97,9 +97,14 @@ EffectView::EffectView( Effect * _model, QWidget * _parent ) :
if( m_controlView )
{
m_subWindow = gui->mainWindow()->addWindowedWidget( m_controlView );
m_subWindow->setSizePolicy( QSizePolicy::Fixed, QSizePolicy::Fixed );
if (m_subWindow->layout()) {
m_subWindow->layout()->setSizeConstraint(QLayout::SetFixedSize);
if ( !m_controlView->isResizable() )
{
m_subWindow->setSizePolicy( QSizePolicy::Fixed, QSizePolicy::Fixed );
if (m_subWindow->layout())
{
m_subWindow->layout()->setSizeConstraint(QLayout::SetFixedSize);
}
}
Qt::WindowFlags flags = m_subWindow->windowFlags();

4
src/gui/widgets/PixmapButton.cpp
View File

@@ -134,11 +134,11 @@ QSize PixmapButton::sizeHint() const
{
if( ( model() != NULL && model()->value() ) || m_pressed )
{
return m_activePixmap.size();
return m_activePixmap.size() / devicePixelRatio();
}
else
{
return m_inactivePixmap.size();
return m_inactivePixmap.size() / devicePixelRatio();
}
}