本發(fā)明涉及溫度測量技術領域，尤其涉及一種溫度場的溫度測量方法及裝置，具體涉及一種結構簡單、占用空間小、測量精度高且對待測溫物體的影響小的溫度場的溫度測量方法及裝置。

背景技術：

溫度是表征物體冷熱程度的物理量，是重要的熱工參數之一，是各種物質的物理化學變化的重要條件，在工農業(yè)生產及生活中均有大量的溫度測量需要。目前，溫度測量的方法已達數十種之多，如膨脹測溫法、壓力溫度測溫法、電學測溫法、磁學測溫法、聲學測溫法和頻率測溫法。電學測溫法中的熱電偶溫度計和熱電阻溫度計是當前工業(yè)上使用極為廣泛的測溫儀器。熱電阻是中低溫區(qū)最常用的一種測溫元件。一般情況下，測量500℃以上的較高溫度時采用熱電偶。

當前工業(yè)上常用的測溫技術，不論是用于測量高溫的熱電偶，還是用于測量中低溫的熱電阻，測量的數值反應的是探頭所在位置的單一點的溫度，無法利用一套測溫裝置同時對多處溫度、有限長度或者平面溫度分布進行測量。若采用多個溫度傳感器測量則存在接線較多且復雜的問題，在工業(yè)應用上，很多設備由于空間限制無法安裝過多的溫度傳感器。除此之外，在測量有限長度或者平面溫度分布時，由于每個傳感器都有鎧甲，多個溫度傳感器同時安裝在設備內部，不但會占用大量的空間，同時對設備的熱特性造成較大的影響，進而降低了測量精度。具體地，圖1為帶端子箱的鎧裝熱電阻，其包括套管2及端子箱4；圖2為圖1中A部件的放大圖，用于測量待測點處的溫度值，具體在套管2內設有填充有氧化鎂1的陶瓷封裝型元件3；圖3為帶引線的鎧裝熱電阻，其包括套管2及引線5；圖4為熱電偶，其包括不銹鋼套管8，不銹鋼套管8的一端設有接線盒11，接線盒11的下端設有引線口10，不銹鋼套管8的另一端設有瓷絕緣套管7及焊點6；熱電偶通過安裝固定件9固定于待測物體上。綜上，現有的熱電偶和熱電阻的測溫部件，均在保護套管的頭部且很小，保護套管和引線占用體積較大，無法使用一個測溫元件實現線溫度場和面溫度場的測量。如果采用多個熱電偶(熱電阻)測量時，存在占用被測空間較大，接線較為復雜的問題。

技術實現要素：

(一)要解決的技術問題

本發(fā)明的目的是：提供一種結構簡單、占用空間小、測量精度高且對待測溫物體的影響小的溫度場的溫度測量方法及裝置，以解決現有的測溫裝置體積大、占用被測空間大且只能測量探頭所在位置的單一點的溫度的問題。

(二)技術方案

為了解決上述技術問題，本發(fā)明提供了一種溫度場的溫度測量方法，利用測溫模塊的電阻測量信息，結合熱流平衡理論以及禁忌搜索算法計算測溫模塊上的溫度分布，具體步驟如下：

步驟1、將所述測溫模塊劃分為m個網格單元，其中，m為不小于3的正整數；

步驟2、采集流經所述測溫模塊上的電流、所述測溫模塊兩端的電壓及每個所述網格單元的位置數據，計算獲得所述測溫模塊的等效電阻R_e及每個所述網格單元的溫度初始值及導熱系數初始值；

步驟3、結合所述測溫模塊的密度值、比熱容值、位置數據、每個所述網格單元的溫度初始值及導熱系數初始值，基于第一類邊界條件，建立所述測溫模塊的熱平衡方程，具體為：

