libdspl-2.0
Библиотека алгоритмов цифровой обработки сигналов
Функции анализа аналоговых и цифровых фильтров

Функции анализа аналоговых и цифровых фильтров Подробнее...

Функции

int filter_freq_resp (double *b, double *a, int ord, double *w, int n, int flag, double *mag, double *phi, double *tau)
 Расчет амплитудно-частотной (АЧХ), фазочастотной характеристик (ФЧХ), а также группового времени запаздывания (ГВЗ) цифрового или аналогового или фильтра. Подробнее...
 
int freqz (double *b, double *a, int ord, double *w, int n, complex_t *h)
 Расчет комплексного коэффициента передачи \( H \left(e^{j \omega} \right)\) цифрового фильтра. Подробнее...
 
int group_delay (double *pb, double *pa, int ord, int flag, double *w, int n, double *tau)
 Расчет группового времени запаздывания цифрового или аналогового фильтра. Подробнее...
 
int phase_delay (double *b, double *a, int ord, int flag, double *w, int n, double *tau)
 Расчет фазовой задержки цифрового или аналогового фильтра. Подробнее...
 
int freqs (double *b, double *a, int ord, double *w, int n, complex_t *h)
 Расчет комплексного коэффициента передачи \( H(j \omega) \) аналогового фильтра. Подробнее...
 

Подробное описание

Функции анализа аналоговых и цифровых фильтров

Функции

◆ filter_freq_resp()

int filter_freq_resp ( double *  b,
double *  a,
int  ord,
double *  w,
int  n,
int  flag,
double *  mag,
double *  phi,
double *  tau 
)

Расчет амплитудно-частотной (АЧХ), фазочастотной характеристик (ФЧХ), а также группового времени запаздывания (ГВЗ) цифрового или аналогового или фильтра.


Функция рассчитывает АЧХ, ФЧХ и ГВЗ аналогового или цифрового фильтра, заданного передаточной характеристикой \(H(s)\), или \(H(z)\) соответственно

Аргументы
[in]bУказатель на вектор коэффициентов числителя передаточной функции \( H(s) \).
Размер вектора [ord+1 x 1].

[in]aУказатель на вектор коэффициентов знаменателя передаточной функции \( H(s) \).
Размер вектора [ord+1 x 1].

[in]ordПорядок фильтра. Количество коэффициентов числителя и знаменателя передаточной функции \( H(s) \) равно ord+1.

[in]wУказатель на вектор значений циклической частоты \( \omega \) (рад/с), для которого будет рассчитаны АЧХ, ФЧХ и ГВЗ аналогового фильтра, если установлен флаг DSPL_FLAG_ANALOG.
В случае если флаг DSPL_FLAG_ANALOG не установлен, то вектор частоты w используется как нормированная частота комплексного коэффициента передачи \( H \left(\mathrm{e}^{j\omega} \right) \) цифрового фильтра.
В этом случае характеристика цифрового фильтра является \( 2\pi \)-периодической, и вектор частоты может содержать произвольные значения, однако целесообразно задавать его от 0 до \( \pi \), а такжет от 0 до \( 2\pi \), или от \( -\pi \) до \( \pi \).
Размер вектора [n x 1].

[in]nРазмер вектора циклической частоты w.

[in]flagКомбинация флагов, которые задают расчет параметров:
DSPL_FLAG_ANALOG  Коэффициенты относятся к аналоговому фильтру
DSPL_FLAG_LOGMAG  АЧХ рассчитывать в логарифмическом масштабе
DSPL_FLAG_UNWRAP  раскрывать периодичность ФЧХ
[out]magУказатель на вектор АЧХ.
Размер вектора [n x 1].
Память должна быть выделена.
Если указатель NULL, то расчет АЧХ не производится.

[out]phiУказатель на вектор ФЧХ.
Размер вектора [n x 1].
Память должна быть выделена.
Если указатель NULL, то расчет ФЧХ не производится.

