求教，非線性測量模型如何進(jìn)行不確定度合成？

測量模型是這樣的Y=(fm-fm0)f1/Mf0,其中只有fm0是個常量。從沒有進(jìn)行這樣的合成，有沒有大神知道？

用模特卡羅方法傳播，此法萬能

hulihutu 發(fā)表于 2016-6-16 19:59
用模特卡羅方法傳播，此法萬能

蒙特卡洛太復(fù)雜了。。。主要是時間太緊張，能不能展開后簡化呢

把(fm-fm0)設(shè)為F，測量模型變?yōu)閅=Ff1/Mf0，這是典型的“黑箱模型”，用相對不確定度評定是非常簡便的。

規(guī)矩灣錦苑 發(fā)表于 2016-6-16 22:12
把(fm-fm0)設(shè)為F，測量模型變?yōu)閅=Ff1/Mf0，這是典型的“黑箱模型”，用相對不確定度評定是非常簡便的。 ...

1、黑箱模型不是無法建立確切表達(dá)式的時候用嗎？
2、怎么用相對不確定度評定呢？
3、如果F和M有相關(guān)性又該怎么辦呢？

啊，我問題太多了，謝謝大神

pennym 發(fā)表于 2016-6-17 09:56
1、黑箱模型不是無法建立確切表達(dá)式的時候用嗎？
2、怎么用相對不確定度評定呢？
3、如果F和M有相關(guān)性又 ...

　　沒人說黑箱模型是無法建立確切表達(dá)式的時候用啊，呵呵。測量模型分為透明箱和黑箱兩種類型，所謂黑箱模型其實就是數(shù)學(xué)中所說的“單項式”，只有乘除沒有加減，黑箱模型和透明箱模型一樣是常用測量模型。黑箱模型的優(yōu)點是不必計算每個輸入量的靈敏系數(shù)，每個輸入量引入的不確定度分量均按相對不確定度處理，省去了微分之苦，化學(xué)分析中的測量模型絕大多數(shù)是黑箱模型。
　　某個輸入量引入的相對不確定度就是其絕對不確定度除以該輸入量全值，計算非常簡單。例如評估得到壓力表某個受檢點的不確定度為ui，壓力表量程為20MPa，則兩者相除就是該壓力表引入的相對不確定度。
　　如果F和M有相關(guān)性，理論上要計算相關(guān)系數(shù)，分量合成時要考慮相關(guān)項。不確定度不是數(shù)學(xué)計算而是估計，實際工作中沒有必要那么“較真”，相關(guān)性不強就仍然按不相關(guān)處理，如果是強相關(guān)，就根據(jù)正相關(guān)還是負(fù)相關(guān)的判斷，將兩者的不確定度相加減。

規(guī)矩灣錦苑 發(fā)表于 2016-6-17 10:25
　　沒人說黑箱模型是無法建立確切表達(dá)式的時候用啊，呵呵。測量模型分為透明箱和黑箱兩種類型，所謂黑箱 ...

我還有幾個問題啊，

1、如果計算出相對合成不確定度，但Y=Ff1/Mf0沒有最佳估計值（因為Y是標(biāo)準(zhǔn)器測量通用計數(shù)器的相對頻偏，這個值每個通用計數(shù)器都不同），那么絕對不確定度是不是求不出了呢？
2、用相對不確定度，能夠進(jìn)行穩(wěn)定性、重復(fù)性考核嗎？因為Y的最佳估計值測不出
3、F和M都是同一個儀器輸出的量，這樣算相關(guān)嗎？
4、F和M都是同一個儀器輸出的量，都有這個儀器的測量不確定度，那么在不確定度合成時，算不算重復(fù)計算了呢？

pennym 發(fā)表于 2016-6-22 09:39
我還有幾個問題啊，

1、如果計算出相對合成不確定度，但Y=Ff1/Mf0沒有最佳估計值（因為Y是標(biāo)準(zhǔn)器測量通 ...

