0

周期数列の母関数と離散フーリエ変換について

65
0

この記事はまだあまり整理されていません.

  1. 周期Nの複素数列全体をCZ/N={xxn=xn+Nfor everynZ}とする.
  2. CZ/N上のラグ作用素LLx=(xn1)nZで定義する.

CZ/Nは加法x+y=(xn+yn)nZ,スカラー乗法λx=(λxn)nZに関するN次元ベクトル空間である.

dimCZ/N=Nのみ示す.δCZ/NNZの指示関数,すなわち
δn={1(n0(modN)),0(otherwise)
とする.このとき,集合S={Lmδ0m<N}CZ/Nの基底である.実際
Lmδn=δnm={1(nm(modN)),0(otherwise)
なので,Sの各元は1次独立であり,任意のxCZ/Nx=m=0N1xmLmδと書ける.

要するに,SCNの標準基底を周期的に拡張した集合である.

ここで
x=m=0N1xmLmδ=(m=0N1xmLm)δ
だから,写像Z:CZ/NC[L]xm=0N1xmLmで定義すれば,(Zx)δ=xが成立する.ただしC[L]={P(L)P(z)C[z]}である.

C[z]/zN1は写像z+zN1LによってC[L]と同型である.

写像ϕ:C[z]C[L]P(z)P(L)で定義する.kerϕ=zN1が示せれば,準同型定理よりC[z]/zN1C[L]がしたがう.また,この同型はz+zN1Lを満たす.

P(z)kerϕを任意にとり,P(z)zN1で割った商をQ(z),あまりをR(z)とする.このときP(z)=Q(z)(zN1)+R(z)LN=1よりϕ(P(z))=R(L)=0である.また,degR(z)<NなのでR(z)=m=0N1rmzmとおける.すると
ϕ(R(z))δ=(m=0N1rmLm)δ=m=0N1rm(Lmδ)=0
であり,SCZ/Nの基底だからrm=0である.つまり,P(z)kerϕならR(z)=0P(z)zN1である.明らかにzN1kerϕも成り立つので,kerϕ=zN1である.

Zは全単射である.

Zˇ:C[L]CZ/NP(L)P(L)δで定義する.Zˇ(Zx)=(Zx)δ=xはすでに示したので,Z(ZˇP(L))=P(L)が示せればZˇ=Z1といえる.

上の同型からdegP(z)<Nを仮定してよい.P(z)=m=0N1pmzmとし,mNで割ったあまりをm%Nと書く,すると
ZˇP(L)=m=0N1pmLmδ=(pm%N)mZ,Z(ZˇP(L))=m=0N1pm%NLm=P(L)
である.

以上により,Zは(加法を保つ)全単射である.よって,CZ/Nの乗法をxy=Z1((Zx)(Zy))で定めれば,3つ組(CZ/N,+,)C[L]と同型な環になる.

xyxyの一般項を用いて表そう.
(Zx)(Zy)=(m=0N1xmLm)(n=0N1ynLn)=m=0N1n=0N1xmynLm+n
であり,k=m+nとおくと
(Zx)(Zy)=m=0N1k=mN1+mxmykmLk
となる.k=mN1+mは1周期にわたる和なので,範囲をk=0N1に変えてよい.よって
(Zx)(Zy)=m=0N1k=0N1xmykmLk=k=0N1(m=0N1xmykm)Lk
であり,z=xyとおくと
z=k=0N1(m=0N1xmykm)Lkδ,zk=m=0N1xmykm
となる.

x,yCZ/Nとする.次式で定義されるxyCZ/Nを,xy巡回畳み込みという.また,xyZ(xy)=(Zx)(Zy)を満たす.
xy=z,zk=m=0N1xmykm

z¯=z+zN1とする.同型写像ψ:C[z]/zN1C[L]P(z¯)P(L)を用いて,imZC[z]/zN1へと変えた写像ψ1ZZ~とおく.つまり
Z~:CZ/NC[z]/zN1,(xn)nZn=0N1xnz¯n
である.

写像f:C[z]CNP(z)(P(ζk))k=0N1ζ=e2πi/N)で定義する.明らかに,fは環準同型である.

kerf=zN1である.

P(z)zN1なら,P(z)=Q(z)(zN1)とおくとf(P(z))=(Q(ζk)(ζkN1))k=0N1=0である.またP(z)kerfなら,因数定理よりP(z)k=0N1(zζk)=zN1で割り切れる.

上の命題から,写像f¯:C[z]/zN1CNP(z¯)(P(ζk))k=0N1はwell-definedな環同型である.よって,写像Ff¯Z~で定義すると,F(CZ/N,+,)からCNへの環同型である.実際にx^=Fxを計算すると
Fx=f(n=0N1xnz¯n),x^k=n=0N1xnζkn
となる.これは(符号の違いを除いて)xの離散フーリエ変換である.

投稿日:20221229
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

数学、特に応用数学が好きです。Mathlogでは主に、数学とプログラミングを絡めたようなことを書けたらいいなと思っています。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中