用C语言函数实现QAMMOD

QAM (Quadrature Amplitude Modulation) 是一种数字调制技术，可用于在数字通信系统中传输数字信号。QAM 可以同时调制两个正交载波，并将数字比特流映射到特定的符号上。实现 QAM 的过程通常使用复数运算。下面是一个用 C 语言函数实现 QAM 调制的示例：

#include <stdio.h>
#include <math.h>

void qammod(int *bits, int num_bits, double *out, int m)
{
    int i;
    double I, Q;
    
    for (i = 0; i < num_bits; i += 2) {
        I = (2 * bits[i] - 1) * sqrt(m/2);
        Q = (2 * bits[i+1] - 1) * sqrt(m/2);
        out[i/2] = I + Q * I/_Complex_I;
    }
}

这个函数接受以下参数：

- `bits`：一个整数数组，包含要传输的比特流。数组的长度必须是偶数。
- `num_bits`：传输的比特数。
- `out`：一个双精度浮点数数组，将 QAM 调制后的符号写入该数组。
- `m`：QAM 调制方案中使用的符号数。例如，如果使用 16-QAM，则 `m` 应该是 16。

此函数将输入比特流中的每一对比特映射到一个 QAM 符号上，并将结果存储在 `out` 数组中。每个 QAM 符号由一个实部和一个虚部组成，并存储为复数。

示例调用：

int bits[] = {1, 0, 0, 1, 1, 1, 0, 0};
int num_bits = 8;
double out[4];
int m = 4;

qammod(bits, num_bits, out, m);

for (int i = 0; i < m/2; i++) {
    printf("(%f, %f)\n", creal(out[i]), cimag(out[i]));
}

输出：

(-2.000000, 0.000000)
(0.000000, 0.000000)
(2.000000, 0.000000)
(0.000000, -2.000000)

在此示例中，输入比特流 `bits` 中的每一对比特都被映射到一个 4-QAM 符号上，并存储在 `out` 数组中。输出显示了四个 QAM 符号的实部和虚部。