　　１從“因為Y是標(biāo)準(zhǔn)器測量通用計數(shù)器的相對頻偏，這個值每個通用計數(shù)器都不同”這句話可以看出，你是在校準(zhǔn)若干個被校對像，用標(biāo)準(zhǔn)器校準(zhǔn)若干個通用計數(shù)器的相對頻偏。對每一個通用計數(shù)器進(jìn)行校準(zhǔn)才是日常的校準(zhǔn)活動，不應(yīng)該把對另一個通用計數(shù)器的校準(zhǔn)與對指定某通用計數(shù)器的校準(zhǔn)相摻和。
　　２相對不確定度和絕對不確定度只是個表達(dá)形式，用絕對不確定度能干的事，用相對不確定度一樣能干。
　　３F和M都是同一個儀器輸出的量，但卻不是同一個量，類似于用萬用表分別測量電壓和電阻，它們應(yīng)屬于不相關(guān)或弱相關(guān)。
　　４F和M雖然都是同一個儀器輸出的量，都有這個儀器的測量不確定度，但它們是不同的兩個“輸入量”，那么在不確定度合成時就不屬于同一個輸入量的重復(fù)，而是兩個不同的輸入量引入的不同不確定度分量，不能取大舍小，而應(yīng)該加以合成。

1、非線性測量模型的不確定度建議使用蒙特卡羅法（MCM），用傳統(tǒng)不確定度法（GUM)有一定的局限性。
2、蒙特卡羅法（MCM）評定測量不確定度要先學(xué)習(xí)數(shù)學(xué)軟件MATLAB、mathematica、mathecad等，學(xué)好其中一數(shù)學(xué)軟件后，你會發(fā)現(xiàn)MCM法評定不確定度不過如此。當(dāng)然，excel2010也能進(jìn)行MCM法評定不確定度，不過操作比較繁瑣。
3、非線性測量模型中輸入量相關(guān)時，可以喬里斯基分解相關(guān)系數(shù)的正定對稱矩陣方法來模擬相關(guān)性隨機(jī)變量。

規(guī)矩灣錦苑 發(fā)表于 2016-6-22 10:46
　　１從“因為Y是標(biāo)準(zhǔn)器測量通用計數(shù)器的相對頻偏，這個值每個通用計數(shù)器都不同”這句話可以看出，你是 ...

謝謝回答，但是我還有些疑問，比如：
1、如果我求出了Y的相對不確定度，但我還想求出絕對不確定度，是否應(yīng)該用相對不確定度乘以Y的最佳估計值？但是Y是通用計數(shù)器的相對頻偏，也就是說每一個通用計數(shù)器都不一樣，如何求出最佳估計值呢？
2、測量模型：Y=(fm-fm0)f1/Mf0，其中f1是測量重復(fù)性帶來的修正因子，請問這樣建立測量模型對嗎？還是說應(yīng)該這樣寫：Y=（(fm-fm0)/Mf0）+f1,f1為測量重復(fù)性帶來的修正項？

pennym 發(fā)表于 2016-6-22 15:29
謝謝回答，但是我還有些疑問，比如：
1、如果我求出了Y的相對不確定度，但我還想求出絕對不確定度，是否 ...

-
       看你的帖子，幾次都看不懂。你在頻標(biāo)比對器的頻差表示式中，加上一個頻率值f1,并不說明含義，誰能知道這是怎么回事？這一貼你才說是“測量重復(fù)性帶來的修正因子”。這種寫法，沒見過。重復(fù)性能修正嗎？不確定度A類評定，評定的就是“重復(fù)性”，哪來的修正因子？
       相對頻差公式Y(jié)=（(fm-fmo)/Mfo）僅僅表示測得的頻差的一個數(shù)。公式中不僅fmo是常數(shù)，倍增次數(shù)M也是常數(shù)，fo是頻率計內(nèi)標(biāo)頻率的標(biāo)稱值，也是常數(shù)。
       例如，用PO7之第四檔測量頻率計的晶振。M=2000，fo=5MHz，頻率計示值差1Hz表示1乘10的-10次方。
       測量20次，得20個Yi，代入貝塞爾公式，計算得到σ，除以根號20，即得不確定度理論的A類不確定度。
       20個頻差的平均值，就是“最佳估計值”。
       B類評定要先寫出誤差函數(shù)公式，再分項評定。參見葉德培先生的樣板評定。載于質(zhì)檢總局計量司編《測量不確定度評定指南》p80。