[out]tauУказатель на вектор ГВЗ.
Размер вектора [n x 1].
Память должна быть выделена.
Если указатель NULL, то расчет ГВЗ не производится.

Возвращает
RES_OK Параметры фильтра рассчитаны успешно.
В противном случае код ошибки.
Пример использования функции filter_freq_resp:
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include "dspl.h"
/* Filter order */
#define ORD 3
/* Frequency response vector size */
#define N 1000
int main(int argc, char* argv[])
{
void* hdspl; /* DSPL handle */
void* hplot; /* GNUPLOT handle */
/* Load DSPL functions */
hdspl = dspl_load();
/* H(s) numerator coefficients vector */
double a[ORD+1];
/* H(s) denominator coefficients vector */
double b[ORD+1];
/* Magnitude ripple from 0 to 1 rad/s */
double Rp = 1.0;
/* Angular frequency (rad/s) */
double w[N];
/* Filter Magnitude (dB) */
double mag[N];
/* Phase response */
double phi[N];
/* Group delay */
double tau[N];
int k;
/* H(s) coefficients calculation */
int res = butter_ap(Rp, ORD, b, a);
if(res != RES_OK)
printf("error code = 0x%8x\n", res);
/* Print H(s) coefficients */
for(k = 0; k < ORD+1; k++)
printf("b[%2d] = %9.3f a[%2d] = %9.3f\n", k, b[k], k, a[k]);
/* Frequency in logarithmic scale from 0.01 to 100 rad/s */
logspace(-2.0, 2.0, N , DSPL_SYMMETRIC, w);
/* Filter frequency parameter calculation */
filter_freq_resp(b, a, ORD, w, N,
DSPL_FLAG_LOGMAG|DSPL_FLAG_UNWRAP|DSPL_FLAG_ANALOG,
mag, phi, tau);
/* Write Magnitude, phase response and group delay to the files */
writetxt(w, mag, N, "dat/butter_ap_test_mag.txt");
writetxt(w, phi, N, "dat/butter_ap_test_phi.txt");
writetxt(w, tau, N, "dat/butter_ap_test_tau.txt");
/* plotting by GNUPLOT */
gnuplot_create(argc, argv, 920, 260, "img/butter_ap_test.png", &hplot);
gnuplot_cmd(hplot, "set logscale x");
gnuplot_cmd(hplot, "unset key");
gnuplot_cmd(hplot, "set grid");
gnuplot_cmd(hplot, "set xlabel 'frequency, rad/s'");
gnuplot_cmd(hplot, "set multiplot layout 1,3 rowsfirst");
gnuplot_cmd(hplot, "set ylabel 'Magnitude, dB'");
gnuplot_cmd(hplot, "set yrange [-100:5]");
gnuplot_cmd(hplot, "plot 'dat/butter_ap_test_mag.txt' with lines");
gnuplot_cmd(hplot, "set ylabel 'Phase response, rad'");
gnuplot_cmd(hplot, "unset yrange");
gnuplot_cmd(hplot, "plot 'dat/butter_ap_test_phi.txt' with lines");
gnuplot_cmd(hplot, "set ylabel 'Groupdelay, sec'");
gnuplot_cmd(hplot, "unset yrange");
gnuplot_cmd(hplot, "plot 'dat/butter_ap_test_tau.txt' with lines");
gnuplot_cmd(hplot, "unset multiplot");
gnuplot_close(hplot);
/* free dspl handle */
dspl_free(hdspl);
return res;
}

Результат работы программы:

b[ 0] =   1.002   a[ 0] =   1.002
b[ 1] =   0.000   a[ 1] =   2.618
b[ 2] =   0.000   a[ 2] =   3.418
b[ 3] =   0.000   a[ 3] =   2.615
b[ 4] =   0.000   a[ 4] =   1.000


В каталоге dat будут созданы три файла:

