MathLabs

解析学

特異積分とカルデロン=ジグムント理論

特異な核を持つ積分作用素を研究する分野で、現代の調和解析や偏微分方程式の中心となる。

直観通常の積分にはなり得ないほど特異な核

写真編集の『シャープ化』フィルタ、音声の位相シフタ、そして弾性材料のひび割れは、いずれも同じ数学的対象を隠し持っている:原点で発散する核、たとえば1次元での 1/x1/x に対する畳み込みである。そのような核は0の近くで絶対可積分ではないため、 ∫K(x−y)f(y) dy\int K(x-y)f(y)\,dy は通常の意味では収束しない——しかし特異点の両側にある核の正の部分と負の部分はほぼ完全に打ち消し合い、適当な ff に対して主値の極限が存在する。特異積分作用素とはまさにこのように構成される作用素であり、カルデロン–ジグムント理論はそれらを意味づけ、期待どおりに良い性質を持つことを証明するための道具立てである。

上半平面上の正則写像を可視化する、ヒルベルト変換に関連したインタラクティブな複素領域ウィジェット。
ヒルベルト変換は上半平面上の正則関数の境界値として現れる。境界付近に拡大したときにこの複素写像がどのように振る舞うか探ってみよ。

大学原型:ヒルベルト変換

定義: ヒルベルト変換

R\mathbb{R} 上の良い関数 ff に対して、そのヒルベルト変換は主値畳み込み Hf(x)=p.v. 1π∫−∞∞f(y)x−y dyHf(x) = \text{p.v.}\ \dfrac{1}{\pi}\int_{-\infty}^{\infty} \dfrac{f(y)}{x-y}\,dy である: y=xy=x における特異点は極限をとる前に対称に取り除かれ、これによってまさに y<xy<x と y>xy>x の間の打ち消し合いが有限な値を生み出すことができる。

Hf(x)=p.v. 1π∫−∞∞f(y)x−y dyHf(x) = \text{p.v.}\ \dfrac{1}{\pi}\int_{-\infty}^{\infty} \dfrac{f(y)}{x-y}\,dy

フーリエ側では、ヒルベルト変換は驚くほど単純である:各周波数に純粋な位相 Hf^(ξ)=−i sgn⁡(ξ) f^(ξ)\widehat{Hf}(\xi) = -i\,\operatorname{sgn}(\xi)\, \hat f(\xi) を掛け、 ξ=0\xi=0 で符号が反転する。すべての周波数は 90∘90^\circ 回転させられる(回転の符号は周波数が正か負かによる)が、大きさは変わらない——これはすでに HH が L2L^2 上の等長写像であるはずだという示唆になっている。

Hf^(ξ)=−i sgn⁡(ξ) f^(ξ)\widehat{Hf}(\xi) = -i\,\operatorname{sgn}(\xi)\, \hat f(\xi)

定義: カルデロン–ジグムント作用素

より一般に、カルデロン–ジグムント作用素とは、 Rn\mathbb{R}^n 上の畳み込み Tf(x)=p.v.∫RnK(x−y)f(y) dyTf(x) = \text{p.v.}\int_{\mathbb{R}^n} K(x-y) f(y)\,dy であって、その核 KK が原点から離れたところでサイズ評価と滑らかさ(ヘルマンダー)評価 ∣K(x)∣≤C∣x∣n,∣∇K(x)∣≤C∣x∣n+1|K(x)| \le \dfrac{C}{|x|^n}, \qquad |\nabla K(x)| \le \dfrac{C}{|x|^{n+1}} を満たし、さらに TT がすでに L2L^2 上で有界であることが分かっているものである。ヒルベルト変換は n=1n=1 、 K(x)=1/(πx)K(x)=1/(\pi x) の場合であり、リース変換 Rjf=p.v. cn∫xj−yj∣x−y∣n+1f(y) dyR_j f = \text{p.v.}\, c_n \int \frac{x_j-y_j}{|x-y|^{n+1}} f(y)\,dy は、以下で楕円型偏微分方程式に用いる高次元版である。

1次元から多次元へ
作用素次元核 K(x)K(x)典型的な用途
ヒルベルト変換n=1n=11/(πx)1/(\pi x)解析信号、位相シフト
リース変換n≥2n\ge 2(xj−yj)/∣x−y∣n+1(x_j-y_j)/|x-y|^{n+1}楕円型偏微分方程式の2階微分評価
一般のカルデロン–ジグムント TTn≥1n\ge 1サイズ+ヘルマンダーの滑らかさ上記すべてに共通する枠組み

研究主要な定理

任意の f∈L2(R)f \in L^2(\mathbb{R}) について ∥Hf∥L2=∥f∥L2\|Hf\|_{L^2} = \|f\|_{L^2} が成り立つ:ヒルベルト変換は単に有界であるだけでなく、 L2L^2 上の等長写像である。