-
       由誤差元可以計算誤差范圍。
       不確定度沒有自己的元素，只好借助于誤差理論的函數(shù)關(guān)系。
       不確定度的模型，實質(zhì)就是誤差范圍對誤差元的函數(shù)關(guān)系。
       你列出的公式，是示值的“指數(shù)”，相當(dāng)于小數(shù)點后多少個零的小數(shù)，不是誤差量的函數(shù)關(guān)系，不能當(dāng)不確定度評定的“模型”，當(dāng)然不能由此公式來評定不確定度。至于頻標(biāo)比對器性能評定，用同頻輸入，看示值尾數(shù)的跳動數(shù)，計算阿侖偏差。
-

pennym 發(fā)表于 2016-6-22 15:29
謝謝回答，但是我還有些疑問，比如：
1、如果我求出了Y的相對不確定度，但我還想求出絕對不確定度，是否 ...

　　1、如果求出了Y的相對不確定度，還想求出絕對不確定度，可以用相對不確定度乘以Y的最佳估計值，也可以乘以名義值。因為最佳估計值和名義值相差并不大，而相對不確定度很小，相乘后得到的絕對不確定度相差會更小，不會影響到顛覆絕對不確定度的地步。也就是說雖然每一個通用計數(shù)器都不一樣，如果名義值是同一個，用各自的最佳估計值乘以相對不確定度得到絕對不確定度不會有區(qū)別。
　　2、測量模型：Y=(fm-fm0)f1/Mf0，其中f1是測量重復(fù)性帶來的修正因子，請問這樣建立測量模型對嗎？還是說應(yīng)該這樣寫：Y=（(fm-fm0)/Mf0）+f1，f1為測量重復(fù)性帶來的修正項？
　　我不是搞時間頻率計量的，時間頻率計量中頻差測量的測量模型，要根據(jù)相關(guān)檢定規(guī)程、校準(zhǔn)規(guī)范或檢測標(biāo)準(zhǔn)規(guī)定的測量方法來寫，一般規(guī)程、規(guī)范或標(biāo)準(zhǔn)中都會給出測量模型的，可以直接引用。史老先生是時間頻率計量的專家，測量模型應(yīng)該是他的拿手好戲。他說“M=2000，fo=5MHz”是常數(shù)，常數(shù)不是輸入量，應(yīng)該將M=2000，fo=5MHz代人公式，測量模型應(yīng)進(jìn)一步推導(dǎo)至：Y=1×10^-4(fm-fmo)f1，若果真如此，測量模型就簡單了，不確定度評定就更容易。

        原式中的M=2000是比對器的倍增“倍數(shù)”，而內(nèi)標(biāo)的標(biāo)稱頻率5MH_z中的“M”是單位詞頭，M=10⁶，故有：

                Y=1×10^-10(f_m-f_mo)                              （1）

        手頭有份不確定度理論著名專家葉德培先生的關(guān)于時頻計量領(lǐng)域評定不確定度的文章，上傳，供參考。我不信不確定度理論，認(rèn)為不確定度的理論、不確定度的評定都是錯誤的。為此，在本欄目已發(fā)抨擊文章三百多篇。先生可以比較、思考。

時間頻率計量中不確定度的評定.pdf (223.7 KB, 下載次數(shù): 8)

2016-6-23 07:57 上傳

點擊文件名下載附件
下載積分: 金幣 -1

規(guī)矩灣錦苑 發(fā)表于 2016-6-23 00:21
　　1、如果求出了Y的相對不確定度，還想求出絕對不確定度，可以用相對不確定度乘以Y的最佳估計值，也可 ...

