Chris@69
|
1 /***********************************************************************
|
Chris@69
|
2 Copyright (c) 2006-2011, Skype Limited. All rights reserved.
|
Chris@69
|
3 Redistribution and use in source and binary forms, with or without
|
Chris@69
|
4 modification, are permitted provided that the following conditions
|
Chris@69
|
5 are met:
|
Chris@69
|
6 - Redistributions of source code must retain the above copyright notice,
|
Chris@69
|
7 this list of conditions and the following disclaimer.
|
Chris@69
|
8 - Redistributions in binary form must reproduce the above copyright
|
Chris@69
|
9 notice, this list of conditions and the following disclaimer in the
|
Chris@69
|
10 documentation and/or other materials provided with the distribution.
|
Chris@69
|
11 - Neither the name of Internet Society, IETF or IETF Trust, nor the
|
Chris@69
|
12 names of specific contributors, may be used to endorse or promote
|
Chris@69
|
13 products derived from this software without specific prior written
|
Chris@69
|
14 permission.
|
Chris@69
|
15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
|
Chris@69
|
16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
Chris@69
|
17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
Chris@69
|
18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
|
Chris@69
|
19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
|
Chris@69
|
20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
|
Chris@69
|
21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
|
Chris@69
|
22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
|
Chris@69
|
23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
|
Chris@69
|
24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
|
Chris@69
|
25 POSSIBILITY OF SUCH DAMAGE.
|
Chris@69
|
26 ***********************************************************************/
|
Chris@69
|
27
|
Chris@69
|
28 #ifdef HAVE_CONFIG_H
|
Chris@69
|
29 #include "config.h"
|
Chris@69
|
30 #endif
|
Chris@69
|
31
|
Chris@69
|
32 /* Filter coefficients for IIR/FIR polyphase resampling *
|
Chris@69
|
33 * Total size: 179 Words (358 Bytes) */
|
Chris@69
|
34
|
Chris@69
|
35 #include "resampler_private.h"
|
Chris@69
|
36
|
Chris@69
|
37 /* Matlab code for the notch filter coefficients: */
|
Chris@69
|
38 /* B = [1, 0.147, 1]; A = [1, 0.107, 0.89]; G = 0.93; freqz(G * B, A, 2^14, 16e3); axis([0, 8000, -10, 1]) */
|
Chris@69
|
39 /* fprintf('\t%6d, %6d, %6d, %6d\n', round(B(2)*2^16), round(-A(2)*2^16), round((1-A(3))*2^16), round(G*2^15)) */
|
Chris@69
|
40 /* const opus_int16 silk_resampler_up2_hq_notch[ 4 ] = { 9634, -7012, 7209, 30474 }; */
|
Chris@69
|
41
|
Chris@69
|
42 /* Tables with IIR and FIR coefficients for fractional downsamplers (123 Words) */
|
Chris@69
|
43 silk_DWORD_ALIGN const opus_int16 silk_Resampler_3_4_COEFS[ 2 + 3 * RESAMPLER_DOWN_ORDER_FIR0 / 2 ] = {
|
Chris@69
|
44 -20694, -13867,
|
Chris@69
|
45 -49, 64, 17, -157, 353, -496, 163, 11047, 22205,
|
Chris@69
|
46 -39, 6, 91, -170, 186, 23, -896, 6336, 19928,
|
Chris@69
|
47 -19, -36, 102, -89, -24, 328, -951, 2568, 15909,
|
Chris@69
|
48 };
|
Chris@69
|
49
|
Chris@69
|
50 silk_DWORD_ALIGN const opus_int16 silk_Resampler_2_3_COEFS[ 2 + 2 * RESAMPLER_DOWN_ORDER_FIR0 / 2 ] = {
|
Chris@69
|
51 -14457, -14019,
|
Chris@69
|
52 64, 128, -122, 36, 310, -768, 584, 9267, 17733,
|
Chris@69
|
53 12, 128, 18, -142, 288, -117, -865, 4123, 14459,
|
Chris@69
|
54 };
|
Chris@69
|
55
|
Chris@69
|
56 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_2_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR1 / 2 ] = {
|
Chris@69
|
57 616, -14323,
|
Chris@69
|
58 -10, 39, 58, -46, -84, 120, 184, -315, -541, 1284, 5380, 9024,
|
Chris@69
|
59 };
|
Chris@69
|
60
|
Chris@69
|
61 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_3_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR2 / 2 ] = {
|
Chris@69
|
62 16102, -15162,
|
Chris@69
|
63 -13, 0, 20, 26, 5, -31, -43, -4, 65, 90, 7, -157, -248, -44, 593, 1583, 2612, 3271,
|
Chris@69
|
64 };
|
Chris@69
|
65
|
Chris@69
|
66 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_4_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR2 / 2 ] = {
|
Chris@69
|
67 22500, -15099,
|
Chris@69
|
68 3, -14, -20, -15, 2, 25, 37, 25, -16, -71, -107, -79, 50, 292, 623, 982, 1288, 1464,
|
Chris@69
|
69 };
|
Chris@69
|
70
|
Chris@69
|
71 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_6_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR2 / 2 ] = {
|
Chris@69
|
72 27540, -15257,
|
Chris@69
|
73 17, 12, 8, 1, -10, -22, -30, -32, -22, 3, 44, 100, 168, 243, 317, 381, 429, 455,
|
Chris@69
|
74 };
|
Chris@69
|
75
|
Chris@69
|
76 silk_DWORD_ALIGN const opus_int16 silk_Resampler_2_3_COEFS_LQ[ 2 + 2 * 2 ] = {
|
Chris@69
|
77 -2797, -6507,
|
Chris@69
|
78 4697, 10739,
|
Chris@69
|
79 1567, 8276,
|
Chris@69
|
80 };
|
Chris@69
|
81
|
Chris@69
|
82 /* Table with interplation fractions of 1/24, 3/24, 5/24, ... , 23/24 : 23/24 (46 Words) */
|
Chris@69
|
83 silk_DWORD_ALIGN const opus_int16 silk_resampler_frac_FIR_12[ 12 ][ RESAMPLER_ORDER_FIR_12 / 2 ] = {
|
Chris@69
|
84 { 189, -600, 617, 30567 },
|
Chris@69
|
85 { 117, -159, -1070, 29704 },
|
Chris@69
|
86 { 52, 221, -2392, 28276 },
|
Chris@69
|
87 { -4, 529, -3350, 26341 },
|
Chris@69
|
88 { -48, 758, -3956, 23973 },
|
Chris@69
|
89 { -80, 905, -4235, 21254 },
|
Chris@69
|
90 { -99, 972, -4222, 18278 },
|
Chris@69
|
91 { -107, 967, -3957, 15143 },
|
Chris@69
|
92 { -103, 896, -3487, 11950 },
|
Chris@69
|
93 { -91, 773, -2865, 8798 },
|
Chris@69
|
94 { -71, 611, -2143, 5784 },
|
Chris@69
|
95 { -46, 425, -1375, 2996 },
|
Chris@69
|
96 };
|