誤差范圍區(qū)間與擴(kuò)展不確定度區(qū)間基于計(jì)算的比較（續(xù)1）

-
                                    誤差范圍區(qū)間與擴(kuò)展不確定度區(qū)間
                                              基于計(jì)算的比較（續(xù)1）
-
                                                                                                        史錦順
-
       本文論題：測量場合下，間接測量的測量結(jié)果區(qū)間的比較。前文《誤差范圍區(qū)間與擴(kuò)展不確定度區(qū)間基于計(jì)算的比較圖》是本文的一個(gè)特例。
       本文計(jì)算題目：一把精致銀米尺，長、寬、厚，形狀規(guī)整。測量該尺的物質(zhì)密度。
      
米尺的基本數(shù)據(jù)（參數(shù)分項(xiàng)測量，各10次）
       測量質(zhì)量M，用普通天平（MPEV 0.03%）。測得值：M_m = 8.900kg
       長度L的規(guī)格：長度的標(biāo)稱值L_o= 1m，MPEV 0.03%
       測量寬度D，用數(shù)顯卡尺（MPEV 0.05%），測得值：D_m =50.00mm。
       測量厚度H，用數(shù)顯千分尺（MPEV 0.03%）。測得值：H_m =20.000mm。

1《史法測量計(jì)量學(xué)》對(duì)誤差范圍的表達(dá)
1.1 測量方程與測得值函數(shù)
       物質(zhì)密度ρ的物理公式
                 ρ = 質(zhì)量/體積=質(zhì)量/ (長×寬×厚)
                    = M/(LDH)                                                                （1）
       物質(zhì)密度ρ的計(jì)值公式（腳標(biāo)o表示標(biāo)稱值；腳標(biāo)m表示測得值）      
                 ρ_m = M_m /(L_oD_mH_m)                                                    （2）
       測量方程為
                 ρ_m / ρ = (M_m/M) / [(L_o/L)( D_m/D)(H_m/H)]                          （3）
       測得值函數(shù)：
                 ρ_m = ρ(M_m/M) (L/L_o)(D/D_m)(H/H_m)                                    （4）
      
       注意，測得值函數(shù)中，自變量是質(zhì)量測得值M_m、長度L、寬度測得值D_m、厚度測得值H_m。
       寬度實(shí)際值D、厚度實(shí)際值H、質(zhì)量實(shí)際值M是物理公式中的值，是真值，是常量。長度標(biāo)稱值L_o是常量(標(biāo)稱值是定義值的一種形式[1]），而長度實(shí)際值L是真值，卻是變量。以往，經(jīng)典誤差理論中的誤差分析，直接從物理公式出發(fā)進(jìn)行微分，常常弄錯(cuò)有標(biāo)稱值的量的誤差元的符號(hào)。如頻率計(jì)誤差分析中的時(shí)標(biāo)誤差。由于誤差量的特點(diǎn)是其“絕對(duì)性”與“上限性”，經(jīng)典誤差理論的結(jié)果分析，沒錯(cuò)。如馬鳳鳴書中直接取絕對(duì)值，不出錯(cuò)。但如果要修正，或者要分析誤差因素的影響，符號(hào)問題就必須講究。誤差理論要精準(zhǔn)，要嚴(yán)格化，不能在符號(hào)上含糊。《史法》的測量方程與測得值函數(shù)，澄清了這一點(diǎn)。關(guān)鍵是辨清常量與變量。微分必須對(duì)變量微分。《史法》的區(qū)分常量與變量的程序，使誤差理論的微分操作，有了嚴(yán)格的數(shù)理邏輯基礎(chǔ)。

1.2 求誤差元
A 差分法
       誤差元是測得值與真值之差。測得值（2）與真值（1）之差就是誤差元。這樣做，物理意義十分清晰。筆者年輕時(shí)讀文獻(xiàn)，不理解“為什么求微分就是求誤差”，只是模仿。臨近退休，才弄明白誤差方程這套道理。取差分，符合誤差（元）的定義。微分不過是求差的數(shù)學(xué)方法。而取微分，必須明確變量。
       差分法求密度誤差的公式為
              Δρ_m =ρ_m - ρ
                = M_m /(L_oD_mH_m) - M/(LDH)                                           （5）   
       對(duì)（5）式做近似計(jì)算，即可求誤差元。（參見《史法測量計(jì)量學(xué)》[^1]第3章。）
B 微分法
       熟悉微積分，用微分法就更方便。
       對(duì)測得值函數(shù)ρ_m（4）
             ρ_m = ρ(M_m/M) (L/L_o)(D/D_m)(H/H_m)   
做全微分，得密度函數(shù)的誤差元      
            dρ_m =ρ (?ρ_m/?M_m)dM_m+ (?ρ_m/?L)dL+ (?ρ_m/?D_m)dD_m+ (?ρ_m/?H_m)dH_m
            dρ_m = ρ(M_m / M) (L / L_o)( D / D_m)(H / H_m)[ dM_m/M_m]  
                 + ρ(M_m / M) (L / L_o)( D / D_m)(H / H_m) [ dL/L]
                 - ρ(M_m / M) (L / L_o)( D / D_m)(H / H_m) [ dD_m/D_m]
                 - ρ(M_m / M) (L / L_o)( D / D_m)(H / H_m) [ dH_m/H_m]
            dρ_m / ρ_m = dM_m/M_m+ dL/L- dD_m/D_m – dH_m/H_m
       標(biāo)成相對(duì)誤差元的形式，有
             rρ_m =r_Mm + r_L - r_Dm - r_Hm                                               （6）   

1.2 誤差合成公式
       由于四項(xiàng)誤差（本例已知分項(xiàng)誤差范圍的指標(biāo)值）,只有一項(xiàng)偏大，其他大小相近，沒有“二、三項(xiàng)大誤差”項(xiàng)，《史法》[1] 誤差合成，取各項(xiàng)誤差范圍的“方和根”，密度測量的誤差范圍（合成誤差元絕對(duì)值99%概率意義上的最大可能值）是：
                 Rρ= √(R_M² + R_L² + R_D² + R_H²)                                       （7）
1.3 本例具體計(jì)算
     尺體質(zhì)量密度的測得值
             ρ_m = M/ L_oD_mH_m
                 =8.900kg/1000.0mm×50.00mm×20.000mm
                 = 8.900×10³ kg/m³
      密度測得值的相對(duì)誤差范圍（7）
             Rρ = √[(0.03%)² + (0.03%)² + (0.05%)² + (0.03%) ²]                 
                =0.072%
      密度測得值的絕對(duì)誤差范圍
             Δ ρ_m =8.9×10³×0.00072 kg/m³               
                 =6.4 kg/m³                                                                 （8）
       米尺質(zhì)量密度的測量結(jié)果：
              ρ=(8900±6) kg/m³                                                          （9）

（未完，待續(xù)）
-

——

（續(xù)前）


2 不確定度體系的計(jì)算
2.1 基本數(shù)據(jù)及分布認(rèn)定
       測量質(zhì)量M：普通天平（MPEV 0.03%，均勻分布）。測得值：M = 8.900kg，u_AM= 0.007%
       長度L的規(guī)格：長度的標(biāo)稱值L_o = 1m，MPEV 0.03%，正態(tài)分布。
       測量寬度D：數(shù)顯卡尺（MPEV 0.05%，均勻分布）。測得值：D_m =50.00mm，u_AD= 0.01%
       測量厚度H：數(shù)顯千分尺（MPEV 0.03%，均勻分布）。測得值：H_m =20.000mm，u_AH= 0.007%

2.2 擴(kuò)展系數(shù)的確定
       用單臺(tái)儀器的直接測量，直接用儀器的誤差范圍的指標(biāo)值MPEV，就夠了。沒有誤差合成問題，也就沒有必要先除以根號(hào)3，再乘以k=1.7。那是將MPEV復(fù)原。似乎道理講得通，但這是沒有實(shí)際意義的折騰。（誤差理論認(rèn)為：基礎(chǔ)測量是常量測量，儀器的隨機(jī)誤差包含在儀器指標(biāo)MPEV中，不另算；不確定度體系認(rèn)為測量對(duì)象可能是常量，也可能是隨機(jī)變量，不確定度評(píng)定要包括測量值的隨機(jī)變化部分。）
       不確定度體系是“方差合成”，一個(gè)關(guān)鍵問題是，多項(xiàng)誤差因素合成為u_C，此后擴(kuò)展因子怎樣取？（《史法》是“范圍合成”不存在這個(gè)問題。）
       不確定度體系的規(guī)范，多處說明，通常要求95%的概率，可取k=2；要求99%的概率，可取k=3。但這樣做的前提是合成后誤差是正態(tài)分布。又說各項(xiàng)誤差是獨(dú)立的，相互之間不發(fā)生作用。取99%概率時(shí)，何以將均勻分布的那一項(xiàng)除以根號(hào)3，卻乘以3，竟毫無道理的擴(kuò)大該項(xiàng)誤差的作用到1.73倍？這是把正態(tài)分布的隨機(jī)誤差的規(guī)律，錯(cuò)誤地用于統(tǒng)計(jì)中保持恒值的系統(tǒng)誤差的惡果。毫無道理，是錯(cuò)誤的。錯(cuò)誤的分布認(rèn)定，錯(cuò)誤的擴(kuò)展因子，只因?yàn)槭前选昂阒档南到y(tǒng)誤差”當(dāng)成“隨機(jī)誤差”來處理。錯(cuò)誤的根源是把測量的“時(shí)域統(tǒng)計(jì)”錯(cuò)誤地當(dāng)成是“臺(tái)域統(tǒng)計(jì)”[2]。
       有人說：乘2乘3是不符合規(guī)矩的錯(cuò)誤操作。其實(shí)不然。且看如下規(guī)范：
      《校準(zhǔn)領(lǐng)域測量不確定度評(píng)估指南》（CNAS-GL09:2008）p8：
       5.2 用實(shí)驗(yàn)的方式通常很難確認(rèn)輸出量的估計(jì)值是否屬于正態(tài)分布。但是，在通常情況下，當(dāng)幾個(gè)不確定度的分量來源于概率分布獨(dú)立的輸入量（即參量個(gè)數(shù)N≥3），如服從于正態(tài)分布或矩形分布時(shí)，則依據(jù)中心極限定理可以假設(shè)輸出量的分布接近于正態(tài)分布。
      （史注：該規(guī)范中的實(shí)例S2、S3、S4、S5、S6、S7都按“輸出量為正態(tài)分布”處理.）
       本題有誤差范圍四項(xiàng)，隨機(jī)誤差三項(xiàng)。常規(guī)的分布認(rèn)定：正態(tài)分布四項(xiàng)，矩形分布（均勻分布）三項(xiàng)，求U₉₉，取k=3，符合不確定度體系作法。至于認(rèn)為各項(xiàng)獨(dú)立，也只能如此。（不確定度體系的“分布論”、“相關(guān)系數(shù)論”本來就是多余的，并無科學(xué)根據(jù)。）