謝謝您，我比較小白，還想問一下，如果寫成Y=(fm-fm0)/Mf0 +f1 這樣的形式，是不是除了蒙特卡洛就沒有更簡便的方法了嗎？

史錦順 發(fā)表于 2016-6-22 18:58
-
       看你的帖子，幾次都看不懂。你在頻標(biāo)比對器的頻差表示式中，加上一個頻率值f1,并不說明含義，誰 ...

謝謝史先生，這樣的測量模型表達(dá)式，主要來源于倪育才先生所著的《實用不確定度評定》一書，書中寫測量模型通式為量=計算公式+隨機(jī)誤差帶來的修正項+系統(tǒng)誤差帶來的修正項。因此我覺得銣原子頻標(biāo)的測量模型為：Y=(fm-fm0)/Mf0+f1+...+fn.
您剛才說M和fo都是常數(shù)，但是M是由頻標(biāo)比對器中的頻差倍增器輸出的，fo是由頻率標(biāo)準(zhǔn)輸出的，是不是也會有誤差，所以也會帶來不確定度呢？

pennym 發(fā)表于 2016-6-23 14:15
謝謝史先生，這樣的測量模型表達(dá)式，主要來源于倪育才先生所著的《實用不確定度評定》一書，書中寫測量模 ...

-
（一）為什么倍增倍數(shù)M與晶振頻率f_o都視為常數(shù)
       這是“微差法”測量技術(shù)帶來的巨大好處。
       頻標(biāo)比對器之計數(shù)器，顯示F變1Hz,代表的是被測晶振頻率f變化（1/Mf_o）Hz.比對器的公式為
             Δf = ΔF/(Mf_o)
       f_o的最大變化不會超過10^-5,引入Δf的誤差不會超過10^-5. Δf是誤差，誤差的相對誤差1/10即可略，10^-5，當(dāng)然可略。M是倍乘的理論值，我做實驗證明至少在10^-6的量級上是準(zhǔn)確的常數(shù)。好，我們把它大大低估，就算千分之一，它引入的Δf的相對誤差，不過是千分之一。Δf是誤差，誤差的相對誤差1/10即可略，千分之一，當(dāng)然可略。以上是對同一路的分析。對兩路來說，可以用“換邊”的辦法，來證明M的常數(shù)性質(zhì)。我做實驗，在10^-14的水平上，證明M是常值。
       由上，把M、f_o視為常數(shù)是沒有問題的。
-
（二）關(guān)于誤差模型
       不確定度計量評定的基本公式（又稱數(shù)學(xué)模型）是
                EX= X―B                                                                                        （1）
                EX(0)+ ΔEX = X(0) + ΔX(分辨)+ ΔX(重復(fù))+ ΔX(其他)―[B(0) +ΔB(標(biāo))]
                ΔEX =ΔX(分辨)+ ΔX(重復(fù))+ ΔX(其他) ―ΔB(標(biāo))                                   （2）
       X是被測量，B是標(biāo)準(zhǔn)量，EX是差值，加（0）表示無誤差時的量。
       ΔEX 是要評定的不確定度（元），ΔX(分辨)表示被檢儀器分辨力因素，ΔX(重復(fù))表示“用測量儀器測量計量標(biāo)準(zhǔn)”時讀數(shù)的重復(fù)性，ΔX(其他)是被檢儀器其他因素的影響；ΔB(標(biāo))是標(biāo)準(zhǔn)的誤差。
       依據(jù)（1）（2）式進(jìn)行不確定度評定，是當(dāng)前計量不確定度評定的常規(guī)。中國的評定如此，歐洲的評定（倪玉才大部分抄錄）也是如此。其本質(zhì)就是GUM的泰勒展開法。
-
       第一步要寫出測得值函數(shù)。第二步進(jìn)行微分，求誤差元。第三步合成為誤差范圍。誤差范圍就是不確定度評定的擴(kuò)展不確定度。
-
       例：數(shù)字式頻率計的分析過程（以下引自《史氏測量計量學(xué)說》第三章，本欄目有）
