数理解析研究所講究録
第 47 巻 1968 年 46-56

Model に於ける相関量のBurgers

特性汎函数による計算方法

航宇技術

\S 1.

細川

厳

統計流体力学の特性汎函数

確率変数としての非圧縮流体の速度場(一般にはベク トルであるが, Burgers$\iota\iota(x)$Model を念頭においてスガラーとして扱う。 ) の分布は, 特性汎函数によって記述されてもよい。

$x$は物理空間の座標 ( $-\infty<x<\infty$ としよう。) 。特牲汎函数は一般に時聞的に変化し,刻に於ける特性汎函数は、初期 $t=0$ に於ける$H\circ pf$

$u$

の確率分布

$P$

。

$(u)$

時

と次のように関係することが.

1) によって証明されている。

$\Phi(t$$\Omega$$t$$(1_{\sim}1)$は考えている函数 $u(x)$ の凡てから成る集合

.

従って,

$f_{\Omega^{\delta P_{0}(u)}}$は測度$P_{o}(x)$による$\Omega$

全体にわたるルベツグ積分を示し, $Tt_{u(x)}$ は. 力学法則、即ち. ナヴイエ. ストークス方程式に従って.
$u(x)$ が$t$

時刻後に発展する場を示す。 従って,$T^{t}u$は, 力学方程式の初期値問題の一般解で陽に与えられてもよい。
(1. 1)の例で理$\ovalbox{\tt\small REJECT}$の汎函数徴分 1) によって.

$t$

時刻に於ける$u$の相関量が容易に与えられることは、次される。

$\frac{\delta^{2}\phi}{i^{g}\delta_{J}(x_{1})\delta y(x_{g})}\{_{y=0}=\int_{\Omega}T^{t}u(x_{1})T^{t}u;(x_{a})\delta P_{o}(u)$$\equiv<u(x_{1})u(x_{2})>_{t}$$-46arrow$(1. 2)

この小論の目的は, 汎函数積分
は 2 点相関であるが,\S 2.

$T^{t}u$の沢山の積を作れば, 容易に多点相関は作り得る。

を計算するためには , 抽象積分の形

操作上有利である。先ず.
-(1. 2)

を直接計算するスキームを与えることである。

汎函数積分(1. 2)

ユ$(1_{-}2)$$x$$\Sigma$$R^{\infty}$

をもっと具体的な形に変えた方がを考え、これによって{と等価な函数集合箇を, 次のようにして作るを変数とする完備正規直交函数系

クリッド空間の無限次元への延長$y(x)arrow$$( \int_{\Omega}\cdot\cdot\delta P_{o})$$si^{(x)\}}$$\circ$(2. 1)$a_{k^{S}}k^{(x)}$

$k=1$

$a$

んは任意の実数。物理的に実現される速度場の成分が凡てこの集合 $-A$ と呼ぶ-の中に含まれるかどうか問題であるが, これは十分い、近似で肯定されると仮定し. そこで $\Omega=A$

元なら $\Omega=A^{3}$ ) とする。 いうまでもなく , A は$A$$R^{\infty}$

(

$u$

が3次

の構造をそのま 1 引継ぐことになる。 従って ,にわたる汎函数積分は, ユークリツド空間の中の体積々分によって具体化される可能性が出てくる。 この場合. 積分される汎函数は変数 $\{a_{k}\}$ に依存する函数と考えてよい。 所で, 空間の次元が。。の時にユークリッドの体積々分の操作が意味を持っには条件があり函数の中に各成分にっき止規化されたガウスの因子.

たとえば, 被積分$(2 \pi\sigma)-1/2_{exp}(-a_{k^{2}}/2\sigma)$が含まれるというようなことが必要である。2) (ガウス測度の汎函数積分。) この他いろいろな形の因子も考えられる。)

然し, 物理学では, 被積分函数を陽に表わさないま 1 で議論を進めることも多く.