$\mathrm{\ρ}c\frac{\∂t}{\∂\mathrm{\τ}}=\frac{\∂}{\∂x}\left(\mathrm{\λ}\frac{\∂t}{\∂x}\right)+\frac{\∂}{\∂y}\left(\mathrm{\λ}\frac{\∂t}{\∂y}\right)$

式中，ρ為所述測溫模塊的密度值；c為所述測溫模塊的比熱容值；λ為導熱系數值；x、y為網格單元的坐標值，t為溫度值；

步驟4、對所述測溫模塊的熱平衡方程進行離散化，采用禁忌搜索算法進行計算得到并計算獲取所述測溫模塊的總電阻值比較所述等效電阻R_e與總電阻值的差值，并當值時，停止計算，進而得到所述測溫模塊上的溫度分布；其中，i指的是沿x軸方向上的第i個網格單元；j指的是沿y軸方向上的第j個網格單元；i×j＝m，且i、j均為正整數；n指的是計算步驟，n≥1。

其中，所述步驟2還包括同時采集所述測溫模塊兩端的溫度值，將所述測溫模塊兩端的溫度值作為所述測溫模塊兩端的網格單元的溫度初始值，并結合所述測溫模塊的密度值、比熱容值、位置數據、剩余每個所述網格單元的溫度初始值及導熱系數初始值，采用禁忌搜索算法對剩余每個所述網格單元的溫度進行計算，直至滿足的條件，以得到所述測溫模塊上的溫度分布。

其中，所述步驟2還包括同時采集所述測溫模塊任意一端的溫度值，將所述測溫模塊任意一端的溫度值作為所述測溫模塊任意一端的網格單元的溫度初始值，并結合所述測溫模塊的密度值、比熱容值、位置數據、剩余每個所述網格單元的溫度初始值及導熱系數初始值，采用禁忌搜索算法對剩余每個所述網格單元的溫度進行計算，直至滿足的條件，以得到所述測溫模塊上的溫度分布。

其中，在步驟4中，的計算公式為：

$t_{i,j}^n=\frac{{\mathrm{\λ}}_{i-1}^{n-1}\·t_{i-1,j}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}\·t_{i+1,j}^{n-1}+{\mathrm{\λ}}_{j-1}^{n-1}\·t_{i,j-1}^{n-1}+{\mathrm{\λ}}_{j+1}^{n-1}\·t_{i,j+1}^{n-1}}{{\mathrm{\λ}}_{i-1}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}+{\mathrm{\λ}}_{j-1}^{n-1}+{\mathrm{\λ}}_{j+1}^{n-1}}$

式中，λ為導熱系數值；t為溫度值；i指的是沿x軸方向上的第i個網格單元；j指的是沿y軸方向上的第j個網格單元；i×j＝m；n指的是計算步驟，n≥1；其中，i、j與n均為正整數。

其中，所述測溫模塊為線型或平面型結構。

本發(fā)明還提供了一種溫度場的溫度測量裝置，其包括：

測溫模塊，設于待測物體上；

電流源模塊，設于所述測溫模塊的兩端，為所述測溫模塊提供電流；

電壓測量模塊，設于所述測溫模塊的兩端，用于測量所述測溫模塊兩端的壓降值；

溫度場轉換模塊，分別與所述測溫模塊、電流源模塊及電壓測量模塊連接。

其中，所述溫度場轉換模塊包括：

第一模塊，用于“將所述測溫模塊劃分為m個網格單元其中，m為不小于3的正整數”；

第二模塊，用于“采集流經所述測溫模塊上的電流、所述測溫模塊兩端的電壓及每個所述網格單元的位置數據，計算獲得所述測溫模塊的等效電阻R_e及每個所述網格單元的溫度初始值及導熱系數初始值”；

第三模塊，用于“結合所述測溫模塊的密度值、比熱容值、位置數據、每個所述網格單元的溫度初始值及導熱系數初始值，基于第一類邊界條件，建立所述測溫模塊的熱平衡方程”；

