• 検索結果がありません。

信号処理とフーリエ変換第 4 回

N/A
N/A
Protected

Academic year: 2021

シェア "信号処理とフーリエ変換第 4 回"

Copied!
48
0
0

読み込み中.... (全文を見る)

全文

(1)

信号処理とフーリエ変換 第 4 回

最短距離

垂直

, Fourier

級数の部分和は直交射影かつ最良近似

かつらだ

桂田 祐史

ま さ し

2020

10

14

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010141 / 21

(2)

目次

1

本日の内容・連絡事項

2 Fourier

級数

Fourier

級数の部分和は直交射影かつ最良近似

垂線の足

(

直交射影

)

は最も近い点

Bessel

の不等式

完全系

, Parseval

の等式

,

内積空間の収束 演習

3

課題

1

について

かつらだまさし

(3)

本日の内容・連絡事項

講義ノート

[1]

§1.4

の部分

(

最短距離

垂直

)

の内容を講義し ます。

次回

(11

21

)

にレポート課題

1

を出します

(

締め切りは

11

11

15:20

の予定

)

。課題

1

のうち今回の授業範囲である部分について

は、

10

14

15:20

に公開します

(

動画

part 7)

。また課題文は

http://nalab.mind.meiji.ac.jp/~mk/lecture/fourier-2020/kadai1.pdf

におくようにします。

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010143 / 21

(4)

1.4 Fourier

級数の部分和は直交射影かつ最良近似 復習と予告

(

復習

)§1.3

で、周期

の区分的

C1

級関数の全体

X

に内積

(f,g) :=

Z π

π

f(x)g(x)dx

を導入し、

Fourier

級数に用いる関数系

(cosmx

sinnx,einx)

の直交性を表した。 例えば

(m,nZ0)m̸=n(cosmx,cosnx) = 0, (m,nN)m̸=n(sinmx,sinnx) = 0, mZ0nN(cosmx,sinnx) = 0.

また一般の内積空間

X

で、直交系

n}

による展開の係数を表す公式を得た。

f =X

n

cnφn cn= (f, φn) n, φn). (

今日の予告

)

次のことが成り立つことを示す。

Fourier

級数の部分和は直交射影と呼ばれるものになっている。

(

一般の内積空間で

)

直交射影は、ノルム

φ=p

(φ, φ)

で測って最良近似である。

(

2

つの結果

) Fourier

級数の部分和は最良近似である。 また、有名な

Bessel

の不等式を示し、ノルム

φ=p

(φ, φ)

についての収束と完全

系の概念を定義する。

かつらだまさし

(5)

1.4 Fourier

級数の部分和は直交射影かつ最良近似 復習と予告

(

復習

)§1.3

で、周期

の区分的

C1

級関数の全体

X

に内積

(f,g) :=

Z π

π

f(x)g(x)dx

を導入し、

Fourier

級数に用いる関数系

(cosmx

sinnx,einx)

の直交性を表した。

例えば

(m,nZ0)m̸=n(cosmx,cosnx) = 0, (m,nN)m̸=n(sinmx,sinnx) = 0, mZ0nN(cosmx,sinnx) = 0.

また一般の内積空間

X

で、直交系

n}

による展開の係数を表す公式を得た。

f =X

n

cnφn cn= (f, φn) n, φn). (

今日の予告

)

次のことが成り立つことを示す。

Fourier

級数の部分和は直交射影と呼ばれるものになっている。

(

一般の内積空間で

)

直交射影は、ノルム

φ=p

(φ, φ)

で測って最良近似である。

(

2

つの結果

) Fourier

級数の部分和は最良近似である。 また、有名な

Bessel

の不等式を示し、ノルム

φ=p

(φ, φ)

についての収束と完全

系の概念を定義する。

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010144 / 21

(6)

1.4 Fourier

級数の部分和は直交射影かつ最良近似 復習と予告

(

復習

)§1.3

で、周期

の区分的

C1

級関数の全体

X

に内積

(f,g) :=

Z π

π

f(x)g(x)dx

を導入し、

Fourier

級数に用いる関数系

(cosmx

sinnx,einx)

の直交性を表した。

例えば

(m,nZ0)m̸=n(cosmx,cosnx) = 0, (m,nN)m̸=n(sinmx,sinnx) = 0, mZ0nN(cosmx,sinnx) = 0.

