dspl/dspl_8c_source.html

/*

* Copyright (c) 2015-2022 Sergey Bakhurin

* Digital Signal Processing Library [http://dsplib.org]

*

* This file is part of libdspl-2.0.

*

* is free software: you can redistribute it and/or modify

* it under the terms of the GNU Lesser  General Public License as published by

* the Free Software Foundation, either version 3 of the License, or

* (at your option) any later version.

*

* DSPL is distributed in the hope that it will be useful,

* but WITHOUT ANY WARRANTY; without even the implied warranty of

* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

* GNU General Public License for more details.

*

* You should have received a copy of the GNU Lesser General Public License

* along with Foobar.  If not, see <http://www.gnu.org/licenses/>.

*/


#ifdef WIN_OS

#include <windows.h>

#endif  /* WIN_OS */


#ifdef LINUX_OS

#include <dlfcn.h>

#endif /* LINUX_OS */


#include <stdio.h>

#include "dspl.h"


#ifndef BUILD_LIB


p_acos_cmplx                            acos_cmplx                    ;

p_addlog                                addlog                        ;

p_array_scale_lin                       array_scale_lin               ;

p_asin_cmplx                            asin_cmplx                    ;


p_butter_ap                             butter_ap                     ;

p_bessel_i0                             bessel_i0                     ;

p_bilinear                              bilinear                      ;

p_butter_ap_zp                          butter_ap_zp                  ;


p_cheby_poly1                           cheby_poly1                   ;

p_cheby_poly2                           cheby_poly2                   ;

p_cheby1_ap                             cheby1_ap                     ;

p_cheby1_ap_zp                          cheby1_ap_zp                  ;

p_cheby2_ap                             cheby2_ap                     ;

p_cheby2_ap_wp1                         cheby2_ap_wp1                 ;

p_cheby2_ap_zp                          cheby2_ap_zp                  ;

p_cmplx2re                              cmplx2re                      ;

p_concat                                concat                        ;

p_conv                                  conv                          ;

p_conv_cmplx                            conv_cmplx                    ;

p_conv_fft                              conv_fft                      ;

p_conv_fft_cmplx                        conv_fft_cmplx                ;

p_cos_cmplx                             cos_cmplx                     ;


p_decimate                              decimate                      ;

p_decimate_cmplx                        decimate_cmplx                ;

p_dft                                   dft                           ;

p_dft_cmplx                             dft_cmplx                     ;

p_dmod                                  dmod                          ;

p_dspl_info                             dspl_info                     ;


p_ellip_acd                             ellip_acd                     ;

p_ellip_acd_cmplx                       ellip_acd_cmplx               ;

p_ellip_ap                              ellip_ap                      ;

p_ellip_ap_zp                           ellip_ap_zp                   ;

p_ellip_asn                             ellip_asn                     ;

p_ellip_asn_cmplx                       ellip_asn_cmplx               ;

p_ellip_cd                              ellip_cd                      ;

p_ellip_cd_cmplx                        ellip_cd_cmplx                ;

p_ellip_landen                          ellip_landen                  ;

p_ellip_modulareq                       ellip_modulareq               ;

p_ellip_rat                             ellip_rat                     ;

p_ellip_sn                              ellip_sn                      ;

p_ellip_sn_cmplx                        ellip_sn_cmplx                ;


p_farrow_lagrange                       farrow_lagrange               ;

p_farrow_spline                         farrow_spline                 ;

p_fft                                   fft                           ;

p_fft_abs                               fft_abs                       ;

p_fft_abs_cmplx                         fft_abs_cmplx                 ;

p_fft_cmplx                             fft_cmplx                     ;

p_fft_create                            fft_create                    ;

p_fft_free                              fft_free                      ;

p_fft_mag                               fft_mag                       ;

p_fft_mag_cmplx                         fft_mag_cmplx                 ;

p_fft_shift                             fft_shift                     ;

p_fft_shift_cmplx                       fft_shift_cmplx               ;

p_filter_freq_resp                      filter_freq_resp              ;

p_filter_iir                            filter_iir                    ;

p_filter_ws1                            filter_ws1                    ;

p_filter_zp2ab                          filter_zp2ab                  ;

p_find_max_abs                          find_max_abs                  ;

p_find_nearest                          find_nearest                  ;

p_fir_linphase                          fir_linphase                  ;

p_flipip                                flipip                        ;

p_flipip_cmplx                          flipip_cmplx                  ;

p_fourier_integral_cmplx                fourier_integral_cmplx        ;

p_fourier_series_dec                    fourier_series_dec            ;

p_fourier_series_dec_cmplx              fourier_series_dec_cmplx      ;

