- メインページ
- モジュール
- ネームスペース
- データ構造
- ファイル

- ファイル一覧
- グローバル

fir.h
説明を見る。
1 //
2 // Copyright (c) 2003-2011, MIST Project, Nagoya University
3 // All rights reserved.
4 //
5 // Redistribution and use in source and binary forms, with or without modification,
6 // are permitted provided that the following conditions are met:
7 //
8 // 1. Redistributions of source code must retain the above copyright notice,
9 // this list of conditions and the following disclaimer.
10 //
11 // 2. Redistributions in binary form must reproduce the above copyright notice,
12 // this list of conditions and the following disclaimer in the documentation
13 // and/or other materials provided with the distribution.
14 //
15 // 3. Neither the name of the Nagoya University nor the names of its contributors
16 // may be used to endorse or promote products derived from this software
17 // without specific prior written permission.
18 //
19 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
20 // IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND
21 // FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
22 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
23 // DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
24 // DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER
25 // IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF
26 // THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
27 //
28
33
34 #ifndef __INCLUDE_FILTER_FIR_FILTER_H__
35 #define __INCLUDE_FILTER_FIR_FILTER_H__
36
37
38 #ifndef __INCLUDE_MIST_H__
39 #include "../mist.h"
40 #endif
41
42 #ifndef __INCLUDE_MIST_TYPE_TRAIT_H__
43 #include "../config/type_trait.h"
44 #endif
45
46 #include <cmath>
47 #include <functional>
48
49 // mist名前空間の始まり
50 _MIST_BEGIN
51
52
53 namespace __fir_filter__
54 {
55
56 template<class T, class Allocator>
57 void fir(
58 const array<T, Allocator>& in,
59 array<T, Allocator>& out,
60 int order,
61 const array<T, Allocator>& coef,
62 {
63 array<double> buf(order, 0.0);
64
65 out.resize(in.size());
66
67 double output;
68 for(int n = 0; n < in.size(); n ++)
69 {
70 for(int k = 0; k < order - 1; k ++)
71 {
72 buf[k] = coef[k] * in[n] + buf[k + 1];
73 }
74 buf[order - 1] = coef[order - 1] * in[n];
75
76 out[n] = buf[0];
77 }
78 }
79
80 }
81
82
83 // mist名前空間の終わり
84 _MIST_END
85
86 #endif // __INCLUDE_FILTER_FIR_FILTER_H__
87