そういう場合には積分を意味あらしめるべく必要な条件は具備されたものとして , 被積分函数を考えておくのが最も便利である。 このような立場から、測度としてはそれ自体では意味をなさない記

写 2}$\delta y=$$\pi$(2. 2)$(da_{k}/\sqrt{2\pi})$$k=1$を導入し, $\int_{A}\delta y$

は

$k=1$

に始まる

$a_{k}\iota_{-}^{}$

ついての積分の無限回繰り返し操作を意味するも$arrow 47-$.
’

のとしよう。 このような積分は, 被積分函数の中に前述した測度を構成するための因子があり.

被積分函数からそれを除いたものが有界であるなら 、確かに存在する。
(1. 2)の場合,$\delta P_{o}=p(u.

と書かれるような汎函数
た(2. 3)

t_{0})\delta u$$p$が与えられ. そして $\Omega=A$ であるならば. これが正規化され$(u\prime t_{0})$(定義上 $f_{\lrcorner}4^{p\delta u=}1$ ) 測度因子となるので. 問題なく積分は存在する。 $p$ は普通確率密度と呼ぶ所のものであるが. これは期待通り$p(u. t_{O})= \int_{A}e$ xp $\{-i\int u(x)y_{0}(x)dx\}\phi(y_{0}. t_{0})\delta y$ 。

によって, その時刻の特性汎函数

べき測度の因子をもつこと.

$\phi$$(y_{\text{。}}.

t_{0})$(2. 4)と関係づけられてもよい。 こ 1 で勿論.

い 1 換えれば $f_{A}$ } $\phi|\delta y_{0}<\infty$が前提されている。 ($p$$\Phi$

はが然る$\Phi$のFourier 逆変換 o2) これに関連して. デルタ函数と同じような性質を持っデルタ汎函数なども扱うことができる。 )(2. 3)$-$$(2.4)$