第四模塊，用于“對所述測溫模塊的熱平衡方程進行離散化，采用禁忌搜索算法進行計算得到并計算獲取所述測溫模塊的總電阻值比較所述等效電阻R_e與總電阻值的差值，并當值時，停止計算，進而得到所述測溫模塊上的溫度分布；其中，i指的是沿x軸方向上的第i個網格單元；j指的是沿y軸方向上的第j個網格單元；i×j＝m，且i、j均為正整數；n指的是計算步驟，n≥1”。

其中，該溫度場的溫度測量裝置還包括一個設于所述測溫模塊任意一端的輔助測溫元件，且所述輔助測溫元件與所述溫度場轉換模塊連接。

其中，該溫度場的溫度測量裝置還包括兩個分別設于所述測溫模塊兩端的輔助測溫元件，且所述兩個輔助測溫元件分別與所述溫度場轉換模塊連接。

其中，該溫度場的溫度測量裝置還包括與所述溫度場轉換模塊連接的環(huán)境溫度測量模塊。

(三)有益效果

本發(fā)明的上述技術方案具有如下優(yōu)點：本發(fā)明提供了一種溫度場的溫度測量方法，利用測溫模塊的電阻測量信息，結合熱流平衡理論以及禁忌搜索算法計算測溫模塊上的溫度分布，具體包括：步驟1、將測溫模塊劃分為m個網格單元；步驟2、采集流經測溫模塊上的電流、測溫模塊兩端的電壓及每個網格單元的位置數據，計算獲得測溫模塊的等效電阻R_e及每個網格單元的溫度初始值及導熱系數初始值；步驟3、結合測溫模塊的密度值、比熱容值、位置數據、每個網格單元的溫度初始值及導熱系數初始值，基于第一類邊界條件，建立測溫模塊的熱平衡方程；步驟4、對測溫模塊的熱平衡方程進行離散化，采用禁忌搜索算法進行計算得到并計算獲取測溫模塊的總電阻值比較等效電阻R_e與總電阻值的差值，并當值時，停止計算，進而得到測溫模塊上的溫度分布。本發(fā)明通過上述測量方法，能一次性地得到待測物體上的溫度分布，測量精度高且既可以獲取線性溫度場，也可獲取平面溫度場的溫度分布；同時，本申請?zhí)峁┑臏y量裝置占用待測物體的空間小，對待測溫度場的影響較小，對待測溫度物體的損傷較小。

附圖說明

圖1是現有技術中帶端子箱的鎧裝熱電阻的結構圖；

圖2為圖1中A部件的結構放大圖；

圖3是現有技術中帶引線的鎧裝熱電阻的結構圖；

圖4是現有技術中熱電偶的結構圖；

圖5是本發(fā)明一種溫度場的溫度測量方法及其裝置實施例的溫度測量方法的主要操作步驟流程圖；

圖6是本發(fā)明一種溫度場的溫度測量方法及其裝置實施例的溫度測量裝置的主要結構圖；

圖7是本發(fā)明一種溫度場的溫度測量方法及其裝置實施例的測溫模塊劃分為五段結構示意圖。

圖中：1：氧化鎂；2：套管；3：陶瓷封裝型元件；4：端子箱；5：引線；6：焊點；7：瓷絕緣套管；8：不銹鋼套管；9：安裝固定件；10：引線口；11：接線盒；A：測溫模塊；B：輔助測溫元件；C：電流源模塊；D：電壓測量模塊；E：溫度場轉換模塊；F：環(huán)境溫度測量模塊。

具體實施方式

為使本發(fā)明實施例的目的、技術方案和優(yōu)點更加清楚，下面將結合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例是本發(fā)明的一部分實施例，而不是全部的實施例?；诒景l(fā)明中的實施例，本領域普通技術人員在沒有做出創(chuàng)造性勞動的前提下所獲得的所有其他實施例，都屬于本發(fā)明保護的范圍。

