999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于PI—PD控制器的四旋翼姿態控制

2014-10-29 13:20:08唐健杰王鑫
電子技術與軟件工程 2014年16期

唐健杰++王鑫

摘 要

四旋翼姿態控制通常選用PID控制器作為主控環節,但PID控制器的參數不易整定,調控結果也存在較大的系統超調量,難以獲得滿意的控制效果。因此設計了PI-PD控制器的四旋翼姿態控制方式,其中PI控制器能夠使系統快速無穩態誤差的收斂,PD控制器可以有效的抑制系統超調量。仿真結果表明:PI-PD控制器能夠很好的抑制系統超調量,縮短系統收斂時間,具有良好的控制效果。

【關鍵詞】四旋翼姿態控制 PI-PD控制器 控制效果 參數整定

近年來,隨著微型系統、微型傳感器、慣導技術以及飛行控制等技術的發展,四旋翼飛行器(以下簡稱四旋翼)引起了人們的廣泛關注。四旋翼是通過改變四個旋翼的轉速來調整其在空中的飛行姿態,包括俯仰角θ、橫滾角φ、偏航角ψ,從而控制機體水平方向上的運動,因此四旋翼的姿態控制是決定其飛行性能的關鍵所在。在工業過程控制和航空航天控制等領域中,PID控制的應用達到80 %以上,不過由于四旋翼系統的強非線性、慣性和延遲,PID控制器對四旋翼姿態的調整效果往往出現較多的系統超調量,或者調整時間較長,控制效果并不令人滿意。因此,設計一種能夠抑制系統超調量,并且保證系統快速收斂的控制器,可以提高四旋翼系統的穩定性和控制性。

1 PID控制器基本原理

PID控制器結構簡單、方便調試,廣泛應用于工業生產中。PID控制器是根據系統輸出的誤差值調節系統輸出的控制形式,包含比例控制(P)、積分控制(I)和微分控制(D),其連續PID控制的結構形式為:

(1)

其中u(t)為系統輸出,Kp、Ki、Kd分別為比例、積分和微分系數,e(t)=y(r)-y(t)為期望值與輸出量的差值,即輸出誤差。而對于數字控制系統,可將PID控制器離散化,得到離散PID的結構形式:

(2)

其中,為所有誤差值累加之和,Δe(t)=e(t)-e(t-1),等效微分運算。當期望值在相鄰的采樣周期保持不變時,y(r)=y(r-1),Δe(t)=-y(t)+y(t-1),Δe(t)即為系統輸出的變化量。

若基于PID控制器來對四旋翼的姿態進行調控,參數整定難度較大,調控效果不佳。鑒于此,本文基于文章[6]提出的PI-PD控制器,設計了基于PI-PD控制器的四旋翼姿態控制方式,用于減小系統超調量,縮短系統收斂時間,提高四旋翼在空中飛行的穩定性和控制性。

2 PI-PD控制器

PID控制器對于高階時滯系統、復雜的模糊系統以及不確定系統而言,控制效果不佳。而在PID控制器基礎上演變而來的PI-PD控制器,對于含有積分、振蕩或不穩定環節的控制對象,可以實現較好的閉環控制。PI-PD控制器其結構圖如圖1所示。

設PI控制器和PD控制器的傳遞函數為:

Gpi(S)=Kp(1+1/tis) (3)

Gpd(S)=Kf(1+tds) (4)

其中Kp、Ti分別PI控制器的比例和積分系數,Kf、Td分別為PD控制器的比例和微分系數。圖中,PI控制器仍處于主控環節上,根據期望值調節輸出量,具有決定系統收斂快慢和消除穩態誤差的作用。而PD控制器成為了反饋環節,具有抑制系統振蕩和超調量的作用,并且只與系統輸出變化量有關,與期望值無關。

為了簡化PI-PD控制器結構,將其進行結構變換,得到圖2所示的等效結構圖。

可得到主控環節PI+PD控制器為:

Gpi(S)+ Gpd(S)=Kp(1+1/tis)+Kf(1+tds) (5)

設PID控制器傳遞函數為:

Gpid*(S)=Kp*(1+1/Ti*s+Td*s) (6)

其中Kp*、Ti*、Td*分別為PID控制器的參數。于是可將式(5)整理成式(6)的類似形式:

(7)

設Kp=βKf,參數β表示Kp與Kf的關系,式(7)可變換為:

(8)