を(1. 1)に入れると , それは特性汎函数の初期値問題の一般解の形になっている。 初期に於いて、速度場の平均量 $U(x)$ と 2 点相関量 $Q$には$(x. X’)$ だけが知れている場合,$\Phi(y_{0}, t_{0})=exp\{i\int U(\alpha)y(x)dxarrow\frac{1}{2}\int\int Q(x_{*}x’)y(x\rangle y(x’)dxdx’\}$とするのが最も簡単な方法である。情報理論によると(2.5), そのような場合、上の如き Gauss の分

布が情報エントロピーを最大にする. 即ち, 最も偏見のない (unbiased) 分布であることがいえる。 以上の結果. われわれは,(1. 2)

を実際の場合にあてはめて計算できる具体的な目途を得たといえるだろう。

$arrow 48-$

\S 3. 確率密度の計算(2.4)(2.5) を仮定し 1
Burgers の乱流では. 平均流は考えないから $U(x)=0$ としたを解析的に計算することを考える。
今.

の代りに $R2N$を取り,$R^{\infty}$

2 $N$$y^{2N_{(X)arrow}}$a$\sum_{i--1}$を作り. 次に. 適当な 2$\{y^{2N}(x\rangle\}$によって, 集合$y^{2N}(x_{m})= \frac{1}{\sqrt{}\pi}$$\sum_{=,j1}^{N}$と置く。 但し.

{$\eta(k_{j})\infty Sk_{j}x_{m}$$N$

個の $x=x_{m}t_{-}^{-}$ っいて

十 $\zeta(k_{j})\sin k_{j}x_{m}$$(0, \infty)$ から抜き出した互いは区間$\triangle k_{j}$(3. 1)$i^{S}i^{(x)}$$t^{-}arrow$

}

重ならない

はその区間の平均座標とする。 2 $N$ 個の方程式 (3. 2) によって

$k_{j}$

に対する一組の$(\eta(k_{1}) , , \eta(k_{N}) \iota\sigma (k_{1}) *\cdots, \sigma (k_{N}))$のヤコビアンを調べて見よう。 この場合. 緩の変数の組を$\sqrt{\triangle N},$ろう$\zeta(k_{1})\sqrt{\triangle^{k_{1}}}$.

$\cdots,$$\sigma$$(k_{N})\sqrt{\triangle^{k_{N}}})$(3. 2)$\triangle k_{j}$

.

$j$

一組の (

$a_{1}$$\ldots.$.a$\angle N^{)}$

が定まるが , この一次変換$(k_{1})\sqrt{\triangle^{k_{1}}}$$($番目の微/J\ 区間で.

$\eta$,.

$(k_{N})$$\eta$とすると , はるかに便利 {こなることが分るであ$\circ$を含む微小区間を.

$x=x_{m}$取るとして.

(3. 1) 及び$\triangle^{x_{m}}$(3. 2)

で示し,$m$の各値に対する$\triangle x_{m}$の両辺に $s_{n}(x_{m})\triangle^{x_{m}}$ をかけては互いに重ならないように$mJ^{-}arrow$

っいて和を取ると$\backslash \text{正^{}f}$規直交函数の性質によって.
$\dot{j}$’a $n^{--} \frac{1}{\sqrt{}\pi_{j}}\sum_{-}N-1\{\eta(k_{j})\int_{-\infty}s_{n}(x)\cos\infty k_{j}xdx+\zeta(k_{j})\underline{\int^{\infty}}s_{n}(\alpha)s\dot{m}k_{j}xdx\}\triangle^{k_{j}}+\epsilon$(3. 3)となる。 こ 1 で

の極限で
これから

$\triangle$$\epsilon$

は.

$\sum_{m}\triangle x_{m}$が全区間$(-\infty.

\infty)$を隙間なく埋めて.

に収束する量である。

.
$-49-$$\triangle^{x_{m}arrow 0}$($2$

亙\rightarrow \infty )$\frac{\partial a}{\partial\{\eta(k_{j})\sqrt{}^{-}\overline{\triangle^{k_{j}}}}=\}|\frac{\triangle^{k_{j}}}{\pi}\ulcorner\underline{f^{\infty}}s_{i}(\alpha)\infty Sk_{j}xd\alpha+O(\epsilon)\equiv$(3. 4)$t_{ij}$sin$k_{j}xdx-+O(\epsilon)\equiv
ti.N+j$(3.4)’

sin(3.5)$\frac{\partial a_{i}}{\partial\{\sigma(k_{j^{)\sqrt{\triangle^{k_{j}}}}}}=\}\sqrt{\frac{\triangle^{k_{j}}}{\pi}}\int^{\infty}s_{i}(x)$

.

更に

$N$

(oos $k_{j}x\infty sk_{j}x’$ +sin$\Sigma$$k_{j^{X}}$$k_{j}x’$) $\triangle k_{j}+O(\epsilon I=\delta ip+O(\epsilon)+\epsilon’$

$J=1$

限で

$0$

に収束する量である。 (3. 4)$\vee\overline{\triangle^{k_{j}}}\}$

が

従って.

$2N$
$\delta y$

$=$

$-$

$(3.5)$

は. 変換 $\{a_{i}\}arrow\{\eta(k_{j})\sqrt{}\triangle^{k_{j}}$.

$C(k_{j})$, 丁度直交変換に近いものになっていることを示している。
(2. 2)に相当して,

2 $N$$(da_{k}\prime 2\pi)$

$\pi$

$k=1$

$N$
$–\sim$$\pi\{d\eta(k_{i})\sqrt{\triangle^{k_{i}/2\pi}}.

d\sigma(k_{i})\sqrt{\triangle k_{i}/2\pi}\}$(3. 6)$i$と書くことができ. この近似関係は$\delta y^{r}arrow\delta y_{\alpha}\Delta N$

そして殆んど凡ての$y(x)= \frac{1}{\sqrt{}\pi}\int_{0}^{\infty}${$x$$\eta(k)\infty Skx+\sigma(x)$$\epsilon$.

$\epsilon’arrow 0$の時に厳密になるであろう。そのような極限では,に対して,Sin $kx$ } $dk$が成り立っ。 そこで. 便宜上.

$-50arrow$(3. 7)

(3. 8)$\delta y=\infty\pi\{d\eta(k\sqrt{dk/2}\Gamma_{4}\cdot d_{\zeta}(k)\sqrt{dk/2\pi}\}$という書き方を許すことにしよう。
$\eta(k)=\frac{z(k)+z(-k)}{\sqrt{2}}$,(3. 7)

1 は.

(3. 9)$\sigma(k)=\frac{i\{z(k)-z(-k)\}}{\sqrt{}\overline{2}}$を導入することによって.
$($$y(x)= \frac{1}{\sqrt{2\pi}}\int_{-\infty}z(k)eikx_{dk}$と書かれてもよい。 但し ,さて,

(2. 5) に戻り

1($z^{*}(k)=z(-k)$ である。

$Q$

$*$

は共役複素数。 )が空間的に一様で,(3. 1 1)$Q(x-x^{l})= \frac{1}{2\pi}\int_{-\infty}\sigma(k)eik(x-x’)_{dk}$とフーリエ分解される時は.

3.1 $0)$$\Phi$の $exp$の肩は次のように簡単な形に帰着する。

$\underline{1}f^{\infty}\int^{\infty}Q(x-x^{l})y_{0}(x)y_{0}(x’)dxdx’=\int\underline{1}$2$2-\infty-\infty$$\sigma(x)z_{0}(k)z_{0}(-k)dk$$-\infty$(3.1 2)$= \frac{1}{2}\int_{o}Re[\sigma(k)]\{\eta_{\theta}^{2}(k)+\sigma_{0^{(k)}}^{2}\}dk$( $Re$ は実数部を示す。 ) 次に$u(x)= \frac{1}{\sqrt{\pi}}\int_{o}^{\infty}\sqrt{Re[\delta(k0}$という, や$\cross${$\Lambda(k)\infty Skx+M(k)$sin $kx$ } $dk$

技巧的なフーリエ分解を行なうと ,

$($

$($

3. 1 2$)_{-}$

を(2. 5)に入れ. これと$3)_{\parallel,:}$$-$(3. 1 4)$\int_{-\infty}y_{0}(x)u(x)dx=\int_{0}^{\infty}\infty\sqrt{\underline{|Re_{1}}[\sigma(k)]}\{\Lambda(k)\eta(k)+M(k)\sigma(k)\}dk$が得られる。

3.1(3. 1 4)$-51-$を (2. 4)に入れて,

(3. 8)のa$p(u. t_{0})=exp$(3. 1 5)$[- \frac{1}{2}\int_{o^{\infty}}\{\Lambda(k\Psi+M(k)^{2}\}dk]\pi\infty(1/\sqrt{Re\sigma(k))}$となる。

以上でこの節の目的は終ったが. 一方,
ると(3. 1 3) で示される $u(x)t_{-}^{-}$ っいて,

(3. 8)

を用

$Aa$

,

$\delta u=\pi\{\sqrt{}\overline{Re[\sigma(k)}]d\Lambda(k)\sqrt{ak/2\pi}$であるから.-(3. 1 5).

$\backslash R\overline{e}[\sigma(k)]dM(k)\sqrt{d:\iota/2\pi}\}$(3.1 6)と合わせて,$p(u. t_{0}) \delta u=exp[-\frac{1}{2}\int_{0}^{\infty}\{2dk]\pi\{\infty d\Lambda(k)\sqrt{dk/2\pi d}M(k)\sqrt{dk/2\pi}\}$(3. 1 7)を得る。 これはまさに変数 $\{\Lambda(k)\sqrt{dk}. M(k)\sqrt{dk}\}$ が標準正規分布に従うことを意味する。

\S 4.

モンテカルロ法の適用(1. 2)の$\delta P_{o}$mmndel の場合は

(u) は

.

$x$

(3. 1 7)を ( $-\infty$$(_{arrow}^{}$

よって具体化されたので. 次に $Ttu$ を考察する。

の中に取った時. Hopf$\epsilon_{\infty)}$

次元型で

.

urger $s’$のと Cole 5) にょって.

$Tt_{u}= \frac{\int_{-\infty}(x-x’)exp\{-\frac{R}{2}f_{o}u(x^{l})dx^{l}-\frac{R(x-x’)}{4^{C}}\}dx’\infty x’a}{t\int_{-}\infty exp\{-\frac{R}{2}f_{0},\infty u(x’)dx^{l}-\frac{R(x-x’)^{2}}{4t}\}dx}$と陽の表現が得られている。 念の為$B$$u(x)$ の時間的発展を律する(4. 1)

Burgers の基礎方程式は無4(4. 2)$\frac{\partial u}{\partial t}---u\frac{\partial u}{-\partial)}+\frac{1\partial^{z}u}{R\partial x^{2}}$と書かれ, $R$ は乱流の初期スケ$-K\downarrow/M$. 初期速度の自乗平均 \mbox{\boldmath$\rho$}( を使ったレイノルズ数である$0$ }$\sqrt{}$ことを注意する。

$-52-$

従って.

の単位は

$u$

$(3. 1^{-}3)\dot{\text{を}}$$\sqrt Q_{\sim}(O)\prime x$$(4.1)$に代の単位は $M$$\nearrow\backslash \text{す^{}-}$

.

ると

.

時間

$\dot{T}^{t_{u}^{\backslash }}$は$t$$\Lambda$の単位は $M/\sqrt Q(0)$ である。
及び$\overline{M}$の-与菟られた汎函数と$\text{な^{}\backslash }$$\dot{\ovalbox{\tt\small REJECT}}^{a}t^{\vee}$

る。 従 っ.
(1. 2)て,は,(3. 1 7)

この場合,ているわけであ F る。

$T^{t}u$が肴界なら積分が春在するのは明らふであるが. 美際の許算ほ.

の複雑さのため. -容坐でない。. そ1)から見られる通り一種のガウス測度の汎函数積分になつ$\iota^{-}$

興味がある。 (そこでは,(4.

.

$\sigma J$

$)$

$c$

で

$\cdot$.‘-.. される ‘k。 特に. なんら、かの近似が要.

$\ovalbox{\tt\small REJECT}$

$R$

(4.

の大きい所

$t_{-}^{-}$から見られるように、非線型効果が相対的に優越する故。)$Y$

高度あ多重積分が本質釣な働き壱しているような現在の問題には. モンテカルロ法が有望のよう$\text{れ^{}-}\text{る^{}-}\sim$$\iota_{-J^{\backslash }}^{-ffl_{\backslash }^{-}}$

わ$A_{\urcorner}$。

$\sqrt dk\}$の$\ovalbox{\tt\small REJECT} 1\Rightarrow$合$tarrow=_{\vee}$

は. モンテカルロ ‘& の:$\dot{\Re}^{-}\dot{4}’F:$は$\dot{g}$た$\dot{\ovalbox{\tt\small REJECT}}\hat{\ddot{\text{単}}}\dot{\text{て}}\cdot$

– —で

あっ

$\prime M(\grave{k})\backslash$\text{て^{}\mathcal{X}}\backslash -\{\Lambda(k)\sqrt{dk},$$X^{-\cdot.\cdot=}:1$の値として. 夫々独立な標準主規乱薮を与 え$\dot{\overline{T}}u\dot{(}x)^{-}T$u-$(X’)^{:}\dot{\text{を}}$:_{-}d$

計 G- し. これを$–$.

多数繰り返えして. その平均値を取れば O. 2) を近似的に計算したことにたる. めである。 実際には.

$\triangle k_{j}$有映で計算する i-1 ら,- このための誤差$(\epsilon i\epsilon’\neq 0)$;にっいては. 適当なチエツクをする必要があろう。

も入ってくる。 従$–$ c-$-$$:^{-}-\backslash \cdot$て. 精度$\dot{;}$-c,$\not\in*$$-\backslash$$j$$\backslash$

今の場合、誤差にっいて概していえることを列記しておこう。
.

(1)

モンテカルロ法の誤差
イ)-

試行回数によるもゐ。 これは$arrow ji^{)_{J}}\cdot\neq\vee\backslash 0$

sampling などいろいろ’ な改善乗もある。
$\mathfrak{o}^{\backslash })a\dot{e}\{I^{\backslash }\mathcal{A}’$Si 数$t_{-}^{\vee}.k\dot{\text{る}}$

もの o

.

$–$

,

る数

..

有 限の課差$-\cdots-$$\nu\backslash -$$\pi\Rightarrow 4ts\iota\ddot{\Re}\prime X\xi E^{:\Leftrightarrow}\dot{\text{し^{}\prime}}\grave{\text{な}}^{-}’-:_{u^{a}\text{の}\dot{\text{て}}\cdot\dot{a}_{\hat{\dagger T}^{t}}\fbox{ }\text{数}\dot{\text{を}}\mathfrak{B}}’--$

は-\tilde 或る範囲で振勅し. そ: れ以上の精 ‘R- は得(2))$\vee^{\backslash }$$\backslash .\backslash _{\wedge^{\backslash }}\sim$

$-$

.

$-$

$t$

$\dot{\text{ら}}*\iota k\dot{b}^{a^{\wedge}}\circ---\cdot$$:.\cdot$

.:

茸

‘

$\mathfrak{t}$$-$,$)_{-}^{\backslash }/1_{1}$上$\ddot{\text{上}}(\neq$て$b\#^{\dot{k}}\backslash -\ovalbox{\tt\small REJECT}/$$j_{\urcorner}^{e}\nu-’$$j|:r_{\grave{1}}-\sim,,$

.

$–$

. .
$t$

$\sim$$\triangle\dot{k}_{j^{e}}$(3. 1 3)の中に$\backslash -($を考慮して$\sigma$$7/$(k) があるが, 通常これは.

. $–$$\bigvee_{-}^{\backslash }$

.

$k$

:
に $n$ し$\dot{\text{て}}$強い減蓑性をもっている切で. これ$-$$.-\vee^{--:}$.

$u(x) \cong\frac{1}{\sqrt{}\pi}J^{N}\sum_{-,-1}\sqrt{Re[\sigma(k_{j})}1f\Lambda(k_{j})^{:}\infty sk_{j}x+\check{M}(k_{j})$$-53-$$\backslash$

.

sm

$.I$

.

.
$\cdot$}$\cdot\triangle^{k}$$k_{j}x\cdot$;:$(\overline{4}^{(}3)$

有限による台形法則の誤蓬よりも小さく守ること瀕できる$\sqrt 7k$

l.M 域 )
$k$$\bullet$

こうすると

9

$T^{t}l4$には $\{/t(k)$} 瀞含まれないことになり. これらの変数にっいての積勇は 1 @,17)$\sqrt ih|k\gg k_{\underline{f}}\psi_{-}^{-}$

によって凡て獄立 (叢霧行できて、常$t_{\approx}^{-}\}_{-}l’arrow$, {@. 17) は案質的儒.

なり$-p \delta u\frac{\approx}{=,-}exp[’\approx\frac{-1}{\overline{g}}\sum_{-}\{\Lambda(k)^{g}+Mtk)^{\S}3\triangle^{k_{j}]}1^{N}-\zeta-.J^{!}\frac{FX}{\frac{arrow---}{-}}\iota\{d_{\Lambda}(k_{j})\frac{-k\overline{/\prime g\pi}}{\overline{\triangle}j}dM(k_{\dot{i}})$$\ovalbox{\tt\small REJECT}_{i}\Lambda;i$と毒えてもよいこと傷なる g

り 1 それは倫)$T^{l}\ddagger\iota$$\Delta k_{j}$$\xi S1\approx$

は有限$\iota\triangle^{k_{j}}$$t_{\approx}^{\vee}$

よる- $g,$$g\nuarrow\neq 0$:{$4_{*}4$の議薫だけが含まれ$\mathfrak{B}$)おの定差 乙よる壷形法則の誤差 ⇔氾戮任△襦

の計算誤瞭

醗究会の講演では, 適当に$Q$を与えた計算の案例を示し. 結果礁ついての討論も行っ- たが 1

$1_{=}^{\geq}$

こでは省略し. 別の機会に詳しく述べることにする。

85. 補

足

以上に述べた方法は.-T$u$の形瀞陽に与えられる時ほは, いつでも硬えるは誇懲汎函数の方法は駄厨かという &, 必ずしも髪うではない。 案瞭初期値$\circ$$\mathfrak{x}-\iota$

そうでない鶏禽を与えて ,$\mathscr{Z}^{\underline{t}}u$を近似的に計算する一般駒畢段泰あれば, 翻的娃達する如雫ある @ たとえ嫁, 初期値間題を解くの [; $u$ の Z 礎方 G 式 ‘; 於いて, $u(x_{*}\ell t--$

$N$

.

$g$

$\iota_{i^{\langle\ell)g_{\dot{i}}(g\}}}$と屑開$L*\xi b_{i}(\ell)3\iota_{-\mathcal{Z}}^{\Phi}$い$\mathfrak{T}$の適$J=:-$

立常微分方程弐を作 っ て, これを$-$$\ddagger i\iota s$

ng@\approx Kuft\S 法で解くというやり方が考えられ嬉。

$t\approx$

種の Galerkin 法といえるかも知れない。 )葡節に関連港せて. 鍵い易い形でいうならば,(講演で述ぺ潅ものと少し異なる潮, )2{$\lambda j=1N^{t}-(k_{\dot{J}^{1}}t)\mathfrak{W}k_{j^{X}}+lt$$u(x_{\iota}t)= \frac{1_{-}}{\sqrt{}\overline{X^{-}}}g$I$k_{j}$

.

$t)$

SZ

と置く。 $t=t$ @ では , これは (4. 3) と一致する。倶し,$\Leftrightarrow g_{-}4\approx-$}$k_{j}\alpha^{-}$(6. 殴$\triangle\tilde{k}_{j}$$N^{l}\backslash >N$

である Q- 従って$\iota$$\lambda(k_{j^{i}}t_{9})-=--\{1_{0}^{\frac{\ddagger\S\vee e[q(k}{j)}}A(k_{j})$.

$j^{j}\sim>N\leqq N$(5, 2)$\cdot$$/-|(k_{\downarrow ’ i}S_{Q})\equiv\{\begin{array}{l}\sqrt{arrow\overline{R}\overline{e_{-}[\sigma(}k_{j})\overline{]}\backslash }M(.k_{j})l/\simeq\prime V\subset O\prime j\geq_{-}N(5\cdot.3\}\end{array}$$y$の基礎方程弍

られる$(S$ ,

をは常$k_{j}$$8_{\wedge}$@$)$$\iota^{-}-arrow-(5,1)$

と A

を入れれ

$\#\ovalbox{\tt\small REJECT}$}

會適当 } 択ぶと\tilde -\tilde \dager$\vee\equiv$$k_{j}$(. 5. @) を籾期条陣とし$\leq$$\lambda,$$5_{:}l$$t^{\underline{\iota}}$

}

を烹配する

$t$

前鯖の奥

$\not\simeq=$

ムは$\not\in$のま$s$

建 $l!t!$ \S $e=$ 瀬 1\ddagger

は, $N’$$t$$\dot{\llcorner}_{-}>$の場合$\iota$瀬屠する笛である 9 始め.

を変えた場合の解の変り方によ.

一$\mathscr{C}-$$!^{\vee}arrow=k\S_{-b}*g$を利用す a9 二れ@ 慈で解いたもの舎逼

とc\acute し\tilde て与えら-,もう一っの誤差の原因としてが余り大きくない限り

a形$\zeta f$$u(x,$ りの申焦高調披威翁が発黛する可能性を窓れで,あることは.

$\xi$)$t\not\leq;^{-}7=$$t_{w}^{\vee}\{Alk_{j}1\prime_{\overline{4}}r_{i}\downarrow M!k_{i}\}\sim\overline{\overline{g}}l_{\overline{j}}3\emptyset_{=}I$

もはや! 非線型の) 常磁倉方程戎’\downarrow 簡巣 ( 得$\iota\veearrow:.$a薫,$T^{\ell}l_{\vee}\iota$$\llcorner\geq--*\iota$$,$-N 有醒であ$N^{j}\geq N_{\sim}\text{と_{}-}$. を\brev ^\acute$\zeta$$9\leq\not\in_{-\}f}$

し\tilde 4 莚のであ登$*_{-}\S$と\sim A$u(x_{\iota}t$$\text{見_{}d\backslash }.\hslash_{J}\mathscr{E}t$

時間き装
$N$

$I\not\cong$が有鰻誓

列認されなければならない 9 精廉の-‘. 判定て-.. 得られる 9組の軍規乱数に対して ‘ Runge $-Kutta$ 法を実行 $*$ る計算時間(4. 1) に与えら- れた -T $tu$ の表現 41 値計算する時間と.. 坊べ得る程度なので , 将米性はあは,ると見てよいだるう。

大

$-!)$

2)
3)

}

$B_{\ell}H_{c}\otimes pf_{:}$$I_{-:}$

J. Rat, Meth. A

Hss\otimes kawa ;$J$$t1_{-}\S 1_{=^{1(19_{-}5_{-}2)@7}}$, Math, Phys, $\underline{8_{-}}(1967)2_{-}21$

K,e. $Fr\dot{:}sdric$ hs et al: I.n. $teg- rati\sigma$)