なぜ正しいのか?

これはあらゆる一般的なカルデロン–ジグムントの定理が拠り所とする基本ケースである: L1L^1 や LpL^p 上で特異積分作用素を制御する前に、まず L2L^2 上でよく振る舞うことが必要であり、そこではプランシュレルの定理によりフーリエ変換を通じた直接的なアクセスが得られる。

証明

ステップ1(正則化を通じて乗数を求める)。核 1/(πx)1/(\pi x) は可積分ではないので近似する: ϵ>0\epsilon>0 に対して kϵ(x)=1πxx2+ϵ2k_\epsilon(x) = \frac{1}{\pi}\frac{x}{x^2+\epsilon^2} と置く。これは奇関数で可積分な近似であり、 ϵ→0\epsilon\to 0 のとき主値核へ収束する。直接計算(留数積分または公式集)により k^ϵ(ξ)=−i sgn⁡(ξ) e−2πϵ∣ξ∣\hat k_\epsilon(\xi) = -i\,\operatorname{sgn}(\xi)\, e^{-2\pi\epsilon|\xi|} が得られる。

ステップ2(極限へ移行する)。 ϵ→0+\epsilon \to 0^+ のとき、固定した各 ξ≠0\xi \ne 0 について e−2πϵ∣ξ∣→1e^{-2\pi\epsilon|\xi|} \to 1 であり、良い ff について kϵ∗f→Hfk_\epsilon * f \to Hf となる。畳み込み定理の中で極限をとると、乗数の等式 Hf^(ξ)=−i sgn⁡(ξ) f^(ξ)\widehat{Hf}(\xi) = -i\,\operatorname{sgn}(\xi)\, \hat f(\xi) が得られる:ヒルベルト変換は周波数側で −i sgn⁡(ξ)-i\,\operatorname{sgn}(\xi) を掛けるという形で作用する。

ステップ3(プランシュレルの定理を適用する)。 ∣−i sgn⁡(ξ)∣=1|-i\,\operatorname{sgn}(\xi)| = 1 が(測度ゼロの一点を除く)すべての ξ≠0\xi \ne 0 で成り立つので、ほとんどすべての ξ\xi について ∣Hf^(ξ)∣2=∣f^(ξ)∣2|\widehat{Hf}(\xi)|^2 = |\hat f(\xi)|^2 となる。両辺を積分してプランシュレルの定理を用いると ∥Hf∥L2=∥f∥L2\|Hf\|_{L^2} = \|f\|_{L^2} が得られる:ヒルベルト変換は L2(R)L^2(\mathbb{R}) 上の等長写像であり、特に有界である。

畳み込み作用素 TT が L2(Rn)L^2(\mathbb{R}^n) 上で有界であり、その核がカルデロン–ジグムント条件 ∣K(x)∣≤C∣x∣n,∣∇K(x)∣≤C∣x∣n+1|K(x)| \le \dfrac{C}{|x|^n}, \qquad |\nabla K(x)| \le \dfrac{C}{|x|^{n+1}} を満たすならば、 TT は弱(1,1)型 ∣{x:∣Tf(x)∣>λ}∣≤Cλ∥f∥L1\big|\{x : |Tf(x)| > \lambda\}\big| \le \dfrac{C}{\lambda}\|f\|_{L^1} であり、したがって(マルツィンキェヴィッチの補間定理と双対性により) TT はすべての 1<p<∞1<p<\infty について Lp(Rn)L^p(\mathbb{R}^n) 上で有界である。

なぜ正しいのか?

この一つの定理が、特異積分——その核が絶対可積分でないにもかかわらず——がほとんどすべての LpL^p 空間上で通常の有界作用素のように振る舞う理由を説明する:唯一の犠牲は端点 p=1p=1 であり、そこでは有界性が強い型から単なる弱い型へと弱まる。

証明

ステップ1(カルデロン–ジグムント分解)。 f∈L1f \in L^1 と λ>0\lambda>0 を固定する。二進立方体上の停止時刻論法(同値に、ハーディ–リトルウッド極大関数 MfMf)を用いて、 f=g+bf = g + b と分解する。ここで『良い』部分は ∥g∥∞≤Cλ\|g\|_\infty \le C\lambda と ∥g∥1≤∥f∥1\|g\|_1 \le \|f\|_1 を満たし、『悪い』部分 b=∑jbjb = \sum_j b_j は互いに素な二進立方体 QjQ_j 上に台を持ち平均ゼロ ∫Qjbj=0\int_{Q_j} b_j = 0 であるピースの和で、総測度は ∑j∣Qj∣≤Cλ∥f∥1\sum_j |Q_j| \le \frac{C}{\lambda}\|f\|_1 を満たす(これはまさに停止立方体に適用した MM の弱(1,1)評価である)。