特別的，本發(fā)明主要采用電阻的溫度效應和熱平衡狀態(tài)下熱流平衡相結合的方式，通過數值計算得到待測溫度場的溫度分布；其中，測溫范圍為-100℃～1800℃。在本實施例中，采用禁忌搜索算法來計算測溫模塊上的溫度分布，還可以采用其它的全局逐步尋優(yōu)算法，如遺傳算法或模擬退火算法等。

如圖5至圖6所示，本發(fā)明實施例提供了一種溫度場的溫度測量方法，利用測溫模塊的電阻測量信息，結合熱流平衡理論以及禁忌搜索算法計算測溫模塊上的溫度分布，具體步驟如下：

步驟1、將測溫模塊A劃分為m個網格單元，其中，m為不小于3的正整數；其中，每個網格單元內的溫度是均勻一致的；

步驟2、采集流經測溫模塊A上的電流、測溫模塊A兩端的電壓及每個網格單元的位置數據，計算獲得測溫模塊A的等效電阻R_e及每個網格單元的溫度初始值及導熱系數初始值；其中，等效電阻R_e＝測溫模塊A兩端的電壓/流經測溫模塊A上的電流；每個網格單元的初始電阻值＝等效電阻R_e/m，根據初始電阻值獲取溫度初始值及導熱系數初始值；

步驟3、結合測溫模塊A的密度值、比熱容值、位置數據、每個網格單元的溫度初始值及導熱系數初始值，測溫模塊A處于熱平衡狀態(tài)，基于第一類邊界條件，假設網格單元的外表面溫度為當前位置的環(huán)境溫度，由于網格單元熱物性及形狀可以認為是薄材，則網格單元的溫度為其所在的環(huán)境溫度，建立測溫模塊A的熱平衡方程，具體為：

$\mathrm{\ρ}c\frac{\∂t}{\∂\mathrm{\τ}}=\frac{\∂}{\∂x}\left(\mathrm{\λ}\frac{\∂t}{\∂x}\right)+\frac{\∂}{\∂y}\left(\mathrm{\λ}\frac{\∂t}{\∂y}\right)$

式中，ρ為測溫模塊A的密度值；c為測溫模塊A的比熱容值；λ為導熱系數值；x、y為網格單元的坐標值，t為溫度值；其中，當測溫模塊A用于線性溫度場時，x≠0、y＝0；當測溫模塊A用于平面溫度場時，x≠0、y≠0；

步驟4、對測溫模塊A的熱平衡方程進行離散化，采用禁忌搜索算法進行計算得到并計算獲取測溫模塊A的總電阻值比較等效電阻R_e與總電阻值的差值，并當值時，停止計算，進而得到測溫模塊A上的溫度分布；其中，i指的是沿x軸方向上的第i個網格單元；j指的是沿y軸方向上的第j個網格單元；i×j＝m，且i、j均為正整數；n指的是計算步驟，n≥1。特別的，在本實施例中，A值為0.001，具體地，根據所需的計算精度來相應確定A值；所需精度越高，則同等條件下所花費的計算時間越久，最終計算得到的溫度分布值更接近于真實值。

步驟2還包括同時采集測溫模塊A兩端的溫度值，將測溫模塊A兩端的溫度值作為測溫模塊A兩端的網格單元的溫度初始值，并結合測溫模塊A的密度值、比熱容值、位置數據、剩余每個網格單元的溫度初始值及導熱系數初始值，采用禁忌搜索算法對剩余每個網格單元的溫度進行計算，直至滿足的條件，以得到測溫模塊A上的溫度分布。

步驟2還包括同時采集測溫模塊A任意一端的溫度值，將測溫模塊A任意一端的溫度值作為測溫模塊A任意一端的網格單元的溫度初始值，并結合測溫模塊A的密度值、比熱容值、位置數據、剩余每個網格單元的溫度初始值及導熱系數初始值，采用禁忌搜索算法對剩余每個網格單元的溫度進行計算，直至滿足的條件，以得到測溫模塊A上的溫度分布。

