python——温度传感器实验代码实现

视频我们就放在这里，B站的链接放在评论区！~

用python实现温度传感器数据处理

温度传感器——RTD

#在这里，我们将使用一段代码来解决温度传感器的数据处理问题
#我们直接copy之前写过的获取表格的函数
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import font_manager

# 设置matplotlib配置参数，使其支持中文显示
plt.rcParams['font.sans-serif'] = ['SimHei']  # 使用黑体
plt.rcParams['axes.unicode_minus'] = False  # 解决保存图像是负号'-'显示为方块的问题
#那么这次，我们对这个函数进行详细的分析
def get_user_input_table():
    n = int(input("请输入表格的行数: "))
    m = int(input("请输入表格的列数: "))
    # 初始化一个空的二维列表（表格）
    table = [[0 for _ in range(m)] for _ in range(n)]
    # 在这句代码里，_ in range(n)是一个外部循环，[0 for _ in range(m)]是一个内循环，它嵌套在外部循环中。
    # 每次外部循环迭代时，内部循环都会执行m次，m是表格的列数，n是表格的行数
    # 内部循环用于生成当前行的每个元素。
    # _ 是一个占位符，它在这里没有实际用途，只是表示我们不关心循环迭代的具体值。
    # 那么0代表的就是给占位符赋值0

    # 循环读取每一行和每一列的数据
    for i in range(n):
        print(f"请输入第{i + 1}行的{m}个元素（用空格分隔）: ")
        # 读取一行数据，并用空格分割成列表
        row_data = input().split()
        # split(): 这是字符串对象的一个方法，用于将字符串按照指定的分隔符（如果没有指定分隔符，则默认使用空白字符，包括空格、换行符、制表符等）分割成一个列表。
        # 组合在一起，input().split() 读取用户输入的一行文本，然后按照空白字符分割这行文本，最终返回一个字符串列表。

        # 确保用户输入了正确数量的元素
        if len(row_data) != m:
            print("错误：输入的元素数量与列数不匹配，请重新输入。")
            # 如果输入错误，重新读取当前行
            i -= 1
            continue

            # 将读取的数据转换为适当的类型（这里假设为整数或浮点数）并存储
        try:# https://www.runoob.com/python/python-exceptions.html这里是一个“异常处理”的学习网站
            table[i] = [float(value) for value in row_data]
            # 如果在列表推导式中有任何 value 不能转换为浮点数，那么 ValueError 异常会被触发。
            # 但是，由于这段代码位于 try 块中，所以不会立即导致程序崩溃。
            # 相反，程序会跳过 try 块中剩余的代码，转而执行与之匹配的 except 块中的代码（如果有的话）。
        except ValueError:
            print("错误：输入的不是有效的数字，请重新输入。")
            # 如果转换失败，重新读取当前行
            i -= 1
            continue
    return table
# 最后返回一个列表
user_table_RTD = get_user_input_table()
temperature = np.array(user_table_RTD[0])
resistance = np.array(user_table_RTD[1])

# np.polyfit(velocity, frequency, 1.0) 是NumPy库中的一个函数，用于进行多项式拟合。
# 这行代码的目的是根据给定的 velocity（速度）和 frequency（频率）数据点来拟合一条直线，并返回该直线的斜率和截距。
# 具体解释如下：
# np.polyfit: 这是NumPy库中的多项式拟合函数。
# velocity: 这是x轴上的数据点，代表速度。
# frequency: 这是y轴上的数据点，代表频率。
# 1.0: 这是多项式的阶数。在这里，1.0 表示我们要拟合一条直线（一阶多项式）。
# 函数返回两个值：
#
# slope: 直线的斜率。
# intercept: 直线的截距。
# 一旦你有了这两个值，你就可以用它们来定义直线方程 y = mx + c，其中 m 是斜率（对应于 slope），c 是截距（对应于 intercept），x 是速度，y 是频率。

