=================================================================== RCS file: /home/cvs/OpenXM/src/asir-contrib/packages/doc/n_wishartd/n_wishartd-ja.texi,v retrieving revision 1.2 retrieving revision 1.3 diff -u -p -r1.2 -r1.3 --- OpenXM/src/asir-contrib/packages/doc/n_wishartd/n_wishartd-ja.texi 2016/08/25 03:13:54 1.2 +++ OpenXM/src/asir-contrib/packages/doc/n_wishartd/n_wishartd-ja.texi 2017/03/30 06:16:36 1.3 @@ -1,5 +1,5 @@ -%comment $OpenXM: OpenXM/src/asir-contrib/packages/doc/n_wishartd/n_wishartd-ja.texi,v 1.1 2016/08/24 07:39:41 noro Exp $ -%comment --- おまじない --- +%comment $OpenXM: OpenXM/src/asir-contrib/packages/doc/n_wishartd/n_wishartd-ja.texi,v 1.2 2016/08/25 03:13:54 noro Exp $ +%comment --- 障 --- \input ../../../../asir-doc/texinfo @iftex @catcode`@#=6 @@ -10,18 +10,18 @@ @overfullrule=0pt @c -*-texinfo-*- @comment %**start of header -@comment --- おまじない終り --- +@comment --- 障腟 --- -@comment --- GNU info ファイルの名前 --- +@comment --- GNU info <ゃ --- @setfilename asir-contrib-n_wishartd -@comment --- タイトル --- +@comment --- 帥ゃ --- @settitle n_wishartd @comment %**end of header @comment %@setchapternewpage odd -@comment --- おまじない --- +@comment --- 障 --- @ifinfo @macro fref{name} @ref{\name\,,@code{\name\}} @@ -33,9 +33,9 @@ @end iftex @titlepage -@comment --- おまじない終り --- +@comment --- 障腟 --- -@comment --- タイトル, バージョン, 著者名, 著作権表示 --- +@comment --- 帥ゃ, 若吾с, , 篏罔茵腓 --- @title n_wishartd @subtitle n_wishartd User's Manual @subtitle Edition 1.0 @@ -48,93 +48,93 @@ Copyright @copyright{} Masayuki Noro 2016. All rights reserved. @end titlepage -@comment --- おまじない --- +@comment --- 障 --- @synindex vr fn -@comment --- おまじない終り --- +@comment --- 障腟 --- -@comment --- @node は GNU info, HTML 用 --- -@comment --- @node の引数は node-name, next, previous, up --- +@comment --- @node GNU info, HTML --- +@comment --- @node 綣違 node-name, next, previous, up --- @node Top,, (dir), (dir) -@comment --- @menu は GNU info, HTML 用 --- -@comment --- chapter 名を正確に並べる --- +@comment --- @menu GNU info, HTML --- +@comment --- chapter 罩g∈筝鴻 --- @menu -* matrix 1F1 の対角領域上への制限パッケージ n_wishartd.rr :: +* matrix 1F1 絲乗筝吾狗宴若 n_wishartd.rr :: * Index:: @end menu -@comment --- chapter の開始 --- -@comment --- 親 chapter 名を正確に --- -@node matrix 1F1 の対角領域上への制限パッケージ n_wishartd.rr ,,, Top -@chapter matrix 1F1 の対角領域上への制限パッケージ n_wishartd.rr -@comment --- section 名を正確に並べる --- +@comment --- chapter 紮 --- +@comment --- 荀 chapter 罩g∈ --- +@node matrix 1F1 絲乗筝吾狗宴若 n_wishartd.rr ,,, Top +@chapter matrix 1F1 絲乗筝吾狗宴若 n_wishartd.rr +@comment --- section 罩g∈筝鴻 --- @menu -* matrix 1F1 の対角領域上への制限:: -* 制限した関数の数値計算:: -* 部分分数係数の微分作用素:: -* Runge-Kutta 法の試験的実装:: +* matrix 1F1 絲乗筝吾狗:: +* 狗∽違医よ膊:: +* 遺違緇篏膣:: +* Runge-Kutta 羈荅薑絎茖:: @end menu -このマニュアルでは, asir-contrib パッケージに収録されている, -matrix 1F1 が対角領域上で満たす微分方程式系を計算するパッケージ -@samp{n_wishartd.rr} について解説する. -このパッケージを使うには, まず @samp{n_wishartd.rr} をロードする. +ャ≪с, asir-contrib 宴若吾蚊, +matrix 1F1 絲乗筝ф緇合綣靗祉荐膊宴若 +@samp{n_wishartd.rr} ゃ茹h. +宴若吾篏帥, 障 @samp{n_wishartd.rr} 若. @example [...] load("n_wishartd.rr"); @end example @noindent -このパッケージの函数を呼び出すには, 全て @code{n_wishartd.} を先頭につける. +宴若吾醇違若喝冴, @code{n_wishartd.} ゃ. -@comment --- section の開始 --- -@comment --- 書体指定について --- -@comment --- @code{} はタイプライタ体表示 --- -@comment --- @var{} は斜字体表示 --- -@comment --- @b{} はボールド表示 --- -@comment --- @samp{} はファイル名などの表示 --- +@comment --- section 紮 --- +@comment --- 娯絎ゃ --- +@comment --- @code{} 帥ゃゃ推茵腓 --- +@comment --- @var{} 絖篏茵腓 --- +@comment --- @b{} 若茵腓 --- +@comment --- @samp{} <ゃ茵腓 --- -@node matrix 1F1 の対角領域上への制限,,, matrix 1F1 の対角領域上への制限パッケージ n_wishartd.rr -@section matrix 1F1 の対角領域上への制限 +@node matrix 1F1 絲乗筝吾狗,,, matrix 1F1 絲乗筝吾狗宴若 n_wishartd.rr +@section matrix 1F1 絲乗筝吾狗 @menu * n_wishartd.diagpf:: * n_wishartd.message:: @end menu -@node n_wishartd.n_wishartd.diagpf,,, matrix 1F1 の対角領域上への制限 +@node n_wishartd.n_wishartd.diagpf,,, matrix 1F1 絲乗筝吾狗 @subsection @code{n_wishartd.diagpf} @findex n_wishartd.diagpf @table @t @item n_wishartd.diagpf(@var{m},@var{blocks}) -@var{m}変数の 1F1 が満たす方程式を, @var{blocks} で指定される -対角領域上に制限した微分方程式系を計算する. +@var{m}紊違 1F1 羣合綣, @var{blocks} ф絎 +絲乗筝狗緇合綣靗祉荐膊. @end table @table @var @item return -@var{[E1,E2,...]} なるリスト, 各 @var{Ei} は -部分分数を係数とする微分作用素で, 制限した 1F1を零化する. +@var{[E1,E2,...]} 鴻, @var{Ei} +違篆違緇篏膣, 狗 1F1九. @item m -自然数 +倶 @item vars -@var{[[s1,e1],[s2,e2],...]} なるリスト. +@var{[[s1,e1],[s2,e2],...]} 鴻. @item options -下の説明参照. +筝茯. @end table @itemize @bullet -@item @var{m}変数の 1F1 が満たす方程式を, @var{blocks} で指定される -対角領域上に制限した微分方程式系を計算する. -@item @var{blocks} の各成分 @var{[s,e]} は @var{ys=y(s+1)=...=ye} を意味する. -このブロックを代表する変数は @var{ye} である. -@item @var{blocks} には全ての変数が現れるように指定する. 特に, 一つの変数からなる -ブロックは @var{[s,s]} と指定する. -@item この関数が常に成功することは証明されていないが, 現在のところ, 各変数ブロックサイズが -36 以下なら成功することは証明されている. -@item 出力される微分作用素のフォーマットに関しては @ref{部分分数を係数とする微分作用素} を -参照. +@item @var{m}紊違 1F1 羣合綣, @var{blocks} ф絎 +絲乗筝狗緇合綣靗祉荐膊. +@item @var{blocks} @var{[s,e]} @var{ys=y(s+1)=...=ye} 潟. +篁h;紊違 @var{ye} с. +@item @var{blocks} 紊違憗絎. 鴻, 筝ゃ紊違 + @var{[s,s]} 絎. +@item ∽違絽吾荐惹, 憜, 紊違泣ゃ冴 +36 篁ヤ荐惹. +@item 阪緇篏膣若≪ @ref{違篆違緇篏膣} +. @end itemize @example @@ -157,27 +157,27 @@ matrix 1F1 が対角領域上で満たす微分方程式系を計算する ] @end example -@node n_wishartd.message,,, matrix 1F1 の対角領域上への制限 +@node n_wishartd.message,,, matrix 1F1 絲乗筝吾狗 @subsection @code{n_wishartd.message} @findex n_wishartd.message @table @t @item n_wishartd.message(@var{onoff}) -計算中のメッセージ出力をon/off する. +荐膊筝<祉若後阪on/off . @end table @table @var @item onoff -@var{onoff=1} のときメッセージを出力し, @var{onoff=0} のときしない. +@var{onoff=1} <祉若吾阪, @var{onoff=0} . @end table @itemize @bullet -@item 計算中のメッセージ出力を on/off する. デフォルトは off である. +@item 荐膊筝<祉若後阪 on/off . off с. @end itemize -@node 制限した関数の計算,,, matrix 1F1 の対角領域上への制限パッケージ n_wishartd.rr -@section 制限した関数の計算 +@node 狗∽違荐膊,,, matrix 1F1 絲乗筝吾狗宴若 n_wishartd.rr +@section 狗∽違荐膊 @menu * n_wishartd.prob_by_hgm:: @@ -185,47 +185,47 @@ matrix 1F1 が対角領域上で満たす微分方程式系を計算する * n_wishartd.ps:: @end menu -@node n_wishartd.prob_by_hgm,,, 制限した関数の計算 +@node n_wishartd.prob_by_hgm,,, 狗∽違荐膊 @subsection @code{n_wishartd.prob_by_hgm} @findex n_wishartd.prob_by_hgm @table @t @item n_wishartd.prob_by_hgm(@var{m},@var{n},@var{[p1,p2,...]},@var{[s1,s2,...]},@var{t}[|@var{options}]) -HGM により重複固有値を持つ共分散行列に対する Wishart 行列の最大固有値の -分布関数の値を計算する. +HGM 茲堺ゃゅ怨h絲障 Wishart 茵紊у堺ゃ +絽∽違ゃ荐膊. @end table @table @var @item return @item m -変数の個数 +紊違 @item n -自由度 +怨墾 @item [p1,p2,...] -重複固有値の個数のリスト +茲堺ゃ違鴻 @item [s1,s2,...] -各重複固有値 +茲堺 @end table @itemize @bullet @item -固有値 @var{si} を @var{pi} 個もつ対角行列を共分散行列とする Wishart 行 -列の最大固有値 @var{l1}の分布関数の値 @var{Pr[l1j} を満たし, さらに因子は -ある一定の順序で整列される. +罸憗 @var{yi0^n0(yi1-yj1)^n1(yi2-yj2)^n2...(yik-yjk)^nk} +@var{[[yi0,n0],[yi1-yj1,n1],...,[yik-yjk,nk]]} 綵≪鴻茵 +. , 絖 @var{yi-yj} @var{i>j} 羣, 絖 +筝絎綺ф翫. @item -@var{f} を上のようなべき積とし, @var{c} を定数とするとき, 単項式にあた -る @var{c/f} は @var{[c,f]} で表現される. @var{f=[]} の場合, 分母が 1 -であることを意味する. +@var{f} 筝鴻腥, @var{c} 絎違, 綣 + @var{c/f} @var{[c,f]} ц;憗. @var{f=[]} 翫, 罸 1 +с潟. @item -最後に, @var{c1/f1+...+ck/fk} は @var{[[c1,f1],...,[ck,fk]]} と表現され -る. ここでも, 各項はある一定の順序で整列される. +緇, @var{c1/f1+...+ck/fk} @var{[[c1,f1],...,[ck,fk]]} 茵憗 +. с, 筝絎綺ф翫. @item -部分分数を通分して簡約した結果, 0 になることもあることに注意する. +違膂∞腟, 0 羈. @end itemize -@node 部分分数係数の微分作用素の表現,,, 部分分数係数の微分作用素 -@subsection 部分分数係数の微分作用素の表現 +@node 遺違緇篏膣茵,,, 遺違緇篏膣 +@subsection 遺違緇篏膣茵 -前節の部分分数を用いて, それらを係数とする微分作用素が表現できる. -@var{f1,...,fk} を部分分数の表現, @var{d1,...,dk} を分散表現単項式 (現 -在設定されている項順序で @var{d1>...>dk}) とするとき, 微分作用素 -@var{f1*d1+...+fk*dk} が@var{[f1,d1],...[fk,dk]]}で表現される. +膀違, 篆違緇篏膣茵憗с. +@var{f1,...,fk} 違茵, @var{d1,...,dk} h;憜綣 ( +荐絎綺 @var{d1>...>dk}) , 緇篏膣 +@var{f1*d1+...+fk*dk} @var{[f1,d1],...[fk,dk]]}ц;憗. -@node 部分分数係数の微分作用素の演算,,, 部分分数係数の微分作用素 -@subsection 部分分数係数の微分作用素の演算 +@node 遺違緇篏膣羲膊,,, 遺違緇篏膣 +@subsection 遺違緇篏膣羲膊 @menu * n_wishartd.wsetup:: @@ -387,7 +387,7 @@ matrix 1F1 が満たす微分方程式の係数は @var{1/yi}, @var * n_wishartd.muldpf:: @end menu -@node n_wishartd.wsetup,,, 部分分数係数の微分作用素の演算 +@node n_wishartd.wsetup,,, 遺違緇篏膣羲膊 @subsubsection @code{n_wishartd.wsetup} @findex n_wishartd.wsetup @@ -397,15 +397,15 @@ matrix 1F1 が満たす微分方程式の係数は @var{1/yi}, @var @table @var @item m -自然数 +倶 @end table @itemize @bullet -@item @var{m} 変数の計算環境をセットする. 変数は @var{y0,y1,...,ym}, @var{dy0,...,dym} -で @var{y0, dy0} は中間結果の計算のためのダミー変数である. +@item @var{m} 紊違荐膊医祉. 紊違 @var{y0,y1,...,ym}, @var{dy0,...,dym} + @var{y0, dy0} 筝腟荐膊弱違с. @end itemize -@node n_wishartd.addpf,,, 部分分数係数の微分作用素の演算 +@node n_wishartd.addpf,,, 遺違緇篏膣羲膊 @subsubsection @code{n_wishartd.addpf} @findex n_wishartd.addpf @table @t @@ -414,16 +414,16 @@ matrix 1F1 が満たす微分方程式の係数は @var{1/yi}, @var @table @var @item return -部分分数係数の微分作用素 +遺違緇篏膣 @item p1, p2 -部分分数係数の微分作用素 +遺違緇篏膣 @end table @itemize @bullet -@item 微分作用素 @var{p1}, @var{p2} の和を求める. +@item 緇篏膣 @var{p1}, @var{p2} 羆. @end itemize -@node n_wishartd.mulcpf,,, 部分分数係数の微分作用素の演算 +@node n_wishartd.mulcpf,,, 遺違緇篏膣羲膊 @subsubsection @code{n_wishartd.mulcpf} @findex n_wishartd.mulcpf @table @t @@ -432,18 +432,18 @@ matrix 1F1 が満たす微分方程式の係数は @var{1/yi}, @var @table @var @item return -部分分数係数の微分作用素 +遺違緇篏膣 @item c -部分分数 + @item p -部分分数係数の微分作用素 +遺違緇篏膣 @end table @itemize @bullet -@item 部分分数 @var{c} と微分作用素 @var{p} の積を計算する. +@item @var{c} 緇篏膣 @var{p} 腥荐膊. @end itemize -@node n_wishartd.mulpf,,, 部分分数係数の微分作用素の演算 +@node n_wishartd.mulpf,,, 遺違緇篏膣羲膊 @subsubsection @code{n_wishartd.mulpf} @findex n_wishartd.mulpf @table @t @@ -452,16 +452,16 @@ matrix 1F1 が満たす微分方程式の係数は @var{1/yi}, @var @table @var @item return -部分分数係数の微分作用素 +遺違緇篏膣 @item p1, p2 -部分分数係数の微分作用素 +遺違緇篏膣 @end table @itemize @bullet -@item 微分作用素 @var{p1}, @var{p2} の積を計算する. +@item 緇篏膣 @var{p1}, @var{p2} 腥荐膊. @end itemize -@node n_wishartd.muldpf,,, 部分分数係数の微分作用素の演算 +@node n_wishartd.muldpf,,, 遺違緇篏膣羲膊 @subsubsection @code{n_wishartd.muldpf} @findex n_wishartd.muldpf @table @t @@ -470,16 +470,16 @@ matrix 1F1 が満たす微分方程式の係数は @var{1/yi}, @var @table @var @item return -部分分数係数の微分作用素 +遺違緇篏膣 @item y -変数 +紊 @item p -部分分数係数の微分作用素 +遺違緇篏膣 @end table @itemize @bullet -@item 変数 @var{y} に対し, 微分作用素 @var{dy} と @var{p} の微分作用素としての -積を計算する. +@item 紊 @var{y} 絲障, 緇篏膣 @var{dy} @var{p} 緇篏膣 +腥荐膊. @end itemize @example @@ -492,67 +492,67 @@ matrix 1F1 が満たす微分方程式の係数は @var{1/yi}, @var ...,[[[a,[[y1,2]]]],(1)*<<0,0,0,0,0>>]] @end example -@node Runge-Kutta 法の試験的実装,,, matrix 1F1 の対角領域上への制限パッケージ n_wishartd.rr -@section Runge-Kutta 法の試験的実装 +@node Runge-Kutta 羈荅薑絎茖,,, matrix 1F1 絲乗筝吾狗宴若 n_wishartd.rr +@section Runge-Kutta 羈荅薑絎茖 @menu * rk_ratmat:: @end menu -@node rk_ratmat,,, Runge-Kutta 法の試験的実装 +@node rk_ratmat,,, Runge-Kutta 羈荅薑絎茖 -@code{n_wishartd.ps_by_hgm} では, パフィアン行列を計算したあと, 与えられたステップ数で -Runge-Kutta 法を実行して近似解の値を計算する組み込み関数 @code{rk_ratmat} を実行している. -この関数を, 値が与えられた精度で安定するまでステップ数を2倍しながら繰り返して実行する. -@code{rk_ratmat} 自体, ある程度汎用性があるので, ここでその使用法を解説する. +@code{n_wishartd.ps_by_hgm} с, c≪活荐膊, 筝鴻違 +Runge-Kutta 羈絎茵菴篌取Вゃ荐膊腟粋昭翠∽ @code{rk_ratmat} 絎茵. +∽違, ゃ筝膕上墾у絎障с鴻違2膵違菴絎茵. +@code{rk_ratmat} 篏, 腮綺羆с, с篏睡羈茹h. @subsection @code{rk_ratmat} @findex rk_ratmat @table @t @item rk_ratmat(@var{rk45},@var{num},@var{den},@var{x0},@var{x1},@var{s},@var{f0}) -有理関数係数のベクトル値一階線形常微分方程式系を Runge-Kutta 法で解く +∽遺違や膩綵√幻緇合綣靗祉 Runge-Kutta 羈цВ @end table @table @var @item return -実数のリスト +絎違鴻 @item rk45 -4 または 5 +4 障 5 @item num -定数行列の配列 +絎域 @item den -多項式 +紊綣 @item x0, x1 -実数 +絎 @item s -自然数 +倶 @item f0 -実ベクトル +絎 @end table @itemize @bullet @item -配列 @var{num} のサイズを @var{k} とするとき, -@var{P(x)=1/den(num[0]+num[1]x+...+num[k-1]x^(k-1))} に対し @var{dF/dx = P(x)F}, @var{F(x0)=f0} を -Runge-Kutta 法で解く. + @var{num} 泣ゃ冴 @var{k} , +@var{P(x)=1/den(num[0]+num[1]x+...+num[k-1]x^(k-1))} 絲障 @var{dF/dx = P(x)F}, @var{F(x0)=f0} +Runge-Kutta 羈цВ. @item -@var{rk45} が 4 のとき 4 次 Runge-Kutta, 5 のとき 5 次 Runge-Kutta アルゴリズムを実行する. -実験的実装のため, adaptive アルゴリズムは実装されていない. +@var{rk45} 4 4 罨 Runge-Kutta, 5 5 罨 Runge-Kutta ≪眼冴絎茵. +絎薑絎茖, adaptive ≪眼冴絎茖. @item -@var{s} はステップ数で, 刻み幅は@var{(x1-x0)/s} である. +@var{s} 鴻違, 祉水@var{(x1-x0)/s} с. @item -@var{f0} がサイズ@var{n} のとき, @var{num} の各成分は @var{n} 次正方行列である. +@var{f0} 泣ゃ@var{n} , @var{num} @var{n} 罨≧f壕с. @item -結果は, 長さ @var{s} の実数リスト @var{[r1,...,rs]} で, @var{ri} は @var{i} ステップ目に計算された -解ベクトルの第0成分である. 次のステップに進む前に解ベクトルを @var{ri} で割るので, 最終的に -解 @var{F(x1)} の第 0 成分が @var{rs*r(s-1)*...*r1} となる. -@item 方程式が線形なので, Runge-Kutta の各ステップも線形となることを利用し, -第0成分を1に正規化することで, 途中の解の成分が倍精度浮動小数の -範囲に収まることを期待している. 初期ベクトル @var{f0} の成分が倍精度浮動小数に収まらない場合 -は, @var{f0} を正規化してから @code{rk_ratmat} を実行し, 前項の結果に @var{f0} の第 0 成分をかければ -よい. +腟, 激 @var{s} 絎違鴻 @var{[r1,...,rs]} , @var{ri} @var{i} 鴻荐膊 +茹c膃0с. 罨<鴻蚊茹c @var{ri} у蚊, 腟 +茹 @var{F(x1)} 膃 0 @var{rs*r(s-1)*...*r1} . +@item 合綣膩綵≪, Runge-Kutta 鴻膩綵≪, +膃01罩h, 筝茹c膕上墾羌絨違 +膀蚊障緇. @var{f0} 膕上墾羌絨違障翫 +, @var{f0} 罩h @code{rk_ratmat} 絎茵, 腟 @var{f0} 膃 0 +. @end itemize @example @@ -573,7 +573,7 @@ Runge-Kutta 法で解く. @end example -@comment --- おまじない --- +@comment --- 障 --- @node Index,,, Top @unnumbered Index @printindex fn @@ -584,5 +584,5 @@ Runge-Kutta 法で解く. @summarycontents @contents @bye -@comment --- おまじない終り --- +@comment --- 障腟 ---