鬼佬大哥大
  • / 12
  • 下載費用:30 金幣  

一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法.pdf

關 鍵 詞:
一種 采樣率 脈沖 UWB 通信 系統 中的 窄帶 干擾 消除 方法
  專利查詢網所有資源均是用戶自行上傳分享,僅供網友學習交流,未經上傳用戶書面授權,請勿作他用。
摘要
申請專利號:

CN201210194499.5

申請日:

2012.06.13

公開號:

CN102710290B

公開日:

2015.01.28

當前法律狀態:

授權

有效性:

有權

法律詳情: 授權|||實質審查的生效IPC(主分類):H04B 1/719申請日:20120613|||公開
IPC分類號: H04B1/719(2011.01)I 主分類號: H04B1/719
申請人: 清華大學深圳研究生院
發明人: 晉本周; 張盛; 潘劍; 林孝康
地址: 518055 廣東省深圳市南山區西麗大學城清華校區A306
優先權:
專利代理機構: 深圳市瑞方達知識產權事務所(普通合伙) 44314 代理人: 張秋紅;張約宗
PDF完整版下載: PDF下載
法律狀態
申請(專利)號:

CN201210194499.5

授權公告號:

102710290B||||||

法律狀態公告日:

2015.01.28|||2013.05.01|||2012.10.03

法律狀態類型:

授權|||實質審查的生效|||公開

摘要

本發明涉及一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法和裝置,該方法包括以下步驟:S1、對接收的第一個未包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,得到采樣向量y1,計算窄帶干擾信號的正交投影矩陣;S2、對后續接收的包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,得到采樣向量yi,基于NBI信號的正交投影矩陣得到消除窄帶干擾信號后的UWB信號。本發明可以有效地消除NBI信號,提高信干比,其系統采樣率遠低于奈奎斯特采樣率,且本發明中NBI估計與消除完全放在數字域進行,無需并行系統結構,也無需產生模擬投影信號,大大降低了系統復雜度和實現難度。

權利要求書

權利要求書
1.  一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法,其特征在于,包括以下步
驟:
S1、對接收的第一個未包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,得到
M×1維采樣向量y1,基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
S2、對后續接收的包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,得到M×1
維采樣向量yi=[yi(0),yi(1),...,yi(M-1)],其中i表示信息包的序號,基于NBI信號的正
交投影矩陣通過 得到消除窄帶干擾信號后的UWB信號。

2.  根據權利要求1所述的方法,其特征在于,所述步驟S1中通過以下步驟基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
S11、對向量y1進行M點的快速傅里葉變換得到其FFT變換向量λ;
S12、設定判決門限Th,通過以下公式得到標識NBI分量位置的向量h,
(j=0,1,...,M-1);
S13、計算窄帶干擾信號的正交投影矩陣正交投影矩陣為以向量f作為第一列、pr作為第一行的托普利茨矩陣;其中向量f為向量h求反傅里葉變換獲得,pr=[f(0),f(M-1),f(M-2),…,f(1)]。

3.  一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除裝置,其特征在于,包括:前端濾波器、低速ADC和數字信號處理器;
所述前端濾波器用于對接收的信息包濾除帶外噪聲;其中,接收的第一個信息包未包含UWB信號,后續接收的信息包包含UWB信號;
所述低速ADC用于對濾除帶外噪聲后的信號進行采樣,得到M×1維采樣向量yi=
[yi(0),yi(1),...,yi(M-1)],其中i表示信息包的序號;
所述數字信號處理器用于基于采樣向量y1計算窄帶干擾信號的正交投影矩陣并基于NBI信號的正交投影矩陣通過得到消除窄帶干擾信號后的UWB信號。

4.  根據權利要求3所述的裝置,其特征在于,所述數字信號處理器進一步包括窄帶干
擾信號估計模塊和窄帶干擾信號消除模塊;所述窄帶干擾信號估計模塊用于基于采樣向量y1計算窄帶干擾信號的正交投影矩陣

