Mirror/aicodix___modem
1
0
Fork 0
mirror of https://github.com/aicodix/modem.git synced 2026-04-27 14:30:34 +00:00
Code Issues Projects Releases Packages Wiki Activity

added support for QAM4096

Browse source
This commit is contained in:
Ahmet Inan 2025-07-18 12:52:01 +02:00
commit c0fedcbc57
4 changed files with 83 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

2
common.hh
View file

@ -12,7 +12,7 @@ Copyright 2025 Ahmet Inan <inan@aicodix.de>
struct Common
{
static const int mod_max = 8;
static const int mod_max = 12;
static const int code_max = 16;
static const int bits_max = 1 << code_max;
static const int data_max = 4096;

6
decode.cc
View file

@ -105,6 +105,8 @@ struct Decoder : Common
return QuadratureAmplitudeModulation<256, cmplx, code_type>::map(b);
case 10:
return QuadratureAmplitudeModulation<1024, cmplx, code_type>::map(b);
case 12:
return QuadratureAmplitudeModulation<4096, cmplx, code_type>::map(b);
}
return 0;
}
@ -125,6 +127,8 @@ struct Decoder : Common
return QuadratureAmplitudeModulation<256, cmplx, code_type>::soft(b, c, precision);
case 10:
return QuadratureAmplitudeModulation<1024, cmplx, code_type>::soft(b, c, precision);
case 12:
return QuadratureAmplitudeModulation<4096, cmplx, code_type>::soft(b, c, precision);
}
}
void demap_hard(code_type *b, cmplx c, int bits)
@ -144,6 +148,8 @@ struct Decoder : Common
return QuadratureAmplitudeModulation<256, cmplx, code_type>::hard(b, c);
case 10:
return QuadratureAmplitudeModulation<1024, cmplx, code_type>::hard(b, c);
case 12:
return QuadratureAmplitudeModulation<4096, cmplx, code_type>::hard(b, c);
}
}
void shuffle(code_type *dest, const code_type *src, int order)

4
encode.cc
View file

@ -200,6 +200,8 @@ struct Encoder : public Common
return QuadratureAmplitudeModulation<256, cmplx, code_type>::map(b);
case 10:
return QuadratureAmplitudeModulation<1024, cmplx, code_type>::map(b);
case 12:
return QuadratureAmplitudeModulation<4096, cmplx, code_type>::map(b);
}
return 0;
}
@ -220,6 +222,8 @@ struct Encoder : public Common
return QuadratureAmplitudeModulation<256, cmplx, code_type>::DIST;
case 10:
return QuadratureAmplitudeModulation<1024, cmplx, code_type>::DIST;
case 12:
return QuadratureAmplitudeModulation<4096, cmplx, code_type>::DIST;
}
return 2;
}

72
qam.hh
View file

@ -258,5 +258,77 @@ struct QuadratureAmplitudeModulation<1024, TYPE, CODE>
}
};
template <typename TYPE, typename CODE>
struct QuadratureAmplitudeModulation<4096, TYPE, CODE>
{
static const int NUM = 4096;
static const int BITS = 12;
typedef TYPE complex_type;
typedef typename TYPE::value_type value_type;
typedef CODE code_type;
static constexpr value_type FAC = 0.8293555858802010;
static constexpr value_type RCP = 63 * FAC;
static constexpr value_type AMP = 1 / RCP;
static constexpr value_type DIST = 2 * AMP;
static constexpr value_type amp(int i)
{
return AMP * i;
}
static code_type quantize(value_type precision, value_type value)
{
value *= DIST * precision;
if (std::is_integral<code_type>::value)
value = std::nearbyint(value);
if (std::is_same<code_type, int8_t>::value)
value = std::min<value_type>(std::max<value_type>(value, -127), 127);
if (std::is_same<code_type, int16_t>::value)
value = std::min<value_type>(std::max<value_type>(value, -32767), 32767);
return value;
}
static void hard(code_type *b, complex_type c)
{
b[ 0] = c.real() < amp(0) ? code_type(-1) : code_type(1);
b[ 1] = c.imag() < amp(0) ? code_type(-1) : code_type(1);
b[ 2] = std::abs(c.real()) < amp(32) ? code_type(-1) : code_type(1);
b[ 3] = std::abs(c.imag()) < amp(32) ? code_type(-1) : code_type(1);
b[ 4] = std::abs(std::abs(c.real())-amp(32)) < amp(16) ? code_type(-1) : code_type(1);
b[ 5] = std::abs(std::abs(c.imag())-amp(32)) < amp(16) ? code_type(-1) : code_type(1);
b[ 6] = std::abs(std::abs(std::abs(c.real())-amp(32))-amp(16)) < amp(8) ? code_type(-1) : code_type(1);
b[ 7] = std::abs(std::abs(std::abs(c.imag())-amp(32))-amp(16)) < amp(8) ? code_type(-1) : code_type(1);
b[ 8] = std::abs(std::abs(std::abs(std::abs(c.real())-amp(32))-amp(16))-amp(8)) < amp(4) ? code_type(-1) : code_type(1);
b[ 9] = std::abs(std::abs(std::abs(std::abs(c.imag())-amp(32))-amp(16))-amp(8)) < amp(4) ? code_type(-1) : code_type(1);
b[10] = std::abs(std::abs(std::abs(std::abs(std::abs(c.real())-amp(32))-amp(16))-amp(8))-amp(4)) < amp(2) ? code_type(-1) : code_type(1);
b[11] = std::abs(std::abs(std::abs(std::abs(std::abs(c.imag())-amp(32))-amp(16))-amp(8))-amp(4)) < amp(2) ? code_type(-1) : code_type(1);
}
static void soft(code_type *b, complex_type c, value_type precision)
{
b[ 0] = quantize(precision, c.real());
b[ 1] = quantize(precision, c.imag());
b[ 2] = quantize(precision, std::abs(c.real())-amp(32));
b[ 3] = quantize(precision, std::abs(c.imag())-amp(32));
b[ 4] = quantize(precision, std::abs(std::abs(c.real())-amp(32))-amp(16));
b[ 5] = quantize(precision, std::abs(std::abs(c.imag())-amp(32))-amp(16));
b[ 6] = quantize(precision, std::abs(std::abs(std::abs(c.real())-amp(32))-amp(16))-amp(8));
b[ 7] = quantize(precision, std::abs(std::abs(std::abs(c.imag())-amp(32))-amp(16))-amp(8));
b[ 8] = quantize(precision, std::abs(std::abs(std::abs(std::abs(c.real())-amp(32))-amp(16))-amp(8))-amp(4));
b[ 9] = quantize(precision, std::abs(std::abs(std::abs(std::abs(c.imag())-amp(32))-amp(16))-amp(8))-amp(4));
b[10] = quantize(precision, std::abs(std::abs(std::abs(std::abs(std::abs(c.real())-amp(32))-amp(16))-amp(8))-amp(4))-amp(2));
b[11] = quantize(precision, std::abs(std::abs(std::abs(std::abs(std::abs(c.imag())-amp(32))-amp(16))-amp(8))-amp(4))-amp(2));
}
static complex_type map(code_type *b)
{
return AMP * complex_type(
b[0]*(b[2]*(b[4]*(b[6]*(b[8]*(b[10]+value_type(2))+value_type(4))+value_type(8))+value_type(16))+value_type(32)),
b[1]*(b[3]*(b[5]*(b[7]*(b[9]*(b[11]+value_type(2))+value_type(4))+value_type(8))+value_type(16))+value_type(32))
);
}
};
#endif