将含变压器噪声生成循环的Matlab代码转换为C语言实现

先帮你理清楚原Matlab代码的核心逻辑：它是通过生成10个不同幅值、不同频率的正弦波分量，再将它们叠加起来得到变压器噪声信号。Matlab靠向量操作简化了循环，但C语言得用显式循环来实现这些逻辑，下面是完整的转换实现：

完整C代码实现 #

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

#define Fs 2000.0    // 采样频率
#define L 4000       // 信号长度

int main() {
    // 计算采样周期
    double T = 1.0 / Fs;
    // 存储噪声信号的数组，初始化为0
    double S[L];
    memset(S, 0, sizeof(S));

    // 噪声幅值数组（对应原Matlab的A）
    double A[] = {16, 18, 28, 33, 38, 41, 59, 43, 40, 58};
    // 噪声分量频率数组（对应原Matlab的forigin）
    double forigin[] = {25, 31.5, 40, 50, 63, 80, 100, 125, 160, 200};
    // 计算分量总数（避免硬编码数字）
    int component_count = sizeof(A) / sizeof(A[0]);

    // 核心循环：叠加每个频率分量到信号S中
    for (int i = 0; i < component_count; i++) {
        double current_amp = A[i];
        double current_freq = forigin[i];
        // 遍历每个时间采样点，计算当前分量的幅值并叠加
        for (int n = 0; n < L; n++) {
            double current_t = n * T;
            // 正弦函数参数用弧度制：2π*频率*时间
            S[n] += current_amp * sin(2 * M_PI * current_freq * current_t);
        }
    }

    // 可添加代码输出信号，比如打印前10个点验证
    printf("前10个噪声信号值：\n");
    for (int n = 0; n < 10; n++) {
        printf("t=%.6f, S=%.4f\n", n*T, S[n]);
    }

    return 0;
}

关键逻辑说明 #

• 循环对应Matlab向量操作：原Matlab里的A(i)*sin(2*pi*forigin(i)*t)是对所有时间点批量计算，在C里我们用嵌套循环实现：外层遍历每个频率分量，内层遍历所有时间采样点，逐个完成叠加计算，和Matlab的逻辑完全一致。
• 内存优化小技巧：我没有单独存储t数组，而是在内层循环里直接计算n*T，这样能节省不少内存空间，如果你需要保留时间点数组，也可以单独定义后提前计算。
• 编译注意事项：因为用到了sin和M_PI，编译时需要链接数学库，比如用GCC编译的话，命令是：gcc transformer_noise.c -o transformer_noise -lm

内容的提问来源于stack exchange，提问作者Kay