-
       計數(shù)式頻率計是最基本最常用的測頻儀器。現(xiàn)行教科書上給出的計數(shù)式頻率計的公式為：
                      f = N/T                                                                                       （3.5）
       式中N為計數(shù)值，T為閘門時間。由于沒有區(qū)分測得值和實際值，用以分析，常常出錯。此式明顯標(biāo)示，頻率與閘門時間成反比。由此，若內(nèi)標(biāo)頻率偏低，則閘門時間長，則頻率值小；其實，恰恰相反：內(nèi)標(biāo)頻率偏低，必有閘門時間值偏大，必定頻率測得值偏大。
       式(3.5)是物理公式，不便直接用于分析測量問題；以往硬這樣做，難免出錯。有些作者看到了這一點，用取絕對值的辦法來避免正負(fù)號的矛盾，這不能算錯，但繞開矛盾，實際上也掩蓋了矛盾。
       要做幾種區(qū)分：區(qū)分頻率的測得值與實際值；區(qū)分閘門時間的標(biāo)稱值與實際值；區(qū)分N的顯示值與實際值。
       計數(shù)式頻率計的計值公式為：
                    f_m = N_r/T_n                                                                                        （3.6）
       式中f_m是測得值（被測頻率的實際值是f），Tn是閘門時間的標(biāo)稱值(閘門實際時間是T)，N_r是計數(shù)器的指示值（N是理論值，等于1/fT），T_n 是閘門時間的標(biāo)稱值，通常為1秒，或1秒的10^±n倍。
       分析測得值，就是分析測得值同實際值（真值）的差別，就是將測得值同實際值相比較。比較的方法之一是二者相除。實際值做除數(shù)，即做標(biāo)準(zhǔn)。
       計值公式(3.6)除以物理公式(3.5)，得測量方程：
                  f_m / f =N_rT/(NT_n)                                                                           （3.7）
       由測量方程，知測得值函數(shù)：
                  f_m=[N_rT/(NT_n)] f                                                                             （3.8）
       注意，我們研究的是測量問題（可設(shè)想是在用幾臺儀器同時測量同一物理量），被測頻率的客觀值f是常量，測得值fm是變量。閘門時間標(biāo)稱值Tn是常量，閘門時間實際值T是變量。理論值N是常量；讀數(shù)Nr是變量。
       （3.7）式是測量方程，（3.8）式是測得值函數(shù)。微分法分析誤差，就是求測得值函數(shù)在常量點上的全微分。
       A 求微分
                  df_m = [N_rf /(NT_n)]dT + [Tf /(NT_n)]dN_r                                               （3.9）
       B 誤差元：變量相對于常量的偏差量
                  Δf_m = [N_rf /(NT_n)] ΔT + [Tf /(NT_n)] ΔN_r                                            （3.10）
       C 相對差
      （3.10）式除以（3.8）式
                   δf_m = ΔT / T + ΔN_r / N_r                                                                    （3.11）
       因閘門時間由內(nèi)標(biāo)（頻率為f_b）分頻而來，有
                   T=K(1/f_b)
                   ΔT/T = - Δf_b/f_b                                                                                   (3.12)
       將（3.12）式代入（3.11）式，得
                   δf_m = - Δf_b/f_b +ΔN_r / N_r  
                   δf_m = - δf_b + δN_r                                                                               (3.13)
      （3.13）式表明，測得值與頻率計內(nèi)標(biāo)頻率成反比，即與時基成正比，是正確的，這糾正了只按物理公式求微分的不當(dāng)認(rèn)識。
       δfb是頻率計內(nèi)晶振引入的誤差項。其中包括：老化率、溫度效應(yīng)、晶振穩(wěn)定度等。δNr包括分辨力，計數(shù)器不穩(wěn)等引入的誤差項。本節(jié)講誤差分析的基本方法，只講主干部分，下續(xù)分析略。
-





補充內(nèi)容 (2016-6-24 06:33):
第18行之“倪玉才”應(yīng)為倪育才。

