-
Notifications
You must be signed in to change notification settings - Fork 5
/
cw-shifter.cpp
52 lines (47 loc) · 1.6 KB
/
cw-shifter.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
#
/*
* Copyright (C) 2020
* Jan van Katwijk (J.vanKatwijk@gmail.com)
* Lazy Chair Computing
*
* This file is part of the SDRuno cw decoder
*
* cw decoder is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* cw decoder 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 General Public License
* along with cw decoder; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*/
#include "cw-shifter.h"
#define _USE_MATH_DEFINES
#include <math.h>
cwShifter::cwShifter (int32_t size):
table (size) {
int i;
this ->tableSize = size;
phase = 0;
for (i = 0; i < size; i ++)
table [i] = std::complex<float>
((float)cos ((float) i / (float)size * 2 * M_PI),
(float)sin ((float) i / (float)size * 2 * M_PI));
}
cwShifter::~cwShifter (void) {
}
std::complex<float>
cwShifter::doShift (std::complex<float> z, int32_t freq) {
std::complex<float> res = z * table [phase];
phase -= freq;
if (phase < 0)
phase += tableSize;
if (phase >= tableSize)
phase -= tableSize;
return res;
}