# 进行线性拟合
slope, intercept = np.polyfit(temperature, resistance, 1.0)
resistance = slope * temperature + intercept
plt.plot(temperature, resistance, '-^k')
plt.xlabel("温度（℃）")
plt.ylabel("电阻（Ω）")
plt.text(0.05, 0.95, f'y = {slope:.4f} * x + {intercept:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
plt.text(0.05, 0.80, f'α = {slope/intercept:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
for a, b in zip(temperature, resistance):  # 添加这个循坏显示坐标
    plt.text(a, b, f'({int(a)}, {int(b)})', ha='center', va='bottom', fontsize=10)
plt.show()

温度传感器——PTC

import numpy as np
import matplotlib.pyplot as plt

from matplotlib import font_manager
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False

def get_user_input_table():
    n = int(input("请输入表格的行数: "))
    m = int(input("请输入表格的列数: "))
    table = [[0 for _ in range(m)] for _ in range(n)]
    for i in range(n):
        print(f"请输入第{i + 1}行的{m}个元素（用空格分隔）: ")
        row_data = input().split()
        if len(row_data) != m:
            print("错误：输入的元素数量与列数不匹配，请重新输入。")
            i -= 1
            continue


        try:
            table[i] = [float(value) for value in row_data]
        except ValueError:
            print("错误：输入的不是有效的数字，请重新输入。")
            i -= 1
            continue
    return table
user_table_PTC = get_user_input_table()
temperature = np.array(user_table_PTC[0])
resistance = np.array(user_table_PTC[1])
height = (9 / 16) * 13.3
figsize = (13.3, height)
# 创建画板
fig = plt.figure(figsize=figsize)
plt.plot(temperature, resistance, '-^k')
plt.xlabel("温度（℃）")
plt.ylabel("电阻（kΩ）")
for a, b in zip(temperature, resistance):  # 添加这个循坏显示坐标
    plt.text(a, b, f'({int(a)}, {b:.3f})', ha='center', va='bottom', fontsize=10)
plt.show()

温度传感器——PN结

import numpy as np
import matplotlib.pyplot as plt

def get_user_input_table():
    n = int(input("请输入表格的行数: "))
    m = int(input("请输入表格的列数: "))
    table = [[0 for _ in range(m)] for _ in range(n)]
    for i in range(n):
        print(f"请输入第{i + 1}行的{m}个元素（用空格分隔）: ")
        row_data = input().split()
        if len(row_data) != m:
            print("错误：输入的元素数量与列数不匹配，请重新输入。")
            i -= 1
            continue


        try:
            table[i] = [float(value) for value in row_data]
        except ValueError:
            print("错误：输入的不是有效的数字，请重新输入。")
            i -= 1
            continue
    return table
# 设置matplotlib配置参数，使其支持中文显示
plt.rcParams['font.sans-serif'] = ['SimHei']  # 使用黑体
plt.rcParams['axes.unicode_minus'] = False  # 解决保存图像是负号'-'显示为方块的问题

user_table_PN = get_user_input_table()
temperature = np.array(user_table_PN[0])
resistance = np.array(user_table_PN[1])
height = (9 / 16) * 13.3
figsize = (13.3, height)
slope, intercept = np.polyfit(temperature, resistance, 1.0)
resistance = slope * temperature + intercept
plt.plot(temperature, resistance, '-^k')
plt.xlabel("温度（℃）")
plt.ylabel("电压（mV）")
plt.text(0.65, 0.95, f'y = {slope:.4f} * x + {intercept:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
plt.text(0.65, 0.80, f'dU/dT = {slope:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
for a, b in zip(temperature, resistance):  # 添加这个循坏显示坐标
    plt.text(a, b, f'({int(a)}, {int(b)})', ha='center', va='bottom', fontsize=10)
plt.show()

温度传感器——NTC

import numpy as np
import matplotlib.pyplot as plt

def get_user_input_table():
    n = int(input("请输入表格的行数: "))
    m = int(input("请输入表格的列数: "))
    table = [[0 for _ in range(m)] for _ in range(n)]
    for i in range(n):
        print(f"请输入第{i + 1}行的{m}个元素（用空格分隔）: ")
        row_data = input().split()
        if len(row_data) != m:
            print("错误：输入的元素数量与列数不匹配，请重新输入。")
            i -= 1
            continue


        try:
            table[i] = [float(value) for value in row_data]
        except ValueError:
            print("错误：输入的不是有效的数字，请重新输入。")
            i -= 1
            continue
    return table
# 设置matplotlib配置参数，使其支持中文显示
plt.rcParams['font.sans-serif'] = ['SimHei']  # 使用黑体
plt.rcParams['axes.unicode_minus'] = False  # 解决保存图像是负号'-'显示为方块的问题

user_table_NTC = get_user_input_table()
user_table_NTC[0] = np.array(user_table_NTC[0])
x = (1.0/(user_table_NTC[0]+273))*1000
resistance = np.array(user_table_NTC[1])
y = np.log(resistance)
height = (9 / 16) * 13.3
figsize = (13.3, height)

slope, intercept = np.polyfit(x, y, 1.0)
y = slope * x + intercept
plt.plot(x, y, '-^k')
plt.xlabel("温度（10-3/K）")
plt.ylabel("电阻（lnΩ）")
plt.text(0.05, 0.95, f'y = {slope:.4f} * x + {intercept:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
plt.text(0.05, 0.80, f'β = {slope:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
for a, b in zip(x, y):  # 添加这个循坏显示坐标
    plt.text(a, b, f'({a:.3f}, {b:.3f})', ha='center', va='bottom', fontsize=10)
plt.show()

温度传感器——LM35

import numpy as np
import matplotlib.pyplot as plt

def get_user_input_table():
    n = int(input("请输入表格的行数: "))
    m = int(input("请输入表格的列数: "))
    table = [[0 for _ in range(m)] for _ in range(n)]
    for i in range(n):
        print(f"请输入第{i + 1}行的{m}个元素（用空格分隔）: ")
        row_data = input().split()
        if len(row_data) != m:
            print("错误：输入的元素数量与列数不匹配，请重新输入。")
            i -= 1
            continue


        try:
            table[i] = [float(value) for value in row_data]
        except ValueError:
            print("错误：输入的不是有效的数字，请重新输入。")
            i -= 1
            continue
    return table
# 设置matplotlib配置参数，使其支持中文显示
plt.rcParams['font.sans-serif'] = ['SimHei']  # 使用黑体
plt.rcParams['axes.unicode_minus'] = False  # 解决保存图像是负号'-'显示为方块的问题

user_table_LM35 = get_user_input_table()
temperature = np.array(user_table_LM35[0])
resistance = np.array(user_table_LM35[1])
height = (9 / 16) * 13.3
figsize = (13.3, height)
slope, intercept = np.polyfit(temperature, resistance, 1.0)
resistance = slope * temperature + intercept
plt.plot(temperature, resistance, '-^k')
plt.xlabel("温度（℃）")
plt.ylabel("电压（mV）")
plt.text(0.05, 0.95, f'y = {slope:.4f} * x + {intercept:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
plt.text(0.05, 0.80, f'dU/dT = {slope:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
for a, b in zip(temperature, resistance):  # 添加这个循坏显示坐标
    plt.text(a, b, f'({int(a)}, {int(b)})', ha='center', va='bottom', fontsize=10)
plt.show()

温度传感器——AD590

import numpy as np
import matplotlib.pyplot as plt

def get_user_input_table():
    n = int(input("请输入表格的行数: "))
    m = int(input("请输入表格的列数: "))
    table = [[0 for _ in range(m)] for _ in range(n)]
    for i in range(n):
        print(f"请输入第{i + 1}行的{m}个元素（用空格分隔）: ")
        row_data = input().split()
        if len(row_data) != m:
            print("错误：输入的元素数量与列数不匹配，请重新输入。")
            i -= 1
            continue


        try:
            table[i] = [float(value) for value in row_data]
        except ValueError:
            print("错误：输入的不是有效的数字，请重新输入。")
            i -= 1
            continue
    return table
# 设置matplotlib配置参数，使其支持中文显示
plt.rcParams['font.sans-serif'] = ['SimHei']  # 使用黑体
plt.rcParams['axes.unicode_minus'] = False  # 解决保存图像是负号'-'显示为方块的问题

user_table_LM35 = get_user_input_table()
temperature = np.array(user_table_LM35[0])
resistance = np.array(user_table_LM35[1])
height = (9 / 16) * 13.3
figsize = (13.3, height)
slope, intercept = np.polyfit(temperature, resistance, 1.0)
resistance = slope * temperature + intercept
plt.plot(temperature, resistance, '-^k')
plt.xlabel("温度（℃）")
plt.ylabel("电压（mV）")
plt.text(0.05, 0.95, f'y = {slope:.4f} * x + {intercept:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
plt.text(0.05, 0.80, f'dU/dT = {slope:.4f}', transform=plt.gca().transAxes, fontsize=36,
         verticalalignment='top', horizontalalignment='left')
for a, b in zip(temperature, resistance):  # 添加这个循坏显示坐标
    plt.text(a, b, f'({int(a)}, {b:.3f})', ha='center', va='bottom', fontsize=10)
plt.show()