ステップ2(良い部分)。 TT は L2L^2 上で有界であり、 ∥g∥22≤∥g∥∞∥g∥1≤Cλ∥f∥1\|g\|_2^2 \le \|g\|_\infty \|g\|_1 \le C\lambda\|f\|_1 なので、チェビシェフの不等式により ∣{∣Tg∣>λ/2}∣≤4λ2∥Tg∥22≤Cλ∥f∥1|\{|Tg| > \lambda/2\}| \le \frac{4}{\lambda^2}\|Tg\|_2^2 \le \frac{C}{\lambda}\|f\|_1 が得られる:良い部分だけですでに弱(1,1)評価を満たす。

ステップ3(悪い部分、2倍立方体の外側)。 Qj∗Q_j^* を QjQ_j と中心が同じで辺の長さが2倍の立方体とする;和集合 ⋃jQj∗\bigcup_j Q_j^* の測度もまた Cλ∥f∥1\frac{C}{\lambda}\|f\|_1 以下である。 Qj∗Q_j^* の外側では、 bjb_j の平均ゼロ条件により核から定数を引くことができる: QjQ_j の中心を cjc_j として Tbj(x)=∫Qj(K(x−y)−K(x−cj))bj(y) dyTb_j(x) = \int_{Q_j} \big(K(x-y)-K(x-c_j)\big) b_j(y)\,dy となり、ヘルマンダーの滑らかさ評価 ∣K(x)∣≤C∣x∣n,∣∇K(x)∣≤C∣x∣n+1|K(x)| \le \dfrac{C}{|x|^n}, \qquad |\nabla K(x)| \le \dfrac{C}{|x|^{n+1}} により、この差は ∫(Qj∗)c∣Tbj(x)∣ dx≤C∥bj∥1\int_{(Q_j^*)^c} |Tb_j(x)|\,dx \le C\|b_j\|_1 となるほど十分小さいことが示される。

ステップ4(まとめる)。ステップ3を jj について足し合わせ、2倍立方体の外側でチェビシェフを適用すると ∣{x∉⋃jQj∗:∣Tb(x)∣>λ/2}∣≤Cλ∑j∥bj∥1≤Cλ∥f∥1\big|\{x \notin \bigcup_j Q_j^* : |Tb(x)|>\lambda/2\}\big| \le \frac{C}{\lambda}\sum_j \|b_j\|_1 \le \frac{C}{\lambda}\|f\|_1 となる。ステップ1の ⋃jQj∗\bigcup_j Q_j^* の測度を加えると、悪い部分もレベル集合に O(∥f∥1/λ)O(\|f\|_1/\lambda) だけ寄与し、ステップ2の良い部分と合わせると完全な評価 ∣{x:∣Tf(x)∣>λ}∣≤Cλ∥f∥L1\big|\{x : |Tf(x)| > \lambda\}\big| \le \dfrac{C}{\lambda}\|f\|_{L^1} が得られる。

大学実世界での応用と具体例

カルデロン–ジグムント理論は、一見まったく異なる2つの対象を同じ道具で制御することで、解析学における中心的な地位を得ている:楕円型偏微分方程式の解の2階微分と、技術者が実世界の波形から瞬時振幅と位相を取り出すために必要とする『直交』信号である。

例: ニュートンポテンシャルの2階微分

楕円型偏微分方程式の理論では、 Rn\mathbb{R}^n 上で Δu=f\Delta u = f を解き、 uu の2階微分を ff によって制御したい。ある 1<p<∞1<p<\infty について f∈Lp(Rn)f \in L^p(\mathbb{R}^n) であるとき、なぜ ∥D2u∥Lp≤Cp∥f∥Lp\|D^2 u\|_{L^p} \le C_p \|f\|_{L^p} が成り立つのか。

解答

解由与ニュートンポテンシャルとの畳み込みで与えられる: u=f∗Nu = f * N 、ここで n≥3n\ge 3 のとき N(x)∼cn∣x∣2−nN(x) \sim c_n|x|^{2-n} である。積分記号の下で形式的に2回微分すると ∂i∂ju=f∗(∂i∂jN)\partial_i\partial_j u = f * (\partial_i\partial_j N) となるが、 ∂i∂jN(x)\partial_i\partial_j N(x) は次数 −n-n の斉次関数であり、原点周りの球面上で平均ゼロを持つ——これはまさにカルデロン–ジグムント核であり、 i=ji=j のときにのみ意味を持つディラック測度の項を除けばそうである。

