[PIC®/AVR®/dsPIC®产品] 采用DeepSeek编写NTC测温度函数

Vcc --- [参考电阻] --- [NTC热敏电阻] --- GND

#include <math.h>



// 常量定义

#define R0      10000.0f  // 25℃时的阻值（10KΩ）

#define T0      298.15f   // 25℃对应的开尔文温度

#define B       3470.0f   // B值



/**

 * [url=home.php?mod=space&uid=247401]@brief[/url] 根据NTC热敏电阻的阻值计算摄氏温度

 * @param resistance: 当前测量的热敏电阻阻值（单位：欧姆）

 * [url=home.php?mod=space&uid=266161]@return[/url] 当前温度（单位：摄氏度）

 */

float CalculateTemperature(float resistance) {

    // 检查阻值是否有效

    if (resistance <= 0) {

        return -273.15f; // 返回绝对零度表示无效输入

    }



    // 计算 ln(R/R0)

    float ln_R_R0 = logf(resistance / R0);



    // 计算 1/T

    float inv_T = (1.0f / T0) + (1.0f / B) * ln_R_R0;



    // 计算开尔文温度 T

    float T_K = 1.0f / inv_T;



    // 转换为摄氏温度

    float T_C = T_K - 273.15f;



    return T_C;

}

#include <stdio.h>



// PID参数

float Kp = 10.0;  // 比例增益

float Ki = 0.1;   // 积分增益

float Kd = 1.0;   // 微分增益



// PID变量

float integral = 0.0;  // 积分项

float previous_error = 0.0;  // 上一次的误差



// 采样时间间隔（单位：秒）

float delta_t = 1.0;



/**

 * @brief PID控制算法

 * @param setpoint: 设定温度

 * @param actual: 实际温度

 * @return 控制输出（如加热膜的功率）

 */

float PID_Control(float setpoint, float actual) {

    // 计算误差

    float error = setpoint - actual;



    // 计算积分项

    integral += error * delta_t;



    // 计算微分项

    float derivative = (error - previous_error) / delta_t;



    // 计算控制输出

    float output = Kp * error + Ki * integral + Kd * derivative;



    // 保存当前误差

    previous_error = error;



    return output;

}



int main() {

    // 示例：设定温度和实际温度

    float setpoint = 37.5;  // 设定温度

    float actual = 36.0;    // 实际温度



    // 模拟PID控制

    for (int i = 0; i < 10; i++) {

        float output = PID_Control(setpoint, actual);

        printf("控制输出: %.2f\n", output);



        // 模拟实际温度变化

        actual += 0.2;  // 假设每次采样后温度上升0.2℃

    }



    return 0;

}

#define K 273.15

#define T0 (25.0+K)

#define R0 10000.0

#define B 3455.0        



/*

电阻值换算到温度

adc_val:NTC当前的阻抗值，单位欧姆

返回值：摄氏度

*/



float Res_to_Tem(float adc_val)

{

        float temp;

        temp=1/((1/T0)+(log(adc_val/R0)/B));

        return temp-K;

}