比較式(6)與式(8),可得出Kp*、Ti*、Td*與Kp、Ti、Kf、Td之間的關系表達式:

Kp=βKp*/(1+β) (9)

Kf=Kp*/(1+β) (10)

Ti=βTi*/(1+β) (11)

Td=(1+β) Td* (12)

根據式(9)和(11)可得:

(13)

根據式(10)和(12)可得:

KfTd=(1+β) Td*=Kf*Td* (14)

由此可以看出,根據PID控制器的Kp*參數以及β值可以計算出PI-PD控制器的Kp和Kf參數,β決定了Kp與Kf的分配比例。式(13)和(14)表明,PI-PD控制器的積分控制和微分控制與PID控制器的參數相同。

因此PI-PD控制器可以根據PID控制參數和β值計算得出Kp、Ti、Kf、Td參數,通過參數再次整定,能使系統在超調量較小、收斂時間較短的情況下平穩收斂,具有良好的調控效果。

3 仿真分析

通過Adams軟件建立四旋翼動力學虛擬樣機,將Adams所建模型與Matlab/Simulink進行聯合仿真,研究控制器對四旋翼姿態控制的調節效果。本文研究的四旋翼參數為:機體質量m=0.67kg,對稱電機軸距l=450mm,旋翼轉速與升力關系8000r/m=9.8N,角度初始值俯仰角θ=0°、俯仰角期望值y(r)=0°,仿真步長t=0.01s。橫滾角φ與偏航角ψ的仿真結果類同,本文不再贅述。

(1)使用PID控制器調節俯仰角θ的角度,整定一組參數Kp=18、Ki=0.2、Kd=260。而根據PID控制器參數,設定不同的β值,計算出PI-PD控制器的Kp、Ti、Kf、Td的參數,仿真結果如圖3所示。

在1s的時刻,期望值y(r)從0°變為1°,俯仰角θ在PID控制器的調節下,收斂時間約為1.8s,系統超調量約為30%。

與PID控制器的仿真結果比較,當β不同時,PI-PD控制器具有不同的調節效果。而當β=3.6時,系統超調量極小,且收斂時間與PID控制器基本相同。由此可以證明,選取合適的β值,PI-PD控制器可以有效的抑制系統超調量,提高了系統的穩定性。

(3)為了使系統收斂時間更短,根據PI-PD控制器的調節特性,重新整定參數Kp、Ki、Kd和β,仿真結果如圖4所示。

仿真結果表明:PID與PI-PD-1的收斂時間約為1.8s,而PI-PD-2的收斂時間約為0.8s,明顯快于前兩種控制效果。由此證明PI-PD控制器對于不同的整定參數,可以在系統無超調量的情況下,縮短系統收斂時間,提高了系統的控制性。

4 結束語

由于四旋翼的非線性和時滯特性,基于PID控制器的四旋翼姿態控制方式的調節效果較難滿足人們需求。本文基于PI-PD控制器,設計了新的四旋翼姿態控制方式。PI-PD控制器方式可以利用已經整定好的PID控制器的參數,根據β值計算出控制參數。仿真結果表明PI-PD控制器能夠使系統穩定收斂,且很好的抑制了系統超調量,獲得了良好的控制效果。在今后的工作中,將繼續研究PI-PD控制器的參數整定、響應時間、魯棒性等問題,獲得更好的四旋翼姿態控制效果。

作者簡介

唐健杰(1986-),男,廣西壯族自治區桂林市人。碩士研究生學歷。研究方向為智能控制與算法應用。

作者單位

桂林電子科技大學計算機科學與工程學院 廣西壯族自治區桂林市 541004endprint

摘 要

四旋翼姿態控制通常選用PID控制器作為主控環節,但PID控制器的參數不易整定,調控結果也存在較大的系統超調量,難以獲得滿意的控制效果。因此設計了PI-PD控制器的四旋翼姿態控制方式,其中PI控制器能夠使系統快速無穩態誤差的收斂,PD控制器可以有效的抑制系統超調量。仿真結果表明:PI-PD控制器能夠很好的抑制系統超調量,縮短系統收斂時間,具有良好的控制效果。

【關鍵詞】四旋翼姿態控制 PI-PD控制器 控制效果 參數整定