所述窄帶干擾信號消除模塊用于基于NBI信號的正交投影矩陣通過 得到消除窄帶干擾信號后的UWB信號。

5.  根據權利要求4所述的裝置,其特征在于,所述數字信號處理器還包括信號檢測模塊,用于將UWB信號轉換為信息符號輸出。

6.  根據權利要求4所述的裝置,其特征在于,所述窄帶干擾信號估計模塊通過以下步
驟基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
對向量y1進行M點的快速傅里葉變換得到其FFT變換向量λ;設定判決門限Th,通過以下公式得到標識NBI分量位置的向量h,
(j=0,1,...,M-1);
計算窄帶干擾信號的正交投影矩陣正交投影矩陣為以向量f作為第一列、pr作為第一行的托普利茨矩陣;其中向量f為向量h求反傅里葉變換獲得,pr=[f(0),f(M-1),f(M-2),…,f(1)]。

7.  一種欠采樣率脈沖UWB通信系統,其特征在于,包括權利要求3-6中任意一項所述的欠采樣率脈沖UWB通信系統中的窄帶干擾消除裝置。

說明書

說明書一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法
技術領域
[0001]本發明涉及無線通信領域,更具體地說,涉及一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法和裝置。
背景技術
[0002] 脈沖超寬帶(UWB)技術作為一種無線通信技術方案,以其低復雜度、低成本、高定位精度等優點,在無線傳感器網絡和高精度定位與導航系統中有著廣泛的應用前景。為了與現有窄帶系統共享頻譜,UWB系統的發射功率譜密極低,遠低于窄帶系統。然而,對UWB通信來說,窄帶信號是一種干擾,并且功率很高,窄帶干擾(NBI)會嚴重影響UWB系統的性能。因此,NBI的估計與消除是UWB通信中非常重要的問題。
[0003]對于傳統數字接收機,NBI信號檢測方案中通常以奈奎斯特率對接收信號r(t)進行采樣,然后基于快速傅里葉變換(FFT)對采樣信號進行分析,估計出NBI信號在頻域的精確位置及譜特征,最后通過設計濾波器完成對NBI信號的消除。然而,由于UWB通信帶寬很大,可達數吉赫茲,根據奈奎斯特采樣定理,這種方案需要極高的采樣頻率,很難由單一模數轉換器(ADC)實現,即使實現,高速ADC也會導致很高的系統功耗。
[0004]近幾年來出現的壓縮感知(CS)理論為NBI信號檢測與消除提供了新的思路。基于CS方案的接收機結構如圖1所示。接收機無需以奈奎斯特率進行采樣,而是采用一組并行的乘積累積器和低速ADC對接收信號進行投影,得到采樣向量y,然后在數字域完成對UWB信號的檢測。其中,投影函數{ψi(t)}的選取主要有如下兩種方法:
[0005] (1)利用NBI信號在UWB信號頻譜范圍內稀疏的特點,基于壓縮感知理論估計出NBI信號子空間,然后基于該子空間設計一個投影矩陣P,使得NBI信號屬于的P零空間,此時,P可將NBI信號映射為零向量。為了系統實現,該方法基于P設計一組投影函數{ψi(t)}將NBI信號映射為零,以此消除窄帶干擾。然而,該方法中NBI信號子空間估計復雜度高,且投影函數波形產生困難。
[0006](2)第二種方法是將{ψi(t)}選取為一組單頻正弦或者余弦信號,其頻率均勻分布在UWB信號頻譜范圍內。通常窄帶干擾的功率遠大于UWB信號,如果在某一個支路所對應的頻點處存在窄帶干擾,采樣值會出現異常,并可以被檢測出來。那么,直接將該支路的采樣值刪除便可抑制NBI。然而,由于UWB信號能量在整個頻譜中分散的特點,采用單頻點投影函數不能有效的獲取UWB信號能量,不利于UWB信號的檢測。
[0007]另外,雖然基于壓縮感知的方案降低了對采樣率的要求,但并行的系統結構使得接收機復雜度非常高。
發明內容
[0008]本發明要解決的技術問題在于,針對現有UWB通信系統中NBI的估計與消除方法復雜度高的缺陷,提供一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法,接收機端直接對濾除帶外噪聲后信號進行低速采樣,NBI信號估計與消除完全在數字域進行,其系統
復雜度低,易于實現。
[0009] 本發明解決其技術問題所采用的技術方案是:構造一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法,包括以下步驟:
[0010] S1、對接收的第一個未包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,
得到M×1維采樣向量y1,基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
[0011] S2、對后續接收的包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,得到
M×1維采樣向量yi=[yi(0),yi(1),...,yi(M-1)],其中i表示信息包的序號,基于NBI信號
的正交投影矩陣通過得到消除窄帶干擾信號后的UWB信號。
[0012] 在根據本發明所述的欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法中,所述步驟S1中通過以下步驟基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
[0013] S11、對向量y1進行M點的快速傅里葉變換得到其FFT變換向量λ;
[0014] S12、設定判決門限Th,通過以下公式得到標識NBI分量位置的向量h,
[0015] (j=0,1,...,M-1);
[0016] S13、計算窄帶干擾信號的正交投影矩陣正交投影矩陣為以向量f作為第一列、pr作為第一行的托普利茨矩陣;其中向量f為向量h求反傅里葉變換獲得,pr=[f(0),f(M-1),f(M-2),…,f(1)]。
[0017] 本發明還提供了一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除裝置,包括:前端濾波器、低速ADC和數字信號處理器;
[0018] 所述前端濾波器用于對接收的信息包濾除帶外噪聲;其中,接收的第一個信息包未包含UWB信號,后續接收的信息包包含UWB信號;
[0019] 所述低速ADC用于對濾除帶外噪聲后的信號進行采樣,得到M×1維采樣向量yi=
[yi(0),yi(1),...,yi(M-1)],其中i表示信息包的序號;
[0020] 所述數字信號處理器用于基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
并基于NBI信號的正交投影矩陣通過得到消除窄帶干擾信號后的UWB信號。[0021]在根據本發明所述的欠采樣率脈沖UWB通信系統中的窄帶干擾消裝置中,所述數字信號處理器進一步包括窄帶干擾信號估計模塊和窄帶干擾信號消除模塊;
[0022] 所述窄帶干擾信號估計模塊用于基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
[0023] 所述窄帶干擾信號消除模塊用于基于NBI信號的正交投影矩陣通過 得到消除窄帶干擾信號后的UWB信號。
[0024] 在根據本發明所述的欠采樣率脈沖UWB通信系統中的窄帶干擾消裝置中,所述數字信號處理器還包括信號檢測模塊,用于將UWB信號轉換為信息符號輸出。
[0025] 在根據本發明所述的欠采樣率脈沖UWB通信系統中的窄帶干擾消裝置中,所述窄帶干擾信號估計模塊通過以下步驟基于采樣向量y1計算窄帶干擾信號的正交投影矩陣

