• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于聲輻射模態(tài)的聲場(chǎng)重建中的測(cè)點(diǎn)優(yōu)化方法

    2017-03-09 08:29:28蘇俊博朱海潮毛榮富蘇常偉
    振動(dòng)與沖擊 2017年3期
    關(guān)鍵詞:階數(shù)聲壓聲場(chǎng)

    蘇俊博, 朱海潮, 毛榮富, 郭 亮, 蘇常偉

    (海軍工程大學(xué) 船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室,武漢 430033)

    基于聲輻射模態(tài)的聲場(chǎng)重建中的測(cè)點(diǎn)優(yōu)化方法

    蘇俊博, 朱海潮, 毛榮富, 郭 亮, 蘇常偉

    (海軍工程大學(xué) 船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室,武漢 430033)

    少測(cè)點(diǎn)條件下,利用聲輻射模態(tài)理論重建聲場(chǎng)時(shí),測(cè)點(diǎn)布置方式是決定聲場(chǎng)重建精度的一個(gè)關(guān)鍵因素。為求得最佳測(cè)點(diǎn)布置方式,提出了一種基于聲輻射模態(tài)的最優(yōu)測(cè)點(diǎn)選擇方法,即基于聲輻射模態(tài)矩陣的奇異值分解,采用循環(huán)迭代的方式,逐次去除對(duì)其最小奇異值最敏感的測(cè)點(diǎn),從而得到了一組使重建方程條件數(shù)最小的測(cè)點(diǎn)。實(shí)驗(yàn)結(jié)果表明,利用文中提出的最優(yōu)測(cè)點(diǎn)選擇方法布置測(cè)點(diǎn),能夠?qū)β晥?chǎng)進(jìn)行有效的重建,重建效果優(yōu)于均勻布置方式,顯著提高了聲場(chǎng)重建精度。

    聲輻射模態(tài);聲場(chǎng)重建;測(cè)點(diǎn)優(yōu)化

    近場(chǎng)聲全息技術(shù)自20世紀(jì)80年代由MAYNARD等[1]提出后,許多學(xué)者進(jìn)行了研究,目前已經(jīng)在噪聲源識(shí)別與定位、聲場(chǎng)可視化等方面取得了良好的效果。但近場(chǎng)聲全息技術(shù)要求在全息面上布置大量測(cè)點(diǎn),特別是應(yīng)用于大型結(jié)構(gòu)時(shí),這限制了近場(chǎng)聲全息技術(shù)在工程實(shí)際中的應(yīng)用。因此如何減少測(cè)點(diǎn)數(shù)目是一個(gè)需要解決的重要問題。HELS(Helmholtz Equation-Least Squares)方法[2-4]是一種所需測(cè)點(diǎn)數(shù)目遠(yuǎn)少于近場(chǎng)聲全息技術(shù)的聲場(chǎng)重建方法,其原理是是將輻射聲場(chǎng)表示為Helmholtz方程一組特解的線性組合,利用聲場(chǎng)測(cè)量數(shù)據(jù)基于最小二乘法求解組合系數(shù),然后再求解聲場(chǎng)中的場(chǎng)點(diǎn)聲壓。HELS方法對(duì)類似球形聲源的聲場(chǎng)能夠?qū)崿F(xiàn)較為精確的重建,但對(duì)長(zhǎng)寬比較大的結(jié)構(gòu)聲源聲場(chǎng)難以進(jìn)行有效的重建。聲輻射模態(tài)理論的提出為解決少測(cè)點(diǎn)條件下任意結(jié)構(gòu)聲源的聲場(chǎng)重建問題提供了新的途徑[5-6]。

    聲輻射模態(tài)理論是BORGIOTTI等[7-8]學(xué)者在九十年代初提出,目前已在結(jié)構(gòu)聲輻射主動(dòng)控制以及輻射聲功率的計(jì)算方面取得了許多成果[9-13]。聲輻射模態(tài)就是輻射體表面一種可能的輻射形式,任何表面法向速度或聲壓都可通過這些聲輻射模態(tài)的線性組合來表示,利用聲輻射模態(tài)進(jìn)行聲場(chǎng)重建的整個(gè)過程可分為求解模態(tài)展開系數(shù),聲場(chǎng)聲壓求解兩個(gè)步驟。由此可見,聲輻射模態(tài)理論與HELS方法具有相似性,但是由于聲輻射模態(tài)中包含了結(jié)構(gòu)表面的幾何形狀信息,因此聲輻射模態(tài)理論適用于任意結(jié)構(gòu)輻射聲場(chǎng)的重建。而在少測(cè)點(diǎn)條件下利用聲輻射模態(tài)進(jìn)行聲場(chǎng)重建時(shí),測(cè)點(diǎn)布置方式是一個(gè)決定聲場(chǎng)重建精度的關(guān)鍵問題,如果測(cè)點(diǎn)的位置分布不當(dāng)會(huì)使各測(cè)點(diǎn)之間有較強(qiáng)的相關(guān)性,在測(cè)點(diǎn)較少的情況下就無法采集到足夠的有效數(shù)據(jù),從而使求解精度降低。此外,在已知測(cè)點(diǎn)布置方式的情況下,聲輻射模態(tài)階數(shù)的選擇也對(duì)求解精度具有重要影響,這是由于利用已知測(cè)點(diǎn)求解模態(tài)展開系數(shù)是一個(gè)數(shù)學(xué)上的逆問題,模態(tài)階數(shù)過少不足以表示,模態(tài)階數(shù)過多又會(huì)使受到噪聲污染的高階模態(tài)參與到求解過程中,從而使求解精度降低。

    宋艷華等[14]研究了利用聲輻射模態(tài)進(jìn)行聲源識(shí)別時(shí)測(cè)點(diǎn)位置對(duì)識(shí)別效果的影響,并得出結(jié)論:測(cè)點(diǎn)位置分布不均勻時(shí),識(shí)別效果較好。但沒有對(duì)測(cè)點(diǎn)分布方式作進(jìn)一步研究。KAMMER[15]提出了一種有效獨(dú)立方法(Effective Independence Method, EFI)。KIM[16]將該方法用于基于邊界元方法的聲場(chǎng)重建中,在聲場(chǎng)中選取了一組線性獨(dú)立的測(cè)點(diǎn),但該方法的局限性在于選取的測(cè)點(diǎn)數(shù)量不能少于結(jié)構(gòu)表面離散的單元數(shù)。針對(duì)EFI方法的這一局限性,ZHANG[17]提出了一種基于奇異值分解的循環(huán)迭代方法。為解決利用聲輻射模態(tài)重建聲場(chǎng)時(shí)所面臨的測(cè)點(diǎn)布置問題,本文提出了一種基于聲輻射模態(tài)的最優(yōu)測(cè)點(diǎn)選擇方法。針對(duì)聲輻射模態(tài)階數(shù)選取問題,將一種基于測(cè)量的迭代方法[18]引入到重建過程中,確定了重建過程中所需選取的聲輻射模態(tài)階數(shù)。文章最后開展了雙層圓柱殼體實(shí)驗(yàn)研究,通過與測(cè)點(diǎn)均勻布置方式下聲場(chǎng)重建效果的比較,驗(yàn)證了本文提出的基于聲輻射模態(tài)的最優(yōu)測(cè)點(diǎn)選擇方法的有效性。

    1 基于聲輻射模態(tài)的聲場(chǎng)重建理論

    將輻射體表面離散成一系列等面積的輻射單元,由文獻(xiàn)[7]可知,輻射聲功率可表示為

    (1)

    R=ΦΛΦH

    (2)

    式中:Φ為N×N維矩陣,Φ的列向量φi(i=1,2,…,N)即為輻射體的聲輻射模態(tài)向量。

    在得到了輻射體的聲輻射模態(tài)以后,輻射體表面的速度向量可表示成如下形式

    v=Φc

    (3)

    式中:c為輻射模態(tài)展開系數(shù)。假設(shè)將輻射體表面離散成N個(gè)輻射單元,并計(jì)算得到了N階聲輻射模態(tài),由于式(3)具有良好的收斂性[14],表面速度可用截?cái)嗟穆曒椛淠B(tài)表示為

    v≈Φ(N×N1)c(N1×1)

    (4)

    式中:矩陣Φ(N×N1)表示只取了前N1階聲輻射模態(tài),當(dāng)輻射體表面速度向量中有N2個(gè)元素是已知的,可以得到由這N2個(gè)元素組成的方程組

    v′=Φ′(N2×N1)c(N1×1)

    (5)

    由式(4)和式(5)可以求得表面速度向量

    (6)

    利用聲輻射模態(tài)由少量測(cè)點(diǎn)聲壓值重建輻射體表面聲壓具有相同的過程。在得到了輻射體表面的聲壓或振速值后,就可利用邊界元等方法方便地計(jì)算出整個(gè)聲場(chǎng)分布。

    2 最優(yōu)測(cè)點(diǎn)選擇方法

    針對(duì)測(cè)點(diǎn)選擇問題,文獻(xiàn)[17]用一種基于奇異值分解的循環(huán)迭代方法,基于使輻射體表面離散點(diǎn)與聲場(chǎng)中測(cè)點(diǎn)之間的傳遞函數(shù)矩陣條件數(shù)最小的準(zhǔn)則,求得了聲場(chǎng)中的一組測(cè)點(diǎn)。本文將其思想應(yīng)用于基于聲輻射模態(tài)的聲場(chǎng)重建中,提出一種基于聲輻射模態(tài)的最優(yōu)測(cè)點(diǎn)選擇方法。假定輻射體表面候選測(cè)點(diǎn)數(shù)量為N,從這N個(gè)候選測(cè)點(diǎn)中選取N2個(gè)作為實(shí)際布置傳感器的位置。

    假設(shè)輻射體表面聲壓向量為p0(N×1),可用聲輻射模態(tài)表示為

    p0(N×1)=Φ(N×N)c(N×1)

    (7)

    對(duì)Φ(N×N)進(jìn)行奇異值分解

    Φ(N×N)=USVH

    (8)

    式中:S為奇異值矩陣,S對(duì)角線元素為聲輻射模態(tài)矩陣Φ(N×N)的奇異值。利用式(7)、式(8)求得的模態(tài)展開系數(shù)重建輻射體表面聲壓pr(N×1)

    pr(N×1)=Φ(N×N)c(N×1)=ΦΦ+p0=

    USVHVS+UHp0(N×1)=USS+UHp0(N×1)

    (9)

    將找到的對(duì)最小奇異值最敏感的測(cè)點(diǎn)從p0(N×1)中去除,并將模態(tài)矩陣Φ(N×N)中與該測(cè)點(diǎn)對(duì)應(yīng)的行去除,并此時(shí)式(7)變?yōu)?/p>

    p0[(N-1)×1)]=Φ[(N-1)×N)c(N×1)

    (10)

    利用式(10)求解模態(tài)展開系數(shù)時(shí),式(10)為一欠定方程,為使方程有解,去除模態(tài)矩陣Φ[(N-1)×N]中的最高階模態(tài)向量,也即模態(tài)矩陣的最后一列,得到新的模態(tài)矩陣Φ1[(N-1)×(N-1)],則式(10)變?yōu)?/p>

    p01[(N-1)×1]=Φ1[(N-1)×(N-1)]c1[(N-1)×1]

    (11)

    然后以式(11)代替式(7)。至此,由式(7)得到式(11)并以式(11)代替式(7)的過程完成了一次循環(huán)迭代,以相同過程迭代N-N2+1次,得到如下方程

    p0(N-N2)(N2×1)=Φ(N-N2)(N2×N2)c(N-N2)(N2×1)

    (12)

    則p0(N-N2)(N2×1)就是我們需要的N2個(gè)測(cè)點(diǎn),式(12)即為由此N2個(gè)測(cè)點(diǎn)組成的重建方程。

    由以上推導(dǎo)過程可知,該方法是基于模態(tài)矩陣的奇異值分解,采用了一種循環(huán)迭代的方式,逐次去除對(duì)最小奇異值最敏感的測(cè)點(diǎn),直至達(dá)到要求的測(cè)點(diǎn)數(shù)目。經(jīng)過此循環(huán)迭代過程直接得到了求解模態(tài)展開系數(shù)的方程,確保了重建方程系數(shù)矩陣的條件數(shù)最小。整個(gè)測(cè)點(diǎn)選擇流程如圖1所示。

    上述介紹的最優(yōu)測(cè)點(diǎn)選擇方法同樣適用于對(duì)某一頻段范圍內(nèi)的噪聲進(jìn)行分析,僅需要在每一次的循環(huán)迭代過程中增加一個(gè)對(duì)角線元素求和的步驟,具體方法如下:根據(jù)式(9)可知,在每一次的循環(huán)迭代過程中,均有下式出現(xiàn):

    pr=Φc=ΦΦ+p0=

    USVHVS+UHp0=USS+UHp0

    (13)

    式(13)是由表面聲壓重建表面聲壓,因此矩陣T=USS+UH為單位矩陣,T的對(duì)角線元素與輻射體表面離散測(cè)點(diǎn)是一一對(duì)應(yīng)的關(guān)系。為得到對(duì)最小奇異值最為敏感的測(cè)點(diǎn),只需將最小奇異值置零。假設(shè)分析的對(duì)象為某個(gè)頻段范圍內(nèi)的噪聲,則在此頻段內(nèi)均勻選取一組離散的頻率點(diǎn)f1,f2, …fn,根據(jù)本文的最優(yōu)測(cè)點(diǎn)選擇方法得到一組矩陣T1,T2, …Tn,將每個(gè)矩陣的對(duì)應(yīng)的S和S+中的最小奇異值置零,并將對(duì)應(yīng)位置的對(duì)角線元素相加,得到一組新的對(duì)角線元素,則將這組新的對(duì)角線元素的最小值所對(duì)應(yīng)的測(cè)點(diǎn)確定為此次迭代過程所要去除的測(cè)點(diǎn)。

    圖1 最優(yōu)測(cè)點(diǎn)選擇流程圖Fig.1 The process ofdetermining the optimum measurement points

    3 實(shí)驗(yàn)研究

    為了驗(yàn)證本文提出的最優(yōu)測(cè)點(diǎn)選擇方法的有效性,在消聲水池中進(jìn)行了雙層圓柱殼體實(shí)驗(yàn)。雙層圓柱殼體高度為2 m,半徑0.5 m,內(nèi)殼厚度0.008 m,外殼厚度0.002 m,近場(chǎng)水聽器陣列布置固定于旋轉(zhuǎn)裝置上,通過旋轉(zhuǎn)近場(chǎng)水聽器陣列達(dá)到測(cè)量四周近場(chǎng)水聲的目的,近場(chǎng)水聽器線陣由11個(gè)水聽器組成,水聽器間距0.22 m,線陣貼近外殼表面,距離外殼表面0.14 m,圓柱殼體內(nèi)部裝有激振器,采用頻率為168 Hz的單頻信號(hào)對(duì)圓柱殼體進(jìn)行激勵(lì)。圖2為實(shí)驗(yàn)裝置簡(jiǎn)圖,圖3為實(shí)驗(yàn)現(xiàn)場(chǎng)圖。

    圖2 實(shí)驗(yàn)裝置示意圖Fig.2 The sketch of experimental device

    圖3 實(shí)驗(yàn)現(xiàn)場(chǎng)Fig.3 The scene of experiment

    當(dāng)近場(chǎng)水聽器陣列以20°為步進(jìn)角度旋轉(zhuǎn)一周時(shí),在圓柱全息面測(cè)得198個(gè)測(cè)點(diǎn)聲壓值,取其中66個(gè)點(diǎn)對(duì)這198測(cè)點(diǎn)進(jìn)行重建,用本文提出的最優(yōu)測(cè)點(diǎn)選擇方法求得這66個(gè)測(cè)點(diǎn)的最優(yōu)分布方式,為便于觀察測(cè)點(diǎn)的分布情況,將圓柱面沿圓周方向展開成平面形式,如圖4所示。

    選取168 Hz為分析頻率,當(dāng)激振器工作時(shí)測(cè)得圓柱全息面聲壓幅值如圖5所示。根據(jù)本文采用的最佳聲輻射模態(tài)階數(shù)選取方法,另外增加一組測(cè)點(diǎn),根據(jù)圓柱的形狀特點(diǎn),另外增加的測(cè)點(diǎn)可選中間的一行和一列,其示意圖如圖6所示。正則化方法對(duì)于逆問題的求解具有非常重要的作用,本文在對(duì)圓柱全息面聲壓值進(jìn)行重建時(shí)均采用了Tikhonov正則化方法。

    圖4 最優(yōu)測(cè)點(diǎn)分布圖Fig.4 The optimum distribution of measurement points

    圖5 圓柱全息面聲壓測(cè)量值Fig.5 Measurement value of the pressure on cylindrical holographic surface

    圖6 增加測(cè)點(diǎn)示意圖Fig.6 The sketch of the increased measurement points

    重建中選取的聲輻射模態(tài)階數(shù)與增加測(cè)點(diǎn)的重建誤差之間的關(guān)系如圖7(a)所示,選取階數(shù)與所有測(cè)點(diǎn)的重建誤差之間的關(guān)系如圖7(b)所示。從圖7(a)與圖7(b)可觀察到二者具有相同的變化趨勢(shì),重建誤差隨著選取的輻射模態(tài)階數(shù)的增加先減小后增大,圖7(a)在28階處重建誤差取得最小值14.57%,圖7(b)在29階處重建誤差取得最小值13.01%,這說明將這種基于測(cè)量的迭代方法引入到基于聲輻射模態(tài)的聲場(chǎng)重建中來確定聲輻射模態(tài)階數(shù)是有效的。此處確定的重建階數(shù)為28,其重建誤差為13.2%,重建效果如圖8所示。

    (a) 增加測(cè)點(diǎn) (b)所有測(cè)點(diǎn)圖7 階數(shù)選取與重建誤差關(guān)系Fig.7 Relationship between number of modes and reconstruction error

    圖8 圓柱全息面聲壓重建值Fig.8 Reconstruction value of the pressure on cylindrical holographic surface

    為驗(yàn)證本文中最優(yōu)測(cè)點(diǎn)選擇方法的有效性,在其他條件不變的情況下,采取下面兩種均勻布置的方式來重建圓柱全息面聲壓值,布置測(cè)點(diǎn)數(shù)目仍為66,其布置方式如圖9、圖10所示。

    圖9 均勻布置方式1Fig.9 The firsteven arrangement

    圖10 均勻布置方式2Fig.10 The second even arrangement

    兩種布置方式條件下選取的階數(shù)與所有測(cè)點(diǎn)的重建誤差之間的關(guān)系如圖11所示,從圖11中可知,均勻布置方式1條件下,重建誤差在46階處取得最小值64.32%,均勻布置方式2條件下,重建誤差在34階處取得最小值23.22%,重建效果如圖12所示。三種布置方式下圓柱全息面聲壓值重建誤差如表1所示。由圖11、12及表1可以看出,測(cè)點(diǎn)分布方式對(duì)最終的重建結(jié)果具有非常重要的影響,當(dāng)測(cè)點(diǎn)按均勻布置方式1布置時(shí),重建誤差達(dá)到64.32%,重建已經(jīng)失效;當(dāng)測(cè)點(diǎn)按本文提出的方法進(jìn)行布置時(shí),重建誤差僅為13.20%,優(yōu)于另外兩種均勻布置方式。對(duì)比結(jié)果表明,利用本文提出的最優(yōu)測(cè)點(diǎn)選擇方法布置測(cè)點(diǎn),能夠?qū)?fù)雜結(jié)構(gòu)聲場(chǎng)進(jìn)行有效的重建,提高了聲場(chǎng)重建精度。

    (a) 均勻布置方式1 (b) 均勻布置方式2圖11 階數(shù)選取與重建誤差關(guān)系Fig.11 Relationship between number of modes and reconstruction error

    (a) 均勻布置方式1 (b) 均勻布置方式2圖12 圓柱全息面聲壓重建值Fig.12 Reconstruction value of the pressure on cylindrical holographic surface

    最優(yōu)布置方式均勻布置方式1均勻布置方式213.20%64.32%23.22%

    上述過程利用全息測(cè)量面部分測(cè)點(diǎn)聲壓值重建了圓柱全息面所有測(cè)點(diǎn)的聲壓值,重建誤差在工程實(shí)際可接受范圍內(nèi),取得了較為理想的效果。求得圓柱全息面聲壓值后進(jìn)而可以求解出聲場(chǎng)中任意一點(diǎn)的聲壓值,此計(jì)算過程本文不再敘述。

    4 結(jié) 論

    研究結(jié)果表明,利用聲輻射模態(tài)理論由輻射體表面少量測(cè)點(diǎn)重建聲場(chǎng)時(shí),測(cè)點(diǎn)布置方式是決定聲場(chǎng)重建精度的一個(gè)關(guān)鍵因素?;谥亟ǚ匠虠l件數(shù)最小準(zhǔn)則,本文提出了一種基于聲輻射模態(tài)的最優(yōu)測(cè)點(diǎn)選擇方法,利用該方法布置測(cè)點(diǎn),能夠?qū)β晥?chǎng)進(jìn)行有效的重建,重建效果優(yōu)于均勻布置方式,聲場(chǎng)重建精度明顯提高,并為實(shí)驗(yàn)驗(yàn)證。

    [ 1 ] MAYNARD J D, WILLIAMS E G, LEE Y. Near field acoustic holography Ⅰ. theory of generalized holography and the development of NAH[J]. J. Acoust. Soc. Am. , 1985, 78(4): 1395-1413.

    [ 2 ] WANG Zhaoxi, WU S F. Helmholtz equation-least-squares method for reconstructing the acoustic pressure field[J]. J. Acoust. Soc. Am. , 1997, 102(4): 2020-2032.

    [ 3 ] WU S F. On reconstruction of acoustic pressure fields using the Helmholtz equation least squares method[J]. J. Acoust. Soc. Am. , 2000, 107(5): 2511-2522.

    [ 4 ] 張海濱, 萬泉, 蔣偉康. HELS法在循環(huán)平穩(wěn)聲場(chǎng)全息重建中的理論與試驗(yàn)研究[J]. 物理學(xué)報(bào), 2009, 58(1):333-340. ZHANG Haibin, WAN Quan, JIANG Weikang. Theorical and experimental research on the HELS method in the cyclo stationary acoustic field[J]. Acta Physica Sinica, 2009, 58(1):333-340.

    [ 5 ] 姜哲. 聲輻射問題中的模態(tài)分析:Ⅲ.聲場(chǎng)重構(gòu)[J]. 聲學(xué)學(xué)報(bào), 2005, 30(3): 242-248. JIANG Zhe. A modal analysis for the acoustic radiation problem: III. Reconstruction of acoustic fields [J].Acta Acustica, 2005, 30(3):242-248.

    [ 6 ] 聶永發(fā), 朱海潮. 利用源強(qiáng)密度聲輻射模態(tài)重建聲場(chǎng)[J]. 物理學(xué)報(bào), 2014, 63(10):104303. NIE Yongfa, ZHU Haichao.Acoustic field reconstruction using source strength density acoustic radiation modes[J]. Acta Physica Sinica, 2014,63(10):104303.

    [ 7 ] BORGIOTTI G V. The power radiated by a vibrating body in an acoustic fluid and its determination from boundary measurements[J]. J. Acoust. Soc. Am. , 1990, 88(4): 1884-1893.

    [ 8 ] CUNEFARE K A. The design sensitivity and control of acoustic power radiated by three-dimensional structures[D]. Ph.D., The Pennsylvania State University, 1990.

    [ 9 ] ELLIOTT S J, JOHNSON M E. Radiation modes and the activecontrol of sound power[J]. J. Acoust. Soc. Am. , 1993, 94(4): 2194-2204.

    [10] PETERS H, KESSISSOGLOU N. Enforcing reciprocity in numerical analysis of acoustic radiation modes and sound power evaluation[J]. Journal of Computational Acoustics, 2012, 20(3): 5-33.

    [11] WU Haijun,JIANG Weikang, ZHANG Yilin. A method to compute the radiated sound power based on mapped acoustic radiation modes[J]. Journal of Computational Acoustics, 2014, 135(2): 679-692.

    [12] 杜向華,朱海潮,毛榮富. 利用聲輻射模態(tài)進(jìn)行聲功率的靈敏度分析[J]. 振動(dòng)與沖擊,2011,30(11):183-185. DU Xianghua, ZHU Haichao, MAO Rongfu. Acoustic sensitivity analysis of sound power using sound radiation modes[J]. Journal of Vibration and Shock, 2011, 30(11):183-185.

    [13] 袁國(guó)清, 姜哲. 基于聲輻射模態(tài)模型求解聲功率靈敏度[J]. 振動(dòng)與沖擊,2009,28(8):109-112. YUAN Guoqing, JIANG Zhe. Method for solving sound power sensitivity based on the theory of acoustic radiationmodes[J]. Journal of Vibration and Shock, 2009, 28(8):109-112.

    [14] 宋艷華, 姜哲. 測(cè)量點(diǎn)位置對(duì)聲源識(shí)別的影響[J]. 振動(dòng)與沖擊,2008,27(3):35-37. SONG Yanhua, JIANG Zhe. Effect of measurement points distribution on sound souce identification [J]. Journal of Vibration and Shock, 2008, 27(3):35-37.

    [15] KAMMER D C. Sensor placement for on-orbit model identification and correlation of large space structures[J]. J. Guid. Control Dyn. 1991,14: 251-259.

    [16] KIM B K, IH J G. On the reconstruction of the vibro-acoustic field over the surface enclosing an interior space using the boundary element method[J]. J. Acoust. Soc. Am. ,1996, 100(5): 3003-3016.

    [17] ZHANG Zhidong, VLAHOPOULOS N. A computational acoustic field reconstruction process based on an indirect boundary element formulation[J]. J. Acoust. Soc. Am. 2000,108 (5):2167-2178.

    [18] 雷宣揚(yáng),陳進(jìn),張桂才,等. 基于Helmholtz方程最小二乘法的聲場(chǎng)重構(gòu)[J]. 上海交通大學(xué)學(xué)報(bào),2006,40(1):129-132. LEI Xuanyang, CHEN Jin, ZHANG Guicai, et al.The reconstruction of sound field using helmholtz equation-least squares method[J]. Journal of Shang Hai Jiao Tong University, 2006, 40(1):129-132.

    Optimization of measurement points in reconstruction of acoustic fieldbased on acoustic radiation modes

    SU Junbo, ZHU Haihao, MAO Rongfu, GUO Liang, SU Changwei

    (National Key Laboratory on Ship Vibration & Noise, Naval University of Engineering, Wuhan 430033, China)

    For an acoustic field reconstruction using the acoustic radiation mode theory, the distribution of measurement points is a key factor, especially, when the number of measurement points is small. In order to determine the optimal arrangement of measurement points, a method based on acoustic radiation modes was proposed. A loop iteration process was adopted in this proposed method. Singular value decomposition (SVD) of the matrix of acoustic radiation modes was employed in each loop iteration, and the measurement point that was the most sensitive to the minimum singular value was removed. Finally, those points that minimize the condition number of reconstruction equations were obtained. Furthermore, compared with the situation when measurement points distribute evenly, the test results showed that the acoustic field of sound sources can be reconstructed effectively and the precision of reconstruction is much better when measurement points distribute according to their optimal arrangement obtained with the proposed method.

    acoustic radiation modes; acoustic field reconstruction; optimal arrangement of measurement points

    國(guó)家自然科學(xué)基金(51305452)

    2015-09-28 修改稿收到日期:2016-01-14

    蘇俊博 男,博士生,1985年生

    朱海潮 男,教授,博士生導(dǎo)師,1963年生

    O423;O429

    A

    10.13465/j.cnki.jvs.2017.03.023

    猜你喜歡
    階數(shù)聲壓聲場(chǎng)
    基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
    關(guān)于無窮小階數(shù)的幾點(diǎn)注記
    確定有限級(jí)數(shù)解的階數(shù)上界的一種n階展開方法
    基于BIM的鐵路車站聲場(chǎng)仿真分析研究
    車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
    汽車工程(2018年12期)2019-01-29 06:46:36
    探尋360°全聲場(chǎng)發(fā)聲門道
    基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
    一種新的多址信道有效階數(shù)估計(jì)算法*
    關(guān)于動(dòng)態(tài)電路階數(shù)的討論
    板結(jié)構(gòu)-聲場(chǎng)耦合分析的FE-LSPIM/FE法
    丰满乱子伦码专区| 你懂的网址亚洲精品在线观看| 免费高清在线观看视频在线观看| 国产高清国产精品国产三级 | 国产一区有黄有色的免费视频| 亚洲最大成人中文| 看非洲黑人一级黄片| 久久人人爽av亚洲精品天堂 | 亚洲精品一二三| 国产午夜福利久久久久久| 国产伦精品一区二区三区视频9| 免费观看性生交大片5| 99热这里只有是精品在线观看| 赤兔流量卡办理| 亚洲欧美一区二区三区黑人 | 3wmmmm亚洲av在线观看| 午夜福利视频精品| 久久久久久久久久久丰满| 99久久九九国产精品国产免费| 欧美亚洲 丝袜 人妻 在线| 久久久久久久午夜电影| 51国产日韩欧美| 在线观看一区二区三区激情| 成人毛片60女人毛片免费| 亚州av有码| 亚洲欧美清纯卡通| 国产高清三级在线| 国内精品美女久久久久久| 97超视频在线观看视频| 国产淫片久久久久久久久| 天堂网av新在线| 免费观看在线日韩| 九色成人免费人妻av| 日韩亚洲欧美综合| 天堂网av新在线| 欧美xxⅹ黑人| 精品午夜福利在线看| 国产亚洲av嫩草精品影院| 老女人水多毛片| 在线天堂最新版资源| 国产成人免费观看mmmm| 久久久久精品性色| 国产精品不卡视频一区二区| 亚洲精品日韩av片在线观看| 天美传媒精品一区二区| 久久精品国产亚洲av天美| 色综合色国产| 少妇丰满av| 91久久精品国产一区二区成人| 成人亚洲精品av一区二区| 国产男女内射视频| 青春草国产在线视频| 欧美老熟妇乱子伦牲交| 欧美日韩视频精品一区| 王馨瑶露胸无遮挡在线观看| 天堂网av新在线| 在线看a的网站| 午夜激情久久久久久久| 日韩 亚洲 欧美在线| 可以在线观看毛片的网站| 中文天堂在线官网| 深爱激情五月婷婷| 欧美高清性xxxxhd video| 看十八女毛片水多多多| 亚洲精品色激情综合| 亚洲最大成人手机在线| 国产精品国产av在线观看| 午夜福利在线在线| 女人十人毛片免费观看3o分钟| 18禁动态无遮挡网站| 大香蕉久久网| 一本色道久久久久久精品综合| 禁无遮挡网站| a级一级毛片免费在线观看| 亚洲精品视频女| 久久久久久国产a免费观看| 丰满少妇做爰视频| 久久这里有精品视频免费| 欧美+日韩+精品| 久久久久久九九精品二区国产| 亚州av有码| 中文乱码字字幕精品一区二区三区| 免费人成在线观看视频色| 纵有疾风起免费观看全集完整版| 亚洲av电影在线观看一区二区三区 | 成人一区二区视频在线观看| 国产黄频视频在线观看| 国产黄片视频在线免费观看| 国产91av在线免费观看| 寂寞人妻少妇视频99o| 国产男女内射视频| 亚洲国产色片| 久久精品国产亚洲av涩爱| 免费观看a级毛片全部| 69人妻影院| 九九爱精品视频在线观看| 成人亚洲欧美一区二区av| 亚洲精品乱久久久久久| 日韩 亚洲 欧美在线| 欧美日本视频| 国产精品秋霞免费鲁丝片| 亚洲va在线va天堂va国产| 国产免费又黄又爽又色| 国产精品一区二区性色av| 97在线视频观看| 美女高潮的动态| 国产精品一区二区三区四区免费观看| 亚洲av福利一区| 亚洲人与动物交配视频| 搞女人的毛片| 亚洲av免费在线观看| 午夜爱爱视频在线播放| 亚洲国产精品999| av国产免费在线观看| 亚洲高清免费不卡视频| 久久久精品免费免费高清| 亚洲精品,欧美精品| 久久97久久精品| av黄色大香蕉| 国产极品天堂在线| 男男h啪啪无遮挡| 国产高清三级在线| 国产精品福利在线免费观看| 亚洲va在线va天堂va国产| 精品视频人人做人人爽| 少妇丰满av| 性色avwww在线观看| 联通29元200g的流量卡| av免费观看日本| 王馨瑶露胸无遮挡在线观看| 国产又色又爽无遮挡免| 搡女人真爽免费视频火全软件| 禁无遮挡网站| 搞女人的毛片| 国产欧美另类精品又又久久亚洲欧美| 免费看av在线观看网站| 高清午夜精品一区二区三区| 国产精品女同一区二区软件| 国产一区有黄有色的免费视频| 国内少妇人妻偷人精品xxx网站| 偷拍熟女少妇极品色| 成年av动漫网址| 亚洲第一区二区三区不卡| 亚洲内射少妇av| 韩国高清视频一区二区三区| 夫妻性生交免费视频一级片| 一级a做视频免费观看| 色视频www国产| 一个人看的www免费观看视频| 各种免费的搞黄视频| 精品人妻熟女av久视频| 亚洲av二区三区四区| 亚洲av福利一区| 国产精品一二三区在线看| 看十八女毛片水多多多| 日韩强制内射视频| 在线观看人妻少妇| 亚洲av国产av综合av卡| 男男h啪啪无遮挡| 精品一区二区三卡| 精品人妻视频免费看| 欧美日韩国产mv在线观看视频 | 天天躁夜夜躁狠狠久久av| 日本爱情动作片www.在线观看| 亚洲天堂国产精品一区在线| 日韩欧美精品v在线| 欧美xxxx黑人xx丫x性爽| 亚洲丝袜综合中文字幕| 一区二区av电影网| 狂野欧美激情性bbbbbb| 插逼视频在线观看| 国产黄片视频在线免费观看| av专区在线播放| av国产免费在线观看| 禁无遮挡网站| 精品酒店卫生间| 自拍欧美九色日韩亚洲蝌蚪91 | 国产成人aa在线观看| 日韩,欧美,国产一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 肉色欧美久久久久久久蜜桃 | 国产男女超爽视频在线观看| 偷拍熟女少妇极品色| 精品人妻一区二区三区麻豆| 国产精品国产三级专区第一集| 一个人观看的视频www高清免费观看| 中文字幕久久专区| 九九久久精品国产亚洲av麻豆| 国产片特级美女逼逼视频| 伦精品一区二区三区| 色吧在线观看| 亚洲av不卡在线观看| 人人妻人人澡人人爽人人夜夜| 五月玫瑰六月丁香| 成人二区视频| 久久久亚洲精品成人影院| 一二三四中文在线观看免费高清| 人妻 亚洲 视频| 免费播放大片免费观看视频在线观看| 麻豆精品久久久久久蜜桃| 成人亚洲欧美一区二区av| 久久ye,这里只有精品| 久久久久久久久久久丰满| 中文字幕av成人在线电影| 亚洲国产精品成人综合色| 美女被艹到高潮喷水动态| 日韩伦理黄色片| 简卡轻食公司| 高清午夜精品一区二区三区| 国产成人免费无遮挡视频| 可以在线观看毛片的网站| 最新中文字幕久久久久| 色视频www国产| 日本免费在线观看一区| 国产黄片视频在线免费观看| 精品视频人人做人人爽| 高清日韩中文字幕在线| 91精品国产九色| 人妻夜夜爽99麻豆av| av国产精品久久久久影院| 大又大粗又爽又黄少妇毛片口| 哪个播放器可以免费观看大片| 热99国产精品久久久久久7| 十八禁网站网址无遮挡 | 日韩制服骚丝袜av| 香蕉精品网在线| 交换朋友夫妻互换小说| 亚洲精品影视一区二区三区av| 国产精品国产av在线观看| 国产精品久久久久久精品电影| 韩国av在线不卡| 高清午夜精品一区二区三区| 国产成人a∨麻豆精品| 国产在线男女| 大片电影免费在线观看免费| 日本三级黄在线观看| 2022亚洲国产成人精品| 久久久久精品性色| 欧美三级亚洲精品| 97热精品久久久久久| 少妇的逼水好多| 免费电影在线观看免费观看| 中文字幕免费在线视频6| 成人毛片a级毛片在线播放| 午夜免费鲁丝| 精品久久久久久久人妻蜜臀av| 亚洲国产精品专区欧美| 晚上一个人看的免费电影| 国产爱豆传媒在线观看| 男人和女人高潮做爰伦理| 岛国毛片在线播放| 久久久精品欧美日韩精品| 丝袜喷水一区| 久久久久久久久大av| 免费看av在线观看网站| 国产精品国产三级专区第一集| 国产亚洲av嫩草精品影院| www.av在线官网国产| 免费观看av网站的网址| 最近中文字幕高清免费大全6| 一级毛片aaaaaa免费看小| 九九久久精品国产亚洲av麻豆| 建设人人有责人人尽责人人享有的 | 国产探花极品一区二区| 青春草亚洲视频在线观看| 又黄又爽又刺激的免费视频.| 美女被艹到高潮喷水动态| 国产伦精品一区二区三区视频9| 建设人人有责人人尽责人人享有的 | 黄色一级大片看看| 日韩成人av中文字幕在线观看| 国产欧美日韩精品一区二区| 国产成年人精品一区二区| 国产欧美日韩一区二区三区在线 | 亚洲美女搞黄在线观看| 哪个播放器可以免费观看大片| 国产黄色视频一区二区在线观看| 三级国产精品欧美在线观看| 亚洲最大成人手机在线| 一级毛片黄色毛片免费观看视频| 美女内射精品一级片tv| 岛国毛片在线播放| av在线播放精品| 97超碰精品成人国产| 国产亚洲91精品色在线| 日本一本二区三区精品| 国产精品国产三级国产av玫瑰| 欧美日韩亚洲高清精品| 成人黄色视频免费在线看| 久久久久国产精品人妻一区二区| 国产爱豆传媒在线观看| 国产精品久久久久久精品电影小说 | av在线蜜桃| 欧美日韩国产mv在线观看视频 | 最新中文字幕久久久久| 久久久久久久大尺度免费视频| 蜜臀久久99精品久久宅男| 在线 av 中文字幕| 日本黄大片高清| 插阴视频在线观看视频| 人体艺术视频欧美日本| 性色avwww在线观看| 精品酒店卫生间| 国产综合精华液| 欧美潮喷喷水| 精品久久国产蜜桃| 午夜免费观看性视频| 黑人高潮一二区| 国产乱人偷精品视频| 黄片wwwwww| 啦啦啦啦在线视频资源| 一区二区三区精品91| 熟女人妻精品中文字幕| 久久热精品热| 亚洲人与动物交配视频| 国产乱人视频| 久久精品夜色国产| 久久精品人妻少妇| 亚洲精品第二区| 国产女主播在线喷水免费视频网站| 亚洲va在线va天堂va国产| 一个人看的www免费观看视频| 最近中文字幕高清免费大全6| 中文资源天堂在线| 国产精品女同一区二区软件| 亚洲,欧美,日韩| 国产高清国产精品国产三级 | 亚洲精品日韩在线中文字幕| 免费黄色在线免费观看| 午夜福利在线在线| 韩国高清视频一区二区三区| av播播在线观看一区| 亚洲精品成人av观看孕妇| 午夜福利网站1000一区二区三区| 麻豆成人午夜福利视频| 国产永久视频网站| 乱系列少妇在线播放| 在线看a的网站| 2021少妇久久久久久久久久久| 欧美区成人在线视频| 亚洲高清免费不卡视频| 99精国产麻豆久久婷婷| 春色校园在线视频观看| 大话2 男鬼变身卡| 在线免费十八禁| 麻豆久久精品国产亚洲av| 国产亚洲av嫩草精品影院| 熟女电影av网| 亚洲最大成人手机在线| 国产免费一级a男人的天堂| 我的老师免费观看完整版| 国产精品爽爽va在线观看网站| 尤物成人国产欧美一区二区三区| 亚洲av成人精品一二三区| 精品99又大又爽又粗少妇毛片| av国产免费在线观看| 91精品国产九色| 成人黄色视频免费在线看| 免费看av在线观看网站| 婷婷色综合大香蕉| 丝瓜视频免费看黄片| 视频区图区小说| 熟妇人妻不卡中文字幕| 久久精品人妻少妇| 精品久久久久久久末码| 小蜜桃在线观看免费完整版高清| 久久99热这里只有精品18| 亚洲美女搞黄在线观看| 日日撸夜夜添| 久久久久久久亚洲中文字幕| xxx大片免费视频| 国产成人精品一,二区| 国产一区亚洲一区在线观看| 久久精品熟女亚洲av麻豆精品| 国产黄a三级三级三级人| 男人添女人高潮全过程视频| 欧美国产精品一级二级三级 | 久久久a久久爽久久v久久| 国产精品国产三级国产av玫瑰| 国产 一区 欧美 日韩| 最近最新中文字幕免费大全7| av在线app专区| 免费播放大片免费观看视频在线观看| 日日啪夜夜爽| av又黄又爽大尺度在线免费看| 天堂网av新在线| 大片电影免费在线观看免费| 少妇的逼水好多| 久久久精品欧美日韩精品| 亚洲国产精品999| 免费大片18禁| 五月开心婷婷网| 亚洲高清免费不卡视频| 亚洲电影在线观看av| 国产精品嫩草影院av在线观看| 国产女主播在线喷水免费视频网站| 在线免费十八禁| 精品人妻偷拍中文字幕| 久久久久性生活片| 婷婷色综合大香蕉| 精品人妻一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 免费人成在线观看视频色| 国产亚洲午夜精品一区二区久久 | 久久韩国三级中文字幕| 老司机影院毛片| 最近手机中文字幕大全| 纵有疾风起免费观看全集完整版| 少妇熟女欧美另类| 亚洲成人中文字幕在线播放| 亚洲精品乱码久久久v下载方式| 免费av毛片视频| 亚洲av在线观看美女高潮| 免费电影在线观看免费观看| 成人一区二区视频在线观看| 国产精品精品国产色婷婷| 久久久久久久久久成人| 男男h啪啪无遮挡| 美女脱内裤让男人舔精品视频| 亚洲四区av| 国内少妇人妻偷人精品xxx网站| 午夜激情福利司机影院| av免费观看日本| 丝瓜视频免费看黄片| 成人国产av品久久久| 大片电影免费在线观看免费| 国产成人午夜福利电影在线观看| 久久99热这里只频精品6学生| 成年版毛片免费区| 亚洲最大成人手机在线| 亚洲av中文av极速乱| 女人被狂操c到高潮| 日韩伦理黄色片| 国产乱人偷精品视频| 最近最新中文字幕免费大全7| 亚洲三级黄色毛片| 一级二级三级毛片免费看| 午夜福利高清视频| 亚洲精品第二区| 亚洲熟女精品中文字幕| 别揉我奶头 嗯啊视频| 97超视频在线观看视频| 一级毛片aaaaaa免费看小| 在线观看国产h片| 看免费成人av毛片| 春色校园在线视频观看| 亚洲精品亚洲一区二区| 亚洲国产成人一精品久久久| 18禁裸乳无遮挡动漫免费视频 | 视频中文字幕在线观看| 老司机影院毛片| videossex国产| 国产免费一级a男人的天堂| 秋霞在线观看毛片| 国内揄拍国产精品人妻在线| 在线观看一区二区三区激情| 激情 狠狠 欧美| 久久精品国产亚洲av涩爱| 18禁在线无遮挡免费观看视频| 亚洲第一区二区三区不卡| 51国产日韩欧美| 日本色播在线视频| 免费黄频网站在线观看国产| 亚洲国产欧美人成| kizo精华| 国产成人免费无遮挡视频| 国产精品久久久久久久电影| 少妇人妻精品综合一区二区| 成人黄色视频免费在线看| 高清在线视频一区二区三区| 久久人人爽人人爽人人片va| 高清视频免费观看一区二区| 日日摸夜夜添夜夜爱| 91精品一卡2卡3卡4卡| 精品一区在线观看国产| 女人久久www免费人成看片| 插阴视频在线观看视频| 亚洲精品日韩在线中文字幕| 精品99又大又爽又粗少妇毛片| 久久久久久久久大av| 中文字幕久久专区| 日本av手机在线免费观看| av天堂中文字幕网| 国产精品99久久99久久久不卡 | 亚洲国产成人一精品久久久| 国产一区亚洲一区在线观看| 国产伦理片在线播放av一区| 欧美成人一区二区免费高清观看| 亚洲av欧美aⅴ国产| 久久精品久久久久久久性| 少妇人妻 视频| 男人添女人高潮全过程视频| 日本av手机在线免费观看| 插逼视频在线观看| 韩国av在线不卡| 少妇人妻 视频| 欧美一级a爱片免费观看看| 噜噜噜噜噜久久久久久91| 亚洲色图av天堂| 国产片特级美女逼逼视频| 日韩欧美精品v在线| 国产日韩欧美亚洲二区| 狂野欧美激情性xxxx在线观看| 国产精品.久久久| 中文乱码字字幕精品一区二区三区| 久久99蜜桃精品久久| 又爽又黄a免费视频| 久久久精品94久久精品| h日本视频在线播放| 亚洲欧美中文字幕日韩二区| 婷婷色综合大香蕉| 男人添女人高潮全过程视频| 最近中文字幕高清免费大全6| 亚洲电影在线观看av| 麻豆精品久久久久久蜜桃| 久久久成人免费电影| 80岁老熟妇乱子伦牲交| 91精品伊人久久大香线蕉| 我的老师免费观看完整版| 一级av片app| 爱豆传媒免费全集在线观看| 亚洲性久久影院| 一本久久精品| 亚洲欧美精品专区久久| 国产精品麻豆人妻色哟哟久久| 一级二级三级毛片免费看| 在线a可以看的网站| 在线免费观看不下载黄p国产| 久久精品国产自在天天线| 大又大粗又爽又黄少妇毛片口| 免费大片18禁| 午夜福利在线观看免费完整高清在| a级毛色黄片| 99久久中文字幕三级久久日本| 欧美老熟妇乱子伦牲交| 亚洲精品国产av成人精品| 自拍偷自拍亚洲精品老妇| 尤物成人国产欧美一区二区三区| 久久久精品欧美日韩精品| 99精国产麻豆久久婷婷| 国产视频内射| 亚洲精品乱久久久久久| 特级一级黄色大片| 最近手机中文字幕大全| 国产一区二区三区av在线| 国产成人a区在线观看| 天美传媒精品一区二区| 麻豆久久精品国产亚洲av| 国产精品蜜桃在线观看| 国产成人精品婷婷| 大香蕉久久网| 51国产日韩欧美| 亚洲伊人久久精品综合| 亚洲无线观看免费| 欧美日韩国产mv在线观看视频 | 亚洲最大成人av| 欧美区成人在线视频| 建设人人有责人人尽责人人享有的 | 97人妻精品一区二区三区麻豆| 黄色欧美视频在线观看| 免费观看无遮挡的男女| 日韩av在线免费看完整版不卡| 国产乱人偷精品视频| 一级二级三级毛片免费看| 免费观看在线日韩| 特大巨黑吊av在线直播| 国国产精品蜜臀av免费| 大陆偷拍与自拍| 热re99久久精品国产66热6| 亚洲自偷自拍三级| 岛国毛片在线播放| 免费黄色在线免费观看| 高清日韩中文字幕在线| 亚洲在久久综合| 简卡轻食公司| 久久99热这里只有精品18| 中国美白少妇内射xxxbb| 国内精品宾馆在线| 久久精品久久久久久噜噜老黄| 国产欧美日韩精品一区二区| 一区二区三区免费毛片| 国产黄频视频在线观看| 中文字幕久久专区| av卡一久久| 亚洲精品,欧美精品| 夜夜爽夜夜爽视频| 男男h啪啪无遮挡| 一区二区三区精品91| 夜夜爽夜夜爽视频| 天天躁日日操中文字幕| 只有这里有精品99| 午夜精品一区二区三区免费看| 卡戴珊不雅视频在线播放| 蜜臀久久99精品久久宅男| 国产日韩欧美在线精品| 卡戴珊不雅视频在线播放| 国产精品国产三级专区第一集| 国产成人精品一,二区| 亚洲国产精品成人综合色| 蜜臀久久99精品久久宅男| 国产日韩欧美在线精品| 在线观看免费高清a一片| 蜜臀久久99精品久久宅男| 精品久久久噜噜| av卡一久久| 黄片无遮挡物在线观看| 国产欧美日韩精品一区二区| 精品久久久精品久久久| 五月伊人婷婷丁香| 国产精品人妻久久久久久| 国产精品熟女久久久久浪| 亚洲欧美成人综合另类久久久| 人妻一区二区av| 黄色视频在线播放观看不卡| 七月丁香在线播放|