近年來,隨著微型系統、微型傳感器、慣導技術以及飛行控制等技術的發展,四旋翼飛行器(以下簡稱四旋翼)引起了人們的廣泛關注。四旋翼是通過改變四個旋翼的轉速來調整其在空中的飛行姿態,包括俯仰角θ、橫滾角φ、偏航角ψ,從而控制機體水平方向上的運動,因此四旋翼的姿態控制是決定其飛行性能的關鍵所在。在工業過程控制和航空航天控制等領域中,PID控制的應用達到80 %以上,不過由于四旋翼系統的強非線性、慣性和延遲,PID控制器對四旋翼姿態的調整效果往往出現較多的系統超調量,或者調整時間較長,控制效果并不令人滿意。因此,設計一種能夠抑制系統超調量,并且保證系統快速收斂的控制器,可以提高四旋翼系統的穩定性和控制性。

1 PID控制器基本原理

PID控制器結構簡單、方便調試,廣泛應用于工業生產中。PID控制器是根據系統輸出的誤差值調節系統輸出的控制形式,包含比例控制(P)、積分控制(I)和微分控制(D),其連續PID控制的結構形式為:

(1)

其中u(t)為系統輸出,Kp、Ki、Kd分別為比例、積分和微分系數,e(t)=y(r)-y(t)為期望值與輸出量的差值,即輸出誤差。而對于數字控制系統,可將PID控制器離散化,得到離散PID的結構形式:

(2)

其中,為所有誤差值累加之和,Δe(t)=e(t)-e(t-1),等效微分運算。當期望值在相鄰的采樣周期保持不變時,y(r)=y(r-1),Δe(t)=-y(t)+y(t-1),Δe(t)即為系統輸出的變化量。

若基于PID控制器來對四旋翼的姿態進行調控,參數整定難度較大,調控效果不佳。鑒于此,本文基于文章[6]提出的PI-PD控制器,設計了基于PI-PD控制器的四旋翼姿態控制方式,用于減小系統超調量,縮短系統收斂時間,提高四旋翼在空中飛行的穩定性和控制性。

2 PI-PD控制器

PID控制器對于高階時滯系統、復雜的模糊系統以及不確定系統而言,控制效果不佳。而在PID控制器基礎上演變而來的PI-PD控制器,對于含有積分、振蕩或不穩定環節的控制對象,可以實現較好的閉環控制。PI-PD控制器其結構圖如圖1所示。

設PI控制器和PD控制器的傳遞函數為:

Gpi(S)=Kp(1+1/tis) (3)

Gpd(S)=Kf(1+tds) (4)

其中Kp、Ti分別PI控制器的比例和積分系數,Kf、Td分別為PD控制器的比例和微分系數。圖中,PI控制器仍處于主控環節上,根據期望值調節輸出量,具有決定系統收斂快慢和消除穩態誤差的作用。而PD控制器成為了反饋環節,具有抑制系統振蕩和超調量的作用,并且只與系統輸出變化量有關,與期望值無關。

為了簡化PI-PD控制器結構,將其進行結構變換,得到圖2所示的等效結構圖。

可得到主控環節PI+PD控制器為:

Gpi(S)+ Gpd(S)=Kp(1+1/tis)+Kf(1+tds) (5)

設PID控制器傳遞函數為:

Gpid*(S)=Kp*(1+1/Ti*s+Td*s) (6)

其中Kp*、Ti*、Td*分別為PID控制器的參數。于是可將式(5)整理成式(6)的類似形式:

(7)

設Kp=βKf,參數β表示Kp與Kf的關系,式(7)可變換為:

(8)

比較式(6)與式(8),可得出Kp*、Ti*、Td*與Kp、Ti、Kf、Td之間的關系表達式:

Kp=βKp*/(1+β) (9)

Kf=Kp*/(1+β) (10)

Ti=βTi*/(1+β) (11)

Td=(1+β) Td* (12)

根據式(9)和(11)可得:

(13)

根據式(10)和(12)可得:

KfTd=(1+β) Td*=Kf*Td* (14)

由此可以看出,根據PID控制器的Kp*參數以及β值可以計算出PI-PD控制器的Kp和Kf參數,β決定了Kp與Kf的分配比例。式(13)和(14)表明,PI-PD控制器的積分控制和微分控制與PID控制器的參數相同。

因此PI-PD控制器可以根據PID控制參數和β值計算得出Kp、Ti、Kf、Td參數,通過參數再次整定,能使系統在超調量較小、收斂時間較短的情況下平穩收斂,具有良好的調控效果。

3 仿真分析