pennym 發(fā)表于 2016-6-23 14:07
謝謝您，我比較小白，還想問一下，如果寫成Y=(fm-fm0)/Mf0 +f1 這樣的形式，是不是除了蒙特卡洛就沒有更 ...

　　并非Y=(fm-fm0)/Mf0 +f1 這樣形式的測量模型，除了蒙特卡洛就沒有更簡便的方法了，GUM法仍然可以評估，特別是在化學(xué)分析領(lǐng)域，類似于這樣的測量模型比比皆是，也都是用GUM法評定的，你也可以參考相關(guān)的不確定度評定案例用GUM法評定。

規(guī)矩灣錦苑 發(fā)表于 2016-6-24 01:15
　　并非Y=(fm-fm0)/Mf0 +f1 這樣形式的測量模型，除了蒙特卡洛就沒有更簡便的方法了，GUM法仍然可以評估 ...

好的，謝謝您

史錦順 發(fā)表于 2016-6-23 20:08
-
（一）為什么倍增倍數(shù)M與晶振頻率f都視為常數(shù)
       這是“微差法”測量技術(shù)帶來的巨大好處。

謝謝耐心回復(fù)。

規(guī)矩灣錦苑 發(fā)表于 2016-6-24 01:15
　　并非Y=(fm-fm0)/Mf0 +f1 這樣形式的測量模型，除了蒙特卡洛就沒有更簡便的方法了，GUM法仍然可以評估 ...

請問有沒有這方面的書推薦呢，我想研究一下

史錦順 發(fā)表于 2016-6-23 20:08
-
（一）為什么倍增倍數(shù)M與晶振頻率f都視為常數(shù)
       這是“微差法”測量技術(shù)帶來的巨大好處。

史先生，我還有一個問題想要請教，可能問題比較小白，但我也一直沒搞清楚。就是在檢定通用計數(shù)器或頻率計時，檢定系統(tǒng)軟件都會提示選擇倍增次數(shù)M，我這里只有兩個選項×100，或×10000，但是研究頻標(biāo)比對器原理的時候，得到fm=fm0+Mfo,也就是前面公式的來由，此公式中的M等于10的n次方乘以N,而N我理解是一開始fx進(jìn)入倍增器的倍增次數(shù)，是隨著輸入頻率也就是fx變化的（因為fx可能為1MHz也可能是5MHz、10MHz），那么M也必然隨著fx大小不同。


  我的問題是：是否此M非彼M呢？也就是說，公式里的M=10的n次方乘以N,隨fx的大小不同而不同，而另一個M=100或者10000？

pennym 發(fā)表于 2016-6-24 11:12
史先生，我還有一個問題想要請教，可能問題比較小白，但我也一直沒搞清楚。就是在檢定通用計數(shù)器或頻率計 ...

       關(guān)于頻標(biāo)比對器的公式，我曾做詳細(xì)推導(dǎo)。現(xiàn)上傳，供參考。弄懂基本原理，處理具體問題就容易了。儀器型號不同，各種應(yīng)用場所，符號意義也不一樣，因此問人家，也不好回答。還是要自己弄明白各種標(biāo)度量之間的關(guān)系。
       請你注意，不確定度評定的模型，必須是儀器示值同誤差因素間的函數(shù)關(guān)系。(Fi-Fo)/Mfo不是誤差函數(shù)關(guān)系，不能用來當(dāng)評定頻率計、晶振、原子頻標(biāo)的模型。
       蒙特卡洛法，對測量計量來說，是“蒙人卡脖”法。它是由輸入量的分布求輸出量的分布。誰知道頻率計、晶振、原子鐘的誤差的分布是什么？假設(shè)一通，算一通，算出來又有什么用？

比對器Microsoft Word 文檔.doc (80.5 KB, 下載次數(shù): 1)

2016-6-24 12:50 上傳

點擊文件名下載附件
下載積分: 金幣 -1

pennym 發(fā)表于 2016-6-24 10:57
請問有沒有這方面的書推薦呢，我想研究一下

可以讀一下化學(xué)成分分析方面的不確定度評定案例。