1 | /* |
---|

2 | Copyright (C) 2003-2009 Paul Brossier <piem@aubio.org> |
---|

3 | |
---|

4 | This file is part of aubio. |
---|

5 | |
---|

6 | aubio is free software: you can redistribute it and/or modify |
---|

7 | it under the terms of the GNU General Public License as published by |
---|

8 | the Free Software Foundation, either version 3 of the License, or |
---|

9 | (at your option) any later version. |
---|

10 | |
---|

11 | aubio is distributed in the hope that it will be useful, |
---|

12 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
---|

13 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|

14 | GNU General Public License for more details. |
---|

15 | |
---|

16 | You should have received a copy of the GNU General Public License |
---|

17 | along with aubio. If not, see <http://www.gnu.org/licenses/>. |
---|

18 | |
---|

19 | */ |
---|

20 | |
---|

21 | #include "aubio_priv.h" |
---|

22 | #include "fvec.h" |
---|

23 | #include "lvec.h" |
---|

24 | #include "temporal/filter.h" |
---|

25 | |
---|

26 | uint_t |
---|

27 | aubio_filter_set_biquad (aubio_filter_t * f, lsmp_t b0, lsmp_t b1, lsmp_t b2, |
---|

28 | lsmp_t a1, lsmp_t a2) |
---|

29 | { |
---|

30 | uint_t order = aubio_filter_get_order (f); |
---|

31 | lvec_t *bs = aubio_filter_get_feedforward (f); |
---|

32 | lvec_t *as = aubio_filter_get_feedback (f); |
---|

33 | |
---|

34 | if (order != 3) { |
---|

35 | AUBIO_ERROR ("order of biquad filter must be 3, not %d\n", order); |
---|

36 | return AUBIO_FAIL; |
---|

37 | } |
---|

38 | bs->data[0] = b0; |
---|

39 | bs->data[1] = b1; |
---|

40 | bs->data[2] = b2; |
---|

41 | as->data[0] = 1.; |
---|

42 | as->data[1] = a1; |
---|

43 | as->data[1] = a2; |
---|

44 | return AUBIO_OK; |
---|

45 | } |
---|

46 | |
---|

47 | aubio_filter_t * |
---|

48 | new_aubio_filter_biquad (lsmp_t b0, lsmp_t b1, lsmp_t b2, lsmp_t a1, lsmp_t a2) |
---|

49 | { |
---|

50 | aubio_filter_t *f = new_aubio_filter (3); |
---|

51 | aubio_filter_set_biquad (f, b0, b1, b2, a1, a2); |
---|

52 | return f; |
---|

53 | } |
---|