通過Adams軟件建立四旋翼動力學虛擬樣機,將Adams所建模型與Matlab/Simulink進行聯合仿真,研究控制器對四旋翼姿態控制的調節效果。本文研究的四旋翼參數為:機體質量m=0.67kg,對稱電機軸距l=450mm,旋翼轉速與升力關系8000r/m=9.8N,角度初始值俯仰角θ=0°、俯仰角期望值y(r)=0°,仿真步長t=0.01s。橫滾角φ與偏航角ψ的仿真結果類同,本文不再贅述。

(1)使用PID控制器調節俯仰角θ的角度,整定一組參數Kp=18、Ki=0.2、Kd=260。而根據PID控制器參數,設定不同的β值,計算出PI-PD控制器的Kp、Ti、Kf、Td的參數,仿真結果如圖3所示。

在1s的時刻,期望值y(r)從0°變為1°,俯仰角θ在PID控制器的調節下,收斂時間約為1.8s,系統超調量約為30%。

與PID控制器的仿真結果比較,當β不同時,PI-PD控制器具有不同的調節效果。而當β=3.6時,系統超調量極小,且收斂時間與PID控制器基本相同。由此可以證明,選取合適的β值,PI-PD控制器可以有效的抑制系統超調量,提高了系統的穩定性。

(3)為了使系統收斂時間更短,根據PI-PD控制器的調節特性,重新整定參數Kp、Ki、Kd和β,仿真結果如圖4所示。

仿真結果表明:PID與PI-PD-1的收斂時間約為1.8s,而PI-PD-2的收斂時間約為0.8s,明顯快于前兩種控制效果。由此證明PI-PD控制器對于不同的整定參數,可以在系統無超調量的情況下,縮短系統收斂時間,提高了系統的控制性。

4 結束語

由于四旋翼的非線性和時滯特性,基于PID控制器的四旋翼姿態控制方式的調節效果較難滿足人們需求。本文基于PI-PD控制器,設計了新的四旋翼姿態控制方式。PI-PD控制器方式可以利用已經整定好的PID控制器的參數,根據β值計算出控制參數。仿真結果表明PI-PD控制器能夠使系統穩定收斂,且很好的抑制了系統超調量,獲得了良好的控制效果。在今后的工作中,將繼續研究PI-PD控制器的參數整定、響應時間、魯棒性等問題,獲得更好的四旋翼姿態控制效果。

作者簡介

唐健杰(1986-),男,廣西壯族自治區桂林市人。碩士研究生學歷。研究方向為智能控制與算法應用。

作者單位

桂林電子科技大學計算機科學與工程學院 廣西壯族自治區桂林市 541004endprint

摘 要

四旋翼姿態控制通常選用PID控制器作為主控環節,但PID控制器的參數不易整定,調控結果也存在較大的系統超調量,難以獲得滿意的控制效果。因此設計了PI-PD控制器的四旋翼姿態控制方式,其中PI控制器能夠使系統快速無穩態誤差的收斂,PD控制器可以有效的抑制系統超調量。仿真結果表明:PI-PD控制器能夠很好的抑制系統超調量,縮短系統收斂時間,具有良好的控制效果。

【關鍵詞】四旋翼姿態控制 PI-PD控制器 控制效果 參數整定

近年來,隨著微型系統、微型傳感器、慣導技術以及飛行控制等技術的發展,四旋翼飛行器(以下簡稱四旋翼)引起了人們的廣泛關注。四旋翼是通過改變四個旋翼的轉速來調整其在空中的飛行姿態,包括俯仰角θ、橫滾角φ、偏航角ψ,從而控制機體水平方向上的運動,因此四旋翼的姿態控制是決定其飛行性能的關鍵所在。在工業過程控制和航空航天控制等領域中,PID控制的應用達到80 %以上,不過由于四旋翼系統的強非線性、慣性和延遲,PID控制器對四旋翼姿態的調整效果往往出現較多的系統超調量,或者調整時間較長,控制效果并不令人滿意。因此,設計一種能夠抑制系統超調量,并且保證系統快速收斂的控制器,可以提高四旋翼系統的穩定性和控制性。

1 PID控制器基本原理

PID控制器結構簡單、方便調試,廣泛應用于工業生產中。PID控制器是根據系統輸出的誤差值調節系統輸出的控制形式,包含比例控制(P)、積分控制(I)和微分控制(D),其連續PID控制的結構形式為:

(1)