在步驟4中，當測溫模塊A用于平面溫度場時，的計算公式為：

$t_{i,j}^n=\frac{{\mathrm{\λ}}_{i-1}^{n-1}\·t_{i-1,j}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}\·t_{i+1,j}^{n-1}+{\mathrm{\λ}}_{j-1}^{n-1}\·t_{i,j-1}^{n-1}+{\mathrm{\λ}}_{j+1}^{n-1}\·t_{i,j+1}^{n-1}}{{\mathrm{\λ}}_{i-1}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}+{\mathrm{\λ}}_{j-1}^{n-1}+{\mathrm{\λ}}_{j+1}^{n-1}}$

式中，λ為導熱系數值；t為溫度值；i指的是沿x軸方向上的第i個網格單元；j指的是沿y軸方向上的第j個網格單元；i×j＝m；n指的是計算步驟，n≥1；其中，i、j與n均為正整數。

在步驟4中，當測溫模塊A用于線性溫度場時，j不取值，則的計算公式為：

$t_i^n=\frac{{\mathrm{\λ}}_{i-1}^{n-1}{t}_{i-1}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}{t}_{i+1}^{n-1}}{{\mathrm{\λ}}_{i-1}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}}$

式中，λ為導熱系數值；t為溫度值；i指的是第i個網格單元，1≤i≤m；n指的是計算步驟，n≥1；其中，i與n均為正整數。

通過上式計算公式得到新一輪的測溫模塊A的溫度分布，然后采用禁忌搜索算法進生成多組測溫模塊A的溫度分布解，即按照禁忌搜索表對每個溫度點進行△t(隨著運算的進行可以調小△t)的加減，生成多組測溫模塊A的溫度分布解，進而再計算得到每組測溫模塊A總電阻值

具體地，測溫模塊A為線型或平面型結構。在本實施例中，當測溫模塊A用于線性溫度場時，采用直徑為0.1-3mm的鉑銠合金導線；當測溫模塊A用于平面溫度場時，采用厚度在0.3-4mm的帶形導體。

優(yōu)選地，測溫模塊A采用導線材質包括不限于：鉑銠合金、純鉑、鎳鉻合金、鎳硅合金、純銅、銅鎳合金、鎳鉻硅、銅鎳和鐵鉻鋁合金等，根據待測溫度選擇材質。其形狀可以根據待測溫度場布置成直線型、L型、U型、S型，測溫模塊A可以直接使用或者放入不銹鋼、剛玉護套內，或者平放到半環(huán)形護套上；平面型測溫模塊A可以按照需要進行折疊。當測量導體(且電導率隨溫度而改變)溫度分布時，可以直接將待測物體作為測溫模塊A使用。

本發(fā)明還提供了一種溫度場的溫度測量裝置，其包括：

測溫模塊A，設于待測物體上；

電流源模塊C，設于測溫模塊A的兩端，為測溫模塊A提供電流；在本實施例中，根據所選用的測溫模塊A的阻值特性，電流源模塊C的所需電流為2-500mA。

電壓測量模塊D，設于測溫模塊A的兩端，用于測量測溫模塊A兩端的壓降值；在本實施例中，電壓測量模塊D的電壓測量范圍為0-12V。

溫度場轉換模塊E，分別與測溫模塊A、電流源模塊C及電壓測量模塊D連接。具體地，根據測溫模塊A的結構參數、輔助測溫元件B的位置和溫度數據、電流源模塊C提供的電流數據、電壓測量模塊D提供的電壓數據，得到測溫模塊A的等效電阻值，采用禁忌搜索算法進行計算并獲取測溫模塊A的總電阻值，比較等效電阻與總電阻值的差值，并當兩者差值小于A值時，停止計算，進而得到測溫模塊A上的溫度分布。特別的，溫度場轉換模塊E也可以采用嵌入式系統。其中，溫度場轉換模塊E為電腦或微處理器或單片機。