また一般の内積空間

X

で、直交系

n}

による展開の係数を表す公式を得た。

f =X

n

cnφn cn= (f, φn) n, φn).

(

今日の予告

)

次のことが成り立つことを示す。

Fourier

級数の部分和は直交射影と呼ばれるものになっている。

(

一般の内積空間で

)

直交射影は、ノルム

φ=p

(φ, φ)

で測って最良近似である。

(

2

つの結果

) Fourier

級数の部分和は最良近似である。 また、有名な

Bessel

の不等式を示し、ノルム

φ=p

(φ, φ)

についての収束と完全

系の概念を定義する。

かつらだまさし

(7)

1.4 Fourier

級数の部分和は直交射影かつ最良近似 復習と予告

(

復習

)§1.3

で、周期

の区分的

C1

級関数の全体

X

に内積

(f,g) :=

Z π

π

f(x)g(x)dx

を導入し、

Fourier

級数に用いる関数系

(cosmx

sinnx,einx)

の直交性を表した。

例えば

(m,nZ0)m̸=n(cosmx,cosnx) = 0, (m,nN)m̸=n(sinmx,sinnx) = 0, mZ0nN(cosmx,sinnx) = 0.

また一般の内積空間

X

で、直交系

n}

による展開の係数を表す公式を得た。

f =X

n

cnφn cn= (f, φn) n, φn). (

今日の予告

)

次のことが成り立つことを示す。

Fourier

級数の部分和は直交射影と呼ばれるものになっている。

(

一般の内積空間で

)

直交射影は、ノルム

φ=p

(φ, φ)

で測って最良近似である。

(

2

つの結果

) Fourier

級数の部分和は最良近似である。

また、有名な

Bessel

の不等式を示し、ノルム

φ=p

(φ, φ)

についての収束と完全

系の概念を定義する。

かつらだ

桂 田 まさし

祐 史 信号処理とフーリエ変換 第4202010144 / 21

(8)

1.4.1 垂線の足 ( 直交射影 ) は最も近い点

直線

=V

とその上にない点

F.

上の動点

G.

平面

π=V

とその上にない点

F. π

上の動点

G.

1:V =

2:V =π

V

上の点

G

をなるべく

F

に近くしたい。FG が最短距離となる

G (それを H

と書く) を求められるか?

F

から

V

に引いた垂線と

V

との交点

H

F

から

V

に下ろした垂線の足

H

」あるいは「

F

V

への直交射影」という。 論理を

(

短く

)

言い切ると「最短

垂直」

かつらだまさし

(9)

1.4.1 垂線の足 ( 直交射影 ) は最も近い点

直線

=V

とその上にない点

F.

上の動点

G.

平面

π=V

とその上にない点

F. π

上の動点

G.

1:V =

2:V =π

V

上の点

G

をなるべく

F

に近くしたい。FG が最短距離となる

G (それを H

と書く) を求められるか?

F

から

V

に引いた垂線と

V

との交点

H

F

から

V

に下ろした垂線の足

H

」あるいは「

F

V

への直交射影」という。

論理を

かつらだ(

短く

)

言い切ると「最短

垂直」

桂 田 まさし

祐 史 信号処理とフーリエ変換 第4202010145 / 21

(10)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 定理にする

定理

4.1 (

垂直

最短

)

K

上の内積空間

X

の部分空間

V

f X,hV

対して

h

f

V

への直交射影

h

が最も

f

に近い

つまり次の

2

つが成り立つ。ただし

∥ ∥

は内積から定まるノルムである。

(1) f hV

となる

hV

があれば、

f h= inf

gVf g.

(2) f h= inf

gVf g

となる

h(

最短距離を達成する

h)

があれば

f hV.

かつらだまさし

(11)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 定理の証明

(1)

の証明

(直角三角形FGH

の図を描く)

g V

に対して

f h2+gh2=f g2 (

ピタゴラスの定理

)

であるから

f h∥ ≤ ∥f g.

(2)

の証明

v:=g h

とおく。関数

f(t) :=f (h+tv)2 (t K).

t= 0

で最小値を与える。これから実は

(f h,v) = 0

が導かれる。 ここでは

K=R

の場合のみ紹介する。