其中u(t)為系統輸出,Kp、Ki、Kd分別為比例、積分和微分系數,e(t)=y(r)-y(t)為期望值與輸出量的差值,即輸出誤差。而對于數字控制系統,可將PID控制器離散化,得到離散PID的結構形式:

(2)

其中,為所有誤差值累加之和,Δe(t)=e(t)-e(t-1),等效微分運算。當期望值在相鄰的采樣周期保持不變時,y(r)=y(r-1),Δe(t)=-y(t)+y(t-1),Δe(t)即為系統輸出的變化量。

若基于PID控制器來對四旋翼的姿態進行調控,參數整定難度較大,調控效果不佳。鑒于此,本文基于文章[6]提出的PI-PD控制器,設計了基于PI-PD控制器的四旋翼姿態控制方式,用于減小系統超調量,縮短系統收斂時間,提高四旋翼在空中飛行的穩定性和控制性。

2 PI-PD控制器

PID控制器對于高階時滯系統、復雜的模糊系統以及不確定系統而言,控制效果不佳。而在PID控制器基礎上演變而來的PI-PD控制器,對于含有積分、振蕩或不穩定環節的控制對象,可以實現較好的閉環控制。PI-PD控制器其結構圖如圖1所示。

設PI控制器和PD控制器的傳遞函數為:

Gpi(S)=Kp(1+1/tis) (3)

Gpd(S)=Kf(1+tds) (4)

其中Kp、Ti分別PI控制器的比例和積分系數,Kf、Td分別為PD控制器的比例和微分系數。圖中,PI控制器仍處于主控環節上,根據期望值調節輸出量,具有決定系統收斂快慢和消除穩態誤差的作用。而PD控制器成為了反饋環節,具有抑制系統振蕩和超調量的作用,并且只與系統輸出變化量有關,與期望值無關。

為了簡化PI-PD控制器結構,將其進行結構變換,得到圖2所示的等效結構圖。

可得到主控環節PI+PD控制器為:

Gpi(S)+ Gpd(S)=Kp(1+1/tis)+Kf(1+tds) (5)

設PID控制器傳遞函數為:

Gpid*(S)=Kp*(1+1/Ti*s+Td*s) (6)

其中Kp*、Ti*、Td*分別為PID控制器的參數。于是可將式(5)整理成式(6)的類似形式:

(7)

設Kp=βKf,參數β表示Kp與Kf的關系,式(7)可變換為:

(8)

比較式(6)與式(8),可得出Kp*、Ti*、Td*與Kp、Ti、Kf、Td之間的關系表達式:

Kp=βKp*/(1+β) (9)

Kf=Kp*/(1+β) (10)

Ti=βTi*/(1+β) (11)

Td=(1+β) Td* (12)

根據式(9)和(11)可得:

(13)

根據式(10)和(12)可得:

KfTd=(1+β) Td*=Kf*Td* (14)

由此可以看出,根據PID控制器的Kp*參數以及β值可以計算出PI-PD控制器的Kp和Kf參數,β決定了Kp與Kf的分配比例。式(13)和(14)表明,PI-PD控制器的積分控制和微分控制與PID控制器的參數相同。

因此PI-PD控制器可以根據PID控制參數和β值計算得出Kp、Ti、Kf、Td參數,通過參數再次整定,能使系統在超調量較小、收斂時間較短的情況下平穩收斂,具有良好的調控效果。

3 仿真分析

通過Adams軟件建立四旋翼動力學虛擬樣機,將Adams所建模型與Matlab/Simulink進行聯合仿真,研究控制器對四旋翼姿態控制的調節效果。本文研究的四旋翼參數為:機體質量m=0.67kg,對稱電機軸距l=450mm,旋翼轉速與升力關系8000r/m=9.8N,角度初始值俯仰角θ=0°、俯仰角期望值y(r)=0°,仿真步長t=0.01s。橫滾角φ與偏航角ψ的仿真結果類同,本文不再贅述。

(1)使用PID控制器調節俯仰角θ的角度,整定一組參數Kp=18、Ki=0.2、Kd=260。而根據PID控制器參數,設定不同的β值,計算出PI-PD控制器的Kp、Ti、Kf、Td的參數,仿真結果如圖3所示。

在1s的時刻,期望值y(r)從0°變為1°,俯仰角θ在PID控制器的調節下,收斂時間約為1.8s,系統超調量約為30%。