-
2.3 不確定度體系的計(jì)算
       尺體密度測得值
              ρ_m = M_m / L_oD_mH_m
                 =8.900kg/1000.0mm×50.00mm×20.000mm
                 = 8.900×10³kg/m³
       合成標(biāo)準(zhǔn)不確定度
              u_C = √ {[(0.03%)/√3]² + [(0.03%)/3]² + [(0.05%)/√3]²+ [(0.03%)/√3]2                 
                    +(0.007%)²+(0.01%)²+(0.007%)²}
                 = √ [(0.0173%)²+ (0.01%)² + (0.0289%)² + (0.0173%) ²                
                    +(0.007%)²+(0.01%)²+(0.007%)²}
                 =0.042%                                                            （10）
       包含概率99%的擴(kuò)展不確定度（相對(duì)值）
              U_99r = 3u_C
                  = 0.13%                                                           （11）
-
2.4 不確定度體系的測量結(jié)果表達(dá)
       測得值
               y = 8.900×10³ kg/m³
       包含概率99%的區(qū)間半寬：               
               U₉₉ = 8.900×10³ kg/m³×0.13%
                  =12 kg/m³                                                           （12）
        測量結(jié)果：
                Y = (8900±12) kg/m³                                                （13）

3 史評(píng)
       拙文《論不確定度體系的公式錯(cuò)誤》[2]中的區(qū)間比較圖以及前文《誤差范圍區(qū)間與擴(kuò)展不確定度區(qū)間基于計(jì)算的比較圖》，畫的是單臺(tái)儀器。有人認(rèn)為擴(kuò)展不確定度的擴(kuò)展系數(shù)的選取是不對(duì)的，因此，比較是不成立的。這樣說似乎有道理，但本質(zhì)上是掩護(hù)了不確定度體系的錯(cuò)誤。
       誤差合成，總是若干量的合成。因此《論不確定度體系的公式錯(cuò)誤》[2]的舉例，以及前文，應(yīng)看成是多個(gè)分量中的一個(gè)。乘3是常規(guī)。
       以（10）式為例，對(duì)整體乘3，本質(zhì)是對(duì)各分項(xiàng)都乘3了。
             u_C = √ {[(0.03%)/3]² +[(0.03%)/√3]² + [(0.05%)/√3]² + [(0.03%)/√3] ²                 
                   +(0.01%)²+(0.02%)²+(0.01%)²}