したがって、その ff 自身の定数倍の項を除けば、 ∂i∂ju=Rijf\partial_i \partial_j u = R_{ij}f が成り立ち、ここで RijR_{ij} は(リース変換型の)カルデロン–ジグムント作用素であり、その核はサイズおよびヘルマンダーの滑らかさ条件を満たす。

カルデロン–ジグムントの定理により、 RijR_{ij} はすべての 1<p<∞1<p<\infty について LpL^p 上で有界であるから、各 i,ji,j について ∥∂i∂ju∥Lp≤C∥f∥Lp\|\partial_i\partial_j u\|_{L^p} \le C\|f\|_{L^p} が成り立ち、これはまさに ∥D2u∥Lp≤Cp∥f∥Lp\|D^2 u\|_{L^p} \le C_p \|f\|_{L^p} である。これは楕円型正則性理論の基礎となる古典的なカルデロン–ジグムント評価であり、まさに p=1p=1 と p=∞p=\infty で破綻するため、偏微分方程式論の研究者は LpL^p( 1<p<∞1<p<\infty )あるいはヘルダー空間で作業せざるを得ない。

例: 無線技術者のための解析信号の構成

通信技術者が実数の搬送波信号 f(t)=cos⁡(2πνt)f(t) = \cos(2\pi\nu t) ( ν>0\nu>0 )を持っており、瞬時振幅と位相を取り出す(単側波帯無線や包絡線検波のような)ために使われる解析信号 f(t)+iHf(t)f(t) + iHf(t) を構成するための90°位相シフトされた『直交』信号を求めたい。 HfHf を計算せよ。

解答

cos⁡(2πνt)=12e2πiνt+12e−2πiνt\cos(2\pi\nu t) = \tfrac12 e^{2\pi i \nu t} + \tfrac12 e^{-2\pi i \nu t} と書くと、そのフーリエ変換は ξ=±ν\xi = \pm\nu における一対のスパイク 12δν+12δ−ν\tfrac12\delta_\nu + \tfrac12\delta_{-\nu} である。

ヒルベルト変換の乗数 Hf^(ξ)=−i sgn⁡(ξ) f^(ξ)\widehat{Hf}(\xi) = -i\,\operatorname{sgn}(\xi)\, \hat f(\xi) を適用する: ξ=ν>0\xi=\nu>0 では −i-i を掛け、 ξ=−ν<0\xi=-\nu<0 では +i+i を掛ける。よって周波数側では、変換された信号は −i2δν+i2δ−ν-\tfrac{i}{2}\delta_\nu + \tfrac{i}{2}\delta_{-\nu} となる。

逆変換すると、これは −i2e2πiνt+i2e−2πiνt=i2(e−2πiνt−e2πiνt)=sin⁡(2πνt)-\tfrac{i}{2}e^{2\pi i\nu t} + \tfrac{i}{2}e^{-2\pi i\nu t} = \tfrac{i}{2}\big(e^{-2\pi i\nu t}-e^{2\pi i\nu t}\big) = \sin(2\pi\nu t) であり( e−iθ−eiθ=−2isin⁡θe^{-i\theta}-e^{i\theta}=-2i\sin\theta を用いた)、 H(cos⁡(2πνt))=sin⁡(2πνt)H(\cos(2\pi\nu t)) = \sin(2\pi\nu t) が得られる。

したがって解析信号は cos⁡(2πνt)+isin⁡(2πνt)=e2πiνt\cos(2\pi\nu t) + i\sin(2\pi\nu t) = e^{2\pi i \nu t} ——単一の正の周波数を持つ回転指数関数であり、その絶対値(ここでは常に1)と位相が技術者の求める瞬時振幅と位相を直接与え、冗長な負の周波数成分は完全に捨てられる。

Hf^(ξ)=−i sgn⁡(ξ) f^(ξ)\widehat{Hf}(\xi) = -i\,\operatorname{sgn}(\xi)\, \hat f(\xi) で定義されるヒルベルト変換のフーリエ乗数は:

カルデロン–ジグムントの定理によれば、適切な核を持ち L2L^2 上で有界な特異積分作用素は、端点 p=1p=1 において:

単側波帯(SSB)無線や包絡線検波において、技術者が『解析信号』 f(t)+iHf(t)f(t) + iHf(t) を構成するのは:

f∈Lp(Rn)f\in L^p(\mathbb{R}^n)( 1<p<∞1<p<\infty )に対するポアソン方程式 Δu=f\Delta u = f について、カルデロン–ジグムント評価が制御するのは:

参考文献

  1. Tuomas P. Hytönen (2012). The sharp weighted bound for general Calderón–Zygmund operators · arXiv:1007.4330
  2. Elias M. Stein (1970). Singular Integrals and Differentiability Properties of Functions