f(t) = ((f h)tv,(f h)tv)

=f h22(f h,v)t+t2v2

t= 0

のとき最小値をとるので、1 次の項の係数

2(f h,v)

0

である。

(K=C

の場合の証明は講義ノートを見よ。)

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010147 / 21

(12)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 定理の証明

(1)

の証明

(直角三角形FGH

の図を描く)  

g V

に対して

f h2+gh2=f g2 (

ピタゴラスの定理

)

であるから

f h∥ ≤ ∥f g.

(2)

の証明

v:=g h

とおく。関数

f(t) :=f (h+tv)2 (t K).

t= 0

で最小値を与える。これから実は

(f h,v) = 0

が導かれる。 ここでは

K=R

の場合のみ紹介する。

f(t) = ((f h)tv,(f h)tv)

=f h22(f h,v)t+t2v2

t= 0

のとき最小値をとるので、1 次の項の係数

2(f h,v)

0

である。

(K=C

の場合の証明は講義ノートを見よ。)

かつらだまさし

(13)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 定理の証明

(1)

の証明

(直角三角形FGH

の図を描く)  

g V

に対して

f h2+gh2=f g2 (

ピタゴラスの定理

)

であるから

f h∥ ≤ ∥f g.

(2)

の証明

v:=g h

とおく。関数

f(t) :=f (h+tv)2 (t K).

t= 0

で最小値を与える。これから実は

(f h,v) = 0

が導かれる。

ここでは

K=R

の場合のみ紹介する。

f(t) = ((f h)tv,(f h)tv)

=f h22(f h,v)t+t2v2

t= 0

のとき最小値をとるので、1 次の項の係数

2(f h,v)

0

である。

(K=C

の場合の証明は講義ノートを見よ。)

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010147 / 21

(14)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 定理の証明

(1)

の証明

(直角三角形FGH

の図を描く)  

g V

に対して

f h2+gh2=f g2 (

ピタゴラスの定理

)

であるから

f h∥ ≤ ∥f g.

(2)

の証明

v:=g h

とおく。関数

f(t) :=f (h+tv)2 (t K).

t= 0

で最小値を与える。これから実は

(f h,v) = 0

が導かれる。

ここでは

K=R

の場合のみ紹介する。

f(t) = ((f h)tv,(f h)tv)

=f h22(f h,v)t+t2v2

t= 0

のとき最小値をとるので、1 次の項の係数

2(f h,v)

0

である。

(K=C

の場合の証明は講義ノートを見よ。)

かつらだまさし

(15)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

4.2 (

直交射影を表す式

)

K

上の内積空間

X

の部分空間

V

が直交系

{φn}Nn=1

で張られている、つまり

V =spanφ1, φ2,· · ·, φNdef.= ( N

X

n=1

cnφn

c1,· · ·,cNK )

であれば、任意の

f X

に対して、f に最も近い

hV

(一意的に存在して)

(1) h=

XN n=1

(f, φn) n, φn)φn.

これは

f

V

への直交射影でもある。

V

の範疇で、

f

を近似するものを選ぶ、と考えると、

(1)

h

は、誤差

f h

を最小にするので、「最良近似」と呼ぶのにふさわしい。

(Fourier

級数の部分和

sN

(1)

h

の形をしている。すなわち、Fourier 級数 の部分和は

(V

の範囲内で) 最も

f

に近い。)

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010148 / 21

(16)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

4.2 (

直交射影を表す式

)

K

上の内積空間

X

の部分空間

V

が直交系

{φn}Nn=1

で張られている、つまり

V =spanφ1, φ2,· · ·, φNdef.= ( N

X

n=1

cnφn

c1,· · ·,cNK )

であれば、任意の

f X

に対して、f に最も近い

hV

(一意的に存在して)

(1) h=

XN n=1

(f, φn) n, φn)φn.

これは

f

V

への直交射影でもある。

V

の範疇で、

f

を近似するものを選ぶ、と考えると、

(1)

h

は、誤差

f h

を最小にするので、「最良近似」と呼ぶのにふさわしい。

(Fourier

級数の部分和

sN

(1)

h

の形をしている。すなわち、Fourier 級数 の部分和は

(V

の範囲内で) 最も

f

に近い。)

かつらだまさし

