計量論壇
標題: 交叉系數決定合成法 [打印本頁]
作者: 史錦順 時間: 2016-5-27 09:30
標題: 交叉系數決定合成法
本帖最后由 史錦順 于 2016-5-27 09:41 編輯
因原符號本網頁不接受,換成另外的表示法,則面目全非。只好先刪掉,另想辦法,再發。
作者: 秦時明月 時間: 2016-5-27 14:14
您可以用截圖、拍照方式把文檔制作成圖片、PDF等格式上傳啊。
作者: 史錦順 時間: 2016-5-27 15:21
本帖最后由 史錦順 于 2016-5-27 15:36 編輯
-
交叉系數決定合成法(1)
-
史錦順
-
引言
誤差,表示測得值與實際值的差距。誤差的概念,有三層意思:誤差元、誤差范圍,或泛指二者。
誤差元定義為測得值減真值。恒值的誤差元,稱為系統誤差;隨機變化的誤差元,稱為隨機誤差。
誤差范圍定義為誤差元的絕對值的一定概率(大于99%)意義上的最大可能值。
測得值與誤差范圍構成測量結果。
誤差范圍又稱為準確度,是測量儀器、計量標準以及測量結果水平的表征量。
誤差合成是由誤差元求誤差范圍。
-
1 誤差合成的原則、途徑與方法
誤差量的特點是其絕對性與上限性。誤差合成的原則是保險性與合理性。保險第一,合理第二;在保險的基礎上追求合理。
保險的含義是確定的誤差范圍值要包括誤差元的最大可能值。合理的含義是確定的誤差范圍值要盡可能接近實際值,就是要利用誤差量之間存在的抵消性。。
-
誤差量要絕對值化,方式有兩種。
第一種方式是直接對誤差元取絕對值。經典誤差理論對系統誤差直接取絕對值,合成取絕對和,保險,但偏于保守。而隨機誤差可正可負,有相互抵消作用,直接取絕對值不能體現隨機誤差的特點。第一種方式不能貫通。
第二種方式是取“方根”。初等數學規定:開平方根取正值。本文提出用“方根法”,可以貫通于隨機誤差與系統誤差。注意保險性與合理性,得出各種使用條件下的誤差合成公式。取“方根”,按交叉系數近于1還是近于零來確定公式,可推導出“絕對和”與“方和根”兩種方法。交叉系數的取值,可以體現誤差量間有無抵消性。
-
誤差合成的途徑也有兩種。
第一種途徑是“方差合成”,其基本條件是隨機性。 不確定度理論合成的途徑是方差合成,其方針是統一采用“方和根法”。對隨機誤差的處理與經典誤差理論相同,沒有問題;但對系統誤差的處理,出現嚴重問題。為實行“方和根法”,產生五項難題:1)認知誤差量的分布規律、2)化系統誤差為隨機誤差、3)假設不相關、4)范圍與方差間的往返折算、5)計算自由度。其中1)很難;2)不可能;3)對系統誤差錯誤;4)與5)都以 1)為基礎,也很難。
第二種途徑是“范圍合成”。本文著眼于范圍,貫通了兩類誤差合成的各種情況。要點是統籌隨機誤差與系統誤差的處理,把隨機誤差元變成是誤差范圍的直接構成單元。為此,用或正或負的恒值β代表系統誤差元;用三倍的隨機誤差元3ξi代表隨機誤差對誤差范圍的貢獻單元。這樣,系統誤差β與隨機誤差元3ξ對誤差范圍的貢獻權重相同。于是,公式推導與合成處理,都簡潔方便。
-
誤差合成新理論的要點與特點如下:
1)體現誤差量的兩大特點:絕對性和上限性。
2)通過取方根,實現誤差量的絕對值化;可以貫通于隨機誤差和各種系統誤差。
3)著眼于“范圍”。進行各誤差元到誤差范圍的合成;進行分項誤差到總誤差范圍的合成。
4)由交叉系數決定合成法的選取。避開有歧義的相關系數概念。
5)合成中,只需辨別誤差的性質(隨機誤差還是系統誤差),大系統誤差還是小系統誤差。不需辨別相關性。與分布無關。
6)依誤差性質、項數的不同,把交叉系數典型化為0或1,由此得到誤差合成的具體方法。
-
誤差合成方法口訣:兩三項大系統誤差,絕對值相加;再與其他項合成,一律方和根。
-
作者: 史錦順 時間: 2016-5-27 15:51
本帖最后由 史錦順 于 2016-5-27 16:00 編輯
-
交叉系數決定合成法(2)
-
史錦順
-
2 隨機誤差元構成的誤差范圍
對隨機誤差序列的處理,誤差理論兩百年前已有“方均根法”,成熟而完美。
測量實踐中,人們易于認識隨機誤差。對常量的重復測量中,測得值的隨機變化量就是隨機誤差。
隨機誤差元可大可小,可正可負。有四個特性:單峰性、對稱性、抵消性、有界性。
按統計理論,隨機誤差是正態分布(在測量次數N不遠大于10時,有t分布成分)。以3σ為半寬的分布區間,包含概率大于99%。
-
對隨機誤差,有如下定義與關系:
1)隨機誤差元等于測得值減測得值的期望值(當無系統誤差時,測得值的期望值是真值)。隨機誤差元的期望值是零。隨機誤差元為:
ξi = Xi - EX (1)
2)標準誤差定義為
σ = √(1/N)∑ξi2
= √(1/N)∑(Xi-EX)2 (2)
3)貝塞爾公式是用測得值的平均值代換(2)式中的期望值,得到:
σ =√{[1/(N-1)]∑[X-X(平)]2} (3)
4)隨機誤差范圍
R(隨) = 3σ=3√(1/N)∑ξi2
=√(1/N)∑(3ξi)2 (4)
5)由公式(4),有:
R(隨) =3σ= σ(3ξ) (5)
如(2)、(3),σ是隨機誤差元標準誤差。
如(5),3σ、σ(3ξ) 是隨機誤差范圍。
-
以3ξ為隨機誤差元,其對誤差范圍的權重為1,與系統誤差元權重相同。 因而以3ξ為隨機誤差元,就可以同系統誤差等權地進行誤差合成。這是方根法的“一從眾”。
-
3 單項系統誤差元構成的誤差范圍
系統誤差元用β表示。β是或正或負的恒值。
單個系統誤差構成的誤差范圍
R(系) =√{(1/N)∑(βi)2} =√(β2)
= |β| (6)
單個系統誤差對誤差范圍的貢獻值是該系統誤差的絕對值。
-
作者: njlyx 時間: 2016-5-27 16:06
本帖最后由 njlyx 于 2016-5-27 16:11 編輯
既然要講“包含概率”,“認知誤差量的分布規律”便是不可回避、必須要辦的事!——辦起來自然有點難,主要靠經驗(包括前輩傳授),有時也可能要“膽識”。
回避“認知誤差量的分布規律”的必然結果——“包含概率”的含糊其辭! 說是“大于99%”,其實也說不清這“大于99%”如何得以保障?
預計的應用效果將是: “愚公”們苦心竭力“認知誤差量的分布規律”,可能會“評估”出一個"誤差(范圍)”值Δ1,承諾“包含概率大于99%”;回避“認知誤差量的分布規律”的“智瘦”們可能會“理直氣壯”的“給”出一個值為3Δ1的“誤差(范圍)”,其“包含概率”是多少呢?也只能說“大于99%”,沒有膽量說“大于99.9%”!
作者: 史錦順 時間: 2016-5-28 07:28
本帖最后由 史錦順 于 2016-5-28 07:51 編輯
-
交叉系數決定合成法(3)
-
史錦順
-
4 誤差合成的理論基礎
函數的改變量,等于函數對各個自變量偏微分的和。就是泰勒展開的一級近似。
f(x,y)= f(xo,yo)+ (?f/?x) (x-xo)+ (?f/?y) (y-yo) (7)
f(x,y) -f(xo,yo) =(?f/?x) Δx+ (?f/?y) Δy (8)
Δf=(?f/?x)Δx + (?f/?y)Δy (9)
公式(9)是偏差關系的普遍形式。對所研究的特定函數來說,?f/?x、?f/?y是常數。
偏差關系用于測量計量領域,x是測得值,xo是真值, Δx是測得值x的誤差元;y是測得值,yo是真值,Δy是測得值y的誤差元;f(x,y)是測量儀器測得值或是間接測量被測量的測得值,簡稱函數值,f(xo,yo) 是函數的真值,Δf=f(x,y)-f(xo,yo) 是函數值的誤差元。
5 交叉系數的一般表達
設函數的誤差由兩項誤差Δx、Δy引起。由此,函數的兩項誤差元為:
Δf(x) =(?f/?x) Δx
Δf(y) =(?f/?y) Δy
把分項誤差作用的靈敏系數與該項誤差歸并,記為:
Δf(x) = ΔX
Δf(y) = ΔY
函數的誤差元式(9)變為:
Δf=ΔX+ΔY (10)
-
對(10)式兩邊平方并統計平均:
(1/N)∑Δf2=(1/N)∑(ΔXi +ΔYi)2
=(1/N)∑ΔXi2 + 2(1/N)∑ΔXiΔYi+(1/N)∑ΔYi2
RΔf2 = RΔX2 + 2(1/N)∑ΔXiΔYi + RΔY2 (11)
(11)式右邊的第一項為ΔX范圍的平方RΔX2;第三項為ΔY范圍的平方RΔY2;第二項是交叉項,是我們研究的重點對象。
-
交叉項為
2(1/N)∑ΔXiΔYi = 2 [(1/N)(∑ΔXiΔYi) / (RΔXRΔY)] × (RΔXRΔY)
= 2 J RΔXRΔY (12)
(12)式中的J為:
J = (1/N)(∑ΔXiΔYi ) / (RΔXRΔY) (13)
稱 J 為交叉系數。
(注:此前,J記為r,稱為相關系數。這和統計理論的相關系數,物理意義有差別。為澄清已有的混淆,本文稱J為交叉系數。)
當交叉系數可略時,誤差范圍的合成公式(11)變為:
RΔf = √ (RΔX2+RΔY2) (14)
(14)式是“方和根”合成公式。
當交叉系數為+1時,誤差范圍的合成公式變為“絕對和”:
RΔf =|ΔX|+|ΔY | = RΔX +RΔY (15)
-
作者: 史錦順 時間: 2016-5-28 10:51
本帖最后由 史錦順 于 2016-5-28 11:04 編輯
-
交叉系數決定合成法(4)
-
史錦順
-
6 隨機誤差間合成的交叉系數
對隨機誤差的合成,ΔX是ξx, 代換為[X-X平];ΔY是ξy,代換為[Y-Y平],有:
J =[1/(N-1)]∑(Xi-X平)(Yi-Y平) / (σΔX σΔY) (16)
由于ξx、ξy是隨機誤差,可正可負,可大可小,有對稱性與有界性,多次測量,是大量的,因此,隨機誤差間的合成的交叉系數為零(或可以忽略)。(15)式是當前不確定度論引用的統計理論的相關系數公式。這個公式對隨機誤差是對的;對系統誤差,不成立。
隨機誤差合成,(14)成立。即隨機誤差的合成公式是“方和根”:
RΔf = √ (RΔX2+RΔY2) (14)
σΔf = √[σΔx2+ σΔy2] (14.1)
7 隨機誤差與系統誤差合成的交叉系數
兩個分項誤差,一個是隨機的,記為ξ,考慮到對誤差范圍的權重,取單元量為3ξ(對應ΔX);一個是系統的(重復測量中不變),記為β(對應ΔY)。
代入公式(13),有
J =(1/N)(∑3ξiβ) / [R(3ξ) R(β)]
系統誤差元β是恒值,可以提出來,有
J =(1/N) (3β∑ξi) / [R(3ξ) R(β)] (17)
大量重復測量(例如N=20,N不得小于10)中,(17)式的∑ξi等于零或可以忽略,因此J近似為0,可以忽略。“方和根法”成立:
R(f) =√[β2+ (3σ)2] (18)
8 系統誤差與系統誤差合成的交叉系數
設(13)式中ΔX為系統誤差βx ,ΔY為系統誤差βy,有
RΔX = √ [(1/N)∑ΔXi2]= |βx| (19)
RΔY= √ [(1/N)∑ΔYi2]= |βy| (20)
則系統誤差的交叉系數為
J = (1/N)(∑βxi βyi) / (|βx||βy|)
= βx βy/ (|βx||βy|)
=±1 (21)
即有
|J|=1 (22)
當βx與βy同號時,系統誤差的交叉系數J為+1;當βx與βy異號時,系統誤差的交叉系數J為-1。
當系統誤差的交叉系數為+1時,(11)式變為:
RΔf2=|βx2|+2|βx||βy| +|βy|2 =(|βx|+|βy|)2
即有
RΔf = |βx|+|βy| (23)
(22)式就是絕對值合成公式。
當系統誤差的交叉因子為-1時,(22)式變為二量差的公式。因為通常只是知道系統誤差之誤差范圍,又鑒于誤差量“上限性”的特點,二量差的公式不能用。
測量儀器的性能指標,給出的都是誤差范圍。該指標值由生產廠家給出,由計量部門公證,測量者按儀器指標應用。直接測量,測量儀器的指標,就可看作是測量的誤差范圍(只要符合儀器使用條件,環境等的影響已包含在儀器的指標中)。間接測量,要按間接測量的函數關系進行誤差合成。測量儀器的誤差范圍指標值因以系統誤差為主,要視其為系統誤差值,按系統誤差處理。
-
作者: 史錦順 時間: 2016-5-28 15:54
本帖最后由 史錦順 于 2016-5-28 16:03 編輯
-
交叉系數決定合成法(5)
-
史錦順
-
9 關于合成方法的主張
誤差合成,統一按“方根法”。對特定的誤差種類,“方根法”分化為“均方根法”、“方和根法”、“絕對和法”、“混合法”。
通常,測量儀器以系統誤差為主。不能無視系統誤差的存在。考慮到系統誤差、隨機誤差都是客觀存在,提出如下主張:
1)隨機誤差序列,用“均方根法”,隨機誤差范圍之間,用“方和根法”;
2)隨機誤差范圍與系統誤差范圍之間,用“方和根法”;
3)有多項中小系統誤差項,僅有一項大系統誤差(或沒有大系統誤差),它們之間的交叉系數,可能是+1,也可能是-1,有相互抵消、或部分抵消的作用,這樣,可以用“方和根法”。
4)直接測量僅有兩三項系統誤差,要用“絕對和法”(適用于研制中確定儀器指標);
5)間接測量,僅有兩三項測量儀器的誤差范圍,要用“絕對和法”;
6)有多項誤差,在兩項或三項大系統誤差之間用“絕對和法”,其余的各種處理,用“方和根法”。總稱謂是“混合法”。
10 間接測量的誤差合成例說
間接測量由若干直接測量構成。各直接測量的誤差,都是間接測量的誤差因素。還加一些綜合性因素。
間接測量,要進行若干項分項誤差的合成。
設函數誤差由以下8項誤差構成:
大系統誤差項β1大、β2大
小系統誤差項β3小、β4小、β5小、β6小、
隨機誤差項ξ7隨、ξ8隨
注:
分項系統誤差的傳遞系數是函數對該自變量的偏微商。
分項隨機誤差的傳遞系數是函數對該自變量的偏微商的3倍(包含概率99%)。
本文中分項誤差項的值,指單項誤差與傳遞系數的乘積。
函數誤差元
Δf =(?f/?x)Δx + (?f/?y)Δy……
Δf =β1大+β2 大
+β3小+β4小+β5小+β6小
+3ξi7隨+3ξi8隨
求“函數誤差元的平方”的統計平均
[(1/N)∑Δfi2]
= (1/N)∑[β1大+β2大]
+β3小+β4小+β5 小+β6小
+3ξi7隨+3βi8隨]2
R2 = (1/N)∑[β1大2+2J大β1大β2大+β2大2
+β3小2+β4小2+β5小2+β6小2
+(3σ7 隨)2+(3σ8隨)2+其他交叉項] (24)
大系統誤差項的交叉系數J大等于+1或-1;因誤差范圍是誤差元的最大可能值,故取+1。由此,大誤差間取絕對和。其他交叉項的交叉因子,凡有隨機誤差項的,交叉因子為零。沒有隨機誤差的,是系統誤差之間的交叉系數,可以是+1,也可以是-1;由于交叉項的數量大,可認為正負項近似抵消,因而其他交叉項之和可略。
合成誤差范圍公式
R =√[(R1大+R2大)2
+R3小2+R4小2+R5小2+R6小2
+(3σ7隨)2+(3σ8隨)2] (25)
二、三項大系統誤差間取“絕對和”;此“絕對和”與所有其他系統誤差、隨機誤差范圍之間,取方和根。
由于測量儀器的誤差范圍,以系統誤差為主,且因誤差范圍是誤差元絕對值的一定概率(99%)意義上的最大可能值,因此某項直接測量的測量儀器誤差范圍指標值,視為間接測量的該項系統誤差。
當分項誤差僅有一項大誤差,或有4項以上大誤差時,考慮交叉項的可能抵消作用,公式(10)變成純“方和根”。
-
全文完。歡迎批評!
-
作者: njlyx 時間: 2016-5-28 17:37
敬佩史先生的執著精神!
可惜本篇長論似乎未現真珠?
若是不計較“包含概率”的確切值【——大于99.xx%就好,不計較‘大于’的代價!】,那么,現有“誤差理論”【譬如費業泰先生著作】對“誤差(范圍)”的“合成”方法已然描述清楚,此篇或未現長處。
長篇難如先生期許效用的“癥結”或在:誤解了所謂“系統(測量)誤差”的本性! 以它與所謂“隨機(測量)誤差”配對,將其歸入了“非隨機”量(確定量?)之列,以致在“推導”中將其【——所謂“系統(測量)誤差”】樣本值(“合成”時未知!)與“范圍”值混為一談,得到有些牽強附會的“結論”。.....誠如本論壇的葉先生所言,所謂的“(未定)系統(測量)誤差”也是一個“不確定量”【——“隨機量”!】,因此才要關注它的“范圍”值!(附言:本人不贊成葉先生全盤否定“誤差分類”效用的觀點!)
若是計較“包含概率”的確切值【——不是大于99.xx%就好,要計較‘大于’的代價!追求“剛好達到”約定的99.xx%】,便應積極響應“不確定度”評估中所提諸法(當然,也包括積極針砭其可能存在的缺點)!
而“相關性”問題,則是新、老處理辦法都不能“回避”的! 可以據“理”簡化處理,不能換個“名字”敷衍行事。
作者: 史錦順 時間: 2016-6-7 07:54
本帖最后由 史錦順 于 2016-6-7 08:09 編輯
-
誤差合成法公式推導中系統誤差恒值的時間要求
-
史錦順
-
誤差是測得值與被測量真值的差距。
誤差元是測得值減真值。
誤差范圍是誤差元絕對值的一定概率(大于99%)意義上的最大可能值。
-
誤差元,有隨機的(在一場測量的N次測量中,大小符號都在變化),也有恒值的。這是誤差量的性質。按性質,誤差被劃分為為隨機誤差和系統誤差。
-
測得值是測量機制中,各種因素共同構成的結果。函數的誤差,取決于各個分項誤差。
推導誤差合成法的公式,要根據各分項誤差的作用機理。誤差合成法必須符合誤差的性質,反映誤差的性質。
現代誤差理論一律地取“方和根”,忽視了誤差性質上的不同。交叉系數理論給出的結果是,合成法取決于誤差的性質或系統誤差的數量,這就反應了誤差性質的不同。
-
測量量值的方式,是測量N次取平均。這個統計平均時間稱“統計時間”。
討論誤差合成公式時,所認定的分項誤差的性質,是指“統計時間”內的性質。
誤差合成法公式推導中系統誤差恒值的時間要求,是統計時間。
-
有人說,系統誤差就長期來說并不是恒值。是的,使用的恰當,物盡其用,系統誤差的變化有接近誤差范圍指標值的。如頻率計準確度指標是5E-8,頻率源是高穩晶振,如果此晶振的老化規律穩定、精確測得其日老化率為+2e-10,,則校頻時,可置晶振的頻率為-4E-8,一年內變化到不超過+4E-8,這樣可使頻率計有5E-8的準確度指標。而計量時把頻率精確地調準到1E-10,而一年內可能達到+7E-8,反而超標了。
這樣的系統誤差變化,是筆者多次面對的。筆者所說的“恒值”,從來沒說過它永久不變。
-
那么在“交叉系數決定合成法”的公式推導中,要求“系統誤差保持恒值”的時段是多大呢?僅僅是統計平均時間,就是一場N次測量所用的時間,大致是幾分鐘到幾小時。在這短短的時段內,“系統誤差保持恒值”這一點,是沒有問題的!
-
系統誤差的長時間后的變化,不影響關于“交叉系數決定合成法”的道理。因為決定誤差大小的,是統計時間內的誤差量的性質。在統計時間內,隨機誤差大小符號都在變化,20個(或100個)隨機誤差元在交叉項中,相互抵消,隨機誤差間合成,隨機誤差與一項系統誤差合成,交叉系數都近于零,誤差可取方和根。兩項系統誤差合成,在統計時間內,兩個系統誤差都是恒值,由它們決定的交叉系數,是+1或-1。就是說,是絕對和,或者是絕對差。誤差范圍的定義是誤差元絕對值的最大可能值,因此交叉系數要取+1,合成方法應是絕對和。
-
作者: csln 時間: 2016-6-7 08:43
如頻率計準確度指標是5E-8,頻率源是高穩晶振,如果此晶振的老化規律穩定、精確測得其日老化率為+2e-10,,則校頻時,可置晶振的頻率為-4E-8,一年內變化到不超過+4E-8,這樣可使頻率計有5E-8的準確度指標。而計量時把頻率精確地調準到1E-10,而一年內可能達到+7E-8,反而超標了。
精確測得晶振日老化率為+2E-10,按檢定規程要求應把晶振頻率準確度校準到優于-2E-9,怎么可以置晶振的頻率為-4E-8,有那一個生產廠會把日化率2E-10的晶振給出4E-8準確度的指標,如此設置生產廠不答應,用戶也不會答應
作者: csln 時間: 2016-6-7 08:51
本帖最后由 csln 于 2016-6-7 09:00 編輯
誤差合成法公式推導中系統誤差恒值的時間要求,是統計時間。
如何知道統計時間內恒值是正、是負、是大、是小、是多少,測量出來了嗎?若測量出來已知恒值是多少,njlyx說了:便可以直接代入{Y=f(X)}折算得到相應的輸出“誤差”(分)量值(不是“范圍”值!),根本沒有取“方和根”還是“絕對和”的問題!!!。若沒有測量出來不知道恒值在多少,就只能知道是在一個區間內,分布、相關性是必須要考慮的。
| 歡迎光臨 計量論壇 (http://m.dy313.com/) |
Powered by Discuz! X3.4 |
久久伊人一区二区_在线日韩av_欧美男男video_国产精品性做久久久久久
99热这里都是精品|
国产很黄免费观看久久|
中文字幕av一区二区三区|
日韩欧美一区在线观看|
欧美一区二区三区免费大片|
91精品婷婷国产综合久久性色|
色香蕉久久蜜桃|
91蝌蚪porny九色|
欧美性感一区二区三区|
欧美亚洲国产一区二区三区|
欧日韩精品视频|
欧美亚洲一区二区三区四区|
欧美日韩在线免费视频|
欧美卡1卡2卡|
欧美成人综合网站|
国产欧美视频在线观看|
亚洲欧美在线aaa|
亚洲精品你懂的|
天堂蜜桃91精品|
激情深爱一区二区|
丰满白嫩尤物一区二区|
91蝌蚪porny九色|
欧美精品v国产精品v日韩精品
|
国产精品自拍在线|
成熟亚洲日本毛茸茸凸凹|
成人毛片在线观看|
欧美在线高清视频|
精品毛片乱码1区2区3区|
国产女同互慰高潮91漫画|
亚洲柠檬福利资源导航|
日韩黄色免费电影|
成人网在线播放|
欧美精品一卡两卡|
欧美激情一区二区三区四区|
亚洲综合一区二区三区|
久久福利视频一区二区|
成人精品国产免费网站|
欧美日韩精品一区二区三区蜜桃|
日韩色在线观看|
亚洲视频一区二区免费在线观看|
蜜臀久久久久久久|
不卡视频在线看|
欧美一级黄色片|
国产精品的网站|
国产在线精品一区二区夜色|
色综合久久久久久久|
欧美成人video|
亚洲靠逼com|
国产精品99久久久久久久女警
|
国产精品久久久久久久浪潮网站|
亚洲免费在线观看视频|
久久99久久久欧美国产|
欧美偷拍一区二区|
中文字幕+乱码+中文字幕一区|
日韩极品在线观看|
色欧美88888久久久久久影院|
精品国产3级a|
日韩精品乱码免费|
91国在线观看|
亚洲欧美aⅴ...|
懂色av一区二区三区免费看|
精品噜噜噜噜久久久久久久久试看|
亚洲精品免费一二三区|
成人av免费在线|
中文字幕乱码久久午夜不卡|
韩国一区二区在线观看|
91精品蜜臀在线一区尤物|
亚洲一区在线视频观看|
99久久精品免费看国产免费软件|
国产欧美一区二区精品性色|
国内精品国产成人国产三级粉色
|
久久99国产精品久久|
欧美日韩精品专区|
亚洲妇女屁股眼交7|
色88888久久久久久影院按摩
|
亚洲一区免费在线观看|
成人免费视频免费观看|
国产精品欧美综合在线|
国产白丝网站精品污在线入口|
久久综合色一综合色88|
国产一区二区三区四区五区美女
|
久久久精品国产免大香伊|
久久国产欧美日韩精品|
日韩亚洲欧美中文三级|
久久成人免费日本黄色|
日韩欧美aaaaaa|
国产综合久久久久久鬼色|
精品国产sm最大网站免费看|
国产伦理精品不卡|
国产精品看片你懂得|
97久久精品人人做人人爽|
亚洲另类春色国产|
4438亚洲最大|
麻豆视频一区二区|
久久日一线二线三线suv|
国产乱码精品一区二区三区五月婷|
久久九九全国免费|
波多野结衣一区二区三区|
亚洲黄色小说网站|
日韩一区二区三区免费观看|
黑人精品欧美一区二区蜜桃|
亚洲国产高清在线|
欧美在线你懂的|
另类人妖一区二区av|
国产精品免费人成网站|
色噜噜夜夜夜综合网|
免费在线观看成人|
国产拍揄自揄精品视频麻豆|
在线观看www91|
国产激情一区二区三区四区|
亚洲人成在线播放网站岛国|
91精品婷婷国产综合久久竹菊|
国产精品一区一区三区|
亚洲精品视频在线|
日韩亚洲欧美高清|
97久久精品人人做人人爽|
免费高清在线视频一区·|
欧美国产成人在线|
制服丝袜中文字幕一区|
成人动漫在线一区|
免费在线观看一区|
亚洲综合网站在线观看|
久久亚洲捆绑美女|
欧美视频完全免费看|
成人一区在线观看|
视频一区国产视频|
椎名由奈av一区二区三区|
久久综合五月天婷婷伊人|
欧美丝袜自拍制服另类|
成人av在线一区二区三区|
免费看日韩a级影片|
一区二区三区四区国产精品|
精品久久久久久久久久久久久久久久久|
www.99精品|
国产在线视频一区二区|
麻豆一区二区三|
天涯成人国产亚洲精品一区av|
18涩涩午夜精品.www|
国产欧美一区二区精品性色超碰|
日韩一二三区不卡|
91精品国产综合久久国产大片|
91免费国产视频网站|
成人精品视频一区二区三区
|
精品国产髙清在线看国产毛片|
欧美视频你懂的|
91久久国产最好的精华液|
成人网页在线观看|
国产成人精品免费|
成人禁用看黄a在线|
国产激情视频一区二区在线观看
|
亚洲欧美国产77777|
国产精品青草久久|
国产精品不卡在线观看|
国产精品色一区二区三区|
亚洲国产精品激情在线观看
|
91亚洲精品久久久蜜桃网站|
懂色一区二区三区免费观看|
国产91丝袜在线观看|
国产 欧美在线|
成人av电影在线|
一本在线高清不卡dvd|
欧美亚洲国产bt|
欧美久久高跟鞋激|
欧美精品丝袜中出|
日韩视频免费观看高清完整版
|
中文一区二区在线观看
|
在线中文字幕一区二区|
在线亚洲免费视频|
6080亚洲精品一区二区|
日韩美一区二区三区|
中文字幕第一区第二区|
中文字幕一区二区三区av|
亚洲黄色性网站|
日韩国产高清在线|
精品午夜一区二区三区在线观看|
国产成人在线视频网站|
91网上在线视频|
91精品国产91久久久久久一区二区|
欧美电视剧在线看免费|
国产日韩成人精品|
亚洲一区在线视频观看|
奇米综合一区二区三区精品视频|
国产精品一二三四|
在线观看不卡一区|
久久众筹精品私拍模特|
亚洲色图在线播放|
日av在线不卡|
caoporen国产精品视频|
欧美精品第1页|
国产精品久久久久一区二区三区
|
国产精品网站在线观看|
一区二区在线电影|
久久综合综合久久综合|
国产成人精品免费看|
欧美日韩国产首页在线观看|
国产亚洲精品中文字幕|
天天综合网 天天综合色|
成人深夜在线观看|
日韩天堂在线观看|
亚洲精品va在线观看|