butter_ap_test_mag.txt    АЧХ фильтра
butter_ap_test_phi.txt    ФЧХ фильтра
butter_ap_test_tau.txt    ГВЗ фильтра

Кроме того программа GNUPLOT произведет построение следующих графиков по сохраненным в файлах данным:

Автор
Бахурин Сергей www.dsplib.org

См. определение в файле filter_freq_resp.c строка 231

Используется в phase_delay().

◆ freqs()

int freqs ( double *  b,
double *  a,
int  ord,
double *  w,
int  n,
complex_t h 
)

Расчет комплексного коэффициента передачи \( H(j \omega) \) аналогового фильтра.


Функция рассчитывает значения комплексного коэффициента передачи \( H(j \omega)\) аналогового фильтра, заданного коэффициентами передаточной функции \( H(s) \):

\[ H(s) = \frac {\sum_{k = 0}^{N} b_k s^k} {\sum_{m = 0}^{N} a_m s^m}, \]

где \( N \) - порядок фильтра (параметр ord).

Комплексный коэффициент передачи рассчитывается путем подстановки \( s = j \omega \).

Аргументы
[in]bУказатель на вектор коэффициентов числителя передаточной функции \( H(s) \).
Размер вектора [ord+1 x 1].

[in]aУказатель на вектор коэффициентов знаменателя передаточной функции \( H(s) \).
Размер вектора [ord+1 x 1].

[in]ordПорядок фильтра. Количество коэффициентов числителя и знаменателя передаточной функции \( H(s) \) равно ord+1.

[in]wУказатель на вектор значений циклической частоты \( \omega \) (рад/с), для которого будет рассчитан комплексный коэффициент передачи \( H(j \omega) \).
Размер вектора [n x 1].

[in]nРазмер вектора циклической частоты w.

[out]hУказатель на вектор комплексного коэффициента передачи \( H(j \omega) \), рассчитанного для циклической частоты w.
Размер вектора [n x 1].
Память должна быть выделена.

Возвращает
RES_OK Комплексный коэффициент передачи рассчитан успешно.
В противном случае код ошибки.
Автор
Бахурин Сергей www.dsplib.org

См. определение в файле freqs.c строка 143

Используется в filter_freq_resp().

◆ freqz()

int freqz ( double *  b,
double *  a,
int  ord,
double *  w,
int  n,
complex_t h 
)

Расчет комплексного коэффициента передачи \( H \left(e^{j \omega} \right)\) цифрового фильтра.


Функция рассчитывает значения комплексного коэффициента передачи \( H \left(e^{j \omega} \right)\) цифрового фильтра, заданного коэффициентами передаточной функции \(H(z)\):

\[ H(z) = \frac {\sum_{k = 0}^{N} b_k z^{-k}} {\sum_{m = 0}^{N} a_m z^{-m}}, \]

где \(N\) — порядок фильтра (параметр ord).
Комплексный коэффициент передачи рассчитывается путем подстановки \(z = e^{j \omega} \).

Аргументы
[in]bУказатель на вектор коэффициентов числителя передаточной функции \(H(z)\).
Размер вектора [ord+1 x 1].

[in]aУказатель на вектор коэффициентов знаменателя передаточной функции \(H(z)\).
Размер вектора [ord+1 x 1].

[in]ordПорядок фильтра. Количество коэффициентов числителя и знаменателя передаточной функции \(H(z)\) равно ord+1.

[in]wУказатель на вектор значений нормированной циклической частоты \(\omega\), для которого будет рассчитан комплексный коэффициент передачи \( H \left(e^{j \omega} \right)\).
Размер вектора [n x 1].

[in]nРазмер вектора нормированной циклической частоты w.

[out]hУказатель на вектор комплексного коэффициента передачи \( H \left(e^{j \omega} \right)\), рассчитанного для циклической частоты w.
Размер вектора [n x 1].
Память должна быть выделена.

