0
大学数学基礎解説
文献あり

実離散フーリエ変換(前編)-基底変換-

59
0

自然数Nは0を含まないものとする。
連続フーリエ変換は連続かつ非周期関数のフーリエ変換を指し、フーリエ変換は広義の、時間の関数から周波数の関数への変換を指すものとする。

はじめに

実離散フーリエ変換とは

離散フーリエ変換は離散時間の関数を離散周波数の関数に変換するものです。
有限のデータを扱うため、音声や画像などのデータ処理において重宝されます。

また、フーリエ変換は三角関数や複素指数関数の基底変換(後述)として捉えることができます。
この場合、連続関数や無限次元ベクトル空間であれば収束条件や積分と極限の入れ替えなどの解析的な議論が必要となるため、無限次元のフーリエ級数展開から連続フーリエ変換を導出するよりも、有限次元の離散フーリエ変換からフーリエ変換を導出した方が分かりやすいと思われます。

一方で、一般的に離散フーリエ変換は複素数を用いて以下のように表現されます。

複素離散フーリエ変換

関数f(x)の複素離散フーリエ変換は
F(m)=n=0N1f(n)exp(imNnτ)
複素離散逆フーリエ変換は
f(n)=1Nm=0N1F(m)exp(imNn)

フーリエ級数展開や連続フーリエ変換が実数から複素数に拡張したように、離散フーリエ変換も実数から考える方が直感的理解につながると思われます。

そこで、このシリーズでは離散フーリエ変換の実数表現、つまり実離散フーリエ変換を、基底変換として導出し、最後にフーリエ級数展開と複素離散フーリエ変換の対応について見ていきます。

このシリーズの道筋

  1. 数ベクトル空間・基底変換とは何かを説明する。(前編)
     正規直交基底により任意のベクトルが展開できることなどを示す。
  2. 離散三角関数の直交性を示す。(中編)
     前編の議論から離散三角関数が正規直交基底になることが示せる。
  3. 正規直交基底の離散三角関数から実離散フーリエ変換を導出し、フーリエ級数展開や複素離散フーリエ変換との対応を調べる。(後編)

数ベクトルと直交

数ベクトル

Rn={[x1x2xn]|x1,x2,,xnR}とする。

a=[a1an], b=[b1bn]Rn, kRnに対して,和と実数倍を

a+b=[a1+b1an+bn], ka=[ka1kan], kb=[kb1kbn]

と定義したとき、Rnを実n次元数ベクトル空間または単に数ベクトル空間といい、その各元を実n次元数ベクトル、または単に数ベクトルという。

0=[00]零ベクトルという

これから説明することは一般のベクトル空間でも定義できますが、情報を減らすために実離散フーリエ変換に関わる実n次元数ベクトル空間に絞って説明します。

一次結合・スパン・一次独立・一次従属

数ベクトルa1,,amRnと実数c1,,cmRに対して
k=1mckak=c1a1+cnan
a1,,am一次結合という。
a1,,amの一次結合で表されるベクトル全体の集合をa1,,amスパンといい、
span(a1,,am)={k=1mckak|a1,,amRn, c1,,cmR}
と表す。また、

k=1mckak=0c1==cn=0
となるとき、a1,,am一次独立であるという。

一次独立の一意性

一次独立な数ベクトルa1,,anに対して
k=1nckak=k=1ndkakck=dk (1kn, kN)
となる。つまり一次独立なら一次結合の表し方は一通りしかない。

仮定より
k=1nckakk=1ndkak=k=1n(ckdk)ak=0
a1,,anは一次独立であるから定義よりckdk=0
よってck=dk

一次独立と零ベクトル

数ベクトルa1,,anが一次独立ならばam0 (0mn)

数ベクトルa1,,anが一次独立であるとする。
aj=0となるjが少なくとも一つ存在すると仮定すると、
k=1nckak=0
としたときにcj0が解の一つになるから、a1,,anが一次独立であることと矛盾する。よって、
am0 (0mn)

内積・ノルム

数ベクトルa=(a1,an), b=(b1,bn)Rnに対して
ab=i=1naibi
ab内積という。また、
a=aa=i=1nai2
aノルムという。
ab=0となるとき、ab直交するという。

直交系・正規直交系