p_fourier_series_rec                    fourier_series_rec            ;

p_freqs                                 freqs                         ;

p_freqs_cmplx                           freqs_cmplx                   ;

p_freqs2time                            freqs2time                    ;

p_freqz                                 freqz                         ;


p_gnuplot_close                         gnuplot_close                 ;

p_gnuplot_cmd                           gnuplot_cmd                   ;

p_gnuplot_create                        gnuplot_create                ;

p_gnuplot_open                          gnuplot_open                  ;

p_goertzel                              goertzel                      ;

p_goertzel_cmplx                        goertzel_cmplx                ;

p_group_delay                           group_delay                   ;


p_histogram                             histogram                     ;

p_histogram_norm                        histogram_norm                ;


p_idft_cmplx                            idft_cmplx                    ;

p_ifft_cmplx                            ifft_cmplx                    ;

p_iir                                   iir                           ;


p_linspace                              linspace                      ;

p_log_cmplx                             log_cmplx                     ;

p_logspace                              logspace                      ;

p_low2bp                                low2bp                        ;

p_low2bs                                low2bs                        ;

p_low2high                              low2high                      ;

p_low2low                               low2low                       ;


p_matrix_eig_cmplx                      matrix_eig_cmplx              ;

p_matrix_eye                            matrix_eye                    ;

p_matrix_eye_cmplx                      matrix_eye_cmplx              ;

p_matrix_mul                            matrix_mul                    ;

p_matrix_pinv                           matrix_pinv                   ;

p_matrix_print                          matrix_print                  ;

p_matrix_print_cmplx                    matrix_print_cmplx            ;

p_matrix_svd                            matrix_svd                    ;

p_matrix_transpose                      matrix_transpose              ;

p_matrix_transpose_cmplx                matrix_transpose_cmplx        ;

p_matrix_transpose_hermite              matrix_transpose_hermite      ;

p_mean                                  mean                          ;

p_mean_cmplx                            mean_cmplx                    ;

p_minmax                                minmax                        ;


p_ones                                  ones                          ;


p_phase_delay                           phase_delay                   ;

p_poly_z2a_cmplx                        poly_z2a_cmplx                ;

p_polyroots                             polyroots                     ;

p_polyval                               polyval                       ;

p_polyval_cmplx                         polyval_cmplx                 ;

p_psd_bartlett                          psd_bartlett                  ;

p_psd_bartlett_cmplx                    psd_bartlett_cmplx            ;

p_psd_periodogram                       psd_periodogram               ;

p_psd_periodogram_cmplx                 psd_periodogram_cmplx         ;

p_psd_welch                             psd_welch                     ;

p_psd_welch_cmplx                       psd_welch_cmplx               ;


p_randb                                 randb                         ;

p_randb2                                randb2                        ;

p_randi                                 randi                         ;

p_randn                                 randn                         ;

p_randn_cmplx                           randn_cmplx                   ;

p_random_init                           random_init                   ;

p_randu                                 randu                         ;

p_ratcompos                             ratcompos                     ;

p_re2cmplx                              re2cmplx                      ;

p_readbin                               readbin                       ;


p_signal_pimp                           signal_pimp                   ;

p_signal_saw                            signal_saw                    ;

p_sin_cmplx                             sin_cmplx                     ;

p_sinc                                  sinc                          ;

p_sine_int                              sine_int                      ;

p_sqrt_cmplx                            sqrt_cmplx                    ;

p_stat_std                              stat_std                      ;

p_stat_std_cmplx                        stat_std_cmplx                ;

p_sum                                   sum                           ;

p_sum_sqr                               sum_sqr                       ;


p_trapint                               trapint                       ;

p_trapint_cmplx                         trapint_cmplx                 ;


p_unwrap                                unwrap                        ;


p_vector_dot                            vector_dot                    ;

p_verif                                 verif                         ;

p_verif_data_gen                        verif_data_gen                ;

p_verif_cmplx                           verif_cmplx                   ;

p_verif_str                             verif_str                     ;

p_verif_str_cmplx                       verif_str_cmplx               ;


p_window                                window                        ;

p_writebin                              writebin                      ;

p_writetxt                              writetxt                      ;

p_writetxt_3d                           writetxt_3d                   ;

p_writetxt_3dline                       writetxt_3dline               ;

p_writetxt_cmplx                        writetxt_cmplx                ;

p_writetxt_cmplx_im                     writetxt_cmplx_im             ;

p_writetxt_cmplx_re                     writetxt_cmplx_re             ;

p_writetxt_int                          writetxt_int                  ;


p_xcorr                                 xcorr                         ;

