// 计算平均化时间常数
float calculate_averaging_time(float cap_value, float input_freq) {
    // 时间常数 τ = R × C，其中R为内部电阻约630kΩ
    const float internal_resistance = 630e3;  // 630kΩ
    float time_constant = internal_resistance * cap_value;
    
    // 对于正弦波输入，建议τ > 10/f
    float min_time_constant = 10.0 / input_freq;
    
    printf("电容值: %.1f μF\n", cap_value * 1e6);
    printf("时间常数: %.2f ms\n", time_constant * 1000);
    printf("最小建议值: %.2f ms\n", min_time_constant * 1000);
    
    return time_constant;
}

// 选择合适的平均化电容
float select_averaging_cap(float input_freq, float max_ripple_percent) {
    // 根据输入频率和纹波要求选择电容
    const float internal_r = 630e3;
    
    // 时间常数需要满足: τ > 10/(2πf × ripple_factor)
    float ripple_factor = max_ripple_percent / 100.0;
    float required_tau = 10.0 / (2 * 3.14159 * input_freq * ripple_factor);
    
    float required_cap = required_tau / internal_r;
    
    // 选择标准电容值
    float std_caps[] = {1e-6, 10e-6, 47e-6, 100e-6};  // 1μF, 10μF, 47μF, 100μF
    
    for(int i = 0; i < 4; i++) {
        if(std_caps[i] >= required_cap) {
            return std_caps[i];
        }
    }
    
    return 100e-6;  // 默认100μF
}