久久伊人一区二区_在线日韩av_欧美男男video_国产精品性做久久久久久
計量論壇
標題:
誤差范圍區間與擴展不確定度區間基于計算的比較(續1)
[打印本頁]
作者:
史錦順
時間:
2018-11-18 08:36
標題:
誤差范圍區間與擴展不確定度區間基于計算的比較(續1)
本帖最后由 史錦順 于 2018-11-18 09:07 編輯
-
誤差范圍區間與擴展不確定度區間
基于計算的比較
(續1)
-
史錦順
-
本文論題:測量場合下,間接測量的測量結果區間的比較。前文《誤差范圍區間與擴展不確定度區間基于計算的比較圖》是本文的一個特例。
本文計算題目:一把精致銀米尺,長、寬、厚,形狀規整。測量該尺的物質密度。
米尺的基本數據(參數分項測量,各10次)
測量質量M,用普通天平(MPEV 0.03%)。測得值:M[sub]m [/sub]= 8.900kg
長度L的規格:長度的標稱值L[sub]o[/sub]= 1m,MPEV 0.03%
測量寬度D,用數顯卡尺(MPEV 0.05%),測得值:D[sub]m[/sub] =50.00mm。
測量厚度H,用數顯千分尺(MPEV 0.03%)。測得值:H[sub]m[/sub] =20.000mm。
1《史法測量計量學》
對誤差范圍的表達
1.1 測量方程與測得值函數
物質密度ρ的物理公式
ρ = 質量/體積=質量/ (長×寬×厚)
= M/(LDH) (1)
物質密度ρ的計值公式(腳標o表示標稱值;腳標m表示測得值)
ρ[sub]m[/sub] = M[sub]m[/sub] /(L[sub]o[/sub]D[sub]m[/sub]H[sub]m[/sub]) (2)
測量方程為
ρ[sub]m [/sub]/ ρ = (M[sub]m[/sub]/M) / [(L[sub]o[/sub]/L)( D[sub]m[/sub]/D)(H[sub]m[/sub]/H)] (3)
測得值函數:
ρ[sub]m [/sub]= ρ(M[sub]m[/sub]/M) (L/L[sub]o[/sub])(D/D[sub]m[/sub])(H/H[sub]m[/sub]) (4)
注意,測得值函數中,自變量是質量測得值M[sub]m[/sub]、長度L、寬度測得值D[sub]m[/sub]、厚度測得值H[sub]m[/sub]。
寬度實際值D、厚度實際值H、質量實際值M是物理公式中的值,是真值,是常量。長度標稱值L[sub]o[/sub]是常量(標稱值是定義值的一種形式[1]),而長度實際值L是真值,卻是變量。以往,經典誤差理論中的誤差分析,直接從物理公式出發進行微分,常常弄錯有標稱值的量的誤差元的符號。如頻率計誤差分析中的時標誤差。由于誤差量的特點是其“絕對性”與“上限性”,經典誤差理論的結果分析,沒錯。如馬鳳鳴書中直接取絕對值,不出錯。但如果要修正,或者要分析誤差因素的影響,符號問題就必須講究。誤差理論要精準,要嚴格化,不能在符號上含糊?!妒贩ā返臏y量方程與測得值函數,澄清了這一點。關鍵是辨清常量與變量。微分必須對變量微分。《史法》的區分常量與變量的程序,使誤差理論的微分操作,有了嚴格的數理邏輯基礎。
1.2 求誤差元
A 差分法
誤差元是測得值與真值之差。測得值(2)與真值(1)之差就是誤差元。這樣做,物理意義十分清晰。筆者年輕時讀文獻,不理解“為什么求微分就是求誤差”,只是模仿。臨近退休,才弄明白誤差方程這套道理。取差分,符合誤差(元)的定義。微分不過是求差的數學方法。而取微分,必須明確變量。
差分法求密度誤差的公式為
Δρ[sub]m[/sub] =ρ[sub]m[/sub] - ρ
= M[sub]m[/sub] /(L[sub]o[/sub]D[sub]m[/sub]H[sub]m[/sub]) - M/(LDH) (5)
對(5)式做近似計算,即可求誤差元。(參見《史法測量計量學》[[sup]1][/sup]第3章。)
B 微分法
熟悉微積分,用微分法就更方便。
對測得值函數ρ[sub]m[/sub](4)
ρ[sub]m[/sub] = ρ(M[sub]m[/sub]/M) (L/L[sub]o[/sub])(D/D[sub]m[/sub])(H/H[sub]m[/sub])
做全微分,得密度函數的誤差元
dρ[sub]m [/sub]=ρ (?ρ[sub]m[/sub]/?M[sub]m[/sub])dM[sub]m[/sub]+ (?ρ[sub]m[/sub]/?L)dL+ (?ρ[sub]m[/sub]/?D[sub]m[/sub])dD[sub]m[/sub]+ (?ρ[sub]m[/sub]/?H[sub]m[/sub])dH[sub]m
[/sub] dρ[sub]m[/sub] = ρ(M[sub]m[/sub] / M) (L / L[sub]o[/sub])( D / D[sub]m[/sub])(H / H[sub]m[/sub])[ dM[sub]m[/sub]/M[sub]m[/sub]]
+ ρ(M[sub]m [/sub]/ M) (L / L[sub]o[/sub])( D / D[sub]m[/sub])(H / H[sub]m[/sub]) [ dL/L]
- ρ(M[sub]m [/sub]/ M) (L / L[sub]o[/sub])( D / D[sub]m[/sub])(H / H[sub]m[/sub]) [ dD[sub]m[/sub]/D[sub]m[/sub]]
- ρ(M[sub]m[/sub] / M) (L / L[sub]o[/sub])( D / D[sub]m[/sub])(H / H[sub]m[/sub]) [ dH[sub]m[/sub]/H[sub]m[/sub]]
dρ[sub]m [/sub]/ ρ[sub]m[/sub] = dM[sub]m[/sub]/M[sub]m[/sub]+ dL/L- dD[sub]m[/sub]/D[sub]m [/sub]– dH[sub]m[/sub]/H[sub]m[/sub]
標成相對誤差元的形式,有
rρ[sub]m[/sub] =r[sub]Mm[/sub] + r[sub]L[/sub] - r[sub]Dm[/sub] - r[sub]Hm [/sub] (6)
1.2 誤差合成公式
由于四項誤差(本例已知分項誤差范圍的指標值),只有一項偏大,其他大小相近,沒有“二、三項大誤差”項,《史法》[1] 誤差合成,取各項誤差范圍的“方和根”,密度測量的誤差范圍(合成誤差元絕對值99%概率意義上的最大可能值)是:
R
ρ
= √(R[sub]M[/sub][sup]2[/sup] + R[sub]L[/sub][sup]2[/sup] + R[sub]D[/sub][sup]2[/sup] + R[sub]H[/sub][sup]2[/sup]) (7)
1.3 本例具體計算
尺體質量密度的測得值
ρ[sub]m[/sub] = M/ L[sub]o[/sub]D[sub]m[/sub]H[sub]m
[/sub] =8.900kg/1000.0mm×50.00mm×20.000mm
= 8.900×10[sup]3 [/sup]kg/m[sup]3
[/sup] 密度測得值的相對誤差范圍(7)
Rρ = √[(0.03%)[sup]2[/sup] + (0.03%)[sup]2 [/sup]+ (0.05%)[sup]2[/sup] + (0.03%) [sup]2[/sup]]
=0.072%
密度測得值的絕對誤差范圍
Δ ρ[sub]m[/sub] =8.9×10[sup]3[/sup]×0.00072 kg/m[sup]3[/sup]
=6.4 kg/m[sup]3[/sup] (8)
米尺質量密度的測量結果:
ρ=(8900±6) kg/m[sup]3[/sup] (9)
(未完,待續)
-
——
作者:
史錦順
時間:
2018-11-18 10:36
本帖最后由 史錦順 于 2018-11-18 10:56 編輯
(續前)
2 不確定度體系的計算
2.1 基本數據及分布認定
測量質量M:普通天平(MPEV 0.03%,均勻分布)。測得值:M = 8.900kg,u[sub]AM[/sub]= 0.007%
長度L的規格:長度的標稱值L[sub]o[/sub] = 1m,MPEV 0.03%,正態分布。
測量寬度D:數顯卡尺(MPEV 0.05%,均勻分布)。測得值:D[sub]m[/sub] =50.00mm,u[sub]AD[/sub]= 0.01%
測量厚度H:數顯千分尺(MPEV 0.03%,均勻分布)。測得值:H[sub]m[/sub] =20.000mm,u[sub]AH[/sub]= 0.007%
2.2 擴展系數的確定
用單臺儀器的直接測量,直接用儀器的誤差范圍的指標值MPEV,就夠了。沒有誤差合成問題,也就沒有必要先除以根號3,再乘以k=1.7。那是將MPEV復原。似乎道理講得通,但這是沒有實際意義的折騰。(誤差理論認為:基礎測量是常量測量,儀器的隨機誤差包含在儀器指標MPEV中,不另算;不確定度體系認為測量對象可能是常量,也可能是隨機變量,不確定度評定要包括測量值的隨機變化部分。)
不確定度體系是“方差合成”,一個關鍵問題是,多項誤差因素合成為u[sub]C[/sub],此后擴展因子怎樣???(《史法》是“范圍合成”不存在這個問題。)
不確定度體系的規范,多處說明,通常要求95%的概率,可取k=2;要求99%的概率,可取k=3。但這樣做的前提是合成后誤差是正態分布。又說各項誤差是獨立的,相互之間不發生作用。取99%概率時,何以將均勻分布的那一項除以根號3,卻乘以3,竟毫無道理的擴大該項誤差的作用到1.73倍?這是把正態分布的隨機誤差的規律,錯誤地用于統計中保持恒值的系統誤差的惡果。毫無道理,是錯誤的。錯誤的分布認定,錯誤的擴展因子,只因為是把“恒值的系統誤差”當成“隨機誤差”來處理。錯誤的根源是把測量的“時域統計”錯誤地當成是“臺域統計”[2]。
有人說:乘2乘3是不符合規矩的錯誤操作。其實不然。且看如下規范:
《校準領域測量不確定度評估指南》(CNAS-GL09:2008)
p8:
5.2 用實驗的方式通常很難確認輸出量的估計值是否屬于正態分布。但是,在通常情況下,當幾個不確定度的分量來源于概率分布獨立的輸入量(即參量個數N≥3),如服從于正態分布或矩形分布時,則依據中心極限定理可以假設輸出量的分布接近于正態分布
。
(史注:該規范中的實例S2、S3、S4、S5、S6、S7都按“輸出量為正態分布”處理.)
本題有誤差范圍四項,隨機誤差三項。常規的分布認定:正態分布四項,矩形分布(均勻分布)三項,求U[sub]99[/sub],取k=3,符合不確定度體系作法。至于認為各項獨立,也只能如此。(不確定度體系的“分布論”、“相關系數論”本來就是多余的,并無科學根據。)
-
2.3 不確定度體系的計算
尺體密度測得值
ρ[sub]m[/sub] = M[sub]m [/sub]/ L[sub]o[/sub]D[sub]m[/sub]H[sub]m
[/sub] =8.900kg/1000.0mm×50.00mm×20.000mm
= 8.900×10[sup]3[/sup]kg/m[sup]3
[/sup] 合成標準不確定度
u[sub]C[/sub] = √ {[(0.03%)/√3][sup]2[/sup] + [(0.03%)/3][sup]2[/sup] + [(0.05%)/√3][sup]2[/sup]+ [(0.03%)/√3]2
+(0.007%)[sup]2[/sup]+(0.01%)[sup]2[/sup]+(0.007%)[sup]2[/sup]}
= √ [(0.0173%)[sup]2[/sup]+ (0.01%)[sup]2 [/sup]+ (0.0289%)[sup]2[/sup] + (0.0173%) [sup]2 [/sup]
+(0.007%)[sup]2[/sup]+(0.01%)[sup]2[/sup]+(0.007%)[sup]2[/sup]}
=0.042% (10)
包含概率99%的擴展不確定度(相對值)
U[sub]99r[/sub] = 3u[sub]C
[/sub] = 0.13% (11)
-
2.4 不確定度體系的測量結果表達
測得值
y = 8.900×10[sup]3[/sup] kg/m[sup]3
[/sup] 包含概率99%的區間半寬:
U[sub]99[/sub] = 8.900×10[sup]3[/sup] kg/m[sup]3[/sup]×0.13%
=12 kg/m[sup]3[/sup] (12)
測量結果:
Y = (8900±12) kg/m[sup]3[/sup] (13)
3 史評
拙文《論不確定度體系的公式錯誤》[2]中的區間比較圖以及前文《誤差范圍區間與擴展不確定度區間基于計算的比較圖》,畫的是單臺儀器。有人認為擴展不確定度的擴展系數的選取是不對的,因此,比較是不成立的。這樣說似乎有道理,但本質上是掩護了不確定度體系的錯誤。
誤差合成,總是若干量的合成。因此《論不確定度體系的公式錯誤》[2]的舉例,以及前文,應看成是多個分量中的一個。乘3是常規。
以(10)式為例,對整體乘3,本質是對各分項都乘3了。
u[sub]C [/sub]= √ {[(0.03%)/3][sup]2 [/sup]+[(0.03%)/√3][sup]2[/sup] + [(0.05%)/√3][sup]2[/sup] + [(0.03%)/√3] [sup]2[/sup]
+(0.01%)[sup]2[/sup]+(0.02%)[sup]2[/sup]+(0.01%)[sup]2[/sup]}
3
u[sub]C[/sub] = √ { [
3×
(0.03%)/3][sup]2[/sup] +[
3×
(0.03%)/√3][sup]2[/sup] + [
3×
(0.05%)/√3][sup]2[/sup]
+ [
3
×(0.03%)/√3] [sup]2[/sup] +[
3×
(0.01%)][sup]2[/sup]+[
3×
(0.02%)[sup]2[/sup]+[
3×
(0.01%)][sup]2[/sup]]
=0.042% (14)
其中的第2項(尺長)、第5項、第6項、第7項,因為單項是正態分布,合成說得過去;但第2項、第3項、第4項都無故擴大1.73倍,是沒有道理的,是錯誤的。因是“方和根”,大項誤差(第3項)較大,自身雖然只擴大到1.73倍,因為平方和凸顯大誤差的作用,且多計了在MPEV中本已包含的隨機誤差項,竟使函數區間半寬擴大到2倍!這算什么玩意兒?不確定度體系,偽科學!
-
本人前幾次的圖解說明,本質是數項誤差合成(N≥3)中單項作用的變化,是通常情況。這是不確定度體系的本質性錯誤之一。本人說明也不夠,易于引起誤解,故另寫兩篇續文在此。在否定不確定度體系這一點上,本人并無錯誤,而表達要盡可能更嚴格,是必要的。故對csln先生表示感謝。
-
不確定度體系的龐大架構:A類不確定度、B類不確定度、合成不確定度、擴展不確定度,這一套的核心是“誤差合成”。直接測量用MPEV就夠了;間接測量必須有誤差合成。就隨機誤差來說,不確定度體系與經典誤差理論并無區別,也與統計理論一致。但系統誤差不行。測量計量的統計,是“時域統計”,而不是“臺域統計”。在單臺儀器的重復測量中,系統誤差是恒值。在間接測量中,每臺儀器各自測量不同的量,各量間,有物理公式所確定的函數關系,但系統誤差卻是各自獨立的。每種儀器的系統誤差對被測函數量的作用,是獨立的。在多次重復測量中,各項系統誤差的量值、各系統誤差間的關系,由量值的函數關系確定,是固定的;即函數的總系統誤差是恒值。函數的系統誤差既不是均勻分布,通常也不可能是標準正態分布(無偏倚值),而一般是有偏正態分布(其偏倚值是系統誤差)。因此,對系統誤差的合成量乘以不等于1的因子k是錯誤的。下文將詳述。
直接進行“范圍合成”,簡單又正確,何不為之?
不確定體系的測量結果區間半寬U[sub]99[/sub]竟是誤差理論測量結果區間R的兩倍(本例計算),孰是孰非,不值得認真研究一番嗎?
不確定度體系錯了,且錯誤嚴重!要它何用?
-
資料及說明
[1]史錦順著《史法測量計量學》,等待出版。質檢出版社在等待國家計量院的鑒別結論。其內容本人已多次在本欄目征求意見。許多網友是知道的。
[2]史錦順《論不確定度體系的公式錯誤》,本欄目中有征求意見稿。宣傳新學說與抨擊不確定度體系的文章,本人已在本欄目貼出有署名的文章513篇(已編為10集)。有意向研究的網友,告訴我你的郵箱號,我可將全部文稿和書稿都寄給你。我已八十一歲半,傳播自己的學術思想,既是義務,也是一種快樂。
-
作者:
liuhuaxing
時間:
2018-11-18 15:38
太有深度了,我得慢慢消化。
作者:
njlyx
時間:
2018-11-18 22:12
關于所謂的"系統(測量)誤差",拎清后的"(測量)誤差理論"(如1990年代后費業泰等先生的著述)已分為"已(確)定"與"未(確)定"兩部分----對于"已(確)定系統誤差","間接測量"的"合成"就是"具體值代入函數運算",不涉及什么"范圍合成",當然也無關"概率分布形式"、"相關性"之類問題; 但對于"未(確)定系統誤差" ,人們只能"估計"它的"可能取值范圍"及相應的"概率分布形式","間接測量"時便涉及"可能取值范圍"的"合成"及"結果"的"概率分布形式"等問題---這些問題在"數學"(概率與統計)上不存在未明的概念問題,"線性化"理論"合成"、"蒙特卡洛"仿真"合成"都是行之有效的方法,實際應用中,各"分量"自身的"概率分布形式"、"分量"之間的"相關性"是不可回避的難題,需要可靠的"經驗"支持!
作者:
csln
時間:
2018-11-19 09:13
本帖最后由 csln 于 2018-11-19 09:16 編輯
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%的擴展不確定度(相對值)
U99r = 3uC
= 0.13% (11)
0.0289/0.042=0.688
有一均勻分布的顯著項占合成標準不確定度三分之二以上,合成不確定度仍然是均勻分布,U99包含因子是1.72或1.71
作者:
規矩灣錦苑
時間:
2018-11-19 11:07
本帖最后由 規矩灣錦苑 于 2018-11-19 11:10 編輯
“誤差區間范圍”與“擴展不確定度”是性質完全不同的兩個概念,無法擺在同一個平臺上相比較。
“誤差區間范圍”是允許的誤差變化范圍的寬度,量化反映了測量過程或測量結果準確性要求,正式術語是最大允許誤差絕對值,寫成英文縮寫就是MPEV。
“測量不確定度”量化反映了測量過程或測量結果的可信性,不是誤差,也不是誤差范圍。擴展不確定度只能用英文單詞的第一個字母U表示,當用小寫字母u時,代表標準不確定度。因此,擴展不確定度不能縮寫為MPEV。
另外,“不確定度”本身就是某個區間(人們估計得到的真值存在的可能區間)的“半寬度”,“區間半寬的區間”非常繞口,也令人費解,因此不存在“擴展不確定度區間”這個概念。
作者:
chownyo
時間:
2018-11-20 10:03
喜歡這種學術性帖子,收藏下來好好消化下
作者:
315046304351501
時間:
2018-12-24 17:04
史錦順 發表于 2018-11-18 10:36
(續前)
380310291@qq.com
作者:
夢回棗鄉
時間:
2021-4-16 14:11
1507933138@qq.com
作者:
史錦順
時間:
2021-4-17 08:48
夢回棗鄉 發表于 2021-4-16 14:11
2018年我的一句承諾,2021年了,我已不再兌現。我已84周歲。年老多病,請諒解。本欄目都有,自己查查吧。
按新版著作權保護法,網上文章、書稿,都在保護之列。可以下載自用,不可復制公開發表,否則即為侵權。作為不同意見的學術討論,可以部分引用,但不可全文另行刊出。
作者:
光風霽月大自在
時間:
2021-4-18 09:22
感謝史老師,健康長壽
歡迎光臨 計量論壇 (http://m.dy313.com/)
Powered by Discuz! X3.4
久久伊人一区二区_在线日韩av_欧美男男video_国产精品性做久久久久久
午夜精品福利在线
|
不卡一区在线观看
|
懂色av一区二区三区蜜臀
|
91精品1区2区
|
精品国产乱码久久久久久夜甘婷婷
|
成人午夜激情视频
|
717成人午夜免费福利电影
|
中文在线一区二区
|
精品一区二区在线免费观看
|
在线视频一区二区三区
|
国产精品久久久一区麻豆最新章节
|
三级在线观看一区二区
|
亚洲天天做日日做天天谢日日欢
|
久久精品国产久精国产爱
|
欧美日韩一卡二卡
|
亚洲综合在线观看视频
|
jlzzjlzz亚洲女人18
|
久久婷婷综合激情
|
激情国产一区二区
|
日韩一区二区三
|
亚洲丶国产丶欧美一区二区三区
|
91蝌蚪国产九色
|
中文无字幕一区二区三区
|
国产麻豆一精品一av一免费
|
日本视频中文字幕一区二区三区
|
91美女精品福利
|
1024精品合集
|
91福利区一区二区三区
|
亚洲天堂网中文字
|
91色porny在线视频
|
国产精品久久久久aaaa樱花
|
国产精品综合一区二区
|
久久久噜噜噜久久中文字幕色伊伊
|
日本成人在线网站
|
欧美日产国产精品
|
日本美女视频一区二区
|
日韩一卡二卡三卡
|
麻豆精品新av中文字幕
|
久久综合久久综合亚洲
|
国模冰冰炮一区二区
|
精品国产第一区二区三区观看体验
|
久久精品国产澳门
|
日韩欧美国产小视频
|
狂野欧美性猛交blacked
|
久久毛片高清国产
|
99久久久久免费精品国产
|
亚洲同性同志一二三专区
|
日本道色综合久久
|
日韩高清在线电影
|
精品国产一区a
|
波多野结衣在线aⅴ中文字幕不卡
|
国产精品福利一区二区三区
|
日本精品一级二级
|
日韩福利视频导航
|
国产区在线观看成人精品
|
国产乱码一区二区三区
|
国产精品视频一区二区三区不卡
|
99久久99久久精品免费观看
|
亚洲成人资源在线
|
精品欧美一区二区三区精品久久
|
国产成人在线视频网址
|
亚洲午夜电影网
|
精品国产精品网麻豆系列
|
色综合天天综合狠狠
|
免费一区二区视频
|
国产精品乱子久久久久
|
制服丝袜av成人在线看
|
国产99久久久国产精品潘金
|
一区二区在线观看不卡
|
337p日本欧洲亚洲大胆色噜噜
|
99久久99久久免费精品蜜臀
|
麻豆成人免费电影
|
中文字幕五月欧美
|
欧美成人午夜电影
|
欧美色图天堂网
|
国产成a人亚洲
|
日本亚洲电影天堂
|
1024亚洲合集
|
精品电影一区二区三区
|
91理论电影在线观看
|
国产一区二区三区久久悠悠色av
|
老司机精品视频导航
|
亚洲欧美另类久久久精品
|
欧美mv和日韩mv的网站
|
日本高清不卡视频
|
成人av在线看
|
国产在线精品不卡
|
免费欧美高清视频
|
午夜成人免费视频
|
亚洲色图制服丝袜
|
国产女人18毛片水真多成人如厕
|
99久久综合精品
|
久久99精品久久久
|
天天操天天色综合
|
亚洲精选视频免费看
|
国产女主播在线一区二区
|
精品国产电影一区二区
|
欧美一区二区三区免费观看视频
|
99久久精品国产毛片
|
成人午夜精品在线
|
国产98色在线|日韩
|
国内一区二区视频
|
麻豆一区二区在线
|
国产精品麻豆网站
|
国产精品乱码一区二三区小蝌蚪
|
国产婷婷色一区二区三区
|
亚洲精品在线三区
|
91精品国产一区二区
|
欧美日韩日日夜夜
|
7777精品久久久大香线蕉
|
欧美二区三区的天堂
|
欧美老人xxxx18
|
欧美一级日韩免费不卡
|
欧美日韩高清一区
|
欧美一区二区三区免费观看视频
|
欧美日韩三级在线
|
欧美一级片在线
|
欧美成人综合网站
|
国产亚洲女人久久久久毛片
|
国产亚洲综合在线
|
中文在线免费一区三区高中清不卡
|
国产亚洲人成网站
|
国产精品日日摸夜夜摸av
|
国产精品视频看
|
亚洲日本中文字幕区
|
一区二区三区四区精品在线视频
|
一区二区三区美女
|
亚洲成人午夜电影
|
日韩不卡一二三区
|
激情综合一区二区三区
|
国产精品一区三区
|
成人午夜短视频
|
色婷婷综合久色
|
欧美一区二区在线播放
|
久久人人爽爽爽人久久久
|
国产精品毛片久久久久久
|
亚洲综合一区二区
|
麻豆极品一区二区三区
|
波波电影院一区二区三区
|
欧美熟乱第一页
|
欧美大片拔萝卜
|
国产精品进线69影院
|
亚洲国产精品影院
|
精品伊人久久久久7777人
|
北岛玲一区二区三区四区
|
欧美午夜免费电影
|
欧美精品一区二区在线播放
|
国产精品久久毛片
|
日韩精品免费专区
|
www.欧美日韩
|
日韩一区二区影院
|
国产精品女同一区二区三区
|
亚洲成在线观看
|
丁香婷婷综合激情五月色
|
欧美日韩国产中文
|
国产精品美女久久久久av爽李琼
|
亚洲一区二区三区在线
|
国产成人精品影视
|
91精品国产综合久久香蕉的特点
|
国产精品入口麻豆九色
|
轻轻草成人在线
|
91黄色免费看
|
国产日韩欧美综合在线
|
亚洲国产成人porn
|
av一二三不卡影片
|
日韩美女视频一区二区在线观看
|
亚洲色图一区二区
|
国产**成人网毛片九色
|
91精品婷婷国产综合久久性色
|
亚洲欧洲日韩综合一区二区
|
精品一区二区三区日韩
|
欧美日韩卡一卡二
|
亚洲欧洲综合另类在线
|
国产精品18久久久久久vr
|
91精品国产综合久久精品
|
亚洲视频在线观看一区
|
懂色av中文一区二区三区
|
日韩免费成人网
|
强制捆绑调教一区二区
|
欧美日韩在线精品一区二区三区激情
|
欧美激情一区二区三区在线
|
另类小说图片综合网
|
51午夜精品国产
|
五月婷婷久久丁香
|
日本道色综合久久
|
一区二区三区精品
|
99精品视频一区二区
|
...av二区三区久久精品
|
成人一区二区三区在线观看
|
国产女人18水真多18精品一级做
|
国产在线播放一区三区四
|
欧美一卡二卡三卡
|
免费的成人av
|
日韩欧美www
|
狠狠色伊人亚洲综合成人
|
精品免费国产一区二区三区四区
|
麻豆freexxxx性91精品
|
精品国产一区二区精华
|
国内久久婷婷综合
|
国产欧美精品一区二区色综合
|