p_xcorr_cmplx                           xcorr_cmplx                   ;


#ifdef WIN_OS

#define LOAD_FUNC(fn) \

        fname = #fn;\

        fn = (p_##fn)GetProcAddress(handle, fname);\

        if(! fn) goto exit_label;

#endif


#ifdef LINUX_OS

#define LOAD_FUNC(fn) \

        fname = #fn;\

        fn = (p_##fn)dlsym(handle, fname);\

        if ((error = dlerror()) != NULL) goto exit_label

#endif


void* dspl_load()

{

    char* fname;

    #ifdef WIN_OS

        HINSTANCE handle;

        handle = LoadLibrary(TEXT("libdspl.dll"));

        if (!handle)

        {

            printf("libdspl.dll loading ERROR!\n");

            return NULL;

        }

    #endif /* WIN_OS */


    #ifdef LINUX_OS

        char* error;

        void *handle;

        /* open the *.so */

        handle = dlopen ("./libdspl.so", RTLD_LAZY);

        if (!handle)

        {

            printf("libdspl.so loading ERROR!\n");

            return NULL;

        }

    #endif  /* LINUX_OS */


    LOAD_FUNC(acos_cmplx);

    LOAD_FUNC(addlog);

    LOAD_FUNC(array_scale_lin);

    LOAD_FUNC(asin_cmplx);


    LOAD_FUNC(bessel_i0);

    LOAD_FUNC(bilinear);

    LOAD_FUNC(butter_ap);

    LOAD_FUNC(butter_ap_zp);


    LOAD_FUNC(cheby_poly1);

    LOAD_FUNC(cheby_poly2);

    LOAD_FUNC(cheby1_ap);

    LOAD_FUNC(cheby1_ap_zp);

    LOAD_FUNC(cheby2_ap);

    LOAD_FUNC(cheby2_ap_wp1);

    LOAD_FUNC(cheby2_ap_zp);

    LOAD_FUNC(cmplx2re);

    LOAD_FUNC(concat);

    LOAD_FUNC(conv);

    LOAD_FUNC(conv_cmplx);

    LOAD_FUNC(conv_fft);

    LOAD_FUNC(conv_fft_cmplx);

    LOAD_FUNC(cos_cmplx);


    LOAD_FUNC(decimate);

    LOAD_FUNC(decimate_cmplx);

    LOAD_FUNC(dft);

    LOAD_FUNC(dft_cmplx);

    LOAD_FUNC(dmod);

    LOAD_FUNC(dspl_info);


    LOAD_FUNC(ellip_acd);

    LOAD_FUNC(ellip_acd_cmplx);

    LOAD_FUNC(ellip_ap);

    LOAD_FUNC(ellip_ap_zp);

    LOAD_FUNC(ellip_asn);

    LOAD_FUNC(ellip_asn_cmplx);

    LOAD_FUNC(ellip_cd);

    LOAD_FUNC(ellip_cd_cmplx);

    LOAD_FUNC(ellip_landen);

    LOAD_FUNC(ellip_modulareq);

    LOAD_FUNC(ellip_rat);

    LOAD_FUNC(ellip_sn);

    LOAD_FUNC(ellip_sn_cmplx);


    LOAD_FUNC(farrow_lagrange);

    LOAD_FUNC(farrow_spline);

    LOAD_FUNC(fft);

    LOAD_FUNC(fft_cmplx);

    LOAD_FUNC(fft_create);

    LOAD_FUNC(fft_free);

    LOAD_FUNC(fft_mag);

    LOAD_FUNC(fft_mag_cmplx);

    LOAD_FUNC(fft_shift);

    LOAD_FUNC(fft_shift_cmplx);

    LOAD_FUNC(filter_freq_resp);

    LOAD_FUNC(filter_iir);

    LOAD_FUNC(filter_ws1);

    LOAD_FUNC(filter_zp2ab);

    LOAD_FUNC(find_max_abs);

    LOAD_FUNC(find_nearest);

    LOAD_FUNC(fir_linphase);

    LOAD_FUNC(flipip);

    LOAD_FUNC(flipip_cmplx);

    LOAD_FUNC(fourier_integral_cmplx);

    LOAD_FUNC(fourier_series_dec);

    LOAD_FUNC(fourier_series_dec_cmplx);

    LOAD_FUNC(fourier_series_rec);

    LOAD_FUNC(freqz);

    LOAD_FUNC(freqs);

    LOAD_FUNC(freqs_cmplx);

    LOAD_FUNC(freqs2time);


    LOAD_FUNC(gnuplot_close);

    LOAD_FUNC(gnuplot_cmd);

    LOAD_FUNC(gnuplot_create);

    LOAD_FUNC(gnuplot_open);

    LOAD_FUNC(goertzel);

    LOAD_FUNC(goertzel_cmplx);

    LOAD_FUNC(group_delay);


    LOAD_FUNC(histogram);

    LOAD_FUNC(histogram_norm);


    LOAD_FUNC(idft_cmplx);

    LOAD_FUNC(ifft_cmplx);

    LOAD_FUNC(iir);


    LOAD_FUNC(linspace);

    LOAD_FUNC(log_cmplx);

    LOAD_FUNC(logspace);

    LOAD_FUNC(low2bp);

    LOAD_FUNC(low2bs);

    LOAD_FUNC(low2high);

    LOAD_FUNC(low2low);


    LOAD_FUNC(matrix_eig_cmplx);

    LOAD_FUNC(matrix_eye);

    LOAD_FUNC(matrix_eye_cmplx);

    LOAD_FUNC(matrix_mul);

    LOAD_FUNC(matrix_pinv);

    LOAD_FUNC(matrix_print);

    LOAD_FUNC(matrix_print_cmplx);

    LOAD_FUNC(matrix_svd);

    LOAD_FUNC(matrix_transpose);

    LOAD_FUNC(matrix_transpose_cmplx);

    LOAD_FUNC(matrix_transpose_hermite);

    LOAD_FUNC(mean);

    LOAD_FUNC(mean_cmplx);

    LOAD_FUNC(minmax);


    LOAD_FUNC(ones);


    LOAD_FUNC(phase_delay);

    LOAD_FUNC(poly_z2a_cmplx);

    LOAD_FUNC(polyroots);

    LOAD_FUNC(polyval);

    LOAD_FUNC(polyval_cmplx);

    LOAD_FUNC(psd_bartlett);

    LOAD_FUNC(psd_bartlett_cmplx);

    LOAD_FUNC(psd_periodogram);

    LOAD_FUNC(psd_periodogram_cmplx);

    LOAD_FUNC(psd_welch);

    LOAD_FUNC(psd_welch_cmplx);


    LOAD_FUNC(randi);

    LOAD_FUNC(randb);

    LOAD_FUNC(randb2);

    LOAD_FUNC(randn);

    LOAD_FUNC(randn_cmplx);

    LOAD_FUNC(random_init);

    LOAD_FUNC(randu);

    LOAD_FUNC(ratcompos);

    LOAD_FUNC(re2cmplx);

    LOAD_FUNC(readbin);


    LOAD_FUNC(signal_pimp);

    LOAD_FUNC(signal_saw);

    LOAD_FUNC(sin_cmplx);

    LOAD_FUNC(sinc);

    LOAD_FUNC(sine_int);

    LOAD_FUNC(sqrt_cmplx);

    LOAD_FUNC(stat_std);

    LOAD_FUNC(stat_std_cmplx);

    LOAD_FUNC(sum);

    LOAD_FUNC(sum_sqr);


    LOAD_FUNC(trapint);

    LOAD_FUNC(trapint_cmplx);


    LOAD_FUNC(unwrap);


    LOAD_FUNC(vector_dot);

    LOAD_FUNC(verif);

    LOAD_FUNC(verif_data_gen);

    LOAD_FUNC(verif_cmplx);

    LOAD_FUNC(verif_str);

    LOAD_FUNC(verif_str_cmplx);


    LOAD_FUNC(window);

    LOAD_FUNC(writebin);

    LOAD_FUNC(writetxt);

    LOAD_FUNC(writetxt_3d);

    LOAD_FUNC(writetxt_3dline);

    LOAD_FUNC(writetxt_cmplx);

    LOAD_FUNC(writetxt_cmplx_im);

    LOAD_FUNC(writetxt_cmplx_re);

    LOAD_FUNC(writetxt_int);


    LOAD_FUNC(xcorr);

    LOAD_FUNC(xcorr_cmplx);


    #ifdef WIN_OS

    return (void*)handle;

    exit_label:

        printf("function %s loading ERROR!\n", fname);

        if(handle)

            FreeLibrary(handle);

        return NULL;

    #endif /* WIN_OS */


    #ifdef LINUX_OS

      return handle;

    exit_label:

        printf("function %s loading ERROR!\n", fname);

        if(handle)

            dlclose(handle);

        return NULL;

    #endif /* LINUX_OS */

}


void dspl_free(void* handle)

{

    #ifdef WIN_OS

        FreeLibrary((HINSTANCE)handle);

    #endif /* WIN_OS */


    #ifdef LINUX_OS

        dlclose(handle);

    #endif /* LINUX_OS */

}


#endif /* BUILD_LIB */

decimate_cmplx
int DSPL_API decimate_cmplx(complex_t *x, int n, int d, complex_t *y, int *cnt)
Complex vector decimation.
Definition: decimate_cmplx.c:144

linspace
int DSPL_API linspace(double x0, double x1, int n, int type, double *x)
Function fills a vector with n linearly spaced elements between x0 and x1.
Definition: linspace.c:169

concat
int DSPL_API concat(void *a, size_t na, void *b, size_t nb, void *c)
Concatenate arrays a and b
Definition: concat.c:146

verif
int DSPL_API verif(double *x, double *y, size_t n, double eps, double *err)
Real arrays verification.
Definition: verif.c:130

ones
int DSPL_API ones(double *x, int n)
Function fills all real vector x by ones values.
Definition: ones.c:103

verif_cmplx
int DSPL_API verif_cmplx(complex_t *x, complex_t *y, size_t n, double eps, double *err)
Complex arrays verification.
Definition: verif_cmplx.c:141

decimate
int DSPL_API decimate(double *x, int n, int d, double *y, int *cnt)
Real vector decimation.
Definition: decimate.c:140

flipip_cmplx
int DSPL_API flipip_cmplx(complex_t *x, int n)
Flip complex vector x in place.
Definition: flipip_cmplx.c:132

array_scale_lin
int DSPL_API array_scale_lin(double *x, int n, double xmin, double xmax, double dx, double h, double *y)
Vector x linear transformation.
Definition: array_scale_lin.c:171

logspace
int DSPL_API logspace(double x0, double x1, int n, int type, double *x)
Function fills a vector with n logarithmically spaced elements between  and .
Definition: logspace.c:178

sum
int DSPL_API sum(double *x, int n, double *s)
Definition: sum.c:43

flipip
int DSPL_API flipip(double *x, int n)
Flip real vector x in place.
Definition: flipip.c:133

fourier_series_dec
int DSPL_API fourier_series_dec(double *t, double *s, int nt, double period, int nw, double *w, complex_t *y)
Fourier series coefficient calculation for periodic signal.
Definition: fourier_series_dec.c:152

idft_cmplx
int DSPL_API idft_cmplx(complex_t *x, int n, complex_t *y)
Inverse discrete Fourier transform of the complex spectrum.
Definition: idft_cmplx.c:163

fft_free
void DSPL_API fft_free(fft_t *pfft)
Free fft_t structure.
Definition: fft_free.c:63

fft_cmplx
int DSPL_API fft_cmplx(complex_t *x, int n, fft_t *pfft, complex_t *y)
Fast Fourier transform for the complex vector.
Definition: fft_cmplx.c:179

fft_create
int DSPL_API fft_create(fft_t *pfft, int n)
Function creates and fill fft_t structure.
Definition: fft_create.c:161

ifft_cmplx
int DSPL_API ifft_cmplx(complex_t *x, int n, fft_t *pfft, complex_t *y)
Inverse fast Fourier transform.
Definition: ifft_cmplx.c:179

fft_shift
int DSPL_API fft_shift(double *x, int n, double *y)
Perform a shift of the vector x, for use with the fft and ifft functions, in order to move the freque...
Definition: fft_shift.c:95

fft
int DSPL_API fft(double *x, int n, fft_t *pfft, complex_t *y)
Fast Fourier transform for the real vector.
Definition: fft.c:173

dft
int DSPL_API dft(double *x, int n, complex_t *y)
Discrete Fourier transform of a real signal.
Definition: dft.c:163

goertzel
int DSPL_API goertzel(double *x, int n, int *ind, int k, complex_t *y)
Goertzel algorithm  individual DFT samples calculation for the real input vector x.
Definition: goertzel.c:125

fourier_series_rec
int DSPL_API fourier_series_rec(double *w, complex_t *s, int nw, double *t, int nt, complex_t *y)
Time signal reconstruction from Fourier series coefficients.
Definition: fourier_series_rec.c:158

dft_cmplx
int DSPL_API dft_cmplx(complex_t *x, int n, complex_t *y)
Discrete Fourier transform of a complex signal.
Definition: dft_cmplx.c:163

goertzel_cmplx
int DSPL_API goertzel_cmplx(complex_t *x, int n, int *ind, int k, complex_t *y)
Goertzel algorithm  individual DFT samples calculation for the complex input vector x.
Definition: goertzel_cmplx.c:127

phase_delay
int DSPL_API phase_delay(double *b, double *a, int ord, int flag, double *w, int n, double *tau)
Phase delay calculation for digital or analog filter corresponds to , or  transfer function.
Definition: phase_delay.c:165

group_delay
int DSPL_API group_delay(double *pb, double *pa, int ord, int flag, double *w, int n, double *tau)
Group delay calculation for digital or analog filter corresponds to , or  transfer function.
Definition: group_delay.c:164

freqs
int DSPL_API freqs(double *b, double *a, int ord, double *w, int n, complex_t *h)
Analog filter frequency response  calculation.
Definition: freqs.c:143

freqz
int DSPL_API freqz(double *b, double *a, int ord, double *w, int n, complex_t *h)
Function calculates the digital filter frequency response  corresponds to transfer function .
Definition: freqz.c:157

filter_freq_resp
int DSPL_API filter_freq_resp(double *b, double *a, int ord, double *w, int n, int flag, double *mag, double *phi, double *tau)
Magnitude, phase response and group delay vectors calculation for digital or analog filter correspond...
Definition: filter_freq_resp.c:235

filter_iir
int DSPL_API filter_iir(double *b, double *a, int ord, double *x, int n, double *y)
Real IIR filtration.
Definition: filter_iir.c:167

conv_cmplx
int DSPL_API conv_cmplx(complex_t *a, int na, complex_t *b, int nb, complex_t *c)
Complex vectors linear convolution.
Definition: conv_cmplx.c:160

conv
int DSPL_API conv(double *a, int na, double *b, int nb, double *c)
Real vectors linear convolution.
Definition: conv.c:159

conv_fft_cmplx
int DSPL_API conv_fft_cmplx(complex_t *a, int na, complex_t *b, int nb, fft_t *pfft, int nfft, complex_t *c)
Complex vectors fast linear convolution by using fast Fourier transform algorithms.
Definition: conv_fft_cmplx.c:192

conv_fft
int DSPL_API conv_fft(double *a, int na, double *b, int nb, fft_t *pfft, int nfft, double *c)
Real vectors fast linear convolution by using fast Fourier transform algorithms.
Definition: conv_fft.c:202

fir_linphase
int DSPL_API fir_linphase(int ord, double w0, double w1, int filter_type, int win_type, double win_param, double *h)
Function calculates linear-phase FIR filter coefficients by window method.
Definition: fir_linphase.c:323

bilinear
int DSPL_API bilinear(double *bs, double *as, int ord, double *bz, double *az)
Transform a s-plane analog filter transfer function  to the digital filter transfer function .
Definition: bilinear.c:210

iir
int DSPL_API iir(double rp, double rs, int ord, double w0, double w1, int type, double *b, double *a)
Digital IIR filter design.
Definition: iir.c:244

cheby2_ap_zp
int DSPL_API cheby2_ap_zp(int ord, double rs, complex_t *z, int *nz, complex_t *p, int *np)
Function calculates arrays of zeros and poles for analog normlized lowpass Chebyshev type 2 filter tr...
Definition: cheby2_ap_zp.c:220

low2high
int DSPL_API low2high(double *b, double *a, int ord, double w0, double w1, double *beta, double *alpha)
Lowpass to highpass filter frequency transform.
Definition: low2high.c:148

cheby1_ap
int DSPL_API cheby1_ap(double rp, int ord, double *b, double *a)
Function calculates the transfer function  coefficients of analog normalized lowpass Chebyshev type 1...
Definition: cheby1_ap.c:181

butter_ap
int DSPL_API butter_ap(double rp, int ord, double *b, double *a)
Function calculates the transfer function  coefficients of analog normalized lowpass Butterworth filt...
Definition: butter_ap.c:173

ellip_ap_zp
int DSPL_API ellip_ap_zp(int ord, double rp, double rs, complex_t *z, int *nz, complex_t *p, int *np)
Function calculates arrays of zeros and poles for analog normlized lowpass elliptic filter transfer f...
Definition: ellip_ap_zp.c:220

cheby1_ap_zp
int DSPL_API cheby1_ap_zp(int ord, double rp, complex_t *z, int *nz, complex_t *p, int *np)
Function calculates arrays of zeros and poles for analog normlized lowpass Chebyshev type 1 filter tr...
Definition: cheby1_ap_zp.c:205

low2low
int DSPL_API low2low(double *b, double *a, int ord, double w0, double w1, double *beta, double *alpha)
Definition: low2low.c:149

ratcompos
int DSPL_API ratcompos(double *b, double *a, int n, double *c, double *d, int p, double *beta, double *alpha)
Rational composition.
Definition: ratcompos.c:184

butter_ap_zp
int DSPL_API butter_ap_zp(int ord, double rp, complex_t *z, int *nz, complex_t *p, int *np)
Function calculates arrays of zeros and poles for analog normlized lowpass Batterworth filter transfe...
Definition: butter_ap_zp.c:211

ellip_ap
int DSPL_API ellip_ap(double rp, double rs, int ord, double *b, double *a)
Function calculates the transfer function  coefficients of analog normalized lowpass elliptic filter ...
Definition: ellip_ap.c:193

cheby2_ap
int DSPL_API cheby2_ap(double rs, int ord, double *b, double *a)
Function calculates the transfer function  coefficients of analog normalized lowpass Chebyshev type 2...
Definition: cheby2_ap.c:184

writebin
int DSPL_API writebin(void *x, int n, int m, int dtype, char *fn)
Save 1D vector to the binary file.
Definition: writebin.c:217

writetxt_3d
int DSPL_API writetxt_3d(double *x, int nx, double *y, int ny, double *z, char *fn)
Save data  to the text file fn for 3D surface plotting.
Definition: writetxt_3d.c:214

readbin
int DSPL_API readbin(char *fn, void **x, int *pn, int *pm, int *dtype)
Read array from the bin-file.
Definition: readbin.c:314

writetxt_int
int DSPL_API writetxt_int(int *x, int *y, int n, char *fn)
Save interger data  to the text file fn.  .
Definition: writetxt_int.c:115

writetxt
int DSPL_API writetxt(double *x, double *y, int n, char *fn)
Save real data  to the text file fn.  .
Definition: writetxt.c:122

gnuplot_open
int DSPL_API gnuplot_open(void **hplot)
Open GNUPLOT program.
Definition: gnuplot_open.c:103

gnuplot_close
void DSPL_API gnuplot_close(void *h)
Close GNUPLOT handle.
Definition: gnuplot_close.c:80

gnuplot_create
int DSPL_API gnuplot_create(int argc, char *argv[], int w, int h, char *fn_png, void **hplot)
Create GNUPLOT chart.
Definition: gnuplot_create.c:202

gnuplot_cmd
void DSPL_API gnuplot_cmd(void *h, char *cmd)
Function sends cmd command to GNUPLOT corresponds to h handle.
Definition: gnuplot_cmd.c:82

sinc
int DSPL_API sinc(double *x, int n, double a, double *y)
Function  for the real vector x.
Definition: sinc.c:88

sqrt_cmplx
int DSPL_API sqrt_cmplx(complex_t *x, int n, complex_t *y)
Square root of the complex vector argguument x.
Definition: sqrt_cmplx.c:136

log_cmplx
int DSPL_API log_cmplx(complex_t *x, int n, complex_t *y)
The logarithm function the complex vector argument x.
Definition: log_cmplx.c:139

ellip_sn_cmplx
int DSPL_API ellip_sn_cmplx(complex_t *u, int n, double k, complex_t *y)
Jacobi elliptic function  of complex vector argument.
Definition: ellip_sn_cmplx.c:105

ellip_cd_cmplx
int DSPL_API ellip_cd_cmplx(complex_t *u, int n, double k, complex_t *y)
Jacobi elliptic function  of complex vector argument.
Definition: ellip_cd_cmplx.c:105

ellip_asn
int DSPL_API ellip_asn(double *w, int n, double k, double *u)
Inverse Jacobi elliptic function  of real vector argument.
Definition: ellip_asn.c:106

ellip_acd_cmplx
int DSPL_API ellip_acd_cmplx(complex_t *w, int n, double k, complex_t *u)
Inverse Jacobi elliptic function  of complex vector argument.
Definition: ellip_acd_cmplx.c:104

ellip_cd
int DSPL_API ellip_cd(double *u, int n, double k, double *y)
Jacobi elliptic function  of real vector argument.
Definition: ellip_cd.c:125

ellip_asn_cmplx
int DSPL_API ellip_asn_cmplx(complex_t *w, int n, double k, complex_t *u)
Inverse Jacobi elliptic function  of complex vector argument.
Definition: ellip_asn_cmplx.c:105

ellip_sn
int DSPL_API ellip_sn(double *u, int n, double k, double *y)
Jacobi elliptic function  of real vector argument.
Definition: ellip_sn.c:127

ellip_landen
int DSPL_API ellip_landen(double k, int n, double *y)
Function calculates complete elliptical integral coefficients .
Definition: ellip_landen.c:175

ellip_acd
int DSPL_API ellip_acd(double *w, int n, double k, double *u)
Inverse Jacobi elliptic function  of the real vector argument.
Definition: ellip_acd.c:102

matrix_eye_cmplx
int DSPL_API matrix_eye_cmplx(complex_t *a, int n, int m)
The complex identity matrix size n x m generation.
Definition: matrix_eye_cmplx.c:81

matrix_eye
int DSPL_API matrix_eye(double *a, int n, int m)
The real identity matrix size n x m generation.
Definition: matrix_eye.c:80

matrix_mul
int DSPL_API matrix_mul(double *a, int na, int ma, double *b, int nb, int mb, double *c)
Matrix multiplication.
Definition: matrix_mul.c:191

matrix_eig_cmplx
int DSPL_API matrix_eig_cmplx(complex_t *a, int n, complex_t *v, int *info)
Eigenvalues calculation of the complex square matrix a.
Definition: matrix_eig_cmplx.c:143

cheby_poly1
int DSPL_API cheby_poly1(double *x, int n, int ord, double *y)
Chebyshev polynomial of the first kind order ord
Definition: cheby_poly1.c:136

cheby_poly2
int DSPL_API cheby_poly2(double *x, int n, int ord, double *y)
Chebyshev polynomial of the second kind order ord
Definition: cheby_poly2.c:133

polyroots
int DSPL_API polyroots(double *a, int ord, complex_t *r, int *info)
Function calculates real polynomial roots.
Definition: polyroots.c:157

randb2
int DSPL_API randb2(double *x, int n, random_t *prnd)
Binary bipolar [-1, 1] pseudorandom vector.
Definition: randb2.c:116

random_init
int DSPL_API random_init(random_t *prnd, int type, void *seed)
Pseudorandom numbers generators initialization.
Definition: random_init.c:120

randb
int DSPL_API randb(double *x, int n, random_t *prnd)
Binary unipolar [0, 1] pseudorandom vector.
Definition: randb.c:118

xcorr_cmplx
int DSPL_API xcorr_cmplx(complex_t *x, int nx, complex_t *y, int ny, int flag, int nr, complex_t *r, double *t)
Estimates the cross-correlation vector for complex discrete-time sequences x and y.
Definition: xcorr.c:449

mean_cmplx
int DSPL_API mean_cmplx(complex_t *x, int n, complex_t *m)
Calculates the mean of the complex input vector x
Definition: mean_cmplx.c:108

mean
int DSPL_API mean(double *x, int n, double *m)
Calculates the mean of the input vector x
Definition: mean.c:105

find_max_abs
int DSPL_API find_max_abs(double *a, int n, double *m, int *ind)
Find maximum absolute value from the real vector a
Definition: find_max_abs.c:121

stat_std_cmplx
int DSPL_API stat_std_cmplx(complex_t *x, int n, double *s)
Calculates the standard deviation of the complex input vector x
Definition: stat_std_cmplx.c:113

stat_std
int DSPL_API stat_std(double *x, int n, double *s)
Calculates the standard deviation of the input vector x
Definition: stat_std.c:116

xcorr
int DSPL_API xcorr(double *x, int nx, double *y, int ny, int flag, int nr, double *r, double *t)
Estimates the cross-correlation vector for real discrete-time sequences x and y.
Definition: xcorr.c:214

bessel_i0
int DSPL_API bessel_i0(double *x, int n, double *y)
Modified Bessel Function of the First Kind –  [1].
Definition: bessel_i0.c:116

sine_int
int DSPL_API sine_int(double *x, int n, double *si)
Sine integral function  for the real vector x.
Definition: sine_int.c:114

cos_cmplx
int DSPL_API cos_cmplx(complex_t *x, int n, complex_t *y)
The cosine function the complex vector argument x.
Definition: cos_cmplx.c:136

acos_cmplx
int DSPL_API acos_cmplx(complex_t *x, int n, complex_t *y)
The inverse of the cosine function the complex vector argument x.
Definition: acos_cmplx.c:141

sin_cmplx
int DSPL_API sin_cmplx(complex_t *x, int n, complex_t *y)
The sine function the complex vector argument x.
Definition: sin_cmplx.c:136

asin_cmplx
int DSPL_API asin_cmplx(complex_t *x, int n, complex_t *y)
The inverse of the sine function the complex vector argument x.
Definition: asin_cmplx.c:138

dspl_load
void * dspl_load()
Perform dynamic linking and load libdspl-2.0 functions.

dspl_free
void dspl_free(void *handle)
Cleans up the previously linked DSPL-2.0 dynamic library.

cmplx2re
int DSPL_API cmplx2re(complex_t *x, int n, double *re, double *im)
Separate complex vector to the real and image vectors.
Definition: cmplx2re.c:130

re2cmplx
int DSPL_API re2cmplx(double *x, int n, complex_t *y)
Convert real array to the complex array.
Definition: re2cmplx.c:120

window
int DSPL_API window(double *w, int n, int win_type, double param)
Window function calculation.
Definition: win.c:329