(17)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

4.2 (

直交射影を表す式

)

K

上の内積空間

X

の部分空間

V

が直交系

{φn}Nn=1

で張られている、つまり

V =spanφ1, φ2,· · ·, φNdef.= ( N

X

n=1

cnφn

c1,· · ·,cNK )

であれば、任意の

f X

に対して、f に最も近い

hV

(一意的に存在して)

(1) h=

XN n=1

(f, φn) n, φn)φn.

これは

f

V

への直交射影でもある。

V

の範疇で、

f

を近似するものを選ぶ、と考えると、

(1)

h

は、誤差

f h

を最小にするので、「最良近似」と呼ぶのにふさわしい。

(Fourier

級数の部分和

sN

(1)

h

の形をしている。すなわち、Fourier 級数 の部分和は

(V

の範囲内で) 最も

f

に近い。)

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010148 / 21

(18)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

証明

hV

であるから、ある

c1,· · · ,cN

が存在して

h= XN n=1

cnφn.

上で述べたことから、f

h

V

と直交するので

(f h, φn) = 0 (n= 1,2,· · · ,N).

ゆえに

(f, φn) = (h, φn)

=

XN

j=1

cjφj, φn

= XN j=1

cjj, φn) =cnn, φn).

(4

つめの等号は前回説明済み) ゆえに

cn= (f, φn) n, φn).

かつらだまさし

(19)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

証明

hV

であるから、ある

c1,· · · ,cN

が存在して

h= XN n=1

cnφn.

上で述べたことから、f

h

V

と直交するので

(f h, φn) = 0 (n= 1,2,· · · ,N).

ゆえに

(f, φn) = (h, φn)

=

XN

j=1

cjφj, φn

= XN j=1

cjj, φn) =cnn, φn).

(4

つめの等号は前回説明済み) ゆえに

cn= (f, φn) n, φn).

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010149 / 21

(20)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

証明

hV

であるから、ある

c1,· · · ,cN

が存在して

h= XN n=1

cnφn.

上で述べたことから、f

h

V

と直交するので

(f h, φn) = 0 (n= 1,2,· · · ,N).

ゆえに

(f, φn) = (h, φn)

=

XN

j=1

cjφj, φn

= XN

j=1

cjj, φn) =cnn, φn).

(4

つめの等号は前回説明済み) ゆえに

cn= (f, φn) n, φn).

かつらだまさし

(21)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

証明

hV

であるから、ある

c1,· · · ,cN

が存在して

h= XN n=1

cnφn.

上で述べたことから、f

h

V

と直交するので

(f h, φn) = 0 (n= 1,2,· · · ,N).

ゆえに

(f, φn) = (h, φn) =

XN

j=1

cjφj, φn

= XN

j=1

cjj, φn) =cnn, φn).

(4

つめの等号は前回説明済み) ゆえに

cn= (f, φn) n, φn).

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010149 / 21

(22)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

証明

hV

であるから、ある

c1,· · · ,cN

が存在して

h= XN n=1

cnφn.

上で述べたことから、f

h

V

と直交するので

(f h, φn) = 0 (n= 1,2,· · · ,N).

ゆえに

(f, φn) = (h, φn) =

XN

j=1

cjφj, φn

= XN

j=1

cjj, φn)

=cnn, φn).

(4

つめの等号は前回説明済み) ゆえに

cn= (f, φn) n, φn).

かつらだまさし

(23)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

証明

hV

であるから、ある

c1,· · · ,cN

が存在して

h= XN n=1

cnφn.

上で述べたことから、f

h

V

と直交するので

(f h, φn) = 0 (n= 1,2,· · · ,N).

ゆえに

(f, φn) = (h, φn) =

XN

j=1

cjφj, φn

= XN

j=1

cjj, φn) =cnn, φn).

(4

つめの等号は前回説明済み)

ゆえに

cn= (f, φn) n, φn).

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第4202010149 / 21

(24)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点 直交射影の式

証明

hV

であるから、ある

c1,· · · ,cN

が存在して

h= XN n=1

cnφn.

上で述べたことから、f

h

V

と直交するので

(f h, φn) = 0 (n= 1,2,· · · ,N).

ゆえに

(f, φn) = (h, φn) =

XN

j=1

cjφj, φn

= XN