Возвращает
RES_OK Комплексный коэффициент передачи рассчитан успешно.
В противном случае код ошибки.
Заметки
Комплексный коэффициент передачи \( H \left(e^{j \omega} \right)\) цифрового фильтра представляет собой \( 2 \pi-\)периодическую функцию нормированной циклической частоты \(\omega\). Поэтому анализ цифровых фильтров целесообразно вести на одном периоде повторения \( H \left(e^{j \omega} \right)\), т.е. в интервале \(\omega\) от 0 до \(2 \pi\), или от \(-\pi\) до \( \pi\).
Кроме того известно, что для фильтра с вещественными коэффициентами \( H \left(e^{j \omega} \right) = H^* \left(e^{-j \omega} \right)\), а значит, анализ цифрового фильтра с вещественными коэффициентами достаточно вести для нормированной частоты \(\omega\) от 0 до \(\pi\).
Автор
Бахурин Сергей www.dsplib.org

См. определение в файле freqz.c строка 155

Используется в filter_freq_resp().

◆ group_delay()

int group_delay ( double *  pb,
double *  pa,
int  ord,
int  flag,
double *  w,
int  n,
double *  tau 
)

Расчет группового времени запаздывания цифрового или аналогового фильтра.


Групповое время запаздывания определяется как:

\[ \tau_g(\omega) = - \frac{d\Phi(\omega)}{d\omega}, \]

где \(\Phi(\omega)\) – ФЧХ фильтра, \(\omega\) циктическая частот в случае аналогового фильтра, или нормированная частота цифрового фильтра.

Аргументы
[in]pbУказатель на вектор коэффициентов числителя передаточной функции аналогового фильтра \( H(s) \) или цифрового фильтра \( H(z) \).
Размер вектора [ord+1 x 1].

[in]paУказатель на вектор коэффициентов числителя передаточной функции аналогового фильтра \( H(s) \) или цифрового фильтра \( H(z) \).
Размер вектора [ord+1 x 1].

Параметр может быть NULL. В этом случае расчет производится для цифрового КИХ-фильтра с коэффициентами, заданными вектором b.

[in]ordПорядок фильтра. Количество коэффициентов числителя и знаменателя передаточной функции \( H(s) \) или \( H(z) \) равно ord+1.

[in]flagФлаг который задает тип фильтра:
DSPL_FLAG_ANALOG  Коэффициенты относятся к аналоговому фильтру
DSPL_FLAG_DIGITAL  Коэффициенты относятся к цифровому фильтру
[in]wУказатель на вектор значений циклической частоты \( \omega \) (рад/с), для которого будет рассчитаны АЧХ, ФЧХ и ГВЗ аналогового фильтра, если установлен флаг DSPL_FLAG_ANALOG.
В случае если флаг DSPL_FLAG_ANALOG не установлен, то вектор частоты w используется как нормированная частота комплексного коэффициента передачи \( H \left(\mathrm{e}^{j\omega} \right) \) цифрового фильтра.
В этом случае характеристика цифрового фильтра является \( 2\pi \)-периодической, и вектор частоты может содержать произвольные значения, однако целесообразно задавать его от 0 до \( \pi \), а такжет от 0 до \( 2\pi \), или от \( -\pi \) до \( \pi \).
Размер вектора [n x 1].

[in]nРазмер вектора циклической частоты w.

[out]tauУказатель на вектор групповой задержки.
Размер вектора [n x 1].
Память должна быть выделена.
Возвращает
RES_OK групповая задержка фильтра рассчитана успешно.
В противном случае код ошибки.
Автор
Бахурин Сергей www.dsplib.org

См. определение в файле group_delay.c строка 161

Используется в filter_freq_resp().

◆ phase_delay()

int phase_delay ( double *  b,
double *  a,
int  ord,
int  flag,
double *  w,
int  n,
double *  tau 
)

Расчет фазовой задержки цифрового или аналогового фильтра.