$n$

献

,

:.

.

of Funetionals-, leeture notes, Inst. Magh. Sci., New則\downarrow v. (1957),

4)

E

5)

J. D. Cole-$H_{-op^{\xi;c_{O\Phi t\mathfrak{n}w\iota!}}.\cdot p_{u_{--}re}}$$jQu_{-a}r\mathfrak{t},$Appl. Math.

$\underline{3_{-}}(1950)2_{-}01$.

AI, $M\S\xi\}_{=}919_{\vee-}=.225-$ ,$\approx\S\ovalbox{\tt\small REJECT}\infty$$Y_{_{-}f}k_{\wedge}$

後記

本論でや
は.

$\backslash$heuristic$(\tau\backslash \wedge’or\Pi)al$というべきか$?$(二扱った汎函数積分の操作について)

無学者による明快な解説がない限り , 納得しない方も多いかも知れない。 モスクワの数理物

理学者 Yagl$\cap \mathfrak{m}$

教授 (大気物理研究所

.

1966)

も, 無限次元ベク トルの分布を与える確率

密度をいきなり出して来る Edwards (J. Fluid

N,lech. l8(1964), $\angle 39$ ) の乱流理論 (二疑問を抱いている一人であった。 これについての, 私の意見は次の Peierls 教授の発言$(_{-}^{-}$

代弁される$\circ$some matheHlaticians as

For many years Dirac$inConSiS\mathfrak{c}ent$$s$

use of the delta function was frowned upon by. So it was, of course, but no more so than Newton’$s$use ofthe coefficient $dy/dx$ with lnfinitely small numerator and denominator. It would have been a correct but unconstructive approach to $teU$ Newton his equations wereconstructive to develop concepts that made them meaningful.
(M. 1. T., 1964).

$arrow 56-$“

meaningless. lt was more

-Analysis in Function Space