j=1

cjj, φn) =cnn, φn).

(4

つめの等号は前回説明済み) ゆえに

cn= (f, φn) n, φn).

かつらだまさし

(25)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点

Fourier

級数

4.3 (Fourier

級数の部分和は直交射影かつ最良近似)

(三角関数を用いた) Fourier

級数の部分和

sN =a0

2 + XN n=1

(ancosnx+bnsinnx)

は、f の

V :=span1,cosx,sinx,cos 2x,sin 2x,· · ·,cosNx,sinNx

への直交射 影であり、f の

V

における最良近似でもある

(f sN

が最も短いという意味)。

(複素指数関数を用いた) Fourier

級数の部分和

sN= XN n=N

cneinx

は、f の

V :=spaneiNx,· · ·,eix,ei0x,eix,· · ·,eiNx

への直交射影であり、f の

V

における最良近似でもある

(やはりf sN

が最も短いという意味)。

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第42020101410 / 21

(26)

1.4.1 直交射影 ( 垂線の足 ) は最も近い点

Fourier

級数

4.3 (Fourier

級数の部分和は直交射影かつ最良近似)

(三角関数を用いた) Fourier

級数の部分和

sN =a0

2 + XN n=1

(ancosnx+bnsinnx)

は、f の

V :=span1,cosx,sinx,cos 2x,sin 2x,· · ·,cosNx,sinNx

への直交射 影であり、f の

V

における最良近似でもある

(f sN

が最も短いという意味)。

(複素指数関数を用いた) Fourier

級数の部分和

sN= XN n=N

cneinx

は、f の

V :=spaneiNx,· · · ,eix,ei0x,eix,· · ·,eiNx

への直交射影であり、f の

V

における最良近似でもある

(やはりf sN

が最も短いという意味)。

かつらだまさし

(27)

1.4.2 Bessel の不等式

有名な

Bessel

の不等式を紹介する。

命題

4.4 (Bessel

の不等式)

内積空間

X

の直交系

{φn}Nn=1

と任意の

f X

に対して

(2)

XN n=1

|(f, φn)|2

φn2 ≤ ∥f2

が成り立つ。無限個の場合は

(極限を取って)

(3)

X n=1

|(f, φn)|2

φn2 ≤ ∥f2 (Bessel

の不等式).

特に

{ψn}

が正規直交系の場合は

(4)

XN n=1

|(f, ψn)|2≤ ∥f2, X n=1

|(f, ψn)|2≤ ∥f2.

かつらだ 桂 田

まさし

祐 史 信号処理とフーリエ変換 第42020101411 / 21

(28)

1.4.2 Bessel の不等式

証明

0,h,f

を頂点とする直角三角形の図を描く。ピタゴラスの定理から

h2+f h2=f2.

ゆえに

h2≤ ∥f2.

左辺はやはりピタゴラスの定理より

h2=

XN n=1

(f, φn) n, φn)φn

2

= XN n=1

(f, φn) n, φn)φn

2

= XN n=1

(f, φn) n, φn)

2φn2= XN n=1

|(f, φn)|2

φn2 .

ゆえに

XN n=1

|(f, φn)|2

φn2 ≤ ∥f2.

かつらだまさし

参照

関連したドキュメント

日頃から製造室内で行っていることを一般衛生管理計画 ①~⑩と重点 管理計画

しかし何かを不思議だと思うことは勉強をする最も良い動機だと思うので,興味を 持たれた方は以下の文献リストなどを参考に各自理解を深められたい.少しだけ案

Maurer )は,ゴルダンと私が以前 に証明した不変式論の有限性定理を,普通の不変式論

Maurer )は,ゴルダンと私が以前 に証明した不変式論の有限性定理を,普通の不変式論

LLVM から Haskell への変換は、各 LLVM 命令をそれと 同等な処理を行う Haskell のプログラムに変換することに より、実現される。

【ご注意点】 ・カタログの中からお好みの商品を1点お 選びいただき、同封のハガキに記載のお

それでは資料 2 ご覧いただきまして、1 の要旨でございます。前回皆様にお集まりいただ きました、昨年 11

Fig.5 The number of pulses of time series for 77 hours in each season in summer, spring and winter finally obtained by using the present image analysis... Fig.6 The number of pulses