一种基于内点罚函数算法的多光谱温度获取方法
【技术领域】
[0001] 本发明涉及多光谱辐射测温技术，尤其是涉及一种基于内点罚函数算法的多光谱 温度获取方法。
【背景技术】
[0002] 多光谱辐射温度测量技术是一种基于探测和分析待测目标热辐射的非接触测温 技术。该项技术以普朗克定律为理论依据，根据待测目标多光谱辐射信息，利用相关算法实 现真实温度和光谱发射率的同时测量。根据多光谱辐射测温理论，利用普朗克公式建立温 度求解的方程组属于欠定方程组（未知数个数大于方程数）。目前，多光谱辐射测温算法需 要事先假设待测对象的光谱发射率与波长之间的函数关系，如果假设的函数关系与实际情 况相吻合，则反演结果较好，反之则偏差很大。但在实际测量过程中，往往难以确定光谱发 射率与波长之间的函数关系，因此需要构建一种无需假设光谱发射率与波长之间函数关系 的多光谱温度获取方法。
[0003] 现有的测温方法为：
[0004] 如果多波长温度计有n个通道，则第T个通道的输出信号￥,可表示为：
[0006] 式中，2?为只与波长有关与温度无关的检定常数，它与该波长下探测器的光谱响 应率、光学元件透过率、几何尺寸以及第二辐射常数有关，e (A"T)为温T时的目标光谱 发射率。AT为工作波长，C2为第二辐射常数
[0007] 将式（1-1)用维恩公式代替普朗克定律改写为：
[0009] 对于有n个通道的多波长温度计来说共有n个方程，却包含n+1个未知量如果不 找出目标真温T和n个光谱发射率e(A"T)之间的关系，问题难以求解。
[0010] 在多波长辐射测温学领域被普遍认可的一种假设是认为光谱发射率随波长的变 化而变化，光谱发射率可用含有n-1个可调参数的波长函数代替，可表达如下：
[0011] Ine(入，T) =a+b入
[0015] e(入，T) =exp[_Cafa1 入）2]
[0016] 将上面的公式带入式（1-2)中就可得到n个方程，因此可以用曲线拟合法求解方 程获取真温和光谱发射率。
[0017]目前对于多光谱辐射温度计的数据处理，所有学者均使用最小二乘法，但是当光 谱发射率与波长函数与实际不相符时，计算得到的结果与实际相差较大。
【发明内容】
[0018] 本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种精度高、适用范 围广、抗噪性好的基于内点罚函数算法的多光谱温度获取方法。
[0019] 本发明的目的可以通过以下技术方案来实现：
[0020] -种基于内点罚函数算法的多光谱温度获取方法，该方法包括以下步骤：
[0021] 1)根据多波长温度计检测的亮度温度，建立基于内点罚函数算法的目标函数；
[0022] 2)建立基于内点罚函数算法的约束条件；
[0023] 3)通过目标函数和约束条件建立基于内点罚函数算法的多光谱温度模型，并采用 内点罚函数算法求解，获取待测量物体的真温。
[0024] 所述的步骤1)中，亮度温度表达式为：
[0026] 其中，L1为第i个通道测得的亮度温度，T为真温，C2为第二辐射常数，Ai为第i 个通道的工作波长，e(AT)为第i个通道的光谱发射率。
[0027] 所述的步骤1)中，基于内点罚函数算法的目标函数表达式为：
[0030]其中，T1为第i个通道获得的真温，E(T)为真温期望，n为多波长温度计的通道数。[0031] 所述的步骤2)中的约束条件为任何材料的发射率都在（0, 1)之间，约束条件的表 达式为：
[0034]X= [ e(入0T)e(入2,T)…e(入n，T)]T
[0035] 其中，A为2nXn阶常数矩阵，b为2n阶常数列向量，X为n阶光谱发射率矩阵，n 为多波长温度计的通道数，e(An，T)为第n个通道的光谱发射率。
[0036] 所述的步骤3)具体包括以下步骤：
[0037] 31)设定初始光谱发射率e(A1,T)、罚因子y>0、缩小系数V> 1、迭代次数k =1以及精度〇 ;
[0038] 32)根据目标函数和约束条件构造增广函数F，F的表达式为：
[0040] 其中，T1为第i个通道获得的真温，E(T)为真温期望，n为多波长温度计的通道数， A为2nXn阶常数矩阵，b为2n阶常数列向量，X为n阶光谱发射率矩阵；
[0041] 33)求解增广函数F，获得最优的光谱发射率和温度真值使得目标函数趋近于0。
[0042] 与现有技术相比，本发明具有以下优点：
[0043] 一、适用范围广：本发明可适用于所有材料光谱发射率和真温的测量，无需事先假 设光谱发射率与其它参量（波长、温度等）的关系，避免了由于设定的真温及光谱发射率之 间的函数关系与实际不符造成的偏差，直接通过内点罚函数算法获取真温及光谱发射率。
[0044] 二、精度高：1800K时真温时，绝对误差小于18K，相对误差小于1%。
[0045] 三、抗噪性好：在输入5%的电压信号随机误差后，1800K真温时，绝对误差小于 30K，相对误差小于1. 6%，显示了良好的抗噪性，对实测数据也反映了良好的测量效果。
【附图说明】
[0046] 图1为本发明的方法流程图。
[0047] 图2为内点法函数算法求解流程图。
[0048]图3为基于本发明的6种材料真温反演结果图，其中，图（3a)为材料A的真温反演 结果图，图（3b)为材料B的真温反演结果图，图（3c)为材料C的真温反演结果图，图（3d) 为材料D的真温反演结果图，图（3e)为材料E的真温反演结果图，图（3f)为材料F的真温 反演结果图。
【具体实施方式】
[0049] 下面结合附图和具体实施例对本发明进行详细说明。
[0050] 实施例：
[0051] 如图1所示，一种基于内点罚函数算法的多光谱温度获取方法，该方法包括以下 步骤：
[0052] 1)根据多波长温度计的检测的亮度温度，建立基于内点罚函数算法的目标函数， 亮度温度表达式为：
[0054] 其中，L1为第i个通道测得的亮度温度，T为真温，C2为第二辐射常数，Ai为第i 个通道的工作波长，e(AT)为第i个通道的光谱发射率，
[0055] 基于内点罚函数算法的目标函数表达式为：
[0058] 其中，T1为第i个通道获得的真温，E(T)为真温期望，n为多波长温度计的通道数；
[0059] 2)建立基于内点罚函数算法的约束条件，约束条件为任何材料的发射率都在 (〇, 1)之间，约束条件的表达式为：
[0062] X = [e (入oT) e (入2，T)…e (入n，T)]T
[0063] 其中，A为2nXn阶常数矩阵，b为2n阶常数列向量，X为n阶光谱发射率矩阵，n 为多波长温度计的通道数，e(An，T)为第n个通道的光谱发射率；
[0064] 3)通过目标函数和约束条件建立基于内点罚函数算法的多光谱温度模型，并采用 内点罚函数算法求解，获取待测量物体的真温，如图2所示，具体包括以下步骤：
[0065] 31)设定初始光谱发射率e(A1,T)、罚因子y>0、缩小系数V> 1、迭代次数k =1以及精度〇 ;
[0066] 32)根据目标函数和约束条件构造增广函数F，F的表达式为：
[0068]其中，T1为第i个通道获得的真温，E(T)为真温期望，n为多波长温度计的通道数，A为2nXn阶常数矩阵，b为2n阶常数列向量，X为n阶光谱发射率矩阵；
[0069] 33)求解增广函数F，获得最优的光谱发射率和温度真值使得目标函数趋近于0。
[0070]本发明选取了 6种材料A-F作为待测对象，每种材料光谱发射率与波长的关系如 表1所示，分别作为代表线性增加，线性降低和波动等关系。
[0071] 表1 6种目标光谱发射率模型
[0072]
[0074] 表2是在真温1800K时，对6种材料的真温反演结果，绝对误差均小于18K，在电压 输入5%随机误差的情况下，绝对误差小于26K，相对误差小于1.5%，图（3a)-(3f)为6种 材料反演的光谱发射率。
[0075] 表2基于本专利的6种材料真温反演结果
[0076]
[0077] 表3是意大利国家计量院提供的4组亮温数据，作为验证本发明的实验数据。
[0078] 表3意大利国家计量院的4组亮温数据
[0079]
[0080] 表4为根据表3的提供数据做出的反演结果。
[0081] 表4根据表3的提供数据做出的反演结果
[0082]
【主权项】
1. 一种基于内点罚函数算法的多光谱温度获取方法，其特征在于，该方法包括以下步 骤： 1) 根据多波长温度计检测的亮度温度，建立基于内点罚函数算法的目标函数； 2) 建立基于内点罚函数算法的约束条件； 3) 通过目标函数和约束条件建立基于内点罚函数算法的多光谱温度模型，并采用内点 罚函数算法求解，获取待测量物体的真温。2. 根据权利要求1所述的一种基于内点罚函数算法的多光谱温度获取方法，其特征在 于，所述的步骤1)中，亮度温度表达式为：其中，L1为第i个通道测得的亮度温度，T为真温，C2为第二辐射常数，λ i为第i个通 道的工作波长，ε ( λ T)为第i个通道的光谱发射率。3. 根据权利要求1所述的一种基于内点罚函数算法的多光谱温度获取方法，其特征在 于，所述的步骤1)中，基于内点罚函数算法的目标函数表达式为：其中，T1为第i个通道获得的真温，E(T)为真温期望，η为多波长温度计的通道数。4. 根据权利要求1所述的一种基于内点罚函数算法的多光谱温度获取方法，其特征在 于，所述的步骤2)中的约束条件为任何材料的发射率都在（0, 1)之间，约束条件的表达式 为：其中，A为2ηXη阶常数矩阵，b为2η阶常数列向量，X为η阶光谱发射率矩阵，η为多 波长温度计的通道数，ε (λη，Τ)为第η个通道的光谱发射率。5. 根据权利要求1所述的一种基于内点罚函数算法的多光谱温度获取方法，其特征在 于，所述的步骤3)具体包括以下步骤： 31)设定初始光谱发射率￡(人1，!')、罚因子以>0、缩小系数￥>1、迭代次数1^=1以 及精度σ ; 32) 根据目标函数和约束条件构造增广函数F，F的表达式为：其中，T1为第i个通道获得的真温，E (T)为真温期望，η为多波长温度计的通道数，A为 2ηXη阶常数矩阵，b为2η阶常数列向量，X为η阶光谱发射率矩阵； 33) 求解增广函数F，获得最优的光谱发射率和温度真值使得目标函数趋近于0。
【专利摘要】本发明涉及一种基于内点罚函数算法的多光谱温度获取方法，该方法包括以下步骤：1)根据多波长温度计检测的亮度温度，建立基于内点罚函数算法的目标函数；2)建立基于内点罚函数算法的约束条件；3)通过目标函数和约束条件建立基于内点罚函数算法的多光谱温度模型，并采用内点罚函数算法求解，获取待测量物体的真温。与现有技术相比，本发明具有精度高、适用范围广、抗噪性好等优点。
【IPC分类】G06F19/00
【公开号】CN105069284
【申请号】CN201510433829
【发明人】邢键, 吴飞, 孔勇, 张玉金
【申请人】上海工程技术大学
【公开日】2015年11月18日
【申请日】2015年7月22日