(57)【特許請求の範囲】
【請求項１】
振幅変調された送信信号を受信し、受信信号の搬送波の上側帯波および下側帯波から、前記上側帯波と前記下側帯波との和信号と、前記上側帯波と前記下側帯波との差信号とを生成する生成手段と、
前記和信号と前記差信号とをそれぞれ周波数領域の値に変換する第１の変換手段と、
周波数領域の値とされた前記和信号と周波数領域の値とされた前記差信号との差分を演算する演算手段と、
前記差分を時間領域の値に変換する第２の変換手段と
を有することを特徴とする受信装置。
【請求項２】
請求項１に記載の受信装置において、
前記生成手段は、受信した信号に搬送波の余弦成分と正弦成分とを乗算することにより、前記和信号と前記差信号とを生成する
ことを特徴とする受信装置。
【請求項３】
請求項１に記載の受信装置において、
周波数領域の値とされた前記和信号と周波数領域の値とされた前記差信号とのそれぞれに２乗平均の演算を適用することにより、前記和信号と前記差信号とを平均化する平均化手段をさらに有し、
前記演算手段は、平均化された前記和信号と、平均化された前記差信号との前記差分を演算する
ことを特徴とする受信装置。
【請求項４】
請求項１に記載の受信装置において、
周波数領域の値とされた前記和信号と、周波数領域の値とされた前記差信号との前記差分のゲインを調整する調整手段をさらに有し、
前記第２の変換手段は、ゲインが調整された前記差分を時間領域の値に変換する
ことを特徴とする受信装置。
【請求項５】
振幅変調された送信信号を受信し、受信信号の搬送波の上側帯波および下側帯波から、前記上側帯波と前記下側帯波との和信号と、前記上側帯波と前記下側帯波との差信号とを生成する生成ステップと、
前記和信号と前記差信号とをそれぞれ周波数領域の値に変換する第１の変換ステップと、
周波数領域の値とされた前記和信号と周波数領域の値とされた前記差信号との差分を演算する演算ステップと、
前記差分を時間領域の値に変換する第２の変換ステップと
を含むことを特徴とする受信方法。
【請求項６】
振幅変調された送信信号を受信し、受信信号の搬送波の上側帯波および下側帯波から、前記上側帯波と前記下側帯波との和信号と、前記上側帯波と前記下側帯波との差信号とを生成する生成ステップと、
前記和信号と前記差信号とをそれぞれ周波数領域の値に変換する第１の変換ステップと、
周波数領域の値とされた前記和信号と周波数領域の値とされた前記差信号との差分を演算する演算ステップと、
上記差分を時間領域の値に変換する第２の変換ステップと
を含む処理をコンピュータに行わせるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は受信装置および受信方法、並びにプログラムに関する。
【背景技術】
【０００２】
AM（Amplitude Modulation（振幅変調））放送波は、USB（Upper Sideband（上側帯波））とLSB（Lower Sideband（下側帯波））に同じ情報を乗せて送信される。USBまたはLSB単独でも通信が成り立つ。USBまたはLSBのいずれか一方を単独で用いる通信は、SSB（Single Sideband）として知られている。従来の受信機は、USBとLSBの和信号を受信信号としていた。
【０００３】
従来、上側波帯信号及び下側波帯信号の少なくとも一方を提供する独立的な側波帯回路と、前記下側波帯信号における望ましくない雑音に対する上側波帯信号における望ましくない雑音を表す品質信号を提供する少なくとも１つの品質検出器と、前記品質信号に応答して、該品質信号によって指定される忠実再生の観点から望ましい低い可聴雑音を有する、前記上側波帯信号のみ、前記下側波帯信号のみ、及び上下両側波帯信号のうちの１つを処理して、雑音の低減した前記変調オーディオ信号の実質上忠実な再生となるオーディオ出力信号を提供する少なくとも１つのセレクタとを設けるようにしたものもある（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第２９９５０７７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、USBとLSBの和信号を受信信号とする場合、USBおよびLSBは同じ信号なので、USBとLSBとの和信号を用いることで信号は２倍になる。一方、USBまたはLSBの雑音は、相関が無いので、２の平方根（以下、２の平方根をsqrt（２）とも記載する）倍にとどまり、USBとLSBとの和信号を用いる場合、USB、LSB単独のものよりsqrt（２）倍S/N比（信号対雑音の比）が改善されるのみである。
【０００６】
そこで、本発明は、上記課題を解決すること、すなわち、雑音をより低減できる受信装置および受信方法、並びにプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明の受信装置の一側面は、振幅変調された送信信号を受信し、受信信号の搬送波の上側帯波および下側帯波から、上側帯波と下側帯波との和信号と、上側帯波と下側帯波との差信号とを生成する生成手段と、和信号と差信号とをそれぞれ周波数領域の値に変換する第１の変換手段と、周波数領域の値とされた和信号と周波数領域の値とされた差信号との差分を演算する演算手段と、差分を時間領域の値に変換する第２の変換手段とを有するものとされている。
【０００９】
さらに、本発明の受信装置の一側面は、上述の構成に加えて、生成手段が、受信した信号に搬送波の余弦成分と正弦成分とを乗算することにより、和信号と差信号とを生成するものとされている。
【００１１】
さらに、本発明の受信装置の一側面は、上述の構成に加えて、周波数領域の値とされた和信号と周波数領域の値とされた差信号とのそれぞれに２乗平均の演算を適用することにより、和信号と差信号とを平均化する平均化手段をさらに設け、演算手段が、平均化された和信号と、平均化された差信号との差分を演算するものとされている。
【００１２】
また、本発明の受信装置の一側面は、上述の構成に加えて、周波数領域の値とされた和信号と、周波数領域の値とされた差信号との差分のゲインを調整する調整手段をさらに設け、第２の変換手段が、ゲインが調整された差分を時間領域の値に変換するものとされている。
【００１３】
また、本発明の受信方法の一側面は、振幅変調された送信信号を受信し、受信信号の搬送波の上側帯波および下側帯波から、上側帯波と下側帯波との和信号と、上側帯波と下側帯波との差信号とを生成する生成ステップと、和信号と差信号とをそれぞれ周波数領域の値に変換する第１の変換ステップと、周波数領域の値とされた和信号と周波数領域の値とされた差信号との差分を演算する演算ステップと、差分を時間領域の値に変換する第２の変換ステップとを含むものとされている。
【００１４】
さらに、本発明のプログラムの一側面は、コンピュータに、振幅変調された送信信号を受信し、受信信号の搬送波の上側帯波および下側帯波から、上側帯波と下側帯波との和信号と、上側帯波と下側帯波との差信号とを生成する生成ステップと、和信号と差信号とをそれぞれ周波数領域の値に変換する第１の変換ステップと、周波数領域の値とされた和信号と周波数領域の値とされた差信号との差分を演算する演算ステップと、差分を時間領域の値に変換する第２の変換ステップとを含む処理を行わせるものとされている。
【発明の効果】
【００１５】
本発明の一側面によれば、雑音をより低減できる受信装置および受信方法、並びにプログラムを提供することができる。
【図面の簡単な説明】
【００１６】
【図１】受信装置の生成部の構成の例を示すブロック図である。
【図２】受信装置の雑音低減部の構成の例を示すブロック図である。
【図３】搬送波、上側帯波USB、および下側帯波LSBを示す図である。
【図４】雑音低減の処理を説明するフローチャートである。
【図５】和信号のスペクトルの例を示す図である。
【図６】差信号のスペクトルの例を示す図である。
【図７】和信号から差信号を引いた結果得られるスペクトルの例を示す図である。
【図８】コンピュータのハードウェアの構成例を示すブロック図である。
【発明を実施するための形態】
【００１７】
以下、本発明の一実施の形態の受信装置について、図１〜図８を参照しながら説明する。
【００１８】
以下の説明において、振幅変調された送信信号を受信する場合を例に説明するが、複数のアンテナによって、同じ信号を２つ以上の経路から受信するダイバーシティ受信装置においても、同様である。
【００１９】
図１は、本発明の一実施の形態の受信装置の構成のうち、和信号および差信号を生成する生成部の構成の例を示すブロック図である。生成部は、搬送波抽出部１１、乗算器１２−１および１２−２、並びにLPF（Low Pass Filter）１３−１および１３−２を含むように構成される。
【００２０】
搬送波抽出部１１は、振幅変調された送信信号（放送波）である、１つの送信信号を受信して得られる中間周波数の信号から搬送波成分を取り出す。搬送波抽出部１１は、搬送波の余弦成分（cos成分）と正弦成分（sin成分）とを取り出す。搬送波抽出部１１は、搬送波の余弦成分（cos成分）を乗算器１２−１に供給し、正弦成分（sin成分）を乗算器１２−２に供給する。
【００２１】
乗算器１２−１は、中間周波数の信号に搬送波の余弦成分（cos成分）を乗算し、その結果をLPF１３−１に供給する。LPF１３−１は、いわゆる低域フィルタであり、乗算器１２−１から供給された信号から所定の周波数より高い成分を除去して、その結果を和信号Saとして出力する。
【００２２】
乗算器１２−２は、中間周波数の信号に搬送波の正弦成分（sin成分）を乗算し、その結果をLPF１３−２に供給する。LPF１３−２は、いわゆる低域フィルタであり、乗算器１２−２から供給された信号から所定の周波数より高い成分を除去して、その結果を差信号Sbとして出力する。
【００２３】
すなわち、生成部において、中間周波数から放送波に含まれている搬送波成分が取り出される。搬送波抽出部１１により搬送波のSIN波成分およびCOS波成分が取り出される。乗算器１２−１および１２−２において、SIN波成分およびCOS波成分にそれぞれ中間周波数成分が乗算され、LPF１３−１および１３−２を通すことでCOS側から和信号Saが得られ、SIN側から差信号Sbが得られる。
【００２４】
図２は、本発明の一実施の形態の受信装置の構成のうち、和信号Saおよび差信号Sbを用いて、雑音を低減する雑音低減部の構成の例を示すブロック図である。雑音低減部は、変換部３１、平均化部３２、スペクトル引き算部３３、ゲイン調整部３４、および逆フーリエ変換部３５を含むように構成される。
【００２５】
変換部３１は、和信号Saおよび差信号Sbとのそれぞれに、フーリエ変換の処理を適用する。変換部３１は、和信号Saにフーリエ変換の処理を適用して得られた値Fa(f)を平均化部３２およびゲイン調整部３４に供給する。また、変換部３１は、差信号Sbにフーリエ変換の処理を適用して得られた値Fb(f)を平均化部３２に供給する。
【００２６】
変換部３１は、フーリエ変換部４１−１および４１−２を含む。フーリエ変換部４１−１は、和信号Saにフーリエ変換の処理を適用する。フーリエ変換部４１−２は、差信号Sbにフーリエ変換の処理を適用する。
【００２７】
平均化部３２は、値Fa(f)および値Fb(f)のそれぞれに、２乗平均の演算を適用することにより、値Fa(f)および値Fb(f)を平均化する。平均化部３２は、値Fa(f)および値Fb(f)のそれぞれの平均化によりそれぞれ得られた値Sag(f)および値Sbg(f)をスペクトル引き算部３３に供給する。
【００２８】
平均化部３２は、スペクトル平均化部４２−１および４２−２を含む。スペクトル平均化部４２−１は、値Fa(f)に２乗平均の演算を適用することにより、値Fa(f)を平均化する。スペクトル平均化部４２−２は、値Fb(f)に２乗平均の演算を適用することにより、値Fb(f)を平均化する。
【００２９】
スペクトル引き算部３３は、値Fa(f)の平均化により得られた値Sag(f)から、値Fb(f)の平均化により得られた値Sbg(f)を引き算し、その結果得られた差分をゲイン調整部３４に供給する。
【００３０】
ゲイン調整部３４は、和信号Saにフーリエ変換の処理を適用して得られた値Fa(f)を参照して、値Sag(f)から値Sbg(f)を引き算して得られた差分のゲインを調整し、ゲインが調整された差分を逆フーリエ変換部３５に供給する。逆フーリエ変換部３５は、ゲインが調整された差分に逆フーリエ変換の処理を適用し、その結果得られた値を出力する。
【００３１】
すなわち、雑音低減部は、和信号Saおよび差信号Sbをそれぞれ、フーリエ変換し、周波数成分に分解する。これにより、差信号Sbの雑音成から雑音のスペクトルが得られる。雑音低減部において、雑音のスペクトルを平均化することで正確な値が得られる。
【００３２】
和信号Saの周波数成分は、信号成分と雑音成分の和と考えられ、平均雑音成分は、和信号Saと差信号Sbとにおいて同じであるので、雑音低減部において、差信号Sbの雑音成分が和信号Saからスペクトル上で引き算されることにより雑音の少ないスペクトルが得られる。さらに、雑音低減部において、スペクトルに逆フーリエ変換を施すことにより、雑音を取り除いた信号が得られる。
【００３３】
受信装置の処理を説明する前に、まず、送信装置から送信されてくる送信信号（以下、放送波とも称し、振幅変調された送信信号をAM放送波とも称する）について説明する。
【００３４】
中間周波数におけるAM放送波の信号を単一周波数としても一般性を失わないので、仮に、周波数fの信号で放送（送信）されたとする。
【００３５】
変調度mの送信側の信号は、式（１）で表される。
mcos(2πft)                                                  ・・・（１）
【００３６】
この場合、送信信号s(t)は、式（２）で表される。
s(t)=A{1+mcos(2πft)}cos(2πωt)                              ・・・（２）
ここで、式（２）において、信号の大きさA、信号周波数f、変調周波数ω、および時間tが用いられる。
【００３７】
さらに計算を進めると、式（３）が得られる。
s(t)=Acos(2πωt)+(Am/2)cos(2π(f+ω)t)+(Am/2)cos(2π(ω-f)t)    ・・・（３）
式（３）において、Acos(2πωt)で表される成分は、搬送波と称され、(Am/2)cos(2π(ω+f)t)で表される成分は、上側帯波USBと称され、(Am/2)cos(2π(ω-f)t)で表される成分は、下側帯波LSBと称される
【００３８】
図３は、搬送波、上側帯波USB、および下側帯波LSBを示す図である。図３において、横軸は、周波数を示し、縦軸は、ゲインを示す。よく知られているように、AM放送波は、図３に示されるように、上側帯波USBと下側帯波LSBとに分かれて送信される。
【００３９】
雑音があると仮定し、ここでは上側帯波USBと下側帯波LSBとが別々の大きさになったとする。すると、上側帯波USBを式（４）のように表すことができ、下側帯波LSBを式（５）のように表すことができる。
(Am/2)cos(2π(ω+f)t)―＞Ucos(2π(ω+f)t)                    ・・・（４）
(Am/2)cos(2π(ω-f)t)―＞Lcos(2π(ω-f)t)                      ・・・（５）
【００４０】
次に、図４のフローチャートを参照して、雑音低減の処理を説明する。
【００４１】
まず、ステップＳ１１において、生成部は、１つの信号から和信号Saと差信号Sbとを生成する。すなわち、例えば、生成部は、式（６）に示されるように、中間周波数の信号に、同期回路の搬送波抽出回路で得られた搬送波成分のCOS成分を乗算し、和信号Saを生成する。例えば、搬送波成分は、搬送波にバンドパスフィルタを通すことで得られる。
s(t)cos(2πωt)=Acos2(2πωt)+Ucos(2π(ω+f)t)cos(2πωt)+Lcos(2π(ω-f)t)cos(2πωt)
=A/2cos(4πωt)+(U/2)cos(2π(2ω+f)t)+(L/2)cos(2π(2ω-f)t)
+A/2+(U/2)cos(2πft)+(L/2)cos(2πft)
・・・（６）
【００４２】
ここでローパスフィルタにより変調周波数成分を取り除くことができる。また直流成分のA/2も取り除くことが出来るので、式（７）に示される和信号Saが得られる。
Sa＝(U/2)cos(2πft)+(L/2)cos(2πft)                          ・・・（７）
【００４３】
ノイズがなければ上側帯波USBと下側帯波LSBの大きさは、同じなので、その大きさをGとすれば、式（８）および式（９）が成り立ち、送信信号と一致する。
(U/2)=(L/2)=(G/2)                                          ・・・（８）
Sa=Gcos(2πft)                                            ・・・（９）
【００４４】
また、生成部は、式（１０）に示されるように、中間周波数の信号に、同期回路の搬送波抽出回路で得られた搬送波成分のSIN成分を乗算し、差信号Sbを生成する。
s(t)sin(2πωt)=Acos(2πωt)sin(2πωt)+Ucos(2π(ω+f)t)sin(2πωt)+Lcos(2π(ω-f)t))sin(2πωt)
=(A/2)sin(4πωt)+(U/2)sin(2π(2ω+f)t)+(L/2)sin(2π(2ω-f)t)
-(U/2)sin(2πft)+(L/2)sin(2πft)
・・・（１０）
【００４５】
ここで、ローパスフィルタで変調周波数成分が取り除くことができるので、式（１１）に示される差信号Sbが得られる。
Sb=-(U/2)sin(2πft)+(L/2)sin(2πft)                        ・・・（１１）
【００４６】
ノイズがなければ上側帯波USBと下側帯波LSBの大きさは、同じなので、その大きさをGとすれば、式（１２）および式（１３）が成り立ち、ノイズのないことが示される。
(U/2)=(L/2)=(G/2)                                        ・・・（１２）
Sb=0                                                  ・・・（１３）
【００４７】
以上のことから、差信号Sbによって、和信号Saに含まれている雑音の大きさがわかる。差信号Sbの雑音は、雑音の大きさのみわかり、雑音の位相がわからないことを理解する必要がある。つまり雑音波形そのものはわからないが、その平均的な大きさはわかる。このため、雑音除去するには信号の大きさのみに着目し、その位相成分を考慮しないで行う。すなわち、以下では、雑音を含んでいる信号から、雑音に相当する大きさを減らすことで、雑音感を少なくする処理を行う。
【００４８】
また、雑音は周波数に依存する量であることを考慮する。このため以下では、信号および雑音を周波数に関係する量、周波数の関数として扱う。
【００４９】
式（７）に含まれるcos(2πft)の大きさは、一定なので、その係数Sagに着目し、さらに周波数の関数であることを明示するため式（７）および式（１１）を書き換える。すなわち、式（７）のSaの大きさは、式（１４）で表すことができる。
Sag(f)=U(f)/2+L(f)/2                                  ・・・（１４）
【００５０】
また、式（１１）の差信号Sbの大きさSbgの係数は、式（１５）のように書き換えることができる。
Sbg(f)=U(f)/2-L(f)/2                                  ・・・（１５）
【００５１】
また、ノイズのあるときは上側帯波USBと下側帯波LSBは違ってくるので、上側帯波USBと下側帯波LSBを信号成分とノイズ成分に分けて考える。
U(f)=S(f)+Nu(f)                                      ・・・（１６）
L(f)=S(f)+Nl(f)                                      ・・・（１７）
【００５２】
信号成分S、上側帯波USBのノイズ成分Nu、および下側帯波LSBのノイズ成分Nlから、式（１８）により、和信号Saの大きさが求められ、式（１９）により、差信号Sbの大きさが求められる。
Sag(f)＝U(f)/2+L(f)/2=(2S(f)+Nu(f)+Nl(f))/2              ・・・（１８）
Sbg(f)=-U(f)/2+L(f)/2=-(Nu(f)-N(f)l)/2                    ・・・（１９）
【００５３】
このように、ステップＳ１１において、１つの信号から和信号Saと差信号Sbとが生成される。
【００５４】
次に、ステップＳ１２において、雑音低減部の変換部３１は、和信号Saと差信号Sbとをそれぞれフーリエ変換する。
【００５５】
フーリエ変換後、和信号Saと差信号Sbとは、各周波数に分解され、周波数ごとの成分に分けられる。一般に周波数ごとの成分は、複素数である。
【００５６】
式（２０）で示されるように、値Fa(f)は、和信号Saがフーリエ変換された結果求められる値である。また、式（２１）で示されるように、値Fb(f)は、差信号Sbがフーリエ変換された結果求められる値である。
Sa  ―＞（フーリエ変換）―＞  Fa(f)              ・・・（２０）
Sb  ―＞（フーリエ変換）―＞  Fb(f)              ・・・（２１）
【００５７】
なお、値Fa(f)および値Fb(f)は、周波数の関数であり、値Fa(f)および値Fb(f)に(f)を付したのは、周波数の関数であることを示すためである。
【００５８】
ちなみにフーリエ変換後の絶対値は、式（２２）および式（２３）で表される。なお差信号に絶対値が付くのはスペクトルが負にならないようにするためである。
｜Fa(f)｜＝Sag(f)＝U(f)/2+L(f)/2=(2S(f)+Nu(f)+Nl(f))/2      ・・・（２２）
｜Fb(f)｜=Sbg(f)=｜-U(f)/2+L(f)/2｜=｜(Nu(f)-N(f)l)｜/2    ・・・（２３）
【００５９】
ステップＳ１３において、平均化部３２は、値Fa(f)および値Fb(f)のそれぞれに、２乗平均の演算を適用することにより、値Fa(f)および値Fb(f)を平均化する。すなわち、スペクトルの大きさを求めるためにそれぞれ２乗平均が行われる。信号と、上側帯波USBの雑音と、下側帯波LSBの雑音とが無相関であるので、平均計算する関数Aveを適用することで、無相関の定義から式（２４）が成り立ち、２乗平均はそれぞれの和になる。
【００６０】
Ave(S(f)*Nu(f))=0
Ave(S(f)*Nl(f))=0
Ave(Nu(f)*Nl(f))=0
・・・（２４）
【００６１】
なお、式（２２）または式（２３）に含まれる２で割る部分は本質的でなく、式が煩雑になるので以下では言及しない。
【００６２】
式（２５）が成り立ち、おのおのの２乗和になる。
Ave((2S(f)+Nu(f)+Nl(f))2)＝Ave(2S(f)2)+Ave(Nu(f)2)+Ave(Nl(f)2)
・・・（２５）
【００６３】
また、差信号Sbにおいて、和信号Saと同様に上側帯波USBと下側帯波LSBとが無相関であることを考慮すれば、差信号Sbでも２乗平均はそれぞれの２乗和になる。
Ave((Nu(f)-Nl(f))2)＝Ave(Nu(f)2)+Ave(Nl(f)2)
・・・（２６）
【００６４】
和信号Saの２乗和に含まれるノイズは、差信号Sbの２乗和に等しくなる。このことで和信号Saに含まれるノイズの大きさを推定できる。
【００６５】
この様子を、図５および図６に示す。図５および図６において、黒丸は、信号のスペクトルを示し、矢印は、雑音（ノイズ）のスペクトルを示す。図７において、同様である。
【００６６】
図５に示されるように、和信号Saに含まれるノイズは、信号とノイズの大きさを加算したものとなる。また図６に示されるように、差信号Sbには、ノイズが現れる。
【００６７】
図７に示されるように、和信号Saの２乗和から、差信号Sbの２乗和を引くことでノイズのない信号の２乗和が得られる。しかしながら和信号Saと差信号Sbの雑音は長時間の平均では一致するものの、短時間では一致しないので適当な長さの平均値が用いられる。
【００６８】
すなわち、ステップＳ１４において、スペクトル引き算部３３は、和信号Saのスペクトルから差信号Sbのスペクトルを引き算することで、ノイズのない信号を得る。
【００６９】
式（２５）に示される、和信号Saに含まれるノイズは、式（２６）に示される差信号Sbから推定できる。
【００７０】
以上のことから、AM変調波のノイズ成分は、上側帯波USBと下側帯波LSBの和信号Saと差信号Sbの２乗和に一致する。AM放送において、和信号Saが聴取される。和信号Saの中のノイズ成分を差信号Sbから求め、周波数軸上でノイズに相当する部分を引き算することにより周波数ごとのノイズが無いと仮定したスペクトルが求まる。
【００７１】
このように、和信号Saと差信号Sbのそれぞれに、フーリエ変換を適用すると、値Fa(f)および値Fa(f)が得られるので、和信号Saと差信号Sbのノイズ成分の２乗和は一致することから、式（２７）に示される計算を全ての周波数で行うことで、ノイズ除去された信号の２乗和が分かり、ノイズ除去された信号A(f)が求められる。
｜A(f)｜2＝｜Fa(f)｜2−｜Fb(f)｜2
・・・（２７）
【００７２】
なお雑音が大きいときに（２７）式が負または極端に小さくなるときがある。このときは０または小さい値に置き換えることができる。また（２７）式の変わりに簡便的に絶対値でも計算できる。すなわち式（２８）を用いることができる。
｜A(f)｜ = ｜Fa(f)｜ − ｜Fb(f)｜
・・・（２８）
【００７３】
このように、｜A(f)｜2の平方根を計算することで、各スペクトルの大きさが分かる。
【００７４】
しかしながら、各スペクトルの大きさは分かるがフーリエ変換後の実数、虚数成分は分からないので、和信号の実数、虚数の大きさに配分する。すなわち実数、虚数成分の比率を変えないで、大きさのみノイズ除去された値A(f)になるように実数、虚数に同じ値を乗算する。
【００７５】
この操作は信号の位相を変えることに相当するが、人間の聴覚特性は位相には鈍感なので、違和感は生じない。すなわち和信号のフーリエ変換を式（２９）で書き換えることで雑音の少ない信号が得られる。
Fa(f)＝Fa(f)｜A(f)｜/｜Fa(f)｜
・・・（２９）
【００７６】
すなわち、ステップＳ１５において、ゲイン調整部３４は、和信号Saのスペクトルから差信号Sbのスペクトルが引き算された結果のゲインを調整する。
【００７７】
ステップＳ１６において、逆フーリエ変換部３５は、式（２９）の演算の結果得られたFa(f)の値を逆フーリエ変換し、時間領域の値に戻すことで雑音が除去された信号を得て、雑音低減の処理は終了する。
【００７８】
このように、AM放送波は上側帯波USBと下側帯波LSBの二つに分かれて送信されている。従来の受信機は上側帯波USBと下側帯波LSBの和信号Saを受信信号としているのに対し、上側帯波USBと下側帯波LSBの差信号Sbから雑音情報を得て、この雑音信号から和信号Saに含まれる雑音成分を推定し、和信号Saに含まれる雑音を低減する。AM放送と同様に、同じ信号を2つの経路から受信されるダイバーシティ受信機においても雑音低減ができる。
【００７９】
以上のように、上側帯波USBと下側帯波LSBの差信号Sbから得られる雑音成分で、和信号Saに含まれる雑音成分を推定し、推定した雑音部分を引くことにより雑音を低減することができる。
【００８０】
元々雑音はランダムであり、AM放送波の下側帯波LSB側に含まれている雑音と上側帯波USB側に含まれている雑音はお互いに無関係な無相関の信号である。無相関の信号の和と差は同じ統計的性質を示す。つまり和信号Saの雑音と差信号Sbの雑音はその２乗和の平均値は同じ値を示す。差信号Sbに含まれている雑音の統計的性質を知ることで、和信号Saの雑音を低減することが出来る。すなわち和信号Saに含まれているノイズの量を差信号Sbから推定し、それを引くことで和信号Saに含まれているノイズを低減することができ、雑音の少ない信号を得ることができる。
【００８１】
これにより単に和信号Saがsqrt（２）倍のS/N改善にとどまるのに対して数倍の改善効果がある。また、AM放送の上側帯波USBと下側帯波LSBと同じように、別々の経路から同じ信号を送信するダイバーシティ受信にも同じように雑音が低減できる。
【００８２】
以上のように、雑音をより低減できる。
【００８３】
上述した一連の処理は、ハードウェアにより実行することもできるし、ソフトウエアにより実行することもできる。一連の処理をソフトウエアにより実行する場合には、そのソフトウエアを構成するプログラムが、専用のハードウェアに組み込まれているコンピュータ、または、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のパーソナルコンピュータなどに、プログラム記録媒体からインストールされる。
【００８４】
図８は、上述した一連の処理をプログラムにより実行するコンピュータのハードウェアの構成例を示すブロック図である。
【００８５】
コンピュータにおいて、CPU（Central Processing Unit）２０１，ROM（Read Only Memory）２０２，RAM（Random Access Memory）２０３は、バス２０４により相互に接続されている。
【００８６】
バス２０４には、さらに、入出力インタフェース２０５が接続されている。入出力インタフェース２０５には、キーボード、マウス、マイクロホンなどよりなる入力部２０６、ディスプレイ、スピーカなどよりなる出力部２０７、ハードディスクや不揮発性のメモリなどよりなる記憶部２０８、ネットワークインタフェースなどよりなる通信部２０９、磁気ディスク、光ディスク、光磁気ディスク、或いは半導体メモリなどのリムーバブルメディア２１１を駆動するドライブ２１０が接続されている。
【００８７】
以上のように構成されるコンピュータでは、CPU２０１が、例えば、記憶部２０８に記憶されているプログラムを、入出力インタフェース２０５及びバス２０４を介して、RAM２０３にロードして実行することにより、上述した一連の処理が行われる。
【００８８】
コンピュータ（CPU２０１）が実行するプログラムは、例えば、磁気ディスク（フレキシブルディスクを含む）、光ディスク（CD-ROM(Compact Disc-Read Only Memory),DVD(Digital Versatile Disc)等）、光磁気ディスク、もしくは半導体メモリなどよりなるパッケージメディアであるリムーバブルメディア２１１に記録して、あるいは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の伝送媒体を介して提供される。
【００８９】
そして、プログラムは、リムーバブルメディア２１１をドライブ２１０に装着することにより、入出力インタフェース２０５を介して、記憶部２０８に記憶することで、コンピュータにインストールすることができる。また、プログラムは、有線または無線の伝送媒体を介して、通信部２０９で受信し、記憶部２０８に記憶することで、コンピュータにインストールすることができる。その他、プログラムは、ROM２０２や記憶部２０８にあらかじめ記憶しておくことで、コンピュータにあらかじめインストールしておくことができる。
【００９０】
なお、コンピュータが実行するプログラムは、本明細書で説明する順序に沿って時系列に処理が行われるプログラムであっても良いし、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで処理が行われるプログラムであっても良い。
【００９１】
また、本発明の実施の形態は、上述した実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において種々の変更が可能である。
【符号の説明】
【００９２】
１１…搬送波抽出部、１２−１および１２−２…乗算器、１３−１および１３−２…LPF、３１…変換部、３２…平均化部、３３…スペクトル引き算部、３４…ゲイン調整部、３５…逆フーリエ変換部、４１−１および４１−２…フーリエ変換部、４２−１および４２−２…スペクトル平均化部、２０１…CPU、２０２…ROM、２０３…RAM、２０８…記憶部、２０９…通信部、２１１…リムーバブルメディア
【図１】
0005658975.tif 000002
【図２】
0005658975.tif 000003
【図３】
0005658975.tif 000004
【図４】
0005658975.tif 000005
【図５】
0005658975.tif 000006
【図６】
0005658975.tif 000007
【図７】
0005658975.tif 000008
【図８】
0005658975.tif 000009