具體地，在本實施例中，溫度場轉換模塊E包括：

第一模塊，用于“將測溫模塊A劃分為m個網格單元其中，m為不小于3的正整數”；

第二模塊，用于“采集流經測溫模塊A上的電流、測溫模塊A兩端的電壓及每個網格單元的位置數據，計算獲得測溫模塊A的等效電阻R_e及每個網格單元的溫度初始值及導熱系數初始值”；

第三模塊，用于“結合測溫模塊A的密度值、比熱容值、位置數據、每個網格單元的溫度初始值及導熱系數初始值，基于第一類邊界條件，建立測溫模塊A的熱平衡方程”；

第四模塊，用于“對所述測溫模塊的熱平衡方程進行離散化，采用禁忌搜索算法進行計算得到并計算獲取所述測溫模塊的總電阻值比較所述等效電阻R_e與總電阻值的差值，并當值時，停止計算，進而得到所述測溫模塊上的溫度分布；其中，i指的是沿x軸方向上的第i個網格單元；j指的是沿y軸方向上的第j個網格單元；i×j＝m，且i、j均為正整數；n指的是計算步驟，n≥1”。

優(yōu)選地，該溫度場的溫度測量裝置還包括一個設于測溫模塊A任意一端的輔助測溫元件B，且輔助測溫元件B與溫度場轉換模塊E連接；或該溫度場的溫度測量裝置還包括兩個分別設于測溫模塊A兩端的輔助測溫元件B，且兩個輔助測溫元件B分別與溫度場轉換模塊E連接；其中，輔助測溫元件B可以為熱電偶或熱電阻或紅外測溫儀等測溫儀器。具體地，根據實際實施條件，輔助測溫元件B可不安裝或只安裝于測溫模塊A的一端或安裝于測溫模塊A的兩端，區(qū)別在于，溫度場轉換模塊E的轉換速度，即當兩端有精確的溫度初始值，利于后續(xù)網格單元的溫度計算值的收斂速度，提高了計算效率。

進一步地，該溫度場的溫度測量裝置還包括與溫度場轉換模塊E連接的環(huán)境溫度測量模塊F，環(huán)境溫度測量模塊F用于測量測溫模塊A附近的環(huán)境溫度，在本實施例中，可將測溫模塊A附近的環(huán)境溫度值作為溫度初始值，以加快溫度場轉換模塊E的轉換速度，有利于提高計算效率及計算精度。

值得說明的是，電流源模塊C和電壓測量模塊D可以由電壓源模塊、電流測量模塊或者電橋等電阻測量模塊取代，用于獲取測量模塊的等效電阻。

特別的，本申請以具有蜂窩體的蓄熱式燃燒器為例加以說明，在測量蜂窩體氣流流動方向的溫度分布時，若采用傳統的熱電偶較為復雜，且對蜂窩體有破壞，而采用蜂窩體作為蓄熱室時，其空間較小，無法布置較多的熱電偶。其中，20萬大卡燒嘴蓄熱室的體積為60cm*50cm*66cm，需要測量蜂窩體在66cm長度內的精確溫度分布。溫度范圍為30℃-1000℃。

具體測溫方法為：采用長度68cm直徑為0.2mm鉑銠合金絲(PtRh30)作為測溫模塊A，將鉑銠合金絲沿氣流方向穿入蜂窩體的孔內(其一般為直徑6mm的六邊形)；在蓄熱室的兩頭貼著蜂窩體分別安裝熱電偶測量，位置設在兩頭的平面上即可。在鉑銠合金絲兩頭沿著蜂窩體安裝電壓測量導線，按照圖7將安裝測量裝置。在溫度場轉換裝置——電腦上設置將66cm的鉑銠合金絲等分為100段，經過計算即可得到該66cm鉑銠合金絲上的100個溫度值，進而得到蓄熱室的溫度分布。