数ベクトルa1,,anが互いに直交するとき、つまり
ijaiaj=0
となるとき、{a1,,an}直交系という。
また、基底a1,,anが互いに直交してそれぞれのノルムが1のとき、つまり
aiaj={0(ij)1(i=j)
となるとき、{a1,,an}正規直交系という。

直交と線形独立

数ベクトルa1,,anが直交系でそれぞれ零ベクトルでなければ線形独立である。

k=1nckak=0とおき、係数ck (1kn)が全て0になることを示す。

一般に任意のm (1mn)に対してam0=0となるから
0=am0=amk=1nckak=k=1nck(amak)

数ベクトルa1,,anが互いに直交するとき、kmならばamak=0であるため、

0=k=1nck(amak)=cm(umum)
仮定よりamは零ベクトルでないため、(amam)0より
ck=0(1kn)
よってa1,,anは線形独立である。

基底変換

基底

数ベクトルa1,,anRn

  1. 一次独立かつ
  2. 任意の数ベクトルxRnをその一次結合で表せる、
    つまりa1,,anのスパンがRnとなるとき

{a1,,an}Rn基底という。

基底と一意性

基底の一次結合による数ベクトルの表し方は一通りしかない。

基底は一次独立であるから補題1より明らか。

基底変換

任意の数ベクトルを、ある基底の一次結合から別の基底の一次結合で表すとき、これを基底変換という。

直交基底・正規直交基底

{a1,,an}が直交系かつ基底となるとき、{a1,,an}直交基底という。
また、{a1,,an}が正規直交系かつ基底となるとき、{a1,,an}正規直交基底という。

基底による展開

任意のベクトルxRnは、Rnの正規直交基底{a1,,an}によって次のように表せる。
x=i=1n(aix)ai
このように与えられたベクトルxを正規直交基底{a1,,an}の一次結合で表すことを、xa1,,anによる展開という。

{a1,,an}Rnの基底であるため、任意のベクトル xRnを次のように線形結合で表せる。
x=i=1nciai(c1,,cmR)
さらに{a1,,an}は正規直交基底であるため、aj (1jn)xの内積は次のようになる。
ajx=aji=1nciai=i=1nci(ajai)=i=1nciδj,i=cj
よってci=aix (1in)となるからxは次のように表せる。
x=i=1nciai=i=1n(aix)ai

取り換え定理

数ベクトルa1,a2,,aMのスパンをVとする。つまり

V=span(a1,a2,,aM)

b1,b2,,bNV (NM)が一次独立のとき、a1,a2,,aMの内、あるN個をb1,b2,,bNに置き換えてもそれらのスパンがVとなる。
つまり、1iN+1<<iMMとするとV=span(b1,,bN,aiN+1,aiM)

数学的帰納法で示す。

まずN=1の場合を考える。
b1が一次独立であるから補題2よりb10
a1,a2,,aMのスパンはVだからb1Vを次のように表せる。
b1=m=1Mcmam

仮にcm=0 (1mM)とするとb10と矛盾するからcm0となるmが存在する。
ck0とすると、k以外の番号を1i2<<iMMとなるように定める。
このとき、a1を次のようにb1aim(2mM)の線形結合で表せる。

b1=ckak+m=2Mcimaimak=1ck(b1m=2Mcimaim)

よって

V=span(a1,a2,,aM)=span(ak,ai2,,aiM)=span(b1,ai2,,aiM,ai2,,aiM)=span(b1,ai2,,aiM)
より、b1akを置き換えられた。

次にN2N1個のときに成り立つとして、N個でも成り立つことを示す。仮定より
bN=span(b1,,bN1,aiN,aiM)=n=1N1(cnbn)+m=NM(dimaim)

dim=0 (NmM)と仮定すると、
bN=n=1N1(cnbn)bNn=1N1(cnbn)=0
これはb1,,bNが一次独立であることと矛盾するため、dk=0となるkが少なくとも一つ存在する。
kを除くiNからiM1jN+1<<jMMとなるように番号を定めると、akを次のようにbn(1nN)ajm(N+1mM)の線形結合で表せる。

bN=n=1N1(cnbn)+dkak+m=N+1M(dimaim)ak=1dk(bNn=1N1(cnbn)m=N+1M(dimaim))
よって
V=span(b1,,bN1,aiN,,aiM)=span(b1,,bN1,ak,ajN+1,,ajM)=span(b1,,bN1,bN,b1,,bN1,ajN+1,,ajM,ajN+1,,ajM)=span(b1,,bN1,bN,ajN+1,,ajM)
よってbNakを置き換えられた。

数学的帰納法より、全ての1NMで示された。

線形独立と基底

N個のN次ベクトルは線形独立ならN次ベクトル空間の基底になる。

定理5の取り換え定理でamN次基底ベクトルとしてbnN個のN次ベクトルとすると、基底になることが分かる。

直交と基底

N個のN次ベクトルは互いに直交するならN次ベクトル空間の基底になる。

補題3と補題7より明らか

中編で離散三角関数の直交性を示し、命題7より基底となることを示します。

参考文献

投稿日:2024818
更新日:2024819
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

sokia
sokia
0
139

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 実離散フーリエ変換とは
  3. このシリーズの道筋
  4. 数ベクトルと直交
  5. 参考文献