與PID控制器的仿真結果比較,當β不同時,PI-PD控制器具有不同的調節效果。而當β=3.6時,系統超調量極小,且收斂時間與PID控制器基本相同。由此可以證明,選取合適的β值,PI-PD控制器可以有效的抑制系統超調量,提高了系統的穩定性。

(3)為了使系統收斂時間更短,根據PI-PD控制器的調節特性,重新整定參數Kp、Ki、Kd和β,仿真結果如圖4所示。

仿真結果表明:PID與PI-PD-1的收斂時間約為1.8s,而PI-PD-2的收斂時間約為0.8s,明顯快于前兩種控制效果。由此證明PI-PD控制器對于不同的整定參數,可以在系統無超調量的情況下,縮短系統收斂時間,提高了系統的控制性。

4 結束語

由于四旋翼的非線性和時滯特性,基于PID控制器的四旋翼姿態控制方式的調節效果較難滿足人們需求。本文基于PI-PD控制器,設計了新的四旋翼姿態控制方式。PI-PD控制器方式可以利用已經整定好的PID控制器的參數,根據β值計算出控制參數。仿真結果表明PI-PD控制器能夠使系統穩定收斂,且很好的抑制了系統超調量,獲得了良好的控制效果。在今后的工作中,將繼續研究PI-PD控制器的參數整定、響應時間、魯棒性等問題,獲得更好的四旋翼姿態控制效果。

作者簡介

唐健杰(1986-),男,廣西壯族自治區桂林市人。碩士研究生學歷。研究方向為智能控制與算法應用。

作者單位

桂林電子科技大學計算機科學與工程學院 廣西壯族自治區桂林市 541004endprint

主站蜘蛛池模板: 国产无码网站在线观看| 亚洲色图狠狠干| 久久人人妻人人爽人人卡片av| 亚洲综合中文字幕国产精品欧美 | 欧美一级黄色影院| 99热国产这里只有精品无卡顿"| 一级爆乳无码av| 六月婷婷激情综合| 欧美精品二区| 国产成人精品男人的天堂| 华人在线亚洲欧美精品| 中国国语毛片免费观看视频| 国产麻豆精品在线观看| 国产理论最新国产精品视频| 亚洲性影院| 亚洲AV无码乱码在线观看代蜜桃| 美女一级毛片无遮挡内谢| 国产裸舞福利在线视频合集| 亚洲国产理论片在线播放| 国产无码高清视频不卡| 色丁丁毛片在线观看| 亚洲一区二区日韩欧美gif| 香港一级毛片免费看| 色欲综合久久中文字幕网| 91丝袜在线观看| 制服丝袜一区| 中文字幕日韩久久综合影院| 欧美精品在线视频观看| 国产免费网址| 18禁色诱爆乳网站| 99免费在线观看视频| 在线精品亚洲国产| 国产精品无码久久久久久| 国产91丝袜在线观看| 国产精品亚洲专区一区| 中文字幕免费视频| 日韩中文字幕亚洲无线码| 亚洲成a人片| 国产在线专区| 最新日韩AV网址在线观看| 狠狠色丁香婷婷| a天堂视频| 综合网天天| 国产黄网永久免费| 亚洲精品波多野结衣| 久久精品中文无码资源站| 欧美黑人欧美精品刺激| 蜜桃臀无码内射一区二区三区| 日本欧美视频在线观看| 3344在线观看无码| 日韩精品一区二区三区中文无码| 亚洲一区网站| 91av成人日本不卡三区| 久久黄色视频影| 性喷潮久久久久久久久| 国产福利小视频在线播放观看| 国产成人久久777777| 无码又爽又刺激的高潮视频| 欧美亚洲一区二区三区导航| 欧美中文字幕第一页线路一| swag国产精品| 国产乱肥老妇精品视频| 国产美女自慰在线观看| 国产精品网址你懂的| 亚洲免费毛片| 午夜少妇精品视频小电影| 久久96热在精品国产高清| 精品国产中文一级毛片在线看| 日本免费福利视频| 在线观看免费人成视频色快速| 午夜丁香婷婷| 国产成熟女人性满足视频| 久久动漫精品| 国产aⅴ无码专区亚洲av综合网| 不卡无码h在线观看| 欧美色综合久久| 精品自窥自偷在线看| 亚洲国产精品日韩av专区| 日韩国产另类| 四虎影院国产| 国产99视频精品免费观看9e| 美女裸体18禁网站|