若以線性溫度場為例，且將測溫模塊A劃分為五段，除首尾兩段長度為l/2，其余各段長度均為l，具體如圖7所示。其中，在測溫模塊A的兩端均設有輔助測溫元件B。

初始化：

T1、T5為已知確定的溫度點，T2、T3、T4為未知溫度點，初始化方法為，由T1、T5根據溫度查出相應的電阻值R1、R5。

由測量值可知：

R1+R5+R2+R3+R4＝U測/I定＝R(U測為測溫元件兩端的壓降，I定給定的流經測溫元件的電流值)。初始化，使R2＝R3＝R4，由此得到R02、R03、R04的初始值。

此時，測溫模塊A處于熱平衡狀態(tài)，使用第一類邊界條件，假設控制單元的外表面溫度為當前位置的環(huán)境溫度，由于控制單元熱物性及形狀可以認為是薄材，則控制單元的溫度為其所在的環(huán)境溫度，建立測溫模塊A的熱平衡方程，具體為：

$\mathrm{\ρ}c\frac{\∂t}{\∂\mathrm{\τ}}=\frac{\∂}{\∂x}\left(\mathrm{\λ}\frac{\∂t}{\∂x}\right)$

由以上假設條件，并離散化可得：

$t_i^n=\frac{{\mathrm{\λ}}_{i-1}^{n-1}{t}_{i-1}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}{t}_{i+1}^{n-1}}{{\mathrm{\λ}}_{i-1}^{n-1}+{\mathrm{\λ}}_{i+1}^{n-1}}$

通過上式計算得到新一輪的測溫模塊A上的溫度分布，然后采用禁忌搜索算法進生成多組測溫模塊A的溫度分布解，即按照禁忌搜索表對每個溫度點進行△t(隨著運算的進行可以調小△t)的加減，生成多組測溫模塊A的溫度分布解；然后計算得到每組測溫模塊A總電阻值，并比較測量值R與總電阻值的差值，選取差值最小值一組數據作為解進入下一次迭代，其中，收斂條件為測量值R與總電阻值的差值小于A(其中A＝0.001，且A值可以根據計算誤差允許范圍取值)。

綜上所述，本發(fā)明提供了一種溫度場的溫度測量方法，利用測溫模塊的電阻測量信息，結合熱流平衡理論及禁忌搜索算法計算測溫模塊上的溫度分布，具體包括：步驟1、將測溫模塊劃分為m個網格單元；步驟2、采集流經測溫模塊上的電流、測溫模塊兩端的電壓及每個網格單元的位置數據，計算獲得測溫模塊的等效電阻R_e及每個網格單元的溫度初始值及導熱系數初始值；步驟3、結合測溫模塊的密度值、比熱容值、位置數據、每個網格單元的溫度初始值及導熱系數初始值，基于第一類邊界條件，建立測溫模塊的熱平衡方程；步驟4、對測溫模塊的熱平衡方程進行離散化，采用禁忌搜索算法進行計算得到并計算獲取測溫模塊的總電阻值比較等效電阻R_e與總電阻值的差值，并當值時，停止計算，進而得到測溫模塊上的溫度分布。本發(fā)明通過上述測量方法，能一次性地得到待測物體上的溫度分布，測量精度高且既可以獲取線性溫度場，也可獲取平面溫度場的溫度分布；同時，本申請?zhí)峁┑臏y量裝置占用待測物體的空間小，對待測溫度場的影響較小，對待測溫度物體的損傷較小。

最后應說明的是：以上實施例僅用以說明本發(fā)明的技術方案，而非對其限制；盡管參照前述實施例對本發(fā)明進行了詳細的說明，本領域的普通技術人員應當理解：其依然可以對前述各實施例所記載的技術方案進行修改，或者對其中部分技術特征進行等同替換；而這些修改或者替換，并不使相應技術方案的本質脫離本發(fā)明各實施例技術方案的精神和范圍。