3 u_C = √ { [3× (0.03%)/3]² +[3× (0.03%)/√3]² + [3× (0.05%)/√3]²
      + [3×(0.03%)/√3] ² +[3× (0.01%)]²+[3× (0.02%)²+[3× (0.01%)]²]
                 =0.042%                                                       （14）
       其中的第2項(xiàng)（尺長）、第5項(xiàng)、第6項(xiàng)、第7項(xiàng)，因?yàn)閱雾?xiàng)是正態(tài)分布，合成說得過去；但第2項(xiàng)、第3項(xiàng)、第4項(xiàng)都無故擴(kuò)大1.73倍，是沒有道理的，是錯(cuò)誤的。因是“方和根”，大項(xiàng)誤差（第3項(xiàng)）較大，自身雖然只擴(kuò)大到1.73倍，因?yàn)槠椒胶屯癸@大誤差的作用，且多計(jì)了在MPEV中本已包含的隨機(jī)誤差項(xiàng)，竟使函數(shù)區(qū)間半寬擴(kuò)大到2倍！這算什么玩意兒？不確定度體系，偽科學(xué)！
-
       本人前幾次的圖解說明，本質(zhì)是數(shù)項(xiàng)誤差合成（N≥3）中單項(xiàng)作用的變化，是通常情況。這是不確定度體系的本質(zhì)性錯(cuò)誤之一。本人說明也不夠，易于引起誤解，故另寫兩篇續(xù)文在此。在否定不確定度體系這一點(diǎn)上，本人并無錯(cuò)誤，而表達(dá)要盡可能更嚴(yán)格，是必要的。故對(duì)csln先生表示感謝。