[0026] 對向量y1進行M點的快速傅里葉變換得到其FFT變換向量λ;設定判決門限Th,通過以下公式得到標識NBI分量位置的向量h,
[0027] (j=0,1,...,M-1);
[0028] 計算窄帶干擾信號的正交投影矩陣正交投影矩陣為以向量f作為第一列、pr作為第一行的托普利茨矩陣;其中向量f為向量h求反傅里葉變換獲得,pr=[f(0),f(M-1),f(M-2),…,f(1)]。
[0029] 本發明還提供了一種欠采樣率脈沖UWB通信系統,包括如上所述的欠采樣率脈沖
UWB通信系統中的窄帶干擾消除裝置。
[0030]實施本發明的欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法和裝置,具有以下有益效果:本發明可以有效地消除NBI信號,提高信干比,其與基于奈奎斯特采樣的方法相比較,系統采樣率可以得到大大降低,降低了ADC的實現難度和功耗;與基于CS的方案相比較,本發明直接對濾波后的接收信號進行低速采樣,NBI估計與消除完全放在數字域進行,無需并行系統結構,也無需產生模擬投影信號,大大降低了系統復雜度和實現難度。
附圖說明
[0031] 下面將結合附圖及實施例對本發明作進一步說明,附圖中:
[0032] 圖1為現有技術中基于CS方案的接收機結構圖;
[0033] 圖2為根據本發明的欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法的流程圖;
[0034] 圖3為根據本發明的方法中計算窄帶干擾信號的正交投影矩陣的流程圖;
[0035] 圖4為根據本發明的欠采樣率脈沖UWB通信系統中的窄帶干擾消除裝置的模塊框圖。
具體實施方式
[0036] 為了使本發明的目的、技術方案及優點更加清楚明白,以下結合附圖及實施例,對本發明進行進一步詳細說明。
[0037] 本發明中,首先假定K個比特組成一個信息包,以包為單位進行信息發送和接收,當前包與下一個包之間的時間間隔足夠大以使得接收信號中無包間干擾。本發明構思了一種欠采樣率脈沖UWB通信系統的窄帶干擾信號消除方法和系統,主要原理在于使第一個信號包不攜帶UWB信號,利用該第一個信號包進行NBI信號估計,利用估計的NBI信號參數值對后續信號包中的NBI信號進行濾除。
[0038] 下面對本發明的NBI信號估計和消除的原理及具體推算過程進行說明。
[0039]在本發明中,欠采樣率脈沖UWB通信系統的接收機首先對接收到的信息包的信號進行前端濾波,濾除系統工作頻帶以外的噪聲,即濾除帶外噪聲,然后由低速ADC采樣得到當前包的M×1維采樣向量yi,最后將yi送入到數字信號處理器進行窄帶干擾的估計或消除。UWB信道估計和信號檢測算法已有相關文獻介紹,本發明中不再論述。
[0040]本發明中欠采樣率脈沖UWB通信系統的結構非常簡單,只包含前端濾波器、低速ADC和數字信號處理器三個模塊。而窄帶干擾的估計或消除主要在數字信號處理器,因此在此對數字信號處理器的原理進行詳細描述。
[0041] 令△t表示低速ADC的采樣時間間隔,Ty表示經過UWB信道后每個包的持續時間,那么當前接收的信息包的低速采樣值可以表示為:
[0042]
[0043] 其中,yi=[yi(0),yi(1),…,yi(M-1)],i表示信息包的序號;uUWB(t)為接收到的UWB信號,其帶寬為B,中心頻率為fc;n(t)為帶限高斯噪聲;nI(t)是帶寬為BI,中心頻率為fI的NBI信號。NBI信號帶寬遠小于UWB信號帶寬,即BI<<B。另外,fc-B/2-BI/2≤fI≤fc+B/2+BI/2,即NBI信號位于UWB通信頻段內。為了衡量系統采樣率降低的程度,定義降采樣系數:
[0044]
[0045]其中,用2B近似表示UWB信號的奈奎斯特采樣頻率。
[0046]先前的NBI消除方案,通常都需要估計出NBI信號在頻域的精確位置,然后通過時域或者頻域濾波的方法去除NBI。然而,本發明指出,在UWB通信中,去除窄帶干擾是最終目的,在這個目標前提下,NBI信號的精確估計不是必須的。
[0047] 本發明的基本思想如下:在欠采樣系統中,雖然通常有Us<<1,但為了獲取足夠多用于信號檢測的UWB信號能量,系統采樣率和2倍NBI信號帶寬相比較仍然足夠大,即滿足
1/Δt>>2BI。因此,除不能反映NBI信號在頻域的精確頻點位置外,采樣向量yi的FFT可以近似完整的反映NBI信號的頻域特征,并且NBI分量在FFT域稀疏。另外,由于NBI信號的能量會遠大于UWB信號和噪聲的能量,那么,在FFT域,NBI分量幅度值也會遠大于其它分量,即容易從FFT域辨別NBI分量所在的位置。NBI信號在頻域的稀疏性和顯著性使得低速采樣下近似實現對NBI信號的檢測和消除成為可能。
[0048]令y1表示第一個包的采樣向量,為了在時域和數字域消除后續包中的NBI信號,
通過y1設計NBI信號的正交投影矩陣那么,后續包中的NBI信號就可以通過如下方法去除:
[0049]
[0050] 其中,i≥2,uUWB和n分別表示UWB信號和噪聲在當前包中的采樣向量。下邊給出矩陣的推導過程。
[0051]由FFT,y1可以重新表示為:
[0052]y1=Wξ+Wκ=Wλ,(4)
[0053] 其中,W∈□M×M 為FFT變換矩陣。ξ∈□M×1和κ∈□M×1分別為n
域表示。
和n的FFT
[0054]由于NBI信號能量遠大于噪聲能量,且1/Δt>>2BI,那么,nI可以由向量λ中的少數幾個絕對值很大的分量近似表示。基于該特點,可以通過對|λ|設定一個門限來檢測出NBI分量在λ中的位置,即,
[0055] (j=0,1,...,M-1), (5)
[0056] 其中,h∈M×1,向量h中元素為零的點表明此處存在NBI分量,否則不存在
NBI。在此,選擇修正的小波(Donoho)閾值作為判決門限,即
[0057] [0058] 其中,α>0,σ=median(|λ|)/0.6745,median表示求中位數。[0059] 那么,NBI信號就可以采用如下方法得到去除:
[0060] [0061] 其中 表示循環卷積,f=IFFT(h).
[0062] 令pr=[f(0),f(M-1),f(M-2),…,f(1)],正交投影矩陣可以設計為:
[0063] (8)
[0064] 是一個以向量f作為第一列,pr作為第一行的托普利茨(Toeplitz)矩陣。在
UWB信號檢測中可以直接作為UWB信息獲取的組成部分。
[0065]請參閱圖2,為根據本發明的欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法的流程圖。下述步驟中選用的參數與前述對NBI信號的估計和消除的原理描述中涉及的參數一致,且原理也一致。如圖2所示,該實施例的方法包括以下步驟:
[0066] 首先,在步驟S1中,對接收的第一個未包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,得到M×1維采樣向量y1,基于采樣向量y1計算窄帶干擾信號的正交投影
矩陣本實施例中可以選取如下參數:B=3GHz,fc=4GHz,BI=10MHz,fI=3.7GHz,Ty=900ns,α=0.5。降采樣系數Us=1/7,即在本例中對應的采樣率fs約為870MHz。第一個接收到的包中不含UWB信號,只有NBI信號和噪聲。在第一個信號包的持續時間內,以采樣率fs對濾除帶外噪聲后的信號進行采樣,得到M×1維采樣向量y1,M的值由Ty和fs共同決定。[0067] 隨后,在步驟S2中,對后續接收的包含UWB信號的信息包濾除帶外噪聲后的信號進行采樣,得到M×1維采樣向量yi=[yi(0),yi(1),...,yi(M-1)],其中i表示信息包的序
號,基于NBI信號的正交投影矩陣通過 得到消除窄帶干擾信號后的UWB信號。本步驟在時域完成對NBI信號的消除。用投影矩陣乘后續包的采樣向量yi,即通過完成對NBI信號的消除。消除NBI信號后,即可根據現有的信號檢測算法完成UWB信號的檢測。[0068] 請參閱圖3為根據本發明的方法中計算窄帶干擾信號的正交投影矩陣的流程圖。
上述步驟S1中通過以下步驟基于采樣向量y1計算窄帶干擾信號的正交投影矩陣
[0069] 首先,在步驟S11中,對向量y1進行M點的快速傅里葉變換得到其FFT變換向量
λ。
[0070] 隨后,在步驟S12中,設定判決門限Th,通過以下公式得到標識NBI分量位置的向量h,
[0071] (j=0,1,...,M-1)。
[0072] 其中,判決門限Th利用前述原理中提到的方法如公式(6)進行計算。
[0073] 最后,在步驟S13中,對向量h求反傅里葉變換得到向量f。計算窄帶干擾信號的
正交投影矩陣正交投影矩陣為以向量f作為第一列、pr作為第一行的托普利茨矩陣;其中pr=[f(0),f(M-1),f(M-2),…,f(1)]。
[0074] 請參閱圖4,為根據本發明的欠采樣率脈沖UWB通信系統中的窄帶干擾消除裝置
的模塊框圖,該窄帶干擾消除裝置包括:前端濾波器10、低速ADC20和數字信號處理器30;[0075] 其中,前端濾波器10用于對接收的信息包濾除帶外噪聲;該欠采樣率脈沖UWB通信系統接收的第一個信息包未包含UWB信號,后續接收的信息包包含UWB信號。
[0076] 低速ADC20用于對濾除帶外噪聲后的信號進行采樣,得到M×1維采樣向量
yi=[yi(0),yi(1),...,yi(M-1)],其中i表示信息包的序號。
[0077] 數字信號處理器30用于基于采樣向量y1計算窄帶干擾信號的正交投影矩陣并基于NBI信號的正交投影矩陣通過 得到消除窄帶干擾信號后的UWB信號。
[0078] 上述數字信號處理器30進一步包括窄帶干擾信號估計模塊31和窄帶干擾信號消
除模塊32,其中窄帶干擾信號估計模塊31用于基于采樣向量y1計算窄帶干擾信號的正交
投影矩陣窄帶干擾信號消除模塊32用于基于NBI信號的正交投影矩陣通過得到消除窄帶干擾信號后的UWB信號。優選地,數字信號處理器還包括信號檢測模塊33,用于將UWB信號轉換為信息符號輸出。
[0079] 其中窄帶干擾信號估計模塊31通過以下步驟基于采樣向量y1計算窄帶干擾信號
的正交投影矩陣
[0080] 首先,對向量y1進行M點的快速傅里葉變換得到其FFT變換向量λ;
[0081] 隨后,設定判決門限Th,通過以下公式得到標識NBI分量位置的向量h,
[0082] (j=0,1,...,M-1);
[0083] 最后,計算窄帶干擾信號的正交投影矩陣正交投影矩陣為以向量f作為第一列、pr作為第一行的托普利茨矩陣;其中向量f為向量h求反傅里葉變換獲得,pr=[f(0),f(M-1),f(M-2),…,f(1)]。
[0084] 本發明還提供了一種欠采樣率脈沖UWB通信系統,包括前述的欠采樣率脈沖UWB
通信系統中的窄帶干擾消除裝置。
[0085]綜上所述,本發明中的NBI信號估計與消除方法和基于奈奎斯特采樣的方法相比較,系統采樣率可以得到大大降低,因此,降低了ADC的實現難度和功耗。另外,與基于CS的方案相比較,本方法直接對濾波后的接收信號進行低速采樣,NBI估計與消除完全放在數字域進行,無需并行系統結構,也無需產生模擬投影信號,大大降低了系統復雜度和實現難度。仿真結果表明,本發明方案可以有效地消除NBI信號,提高信干比。
[0086] 本發明是根據特定實施例進行描述的,但本領域的技術人員應明白在不脫離本發明范圍時,可進行各種變化和等同替換。此外,為適應本發明技術的特定場合或材料,可對
本發明進行諸多修改而不脫離其保護范圍。因此,本發明并不限于在此公開的特定實施例,而包括所有落入到權利要求保護范圍的實施例。

關于本文
本文標題:一種欠采樣率脈沖UWB通信系統中的窄帶干擾消除方法.pdf
鏈接地址:http://www.wwszu.club/p-6420776.html
關于我們 - 網站聲明 - 網站地圖 - 資源地圖 - 友情鏈接 - 網站客服 - 聯系我們

[email protected] 2017-2018 zhuanlichaxun.net網站版權所有
經營許可證編號:粵ICP備17046363號-1 
 


收起
展開
鬼佬大哥大