JP6038562B2 - 受信装置、雑音抑圧方法及び雑音抑圧プログラム - Google Patents
受信装置、雑音抑圧方法及び雑音抑圧プログラム Download PDFInfo
- Publication number
- JP6038562B2 JP6038562B2 JP2012204504A JP2012204504A JP6038562B2 JP 6038562 B2 JP6038562 B2 JP 6038562B2 JP 2012204504 A JP2012204504 A JP 2012204504A JP 2012204504 A JP2012204504 A JP 2012204504A JP 6038562 B2 JP6038562 B2 JP 6038562B2
- Authority
- JP
- Japan
- Prior art keywords
- correlation matrix
- unit
- matrix
- noise
- eigenvalue
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Description
X=[x1(mT) x2(mT)…xN(mT)]T (1)
と定義する。ここで、mはサンプル番号を示し、Tはサンプルの時間間隔を示し、[・]Tは転置の操作を示す。そうすると、第1の相関行列Rxx Aは、
Rxx A=〈X・XH〉 (2)
と表せる。ここで、〈・〉はアンサンブル平均を示し、[・]Hは共役複素転置を示す。なお、第1の相関行列Rxx Aは、
Rxx A=ASSHAH+NNH (3)
と分解される。ここで、Aは信号の方向の情報を含むステアリング行列を示し、Sは信号成分を含む行列を示し、Nは雑音成分の平均電力を含む行列を示す。
〈NNH〉=σ2I (4)
と表される。ここで、σ2は単一のアンテナで受信された場合の平均雑音電力を示し、Iは単位行列を示す。この場合、相関行列を固有値分解すると、取得される固有値は、信号成分についての固有値に一定の雑音成分が付加されたものとなる。つまり、本質的には、雑音成分は、固有値分解の結果に影響を与えないこととなる。
〈NNH〉=diag[σ1 2 σ2 2…σN 2] (5)
と表される。この場合、相関行列を固有値分解すると、固有値毎に付加される雑音成分の大きさが異なることになる。
Z=[1 k k2…kN−1] (6)
と設定する。kは実数であり、0<k<1である。第2の相関行列演算部42は、変換行列Zにおけるkの値を、測定部47により測定されるC/N(Carrier-Noise Ratio)に応じて変化させる。すなわち、第2の相関行列演算部42は、C/Nが所定の値よりも高い場合にはkを1又は1に近い値に設定し、低い場合にはkを0又は0に近い値に設定する。なお、kの値は、受信装置の操作者がC/Nに応じて手動で変更するようにしてもかまわない。
Rxx B=ZRxx AZ−1 (7)
の計算を行い、第2の相関行列Rxx Bを算出する。このとき、〈N・NH〉の相似変換は、
Z〈N・NH〉Z−1=〈N・NH〉 (8)
となる。
Rxx=〈Rxx A〉−〈Rxx B〉 (9)
を計算することにより、第3の相関行列Rxxを算出する。これにより、第1の相関行列Rxx Aに含まれる雑音成分が抑圧される。第3の相関行列Rxxは、雑音成分が解消されているため、
Rxx=EΛEH (11)
と表される。ここで、Λは固有値を示し、Eは固有ベクトルを示す。
したがって、本実施形態に係る受信装置によれば、受信信号に含まれる雑音が複数のアンテナ間で不均一であっても、正確な結果を取得できる。
Claims (9)
- 雑音成分を含む電波を受信する複数のアンテナと、
前記複数のアンテナで受信された複数の電波を予め設定された周波数帯へ変換するRF変換部と、
前記周波数変換された複数の信号をデジタル信号へ変換するアナログ−デジタル変換部と、
前記複数のデジタル信号についての第1の相関行列を算出する第1の演算部と、
前記複数のデジタル信号についての前記第1の相関行列を算出し、前記算出した第1の相関行列を、変換行列を用いて第2の相関行列に相似変換する第2の演算部と、
前記第1の相関行列と前記第2の相関行列との差を取ることで、前記雑音成分を抑圧した第3の相関行列を作成する雑音抑圧部と、
前記第3の相関行列に対して固有値分解を施し、固有値と固有ベクトルとを算出する固有値分解部と
を具備する受信装置。 - 前記固有値に基づいて波数を推定する波数推定部をさらに具備する請求項1記載の受信装置。
- 前記固有ベクトルと前記推定した波数とを参照し、前記電波の到来方位に関する情報を算出する方位算出部をさらに具備する請求項2記載の受信装置。
- 前記固有ベクトルと前記推定した波数とを用いて、前記電波の到来方位を示す方位スペクトラムを算出する方位算出部と、
前記方位スペクトラムを累積する累積部と
をさらに具備する請求項2記載の受信装置。 - 前記受信した電波のC/N(Carrier-Noise Ratio)を測定する測定部をさらに具備し、
前記第2の演算部は、前記C/Nの値に応じて前記変換行列の値を変化させる請求項1乃至4のいずれかに記載の受信装置。 - 雑音成分を含む電波がデジタル形式に変換された複数のデジタル信号を受信し、
前記複数のデジタル信号についての第1の相関行列を算出し、
前記複数のデジタル信号についての前記第1の相関行列を算出し、前記算出した第1の相関行列を、変換行列を用いて第2の相関行列に相似変換し、
前記第1の相関行列と前記第2の相関行列との差を取ることで、前記雑音成分を抑圧した第3の相関行列を作成し、
前記第3の相関行列に対して固有値分解を施し、固有値と固有ベクトルとを算出する雑音抑圧方法。 - 前記受信した電波のC/N(Carrier-Noise Ratio)を測定し、
前記C/Nの値に応じて前記変換行列の値を変化させる請求項6記載の雑音抑圧方法。 - 雑音成分を含む電波を複数のアンテナにより受信し、前記受信した電波をデジタル形式の複数のデジタル信号へ変換する受信装置で用いられる雑音抑圧プログラムであって、
前記複数のデジタル信号についての第1の相関行列を算出する処理と、
前記複数のデジタル信号についての前記第1の相関行列を算出し、前記算出した第1の相関行列を、変換行列を用いて第2の相関行列に相似変換する処理と、
前記第1の相関行列と前記第2の相関行列との差を取ることで、前記雑音成分を抑圧した第3の相関行列を作成する処理と、
前記第3の相関行列に対して固有値分解を施し、固有値と固有ベクトルとを算出する処理と
を前記受信装置のコンピュータに実行させる雑音抑圧プログラム。 - 前記受信した電波のC/N(Carrier-Noise Ratio)を測定する処理と、
前記C/Nの値に応じて前記変換行列の値を変化させる処理と
を前記コンピュータにさらに実行させる請求項8記載の雑音抑圧プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012204504A JP6038562B2 (ja) | 2012-09-18 | 2012-09-18 | 受信装置、雑音抑圧方法及び雑音抑圧プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012204504A JP6038562B2 (ja) | 2012-09-18 | 2012-09-18 | 受信装置、雑音抑圧方法及び雑音抑圧プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014059225A JP2014059225A (ja) | 2014-04-03 |
JP6038562B2 true JP6038562B2 (ja) | 2016-12-07 |
Family
ID=50615812
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012204504A Expired - Fee Related JP6038562B2 (ja) | 2012-09-18 | 2012-09-18 | 受信装置、雑音抑圧方法及び雑音抑圧プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6038562B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109309506B (zh) * | 2018-10-16 | 2019-12-17 | 南京邮电大学 | 一种基于信号空间特征值分解的无线信号噪声抑制方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3477132B2 (ja) * | 1999-12-27 | 2003-12-10 | 株式会社東芝 | 到来波の入射方位角測定装置 |
JP3596517B2 (ja) * | 2000-12-12 | 2004-12-02 | 松下電器産業株式会社 | 電波到来方向推定装置、電波到来方向推定方法及び指向性可変送受信装置 |
JP3966860B2 (ja) * | 2004-02-12 | 2007-08-29 | 株式会社国際電気通信基礎技術研究所 | 電波到来方向探知装置及び方法 |
JP2006258529A (ja) * | 2005-03-16 | 2006-09-28 | Fujitsu Ten Ltd | 電波到来方向推定装置及び方法 |
JP4593508B2 (ja) * | 2006-03-31 | 2010-12-08 | Kddi株式会社 | 空間角度広がり推定方法および受信装置 |
JP4977849B2 (ja) * | 2006-05-19 | 2012-07-18 | 三菱電機株式会社 | 電波到来方向探知装置 |
JP2010121975A (ja) * | 2008-11-17 | 2010-06-03 | Advanced Telecommunication Research Institute International | 音源定位装置 |
JP5675285B2 (ja) * | 2010-11-10 | 2015-02-25 | 富士通テン株式会社 | レーダ装置 |
-
2012
- 2012-09-18 JP JP2012204504A patent/JP6038562B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2014059225A (ja) | 2014-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8879672B2 (en) | Signal compression apparatus and method for dynamic compressive sensing using sparsity information of input signal | |
JP5659472B2 (ja) | 到来方向推定装置及び方法 | |
JP4962422B2 (ja) | 到来電波方位測定装置、到来電波方位測定方法及び到来電波方位測定プログラム | |
JP2017036990A (ja) | 到来方向推定装置 | |
Park et al. | Online compressive covariance sensing | |
CN109923430B (zh) | 用于进行相位差展开的装置及方法 | |
JP6038562B2 (ja) | 受信装置、雑音抑圧方法及び雑音抑圧プログラム | |
JP5047002B2 (ja) | 波数推定装置 | |
JP4977849B2 (ja) | 電波到来方向探知装置 | |
WO2015173861A1 (ja) | キャリブレーション装置 | |
JP5705066B2 (ja) | パッシブレーダ装置 | |
JP4015989B2 (ja) | 波数推定装置および波数推定方法 | |
JP5289467B2 (ja) | 校正装置 | |
JP2013152112A (ja) | 時間差方位探知装置 | |
JP2006153579A (ja) | 到来波数推定装置および到来波数推定方法 | |
JP5078717B2 (ja) | 入射波数推定装置及び入射波数推定方法 | |
JP2014153206A (ja) | 信号処理装置及び信号処理方法 | |
CN112363109A (zh) | 一种基于宽带强弱临近信号的加权波达方向估计方法 | |
JP2014196957A (ja) | センサ装置、ターゲット応答推定方法、及びセンサ装置用ターゲット応答推定プログラム | |
JP6030012B2 (ja) | 方位測定装置、方位測定プログラム及び方位測定方法 | |
JP5740218B2 (ja) | 受信装置及び背景雑音推定方法 | |
JP5724548B2 (ja) | 方位測定装置、方位測定方法及び方位測定プログラム | |
JP4391771B2 (ja) | 測角装置 | |
JP2014137241A (ja) | 信号処理装置及び信号処理方法 | |
JP4959228B2 (ja) | 到来波数推定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20131219 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20131226 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20140109 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150727 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20160427 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160531 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160726 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20161004 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161102 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6038562 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
LAPS | Cancellation because of no payment of annual fees |