-
       不確定度體系的龐大架構(gòu)：A類不確定度、B類不確定度、合成不確定度、擴(kuò)展不確定度，這一套的核心是“誤差合成”。直接測量用MPEV就夠了；間接測量必須有誤差合成。就隨機(jī)誤差來說，不確定度體系與經(jīng)典誤差理論并無區(qū)別，也與統(tǒng)計(jì)理論一致。但系統(tǒng)誤差不行。測量計(jì)量的統(tǒng)計(jì)，是“時(shí)域統(tǒng)計(jì)”，而不是“臺(tái)域統(tǒng)計(jì)”。在單臺(tái)儀器的重復(fù)測量中，系統(tǒng)誤差是恒值。在間接測量中，每臺(tái)儀器各自測量不同的量，各量間，有物理公式所確定的函數(shù)關(guān)系，但系統(tǒng)誤差卻是各自獨(dú)立的。每種儀器的系統(tǒng)誤差對(duì)被測函數(shù)量的作用，是獨(dú)立的。在多次重復(fù)測量中，各項(xiàng)系統(tǒng)誤差的量值、各系統(tǒng)誤差間的關(guān)系，由量值的函數(shù)關(guān)系確定，是固定的；即函數(shù)的總系統(tǒng)誤差是恒值。函數(shù)的系統(tǒng)誤差既不是均勻分布，通常也不可能是標(biāo)準(zhǔn)正態(tài)分布（無偏倚值），而一般是有偏正態(tài)分布（其偏倚值是系統(tǒng)誤差）。因此，對(duì)系統(tǒng)誤差的合成量乘以不等于1的因子k是錯(cuò)誤的。下文將詳述。
       直接進(jìn)行“范圍合成”，簡單又正確，何不為之？

       不確定體系的測量結(jié)果區(qū)間半寬U₉₉竟是誤差理論測量結(jié)果區(qū)間R的兩倍（本例計(jì)算），孰是孰非，不值得認(rèn)真研究一番嗎？
       不確定度體系錯(cuò)了，且錯(cuò)誤嚴(yán)重！要它何用？
-
資料及說明
[1]史錦順著《史法測量計(jì)量學(xué)》，等待出版。質(zhì)檢出版社在等待國家計(jì)量院的鑒別結(jié)論。其內(nèi)容本人已多次在本欄目征求意見。許多網(wǎng)友是知道的。
[2]史錦順《論不確定度體系的公式錯(cuò)誤》，本欄目中有征求意見稿。宣傳新學(xué)說與抨擊不確定度體系的文章，本人已在本欄目貼出有署名的文章513篇（已編為10集）。有意向研究的網(wǎng)友，告訴我你的郵箱號(hào)，我可將全部文稿和書稿都寄給你。我已八十一歲半，傳播自己的學(xué)術(shù)思想，既是義務(wù)，也是一種快樂。
-