Фазовая задержка определяется как:

\[ \tau_{\varphi}(\omega) = - \frac{\Phi(\omega)}{\omega}, \]

где \(\Phi(\omega)\) – ФЧХ фильтра, \(\omega\) циктическая частот в случае аналогового фильтра, или нормированная частота цифрового фильтра.

Аргументы
[in]bУказатель на вектор коэффициентов числителя передаточной функции аналогового фильтра \( H(s) \) или цифрового фильтра \( H(z) \).
Размер вектора [ord+1 x 1].

[in]aУказатель на вектор коэффициентов числителя передаточной функции аналогового фильтра \( H(s) \) или цифрового фильтра \( H(z) \).
Размер вектора [ord+1 x 1].

Параметр может быть NULL. В этом случае расчет производится для цифрового КИХ-фильтра с коэффициентами, заданными вектором b.

[in]ordПорядок фильтра. Количество коэффициентов числителя и знаменателя передаточной функции \( H(s) \) или \( H(z) \) равно ord+1.

[in]flagФлаг который задает тип фильтра:
DSPL_FLAG_ANALOG  Коэффициенты относятся к аналоговому фильтру
DSPL_FLAG_DIGITAL  Коэффициенты относятся к цифровому фильтру
[in]wУказатель на вектор значений циклической частоты \( \omega \) (рад/с), для которого будет рассчитаны АЧХ, ФЧХ и ГВЗ аналогового фильтра, если установлен флаг DSPL_FLAG_ANALOG.
В случае если флаг DSPL_FLAG_ANALOG не установлен, то вектор частоты w используется как нормированная частота комплексного коэффициента передачи \( H \left(\mathrm{e}^{j\omega} \right) \) цифрового фильтра.
В этом случае характеристика цифрового фильтра является \( 2\pi \)-периодической, и вектор частоты может содержать произвольные значения, однако целесообразно задавать его от 0 до \( \pi \), а такжет от 0 до \( 2\pi \), или от \( -\pi \) до \( \pi \).
Размер вектора [n x 1].

[in]nРазмер вектора циклической частоты w.

[out]tauУказатель на вектор фазовой задержки.
Размер вектора [n x 1].
Память должна быть выделена.
Возвращает
RES_OK фазовая задержка фильтра рассчитана успешно.
В противном случае код ошибки.
Автор
Бахурин Сергей www.dsplib.org

См. определение в файле phase_delay.c строка 162

void gnuplot_close(void *h)
Закрыть хэндл GNUPLOT.
Definition: gnuplot_close.c:80
int logspace(double x0, double x1, int n, int type, double *x)
Функция заполняет массив значениями логарифмической шкале
Definition: logspace.c:178
void * dspl_load()
Произвести динамическую линковку и загрузить функции libdspl-2.0.
void gnuplot_cmd(void *h, char *cmd)
Функция посылает команду cmd пакету GNUPLOT, для построения или оформления графика,...
Definition: gnuplot_cmd.c:82
int butter_ap(double rp, int ord, double *b, double *a)
Расчет передаточной характеристики аналогового нормированного ФНЧ Баттерворта.
Definition: butter_ap.c:172
int writetxt(double *x, double *y, int n, char *fn)
Сохранить вещественные данные в текстовый файл
Definition: writetxt.c:122
int filter_freq_resp(double *b, double *a, int ord, double *w, int n, int flag, double *mag, double *phi, double *tau)
Расчет амплитудно-частотной (АЧХ), фазочастотной характеристик (ФЧХ), а также группового времени запа...
void dspl_free(void *handle)
Очищает связанную ранее динамическую библиотеку DSPL-2.0.
#define RES_OK
Функция завершилась корректно. Ошибки отсутствуют.
Definition: dspl.h:558
int gnuplot_create(int argc, char *argv[], int w, int h, char *fn_png, void **hplot)
Создать график GNUPLOT.