太有深度了，我得慢慢消化。

關(guān)于所謂的"系統(tǒng)(測量)誤差"，拎清后的"(測量)誤差理論"(如1990年代后費(fèi)業(yè)泰等先生的著述)已分為"已(確)定"與"未(確)定"兩部分----對(duì)于"已(確)定系統(tǒng)誤差"，"間接測量"的"合成"就是"具體值代入函數(shù)運(yùn)算"，不涉及什么"范圍合成"，當(dāng)然也無關(guān)"概率分布形式"、"相關(guān)性"之類問題;  但對(duì)于"未(確)定系統(tǒng)誤差" ，人們只能"估計(jì)"它的"可能取值范圍"及相應(yīng)的"概率分布形式"，"間接測量"時(shí)便涉及"可能取值范圍"的"合成"及"結(jié)果"的"概率分布形式"等問題---這些問題在"數(shù)學(xué)"(概率與統(tǒng)計(jì))上不存在未明的概念問題，"線性化"理論"合成"、"蒙特卡洛"仿真"合成"都是行之有效的方法，實(shí)際應(yīng)用中，各"分量"自身的"概率分布形式"、"分量"之間的"相關(guān)性"是不可回避的難題，需要可靠的"經(jīng)驗(yàn)"支持！

uC = √ {[(0.03%)/√3]2 + [(0.03%)/3]2 + [(0.05%)/√3]2+ [(0.03%)/√3]2                 
                    +(0.007%)2+(0.01%)2+(0.007%)2}
                 = √ [(0.0173%)2+ (0.01%)2 + (0.0289%)2 + (0.0173%) 2                 
                    +(0.007%)2+(0.01%)2+(0.007%)2}
                 =0.042%                                                            （10）
       包含概率99%的擴(kuò)展不確定度（相對(duì)值）
              U99r = 3uC
                  = 0.13%                                                           （11）

0.0289/0.042=0.688

有一均勻分布的顯著項(xiàng)占合成標(biāo)準(zhǔn)不確定度三分之二以上，合成不確定度仍然是均勻分布，U99包含因子是1.72或1.71

　　“誤差區(qū)間范圍”與“擴(kuò)展不確定度”是性質(zhì)完全不同的兩個(gè)概念，無法擺在同一個(gè)平臺(tái)上相比較。
　　“誤差區(qū)間范圍”是允許的誤差變化范圍的寬度，量化反映了測量過程或測量結(jié)果準(zhǔn)確性要求，正式術(shù)語是最大允許誤差絕對(duì)值，寫成英文縮寫就是MPEV。
　　“測量不確定度”量化反映了測量過程或測量結(jié)果的可信性，不是誤差，也不是誤差范圍。擴(kuò)展不確定度只能用英文單詞的第一個(gè)字母U表示，當(dāng)用小寫字母u時(shí)，代表標(biāo)準(zhǔn)不確定度。因此，擴(kuò)展不確定度不能縮寫為MPEV。
　　另外，“不確定度”本身就是某個(gè)區(qū)間（人們估計(jì)得到的真值存在的可能區(qū)間）的“半寬度”，“區(qū)間半寬的區(qū)間”非常繞口，也令人費(fèi)解，因此不存在“擴(kuò)展不確定度區(qū)間”這個(gè)概念。

喜歡這種學(xué)術(shù)性帖子，收藏下來好好消化下

史錦順 發(fā)表于 2018-11-18 10:36
（續(xù)前）

380310291@qq.com

1507933138@qq.com

夢回棗鄉(xiāng) 發(fā)表于 2021-4-16 14:11

       2018年我的一句承諾，2021年了，我已不再兌現(xiàn)。我已84周歲。年老多病，請(qǐng)諒解。本欄目都有，自己查查吧。

       按新版著作權(quán)保護(hù)法，網(wǎng)上文章、書稿，都在保護(hù)之列。可以下載自用，不可復(fù)制公開發(fā)表，否則即為侵權(quán)。作為不同意見的學(xué)術(shù)討論，可以部分引用，但不可全文另行刊出。

感謝史老師，健康長壽