福島事故原発の取り壊し方法を考えるスレδ

週刊□福島廃炉
α=1486207162
β=1584849320
γ=1655045111

率直に言ってプリンキピアは見えて来たから、今年いっぱいフルにするんだけど
少しゆとりを持とう。入手している人は、1編の3分の1くらいまで円錐曲線の
作図法を与えよの問題、その後に変なわけのわからない図がいくつか並んでいる
のを見ていると思う。文脈を押さえるとそうでもなくて、右側は2つの物理量の
普通のグラフで左側は解釈に使う同心円とその中の曲線。こんな風に読める。
解釈の方法はケプラー方程式からのホドグラフをニュートンは普通使っている。

だからわけのわからない物は無い。残っているところ（まだ８割ほどは残ってる？）
を攻略して解釈インフラとして社会に提出てか提供しようと思う。
第2編などではなぜこんな言いたいことが入ってこないような書き方をという
印象をみなさん持ち、同じである。そこの分析も少し批判視点も交えてしてみる。

読み方が伝わり、読みにくさの構造性がきれいに解きほぐされると、哲学書の
カントやヘーゲルのように、一般的に解説注釈書が出来ていくものになるだろう。

さらに、なぜこの結論を出せてる？という突っ込み方。
r = sinθの時に逆5乗則だと自信を持って書いている。我々は微分の記号を自由に
使えている時代だからそうだね、と共通推論の形式で人と人との同意に辿りつけるが、
この体系運用して自信が持てないことにはならなかったのかなあ、とか。
よくどんどん進めれるなあと。そういうことを解いて行ってみたいと思う。

今日3月8日からしばらくはバイオをする。

その後でまた戻り系の隅々まで押さえたり、上級の構造化へ向かう。
まず円錐曲線の幾何と二重化法による物理運動の初等幾何、その論証化。
その後で少し話題を変え(PCや建築)、また戻り
円錐曲線を楕円曲線の次元落としとコホモロジーでそれらをつかむ（今狙って準備中）。
その後でまた行って戻りで
モデル理論と代数解析K理論でそれらをつかむ。
どこまで出来るか。最後の行の以外は現実プラン。

昨年最終の予告が臓器移植だった。これをする(2025-12)。
気分的にバイオをいっぱい時間取りたくて、今日のノーカンにしようかな。
もちろん被曝した危ない人を臓器移植で救出することはありうる。
そのためにスレのトピとなるが、しかし他人の本来の所有物を犠牲にしている。
移植はこの問題を離れることはずっと出来ない。

我が国では移植件数は国際的にかなり少なくて、逆に米国やスペインが多いそう。
いいのかなと言う問題もいつまでも付きまとうから、どっちの国も正しいのだろう。
技術の質の向上は目指すが、必ずしも施行症例を量として増やすとは
少なくとも私は思っていなくて、だから日本はダメなんだと言われれれば一面ではそうで。

また次第に技術が向上するから安全性も高まってくる。
そうすると需要が増え、多くの人の本来的所有権を犠牲にする危なさが近づく。
話題もほどほどでいいのだろう。
　
　
さて法律も定まり、心臓、肺、肝臓、腎臓、膵臓、小腸が保険にまで掛かる。
中でも腎臓が多く透析からの脱出意図で登録され10年も待つほどの状況だそう。
他のは2-3年とか。登録してからリスト上で点数評価され順位が上がって行く。
禁忌は大量飲酒癖などコンプライアンスの問題の人と全身病で適応から外れる人。

何らかでリスト上位者に声が掛かる時、ABO赤血球血液型とHLA(白血球血液型)が
重視される。公共の所からではなく知り合いでする時、ABOは乗り越えられ
HLAの抗体が反応し合うかが最も重視される。ここに問題がある時は脱感作療法と
いうアレルギーの時にする方法をして抗体を反応しないように学習させる。
アレルギーの究極はアナフィラキシーだからそこを抑え込んでするとは
かなりの医療技術になったのだとはわかるだろう。免疫抑制薬も用いる。
そうして大抵進む。拒絶反応も時々ある。
まだ初期段階技術だからなのか対症療法(有り合わせの手段で何とかするだけ)である。
だがABO異種の人へは現代では余裕でできる。

元々は外科では腫瘍の外科手術が大量にあり臓器技術が相当に向上している。
移植はその意味では一歩進むだけらしい。心臓以外は腫瘍手術がよくあり
そこのよくできる人が移植にまで進む。血管をクリップで止めて切って
必要な場所において縫合して、後は患者の不快な思いを最小にするように精一杯配慮するだけ。

不思議だよね。不思議だ。何が？どうしてそれが置かれてレシピエント(受け手)
の脳の指示を受けて動いているのか。或る程度疎遠な関係の生物種同士なら
縫合しようが全然関係ないと扱われるだろう。しかし実験して神経とかは
つなごうと特別にしなくても動くのである。人間の体のつくりがそうなっていた。
小児のレシピエントも多いが移植臓器もちょうどのサイズで成長を共にしてくれる。
　
　
我々はそういう手段が可能であるという既に一般化した実践技術からもっと
学ぶことが出来る。手段として持っておいて様々な危ない状態の人への適用を
時に検討するという本来のとは別に前線を広げていくことも。

胃はどうか。歯は。それぞれの骨は。眼球はあるけれど視力には関係がなく。
しかし聴力系はどうか。脳だって出来そうに思う。生殖器も。腕一本とか。
今回10冊以上文献を見たがどうして動くの？と言うことは書いていなかった。
なぜなんだろう。そこの仕組みは本質的だと思う。歯や骨ぐらい簡単そう。違う？
　
　
さて移植を受けた人は免疫抑制薬は一生続く。カルシニューリン阻害薬というのが多い。
カルシニューリンだけ検索して調べておいてくれればいいと思う。
血中濃度を見ながら投与量を調整して決定していく。すなわちずっとその薬が
体内にあるように管理されている。

拒絶反応の再発もある。アレルギーのように免疫の一つの形である。
リンパ球増殖の悪性腫瘍が多く、皮膚とその臓器に一般人よりも多く発生する。
その理由に免疫抑制薬もあるが因子として元が異物だからや赤白血液型の強行。
大抵は臓器の末期の人が適応と許可が出る。余命診断までされるような人である。
倫理問題があるから。肝臓も膵臓糖尿も。骨髄少々でいい白血病は標準治療。

原子炉直下に想定外「消えたコンクリート」　福島第一原発、今も残る謎　#知り続ける（朝日新聞） https://news.yahoo.co.jp/articles/6f6855ab79113f81c7be0229779b5eed0c60ce23

>>125
コンクリートだけが酸に侵されたんじゃないかな
鉄筋は残ってるわけだし

バーコードやQRコードがあり様々な情報を伝えるのに使われている。
この手法で問題の解法を適切に伝えれるか？
積極的開発をする話題。

思うに理数書も散文ばかりである。論理にしようとすると、
言う人は居るものの現実にはアセンブラ言語（の記述的手間さと書かれたもの
が意味を読み取れなくなる危惧）にも思われて中々進んで来ない。
この中間の形があってPC用にするの可能性。
　
　
さてAIが何でも知っているのではなく高校3年の標準ぐらいの学力としよう。
とすると大抵の問題は中々解けない。段階があるが小5で中2で高2で大1で、
昭和のコイン転がし落とし駄菓子屋ゲームのようにコースを外れて多数脱落する。
まあ受験勉強途中の偏差値55ぐらいの学生を想像すればいいだろう。

ここにバーコードやQRコード類似の限られた範囲のデジタルを渡す装置がある。
ヒューマン的な有限さに限定している情報渡し装置である。
（対立概念は理論的な有限さ。指数関数で上限を押さえれるとか言っても仕方なく）
　
　
そのコードマークを読み取ると或る問題が解けるようになる。
関連する近い問題には類推力が働き学力がその影響範囲で付く。
これが学生が学ぶ時にそっくりになり、関連する問題への応用がいかにも人間的で、
教育学習を模しているような状況。

(1)こういう状況をAI-PC体系の土台に於いて作ることは可能か？
(2)コードはどんな言語なら出来るのか？
(3)問題や話題への特異度を上げて行く（それだけが解けるようになる）

廃炉はこうだよとQRコード類似物1つ記号に収まるようになって仕上がるといいな。
幾何の証明なども論理に仕上がるのではなくQRコードになって伝え合い
AI同士はそれで概念発見等と通知し合い、ユークリッドもそう再構成される。

今日は雑談。
3/22は血管・血液・輸血を題材に分子主義でストーリー仕立て(2026-1)
3/29植物(2)、4/5女子化粧(3)、4/12-26理数
5/3本格派分子医学(4)、5/10生物統計(5)、5/17建築、5/24-31パソコン論

ゴールデンウィークを充当したくなる筈と予測して食いこませてる部分あり。
統計はこのスレで出て来ては消えるが次はする。
ファイナンスにマウントするという次の目標が出来ているためである。
実践介護、検査技術、高分子数値計算、テキスト処理型での健康談義集成
みたいなのをまた期間置いた後に。
どれも原子力のおまけで置いておいていいと思う。
　
　
トピは大きく2つ思いついて、輸血をすると抗体が出来る。
手術を受ける時に輸血をしたことがあるか聞かれると思う。
体の中に履歴が残るからである。さらではなく免疫的に違う状態になってる。
何かの時にそのことを勘案して対処を探るからという情報取り。

そのような事情のため今では輸血を積極的には用いない。
輸液や点滴というのは全く別である。が外科と産科でまたICU等でそれぞれある。

別方面のトピとは、血管や血液の周辺で間葉系細胞への細胞転換があるらしい。
再生医療で細胞を操作して変える技術があり、ストレスがそれを起こすという
話題すらもあったはず。低酸素がそれを起こすということになっていて
これは違う動作性を起こし、つまり言ってることは臓器の低酸素で直接腫瘍が発生
することはある。実験で言っている。DNAの傷とはまた別である。
エピジェネ指向の腫瘍で、既存のが悪性化する場面としても注目されている。
我々原子力でそういう機序で起きてる場面が無いかはまた見てみよう。

２０５１年の廃炉、極めて厳しく　デブリ本格回収、見通せず―東電福島第１原発―東日本大震災１５年：時事ドットコム
https://www.jiji.com/jc/article?k=2026031100779

２０５１年の廃炉、極めて厳しく　デブリ本格回収、見通せず―東電福島第１原発―東日本大震災１５年：時事ドットコム
https://www.jiji.com/jc/article?k=2026031100779

今日はインプット偏重で何を書こう。
ほぐれて来たらそっちのが得る物多いからという算段でですね(インプットに比重を置いたの)。

分子生物学を物質分子のキャラを立たせてストーリーで、と予告したが
それが出来たらもう結構なものなわけで、学習用の漫画にしたら
素人が楽しみ食いついて読んで、数百個の知識が入って来るようなものになるのだが。
そのうち出来ると思うが、今日その内容は出来ない。とりあえず夜半までは粘る。
　
　
感想を言うとバイオの一般の教科書は平板な記述が多い。
何がどうと言うような内容が3行に2つぐらいあるような感じで延々と続く。
それぞれが報告と時間をかけた賛否によって決まって行った言述で
なんて大量にあるんだと思ってしまうことも。
どうこれを入り込んだ把握できるかという問題になる。

その問題を読者と一緒に解こうと思う。
思えば趣味の世界もそんな入り方をしていたはずで、ボードゲームの人とかは
すごい細かい分類でそこを理解しているが、外部から見ると、それが相当。

その問題の解決案は、結局は体系をつかむ方法で、或る程度は大量にして
単語や概念は詰めて、収穫曲線も意識して、収穫が減り始めることを
感じれるぐらいが大量の目安で、有機的つながりを自分の方から指摘できる
ように段々なるとホーム感。時間にあくせくはしないで、はまるほどさらに
時間が費やされて、なんとなくそこの人になる。
　
　
忙し過ぎてそんな一つの世界に入り込むような気構えや何やを構成し得ない
と言う人が普通で、だから最短コースで入り込める提供法を思っている。
そうすると、これを自分の研究にしようという人も出てその人の生産で豊かになるから、
戦略としては妥当である。
再生医療には発生学つまり分子生物学が必要で、そこを少し面白く見せられれば。

和算の体系を学んでみよう。

原子核のトピで和算て。
あれこれ調べる→書かないともったいなくなる→書く。
こんな感じで読者に知識が。
あとひそかにあの丸いのを陽子とかに見立てたいと思ってる。
　
　
(1)二円の共通接線に関し
大きさの違う重なりは無く離れる2つの円について、だがほぼ近い大きさでイメージしてもらえば
外側で平行に近く包むようなのと、間の領域でクロスするようなのと
全部で4つの共通接線が描かれることがわかる。

線対称性を考えればそれぞれ2つの共通外接線、2つの共通内接線(そういう名前を付けるとして)
はそれぞれ長さが同じ。さらに内接線と外接線の交点が4つ現れる。

接点と交点との距離、が双方の円について等しいというのが定理である。

証明。円をO1とO2、共通外接線をACとBDとする。AC=BD。
AとBはO1上、CとDはO2上、でそれぞれ接点で表している。
クロスして来て、Aの近くに交点を作る共通内接線がある。これを直線IEHLとする。
EはO1上、HはO2上、Iは交点でAC上、Lも交点でBD上。

円O1から見た線対称性でAI=IE。 同じくHL=LD。
円O2から見た線対称性でIC=IH=IE+EH。 同じくEH+HL=EL=BL。

AC = AI + IC = 2 IE + EH
BD = BL + LD = EH + 2 HL
これより IE = HL。

(2)余弦定理
△の頂点を上A、左下B、右下C、AからBCに下した垂線の足Hと名付ける。

AH^2 = AB^2 - BH^2 = AC^2 - CH^2 = AC^2 - (BC - BH)^2
すると AB^2 = AC^2 - BC^2 + 2 BC･BH
AC^2 = AB^2 + BC^2 - 2 BC･BH

BH = AB cos∠B なのでこう書くと三角関数型の余弦定理である。
ともかく角Bを主人公な場所とした形の余弦定理を証明した。

垂線の足までの長さとの積をおまけ項に、という形で使われる。
プリンキピアにもこういう形で辺の長さの積が導入されてる所がある。
　
　
(3)ヘロンの公式（ギリシャ型証明）
(AH, BH, BC, CA, AB) = (h, d, a, b, c) と書き換える。
余弦定理の稿、するとの行は、c^2 = b^2 - a^2 + 2 a･d
変形して
d = (a^2 + c^2 - b^2)/(2a)
直下の第3等号ではcではなくd^2の方を消去したいがためにこういうくくり出しをしている。

面積の普通の式を書き換えて行く。
△の面積S = 1/2 a h = 1/2 a √(c^2 - d^2)
　= 1/2 a √{c^2 - (a^2 + c^2 - b^2)^2 / (2a)^2}
　= 1/2 √{a^2 c^2 - (a^2 + c^2 - b^2)^2/4}

16 S^2 = 4 a^2 c^2 - (a^2 + c^2 - b^2)^2
　= (2 a c + a^2 + c^2 - b^2) (2 a c - a^2 - c^2 + b^2)
　= {(a+c)^2 - b^2} {-(a-c)^2 + b^2}
　= (a+c+b)(a+c-b)(b+a-c)(b-a+c) = 16 s(s-b)(s-c)(s-a)

但し s = (a+b+c)/2。

内心と各頂点A,B,Cを結んだ線分で△を3つに分けることにより
△の内接円O1の半径をrとして、S = 1/2 r (a + b + c)が言える。
これと(3)の式でrも求められる。
一方、r(a+b+c)を別の方法を使ってa,b,cで表せれば、それが和算型のヘロンの公式とも言える。

(4)ヘロンの公式（和算型の証明）
辺の呼び方の名づけを変える。
Aから内接円接点までの長さをa(2つあるがどちらも長さ同じ)、
Bからのをb、Cからのをcとする。AB=a+b、BC=b+c、CA=c+a。

ABの左外側に△の傍接円(O2、半径はR)を描いてみる。
CB,CAは内接円と傍接円の共通外接線になっている。
ABは内接円と傍接円の共通内接線になっている。

O1からBCに下した垂線の足をD、O2からBCに下した垂線の足をEとする。
Dは△の内接円の接点、Eは△の傍接円の接点である。

(1)の定理よりEB=a、またBD=b、DC=c。

BからO1への2つの接線はBCとBAであり、
BからO2への2つの接線はBAとBEであり、それぞれ線対称性を持っていて
O1BO2は直角となる。

これより角の等しい直角三角形の相似で、R:a = b:r

また相似の関係からR:r = EC:DC = (a+b+c):c

R = a b / r = r (a+b+c) /c

ゆえに r^2 (a+b+c)^2 = (a+b+c) a b c
辺の名替えは一次式なので直してくれれば。

もっと多くしようと思ったんだが出来上がらなくて。これだけ。
楕円の幾何学的中心から近日点を向いて、長さが離心率eに比例するような
レンツベクトルLnという物を作ると、運動のこれまでの考察からは保存量のはずである。
原子の電子に使えるだろうし、もっと難しい系の場合の参考用にしよう。

Ln = G M m 離心率e
Ln = m v×(r×v) - G M m rhat
rhat = r/|r| でベクトルの方向だけを取得したもの
単位次元は kg m^3 /s^2である。
rやvは無言及ではベクトルで、ベクトルの外積と内積は×と･で明記する。|r|^2 = r^2

・実際に保存すること
・離心率とその関係であること
・その方向が中心→近日点であること
　
　
m r×vは角運動量だから既に保存する。
万有引力は dv/dt = - G M |r|^-3 r
d|r|/dt = (v･r)/|r|
右辺はvのr方向への射影という意味。

ベクトルについて a×(a×b) = (ay (ax by - ay bx) - az (az bx - ax bz), …)
　= (- (ax^2 + ay^2 + az^2) bx + (ax bx + ay by + az bz) ax, …) = - (a･a) b + (b･a) a

d(r/|r|)/dt = |r|^-1 dr/dt - |r|^-2 d|r|/dt r = |r|^-3 {r^2 v - (v･r) r} = - |r|^-3 r×(r×v)

d(Ln/(G M m))/dt = - |r|^-3 r×(r×v) + |r|^-3 r×(r×v) = 0

Lnとr×vの内積を求めると、定義右辺1項はr×vと垂直のベクトルのはずで
右辺2項はrの方向でこれもr×vとは垂直のベクトルで0を得る。
角運動量ベクトルと垂直なのでLnは軌道平面内にある。

rとLnの為す角θを真近点離角と言う。一番普通の軌道平面内での近日点方向から測った角なのだが、
代数的定義が実際にそう幾何的なものであることを見よう。

Ln/(G M m)とrの内積を作る。
rhatとrの内積は|r|。
(v×(r×v))・r = (r×v)・(r×v) はr,v,r×vの3つのベクトルについての平行六面体公式。

すると Ln/(G M m)・r = (r×v)^2/(G M) - |r|

|r| {1 + |Ln/(G M m)| cosθ} = (r×v)^2/(G M)
極座標での円錐曲線の方程式で、|Ln/(G M m)|を離心率eと解するのが良いことがわかった。
似た形式が出ただけでてんで違う物の可能性は？とか突っ込む人がいるが、
それがそうなら双対とかさらに豊かになるので話をまとめて指摘してもらいたい。

4行上の方程式で右辺はr×vが保存量。
左辺はθが0のとき{}が最大で|r|が最小となる。
rとLnが角度0のときrが近日点方向ということだから
ベクトルLnの方向が幾何学中心→焦点→近日点を表わしている。

天体摂動や宇宙航行でもこの動きを見るといい場合がある。
原子力のエネルギー問題絡みで宇宙輸送をする手筈なのだから関係する。

運動法則と万有引力ポテンシャル→しかる保存則対称性が示されたが
万有引力ポテンシャルと保存則対称性を人為に初めに自由に与えて運動法則を求める話。

求心力を求める例をもっとやってみる。2/15の内容で、
軌道式 f(r,θ)=0の微分を求めて、dθ/dt=r^-2 (角運動量保存式)と軌道式自体を使いながら、
θとdθ/dtを消去し続けて r^-3 - d(dr/dt)/dt が求心力 F(r)になることがわかっている。
数学的な形状だけ見て、定数倍などの係数はあまり扱わない。必要な人はそれぞれで。

円錐曲線 f(r,θ) = r (1 + e cosθ) - 1 = 0をやってみる。
1 + e cosθ = r^-1
df = dr (1 + e cosθ) + r (- e sinθ dθ) = 0
dr/dθ = r e sinθ/(1 + e cosθ) = r^2 e sinθ
dr/dt = dr/dθ dθ/dt = e sinθ
d(dr/dt)/dθ = e cosθ
d/(dr/dt)/dt = e cosθ r^-2 = r^-2 (r^-1 - 1)
F(r) = r^-3 - (r^-3 - r^-2) = r^-2

軌道が円錐曲線ならそれを求心力が起こしているのなら逆二乗の形と証明されたのだが、
e=0とする退化を辿ってもらいたい。dr/dθが消えてしまい以後が辿れなくなる。
即ちこの方法では円軌道についても楕円軌道のe=0極限と見ることで力の形状を得る流れとなる。
　
　
r = e^θの例をやってみる。
dr/dt = dr/dθ dθ/dt = e^θ r^-2 = r^-1
d(dr/dt)/dt = d(r^-1)/dr dr/dt = - r^-2 dr/dt = - r^-3
F(r) = r^-3 - (- r^-3) = 2/r^3

これは等角らせんである。∵)角θ測る用のθ原点を内的に変えると係数Aに入り r = A e^θ。
シフト角分だけ回転させて縮尺を掛ける。変化物は元の物にぴたり重なろう。
回転と等倍で角度は不変。どこを回転縮尺で注視点選びにしても同じ形状で重なりらせんは等角性を持つ。

角運動量の正指定で反時計回り。r=e^-θはr=e^θと同じ結果を出す。
というのはらせんを描く時縮もうが広がろうが引かれて中心に向いて曲がっている事情は同じ。
もし斥力なら外側へカーブするので内側にカーブするなら求心力。

回転円軌道上で宇宙機同士がランデブーする際のヒルの方程式を基礎から学ぶ。
遠心力・コリオリ力・加減速力・回転面の上下振動が見られる。
上下振動は本当は少し斜めなのだと思うと周期に同期して上下しているとも見れるが解釈。
地球の赤道上空を公転する例をする。zが北極星方向、xが半径外向き、yが進行方向。

さてdA/dt = δA/δt + ω×A。左辺は静止系。δ/δtは回転座標系で見る時間微分。ωは角速度。
Aが回転座標系での外向き方向ベクトルと決め打ちすると、回転に連動して静止でのy成分を獲得して向きが回って行く。
y方向ベクトルなら-x成分を獲得。少なくとも考察微小の時間範囲では右辺第2項の置き方形式で良さそうである。

d(dA/dt)/dt = δ(dA/dt)/δt + ω×dA/dt
　= δ{δA/δt + ω×A}/δt + ω×{δA/δt + ω×A}
　= δ(δA/δt)/δt + δω/δt×A + ω×δA/δt + ω×δA/δt + ω×(ω×A)

二階微分まで片付けたから以後A=rの場合だけを扱う。m掛ける左辺は力。mで割ったまま考察。
… = r'' + ω'×r + 2 ω×r' + ω×(ω×r)
運動系での動作が静止系での力をも受けてどうなっているかの読む方程式。
dω/dt = δω/δt も(Aにω入れ)わかる。
　
　
回転座標系でのr=(x,y,z)と書く。宇宙機をその原点に置いてる。遠い所に原点がある静止座標系ではない。
ちなみにその場合はδrとrが同一物になる。角速度ベクトルω=(0,0,n)としよう。

ベクトル外積計算で
ω×r' = (0,0,n)×(x',y',z') = (-n y', n x', 0)
ω×(ω×r) = (0,0,n)×(-n y, n x, 0) = (-n^2 x, -n^2 y, 0)
加減速はnの変化ではなくxyzで表示されnは定数とする。ω'×rの項は0。

ところでここまでの考察では座標を回転系にしただけで重力を思っていない。重力を入れる。
遠心力 m r n^2 = G M m / r^2 の平衡式からは n = √(G M /r^3) の関係式が付く。
重力は左辺の…部に入る。しかも潮汐力だけである。主要部は前行型で系作りに今使用された。

スラスタ等制御力加速度 + 重力潮汐型加速度 = (x''- 2n y' - n^2 x,　y''+ 2n x' - n^2 y,　z'')

重力ベクトルの位置による変化を探る。
そのために重力ベクトル g = - G M r /|r|^3 = - n^2 r を各方向に偏微分する。
その9成分量と位置変位ベクトルとの内積で、位置による変化がわかる。
n^2とはまだ置いてはいけない。円軌道に固定せずもっと自由な動作を見る段階なので|r|が固定していない。

∂(|r|^3)/∂x = ∂{(r^2)^3/2}/∂x = 3/2 (r^2)^1/2 2x = 3x |r|

∂{g_x/(- G M)}/∂x = ∂(x/|r|^3)/∂x = |r|^-6 {|r|^3 - x ∂(|r|^3)/∂x} = |r|^-5 (r^2 - 3x^2)
∂{g_y/(- G M)}/∂x = ∂(y/|r|^3)/∂x = |r|^-6 {- y ∂(|r|^3)/∂x} = |r|^-5 (-3x y)

結果9成分行列 ∂g/∂r = - G M |r|^-5 (r^2 I - 3 [r,r]) = - G M |r|^-3 (I - 3 [rhat,rhat])
[r,r]等はベクトル同士の全成分積のつもりだがその内容はあらわに書くとその上の物。
　
　
さて微分計算は終わったのでnを使っても良く、またr/|r| = (1,0,0)。
∂g/∂r = - n^2 (I - 3[但xx成分のみ])
これは対角行列である。だから3成分ベクトルに掛ける時も簡単。
重力それ自体は回転座標系の構成に含まれていて、位置による小さな変化が潮汐力的に力の方程式に入るのである。
それは ∂g/∂r・δr = - n^2 (x - 3x, y, z)

前リプ最後の式左辺第2項に入れよう。- n^2 x 的な物が良く打ち消し合い
(ax,ay,az) = (x'' - 2n y' - 3n^2 x,　y'' + 2n x',　z'' + n^2 z)

これがヒルの方程式である。
もちろん夜でも成り立つ。六本木夕方ズとか芸人にいそう。
-2n y'と2n x'はコリオリ力の項。
n^2 zは単振動に見せる項。但し周期が回転周期なので単に軌道が傾いてる。
またx'' = 3n^2 x + …、と書くと半径方向に広がって行く不安定さがある。(本当？)

1リプ半でケプラー方程式を導出する。準備1で次の半で本導出。
長半径a、短半径b、中心焦点距離c、離心率eの楕円を考える。
c = a e のこと。∵)
r = l/(1 + e cosθ) においてlの定数倍時にrは定数倍される。当面l=1で略。
この方程式が楕円であることはr cosθ = xと置いてr = l - e xを2乗すれば見える。
θ=0でr=1/(1+e)、θ=πでr=1/(1-e)、θ=π/2でr=1。

楕円の方程式はx^2/a^2 + y^2/b^2 = 1である。
θは中心ではなく焦点から測る角度であるため上の結果からb=1はない。
θ=π/2でのrは短軸とは並行する別の場所にある物であり半通径と呼ばれる。
しかし 2 a = 1/(1+e) + 1/(1-e) = 2/(1-e^2) であり
対して 2 c = 1/(1-e) - 1/(1+e) = 2e/(1-e^2)。
より所期の件については証明される。

極座標では焦点を、直交座標では中心を原点に置く流儀が普通。
bについては(c,0)と(0,b)の距離がaから決める。b = √(a^2-c^2) = a√(1-e^2) = 1/√(1-e^2)。
ところでl=(1-e^2)と置いてみればスケールがl倍され、x^2 + y^2/(1-e^2) = 1 で焦点は(±e,0)の半通径は1-e^2。
　
　
系の記述を中心から測る角度αで表したいのがケプラーの発想。
しかし例えば縦方向に2倍に延ばすだけで相当点への角度が変わってしまう。
これを標準化し一意感の物にするために単位円周への角度、それを横a倍、縦b倍。
このように2段階に分けて表して行くことにする。
すると楕円周のパラメータ表示は (a cosα, b sinα)。かえってきれいに。

今一般点(a cosα, b sinα)と焦点(c,0)の距離rを求める。
r^2 = (a cosα - a e)^2 + (b sinα)^2
　= a^2 {(cosα)^2 - 2 e cosα + e^2} + a^2 (1 - e^2) {1 - (cosα)^2}
　= a^2 {- 2 e cosα + 1 + e^2 (cosα)^2}
　= a^2 (1 - e cosα)^2
これr = a(1-e cosα)は次々ランベルトOP,OQでも出てる。

ケプラー方程式
少し丁寧に焦点(c,0)に太陽があって楕円軌道上を地球が周回している様子を思い浮かべる。
楕円のaとbはわかっていて、面積速度則(=角運動量保存則)がある。
近日点(x軸上正位置の点)から計測を開始し、今第2象限(左上方面)に来ているとしよう。

中心と焦点が楕円において違う点であることが問題になる。
これまでに掃いた面積を中心から測る角度αで書く式を求めよう。
ナイーブに短軸方向をa/b倍して円にしてわかりやすくする。

すると弧度法面積で1/2 α a^2。それより中心・現在点・焦点で作る三角形の面積を引いた物。
これがこれまでに掃いた面積。それは時間に比例してh tなどとも書ける。

中心焦点距離 = a e。三角形の高さ = a sinα。
これまでに掃いた面積 h t = 1/2 α a^2 - 1/2 a e ･ a sinα = a^2/2 (α - e sinα)
nを適当な定義としてケプラー方程式 n t = α - e sinαが得られた。
　
　
ケプラー方程式からランベルト定理が導ける。
ケプラー用の今のα(離心近点離角)角をAやB化。
三角関数の加法定理
sin(a±b) = sin(a)cos(b)±cos(a)sin(b)
cos(a±b) = cos(a)cos(b)-±sin(a)sin(b)
複号でcos同士は足しsin同士は引くことが多いとするとcos(a)が残る。

a=(B+A)/2=C、b=(B-A)/2=D を入れる。するとB=a+b、A=a-b。(Bのが未来なのでこんな書き方)
a,bを捨てて式の略記という思想でC,Dを使ってる。(以後は楕円の長半径短半径的な別の意味でab)
cosB + cosA = 2 cosC cosD
cosB - cosA = -2 sinC sinD
sinB - sinA = 2 cosC sinD
cosγ = e cosC と定義し、2倍角の公式(実際には使わない) 1 - cos(2A) = 2 sinA^2

軌道上のP→Q移動時間はPとQのケプラー角A,Bを用い n t = (B - e sinB) - (A - e sinA) と書けるが
直上のγとDを使い、n t = (γ+D - sin(γ+D)) - (γ-D - sin(γ-D)) とも書ける☆。という内容。
離心率eが1化して仮想的な抽象放物線上のケプラー運動で問題が表示されるという。

∵) 証明は3つの場面でe cosCをcosγで置換する方が向上することを発見したことから。

n t = (B - e sinB) - (A - e sinA) = (B - A) - 2 e cosC sinD = (B - A) - 2 cosγ sinD

求心力中心OからPへの距離をp、Qへはqとする。PQ間距離をsとする。
p = a (1 - e cosA)、 q = a (1 - e cosB) から
(p + q)/a = 2 - 2 e cosC cosD = 2 - 2 cosγ cosD

Pの座標は(a cosA, b sinA)、 Qの座標は(a cosB, b sinB)。また b^2 = a^2 (1 - e^2)
s^2 = (a cosB - a cosA)^2 + (b sinB - b sinA)^2
s^2/a^2 = (cosB - cosA)^2 + (1 - e^2) (sinB - sinA)^2
　= 4 sinC^2 sinD^2 + (1 - e^2) 4 cosC^2 sinD^2
　= 4 sinD^2 (1 - e^2 cosC^2) = 4 sinD^2 (1 - cosγ^2) = 4 sinD^2 sinγ^2
s/a = 2 sinγ sinD

直近2段落ではcosの加法定理が見えている。
(p+q+s)/(2a) = 1 - cos(γ+D)
(p+q-s)/(2a) = 1 - cos(γ-D)
これは連立方程式だからOP+OQとPQとaからγとDが求まる。

(γ+D - sin(γ+D)) - (γ-D - sin(γ-D)) = 2D - {sin(γ+D) - sin(γ-D)} = B-A - 2 cosγ sinD
この時点で15行上の形になっていて主張☆が証明された。

ベクトル外積についてトリビアを学ぼう。

・a×(a×c)は0ではない。2回目のとこ内積なら0だけど？外積は垂直同士はフルになる演算。
・外積に結合則は成り立たない。
・3つの外積は先に掛け合わせた方のベクトルの線形和。係数は残りの内積で真ん中ベクトルが正。

a×bはεijk aj bkと書ける。εは123かその偶置換で1、奇置換で-1、他では0の、
6パターンでだけ非0を出す記号。暗黙にΣj Σkが略されていて残ったi添え字が方向値。
例えばi=1なら定義からa2 b3 - a3 b2がこの値。

a×(b×c) = εijk aj (ε b c)k = εijk aj (εklm bl cm)

i=1の時、jkは23か32、jk=23ならε123=1、a2 (b1 c2 - b2 c1)
jk=32ならε132=-1、a3 (b3 c1 - b1 c3)

よって、第1成分は、a2 (b1 c2 - b2 c1) - a3 (b3 c1 - b1 c3)
= a1 b1 c1 + a2 b1 c2 + a3 b1 c3 - a1 b1 c1 - a2 b2 c1 - a3 b3 c1
= [(a･c) b - (a･b) c]1
対称性による変数回転でどの成分にも同じ。
a×(b×c) = (a･c) b - (a･b) c　☆

まず単純に置き換えて変形。
c×(a×b) = (c･b) a - (c･a) b
(a×b)×c = (c･a) b - (c･b) a

明らかに結合則は成り立ってない。主張の2と3は示された。
その具体的な例を、a=bとすると、
a×(a×c) = (a･c) a - (a･a) c
(a×a)×c = (c･a) a - (c･a) a

成分を見ないで幾何学的な意味を1つはつかもう。
☆のような一般は言いにくいが出来る人居たら頼む。

nを単位ベクトルとして、
n×(v×n) = (n･n) v - (n･v) n = v - (v,n) n
右辺第2項はvのn軸への射影。
ということは、左辺はn軸からvの先端へ垂線を表示しているベクトル。
一端それを理解して左辺をもう一度解釈すると
｢vn平面に垂直なベクトル｣とn に垂直なベクトルは、n軸からvの先端へ向かう垂線。

六面体公式
abc 3つのベクトルに対し、a×bはab平面に垂直で大きさが|a||b|sinθ
だからちょうどabで張られる平行四辺形の面積を大きさとして持っている。
cはそれと斜めの関係にあろうけれど内積を取るとそのab平面垂直への射影で
(実効高さだけがほしかったのが与えられて)結局平行六面体の体積を得る。
体積という具体な物に対応したので変数回転の公式は従属する。
成分では(εijk aj bk) ci これはさらに自明的に見えているはず。

以上が要点。4元数にこれを使い、8元数をこの体系の7重マルチとして見ようと思う。
8元数の7もこういうεに類似の。(x×x)×y = 0 と x×(x×y) = -y。
結合則の成り立たなさを幾何学的意味を外して抽象公理の当てはめ様子でつかむ。

空間座標で回転は縦ベクトルに回転行列を左から掛ける感じだった。
固有座標で回転は代表基底としてのベクトルらを混ぜる物だった。
実は後者は右から回転行列を掛けている。それを見よう。
だから[Z][Y][1] = [1][z][y]なのである。大文字と小文字は同じ行列でもある。

ベクトルの方を混ぜるような操作は行列になるのか？
|a1 b1 c1| |cθ　s' 0|　 |a1 cθ + b1 s　 a1 s' + b1 cθ　 c1|
|a2 b2 c2| | s　cθ 0| = |a2 cθ + b2 s　 a2 s' + b2 cθ　 c2|
|a3 b3 c3| | 0　 0　1|　 |a3 cθ + b3 s　 a3 s' + b3 cθ　 c3|

見るとaベクトルとbベクトルを混ぜているよね。
左に掛けるとベクトルの基底を変化させている(xをcx+syで置換等)。
右に掛けるとベクトルを係数付けて混ぜている。
回転操作表示行列を左から掛ける時、元のxの係数に、その左上･x+左中･y+･が掛かる状態になる。
右から掛ける時、その左上･第1列+左中･第2列+、が第1列になる結果になる。
前リプはベクトルを横にしているので転置にして本リプと合わせ。s符号まで同じ物。
常に物体固有の軸で回す方が素直である。3行下のθ'係数(前リプ下から2行目最左)でも。
　
　
角速度ベクトルωをオイラー角体系で表示すること。時間微分('で表す)が絡む内容へ進む。
計算方法は瞬間回転軸とφ'等の積和。
(ωX,ωY,ωZ) = (0,0,1)φ' + (-sφ,cφ,0)θ' + (sθcφ, sθsφ, cθ)ψ'
= (-sφθ'+ sθcφψ',　cφθ'+ sθsφψ', φ'+ cθψ')

これに(x,y,z)=Ｃ (X,Y,Z)を表示していた上の結果行列を掛けて(ωx,ωy,ωz)を得る。
そこ(下の3行)は高校1年生の計算過ぎるので自分で。

ωz = (sθcφ)(-sφθ'+ sθcφψ') + (sθsφ)(cφθ'+ sθsφψ') + (cθ)(φ'+ cθψ') = cθφ'+ ψ'
ωx = ( cθcφcψ-sφsψ)(-sφθ' + sθcφψ') + (cθsφcψ+cφsψ)(cφθ' + sθsφψ') + (-sθcψ)(φ' + cθψ') = sψθ' - sθcψφ'
ωy = (-cθcφsψ-sφcψ)(-sφθ' + sθcφψ') + (-cθsφsψ+cφcψ)(cφθ' + sθsφψ') + (sθsψ)(φ' + cθψ') = cψθ' + sθsψφ'
特に
ωx^2 +ωy^2 = θ'^2 + sθ^2 φ'^2

コワレフスカヤのコマ。
最近やったレンツベクトルにしろ宇宙機ラグランジュ点にしろ組合せで新数理を見出して来た事情があった。
コワレフスカヤの模型はI1=I2 = 2 I3とし重心がずれている状況を作り同様に新性質に至る。
主慣性軸の方向にコマ座標を取って慣性モーメントをI1,I2,I3とする。x軸座標aに重心がある。重力は-Z方向に働く。

コマの運動エネルギー T = 1/2 (I1 ωx^2 + I2 ωy^2 + I3 ωz^3) = I3 (θ'^2 + sθ^2 φ'^2) + I3/2 (cθφ'+ ψ')^2
コマの位置エネルギー U = m g h = m g a (-sθcψ) = - m g a sθcψ
Lagrangean = T - U。 I3で割って m g a /I3 = kと置いておこう。
L = θ'^2 + sθ^2 φ'^2 + 1/2 (cθφ'+ ψ')^2 + k sθcψ　☆

力学変数はθφψであるが、d/dt(∂L/∂φ') = ∂L/∂φ が運動方程式であるから
Lがφを陽変数として含まないなら↑は右辺0につき保存量。
qで全変数を代表させてL(q,q')から L'= q' ∂L/∂q + q'' ∂L/∂q' = q' d/dt(∂L/∂q') + q'' ∂L/∂q' = d/dt (q' ∂L/∂q')
これよりq' ∂L/∂q' - L は保存量でエネルギーと名付けられる。
コマの解とはもう1つ方程式を求めて、力学は終わりで3方程式を数学的に解いてどうぞとすること。だからあと1つ。
　
　
では解くがまず運動方程式。保存則への志向があるから∂L/∂q'まででd/dtをしない流儀する。
ψ: {cθφ' +ψ'}' = - k sθsψ
φ: {2 sθ^2 φ' + cθ(cθφ'+ ψ')}' = 0
θ: {2 θ'}' = 2 sθcθ φ'^2 + (-sθφ')(cθφ'+ ψ') + k cθcψ　= (cθφ'- ψ') sθφ' + k cθcψ

φの左辺 = 2 sθ^2 φ'' + 4 sθcθθ'φ' - sθθ'(cθφ'+ ψ') + cθ(- k sθsψ)
ψの式を使って代入してる。sθで割る。
0 = 左辺/sθ = 2 sθφ'' + 3 cθθ'φ' - θ'ψ' - k cθsψ
= 2 {sθ φ'}' + (cθφ' - ψ')θ' - k cθsψ

θの式にi(虚数単位)を掛けsθφ'のに加える。
2 {sθ φ' + i θ'}'
= - (cθφ'- ψ')θ' + k cθsψ + i {(cθφ'- ψ') sθφ' + k cθcψ}
= i (cθφ'- ψ') (sθφ' + i θ') + i k cθ e^-iψ

この形になると常微分方程式でするような工夫をしたくなる。
適切な形を見つけるまでパラメータを置いたり係数を見たり紆余曲折あろうが下のFがかなめを与える。
A = cθφ'- ψ'、 B = sθφ' + i θ'、 C = k e^-iψ を定義。
c = cosθ、s = sinθ、C' = - i ψ' C。 s' = c θ'等。
2 B' = i (A B + c C)

F = B^2 + s C と置く。
F' = i A B^2 + i c C B + (c θ' + s (- i ψ')) C
= i A B^2 + {i c (s φ' + i θ') + (c θ' - i s ψ')} C
= i A B^2 + i s (c φ' - ψ') C　= i A B^2 + i s A C　= i A F

これより (logF)' = F'/F = i A。
Aは実だからFの複素共役をGと置くと微分操作自体は実性で (logG)' = - i A のはず。
即ち足して (log(|F|^2))' = 0。
log(|F|^2)=一定。|F|^2=一定。これがコワレフスカヤ積分で第3の方程式。
　
　
sθ φ'+ i θ' や B^2 + s C が以前の例題達に似てる感を感じさせると思う。
可積分系と言われ、他に戸田格子の例があり、
うまく見つけなくても演繹していけるようなシステムをまた作る。
この問題らの解存在の概念をもっと粒度を細かくしている理論があるのである。
それを工学や物理に様々な箇所から側面支援するのに使える。

コマに関しては方程式が定まるまでとその解釈とをしっかり段階分離。
Kovalevskayaの解釈も解釈中身はこれからだが定性的なことは色々あるのだろう。
眠りゴマをシームレスにコワレフスカヤに移した時に条件を見たり慣性楕円体使う等。

ラグランジアンから開始したが、オイラーの剛体運動方程式。
I3 ωz' = (I1 - I2) ωx ωy + Nz。 Nは外部入力トルクこと、重力や接地点静摩擦力或いは磁界。
は解析力学的な最小作用運動方程式でもある。式を立てて同じ結果に至るはず。

先週は {x,y,z} = Ｃ{X,Y,Z} (|x>=Ｃ|X>)と動軸ベクトルxなどを静成分で表わせて、
{ωx,ωy,ωz} = Ｃ{ωX,ωY,ωZ} ともしたのだった。
{ωx,ωy,ωz} の別の求め方を学ぶ。というのは瞬間回転軸とθ'等との積和
というのが四元数と一回回転表現などに行くと意味を取れなくなるから。
より一般の方法の結果をθ'等ごとの線型表現にしたら前の求め方になっているとまとまる。

ベクトルは波括弧{}、行列は角括弧[](略時有)。今週は縦ベクトル|>と横ベクトル<|も。
小文字の動座標系と大文字の静座標系。基本的には動座標系の量を静座標系で表す話。

回転系の考察。任意の空間ベクトルをLとする。Lは回転系の上にある量である。
(Lの時間微分の静系での成分表示) = (Lの時間微分の動系での成分表示) + ω×L　= () + [ω] L
は承認する。
ω=z=Z軸、L=x(動系では定数)とすると合ってはいる。()の部分は今興味ない。
外積の行列化 [ω] = {{0,-ω3,ω2},{ω3,0,-ω1},{-ω2,ω1,0}} を使ってる。
　
　
座標に関するよりセンシティブな方法を導入する。
5行上の式は d<X|[CT]|L>/dt = <x|dL/dt> + <x|[ω]|L>　☆と変わる。
説明しよう。

x座標基底で測定しているという意味を込めている。
右辺は動系が自己の動きを感知しない量で、dL/dtをxで測定することになる。
Lは(回転系の上にある量とはしているが)いささか抽象的な量で測定するまでその記述法も
定まらないような、基底を定めてから記述法が定まるような物。具体計算は3D内で
(<x|,<y|,<z|)があり、微分量である|dL/dt>と内積し座標数3つ組を与える、が右辺第1項。

左辺は、そもそもx系とX系は相当に方角がずれている。
それを回してから測定(実質は基底ベクトル並びと内積を取る形)をする。
そのために上のような式。Tは転置でもあり逆行列でもある(両者は同一)。

d<x|L>/dt = d<X|[CT]|L>/dt = <X'|[CT]|L> + <X|[CT]|L'> + <X|[CT]'|L> = 0 + <x|L'> + <x|[C][CT]'|L>

第1の等号は|x>=C|X> の転置で <x|=<X|CT。また|X>=CT|x> から最右では<X|=<x|C。
上の第1の等号の右は C^-1 |L>と座標を回して静系で見る量にしてからXと内積しようとしている。
Lを回してからX系基底と内積、Lをx系基底と内積、は同じ物。
時間性が無い範囲ではそれでいい。時間が入ると項が少し出る。Xは静止系だからX'=0。

そのLに付いて回す用の行列は<x|の変換から結合法則も適用にて供給されている。
これより前リプ下から10行目☆と比較し [ω] = [C][CT]'　がω求め方の結論。

ところでオイラー角行列Cについて転置が逆行列は3度の単純回転の合成であるということから或る意味自明。
Cの表記だけにした物から転置との積を計算もミスさえなければ単位行列に結果する筈だけれど。
先週の角速度ベクトルωがこちらの方法で再現されることは未確認だが略。
　
　
以上でコマでもトルクのある物でも四元数や一回回転表現でオイラー角不使用で
式を立て解を得ることまで出来るはずである。
論理を追ってやれば出来るはず。

まずはその角度表記法によるCまでを求め、Ｃからωを求め、
使われ方は、ラグランジアン・オイラー剛体方程式・キネマティック方程式の3通り。
ωまで分かっていればどれも書けるから。
宇宙工学と原子力解析(原子核の回転とMRIなど)用途、また地球の動きで地震解析にするから、
ここでなるべく具体を準備して来る。また来週ね。特殊関数バリバリの地震学を近日中。

さてオイラー角以外の候補は、一回回転表現と四元数がある。
なんだか四元数が定番ということになっていて、大勢敷居をまたげない状況になってるのは
どうなんだ。一回回転表現が合成は不便だが私はいいと思うが。
四元数はどこか隔絶感があってすぐ意味を見失って道に迷った感になってしまい易いのよね。

ところでここでまったり八元数トピをしよう。四元数なんか初等だと強がって学びやすくするために。
八元数の積は247, 375と覚える。他は自然収納。
実数同士、実数と虚数、同じ虚数同士の積は初等的。虚数の積は
ε(123)ijk + ε(145)ijk + ε(167)ijk + ε(247)ijk + ε(256)ijk + ε(346)ijk + ε(375)ijk

外部に出るべき添字のiかkかは巡回することで等しくなり同値。
7つの虚数。対角線を除き半分三角(7×7-7)/2で21個の積パターンがある。3巡回7つの合わせである。
145,167は123の並べ続きで、247と375を覚え他は昇順で正とすると全譜出せる。

結合法則の崩壊を見る。3つ組でない虚数を3つ取ると。
(1 4) 6 = 5 6 = 2
1 (4 6) = 1 3 = -2
この計算には1 6 = 7と7の関わる2つ以外の4つのεを使ってる。
　
　
戻ろう。キネマティック方程式と言うのがある。左辺がθ'等で右辺は時間微分を持たない項
で書かれる式。右辺はωx等は含む。
角運動量は保存し角速度ωは時間保存しないからωは時間変化する。
これはθφψから今しがたまでやっていた方法で再計算されて新しい時間ステップの値になる。

式の形を聞くだけでシミュレーション向けだもわかるね。
角運動量と角速度が本質的に一致していないは力学の有名な話。
特に宇宙機の異常回転や航空機のダッチロールなど、どんどん回っていくのを追いかけれる式。
オイラー角が変化して行くのを何秒後にθφψがどうなっているかわかる。
それを予測して噴射して立て直しの試み。

オイラー角についてはωx = f(θ',…)としたのを先週したから、その三元方程式をθ'等について
陽表示に解き直して得れる。
他の角表示法流儀の場合はまたコメントする。

四元数の基本的な関係式。四元数をp'=[p0,p]と書く。
[a0,a] [b0,b] = [a0 b0 - a･b, a0 b + b0 a + a×b]
a×bの所はi×j=kで全ての組合せがあってその線形和で得れるから。
　
　
四元数の基本的な計算。球関数相当特殊関数を四元数内に定めたいと思う。
([p0,p] [0,v]) [p0,-p]
= [- p･v, p0 v + p×v] [p0, -p]
= [(-p･v) p0 - (p0v+p×v)･(-p), (-p･v) (-p) + p0 (p0v+p×v) + p×(p0v+p×v)]
= [0, (p･v)p + p0^2 v + p0(p×v) + p0(p×v) + p×(p×v)]

p×(p×v) = (p･v)p - (p･p)v = (p･v)p + (p0^2-1) v を入れ
= [0, 2 (p･v) p + (2 p0^2 - 1) v + 2 p0 (p×v)]

回転用四元数が [p0, p] = [cosα, sinα n](α=χ/2) という式を入れる。
= [0, 2 sα^2 (n･v) n + (2 cα^2 - 1) v + 2 cα sα (n×v)]
= [0, (1 - cosχ) (n･v) n + cosχ v + sinχ (n×v)]
直交ではない。(n･v)n はvのn方向への射影。重要なのはここまで。
　
　
外積の時結合法則が成り立たなかった。四元数としては成立することを見る。
四元数の結合則を記号だけで迅速に言う、a'b'= d'= [d0,d] と置こう。

(a'b') c' = [d0c0-d･c, d0c+c0d+d×c]
= [(a0b0-a･b)c0 - (a0b+b0a+a×b)･c,　(a0b0-a･b)c + c0(a0b+b0a+a×b) + (a0b+b0a+a×b)×c]
= [a0b0c0-{a0b･c+b0c･a+c0a･b}-{abc}, {a0b0c+b0c0a+c0a0b}+{c0(a×b)+a0(b×c)+b0(a×c)} - (a･b)c + (a×b)×c]

上ではb0(a×c)以後が巡回対称ではないが、
そもそも巡回対称は不要でaとcについて左右対称なら演算の左右対称はあり結合則の成立である。
(地道にやっても同じ式が出るという意味)。
(a×b)×c=(a･c)b-(b･c)a、-+-と中央ベクトルを残すのだけ正。つまり-(a･b)cの付加で結合則は回復してる。

四元(数)解析って知っているだろうか。複素(数)解析のもじりである。
どうも四元数については-1の平方根が無限個あるという水準の話が多く
それ以上の構造は何？と思われている。ここでそれを述べよう。
これから話を作るようなことばかりで今日は本質的なこと。

三角関数という円を回るのに適した関数が既にあって複素解析で活躍した。
たまたま性質が一致していただけで本当に複素数の住人なのか？というのは
どうなんだろうと今も再検討されていいことである。
電気の振動だとか地震の地殻揺れとかは複素数の実存を表わしているのか否か。

ルジャンドル球面調和関数というのがある。四元数解析はまさにそのまま。
もう優秀な人ならこの一言から一つのゲージ理論クラスの理論作って来れるのでは。
　
　
球面調和関数は、学校で一番最初に触れるのは2s, 2p, 3dなどである。
これは確かに方向を峻別し虚数の間に違う役目を分配することが出来ている。
球面調和関数は最初に見える超幾何関数でもある。
それは3DのSphereSurfaceに制限されたラプラス作用素の固有関数系列。

四元数は4Dの超球面SuperSphereSurfaceに超幾何関数をはべらせて、
はべらせることが出来れば大学専門課程程度の数学にはすぐ成りそう。
超幾何関数の級数定義はΣ(k=0,∞) {(a+k)!(b+k)!x^k}/{(c+k)!k!}。
分子分母(2,1)型で、分子分母(m,n)型を一般超幾何関数という。

結局4次元で動く作用素を持って来てその固有関数を使って方向性解像度を上げる。
そこには特異点が動いて行ったり合流したり無限遠点やリー局所構造や
収束級数になりそうなのにニアミスを起こす漸近級数などが現れ、
力学のベッセル関数も、球面調和→超幾何→合流型→制限でベッセル
という形で抽象的に方向虚数に乗る。これらの現象を全部見ればできるわけ。
微分作用素は代数解析で代数幾何の特異点論もそこにつなぐ。

ニュートン力学→一般相対論→弦理論、と空間処理場の方法が進むと思う。
一般相対論を丁寧にやってみたい。するとあまり載っていない、量子的でない
古典力学型の弦理論が成立するのか、など判定が出来るようになりそう。

原子力は、相手はあまり動かない物であるし、こういうことをやって
いるうちに何とかなるだろう、と。
　
　
2つの項の和がうまく解を与える現象をここで最近よくやって来た。
そういう物はLとBという2つのx微分の微分演算子のラックス方程式
というのが式そのものを与えている。

一般相対論にはカー解以外にゲーデル、冨松佐藤などやはり特殊な解があり
系列を作っていることがある。これが可積分系なのだろうと、だが
そこまでは作れなかったのようなことが書いてある。
　
　
ニュートン力学と一般相対論でそれぞれの可積分系は対応するのだろうか。
KdV方程式を持ち上げた物は何だろうか。そもそもカー解をチェックしたことがない。

ということでこういうことをしばらくやって、一般相対論から落した物として
プリンキピアを見れるような視点を得よう。
おそらくこの方向は宇宙解も新しい物を出すと思う。
またラグランジュ点安定の逆で、共振とはまた別の不安定現象をプラズマに出すかもしれない。
　
　
それと、一般相対論で空間が曲がるけれど、どういう曲がり方をするかが導けると思う。
等価原理と大域的ニュートン性は、それを決めるだけの制約がある。
とするとその論理を明示に書き出すことはすべきこと。

ややこしいアインシュタイン方程式（隠ぺいされた変数でなくテンソル実体）は
複雑な解みたいな計算結果だ、とできれば話がまた進む。

制御工学で言うコントローラは、ビークル即ち、航空・電気鉄道・自動車また船の
振動の低減や破壊防ぎに数理的に直接証明されて有用である。橋や歩道橋(駅のコンコース
やビル間の連絡通路等)も乗り物のようなもの。
無造作に作った物は商品として質が低く、ビルの耐震や、イヤホンの音消しのような
リアルタイムで観測して、対抗制御を投入して、乗り心地を専用に提供物或いは
商品価値に仕上げる技術は現代で普通である。

外乱に有効な対抗制御を入れられる事実に、そうだなるほどと納得できる例示を
ここでしたいし、知識を持っていなかった者が学ぶと、直ぐに現物に応用して
自分の商品を向上させれるだろう。そういうまとめをしてみたい。

そこから建築機械やちょっとした人型ロボットがたわんでしなったりなど、幻滅して
しまうような機械っぽさを、たわんでしなったりがまるでないような、であって
使い易い形の動作物に換えることも出来る。

幻滅さの部分も、制御で対抗入力を入れる制御対象である。
建設機械は廃炉になるし。だからこだわって満足出来るまで作り込むことが有意味。
　
　
機械機器からエンジンに移ると、ロケットの現代化につながる内容があるだろう。
エンジン内には流体(反応性流体)が現れる。
流体にソリトンがあり、界面近くでは際立ち傾向になる。

当然、反応性流体でもソリトンを見て、有り得る形を見ておくことは必要で
そこの数理には可積分系がある。圧縮性・反応性・電磁・粉体・(フェルミ)の5流体シリーズ。
そこでのソリトン構成はとても思いつかないが、先に抽象形式を作ると自動で出る。

反応性流体をセンサと推定から管理する仕組みを作り、半世紀前のロケットから改善しよう。
ところで現時点の最適方法を仕上げ、流体性を落とせば建築のまた水力発電等の管理に。
逆に電磁気的なことを入れることで送電配電
というまた新たにするべきテーマ(流体場から電磁場へ)。

バイオもすべきなのだけれど(1月まで済で今は2-5月残)、
各分野の中級上級をつなぐことに今こだわっている。
数学の一級の研究結果もそういう手合いのことが多いから、上手く行けば何かは。

系統的な方法をつかむと新しい解セットが手に入って来る、というのは
このスレのここまで(に限らず一般の科学の文脈で)でも感じ取られていると思う。
ガウスの正17角形もそうだったね。
その系統的な方法を、可能な上限まで押さえに行きたいのである。

毎週新視点仕上がりなどとは到底行かないだろうが、現在進行形のことを
雑然書きにでもして他の同業研究者推進用にもしてみよう。
だから今日は以後は雑書き。感想や内心やあすなろパートが多くなろうの。
　
　
来週に航空カルマンフィルタをしてみようと予定している(早速あすなろ)。
理学系の人が工学系の文献の記述を読みこなせずついていけなくなることも多く、
そこのほぐしは重要と思っている。それを狙いつつするね。

教科書的にはリャプノフ・リッカチ・逆行列定理という制御工学の行列性質を使い
モデル(数理モデル)も不明な機械対象を、誤差とその共分散を想定して推定して
制御入力を決める。モデルの既知未知で様子が変わるのは統計学で知っていると思う。
分散と不偏分散があって、モデルが既知なら正規分布の分散で、未知ならt分布の不偏分散
というような話。制御の未知部の多い機械での話。

外乱や制御入力と観測と、何通りも誤差を入れ2次形式でとにかくこうコントロール
するんだとセンサ値に対する入力内容を決める式。
機械に数理モデルはわかりきらない。

カルマンという人は有名学者で航空本には必ずこの章があるような感じになっているが、
それを学習しプラント系に。読んで読み取れない本が多いが、難所解体して
どの本にもアプローチできるぜというような提供物になるといいな。

現代制御の可観測・可制御・可検出・可安定は基本性質である。
物理学の端の方で可観測・可制御・可検出・可安定はどうなっているか？
宇宙地平線の向こうや、アプローチに使うエネルギーが質量エネルギーを超えてしまう
クォーク以下の世界。ひもやブレーンではアプローチエネルギーが質量の何京倍。

見ることもアプローチも出来ないのに、存在していることは可能と
数理的なことが制御に似た感覚から言えるかもしれない。
しかし実は制御は線形代数に終始していて議論の水準が低い。物理の可観測を持ち込んで
可制御は双対操作と定義して、するとほしかった結論を得れるのかも。
物理の上限は不可観測で閉じているのかもしれない。温度上限なども。個人的には無関係派だが。
可観測の弱化である可検出は、制御から持ち出すと何だろう。
　
　
制御本ではブロック図や複素平面の軌跡など、どの本も大学の講義ノートの出版みたいで
重複していて興を惹く良トピックという点で少し弱いと思う。そこは構成を変えるべきで
抽象ではない具体的な機械構成の描写で、中における機構の現れ方を例示して学ぶ。

レギュレータで極を配置すると言う。共振周波数を変えて不安定入力だったものが安定入力
の一つでしか無いように安定体系に変えられる重要トピックなのに、
伝達関数の数理の話しかない。機械や建築でつまりこういうことだと数理の方を逆に隠す。

力学に戻りブーメランに付加機器を付けてマジックめいた動作をさせてみるという課題。
ドローンはやってる人が多いがこちらは居なくて面白いだろう。
人工衛星・人工惑星が希望軌道にあることを、制御問題と捉えて数理に当てはめる概念調整
をすることで、航行力学がもう一度揉まれよう。のんびり時間がかかるだけで制御である。
　
　
さて制御には行列の固有値から特異値、四元数、確率統計の中級理論援用、関数解析、流体
などの概念を直ぐに投入したくなる。そのまとめを押さえて定番扱い化に。

高階微分方程式を連立微分方程式に変えれる。力学のハミルトン形式が例。
するとやや難しい概念の連立微分方程式表現があるはず。文脈違うが2階微分方程式の確定特異点。

現場では残留物の取り出しが着々と(かどうか？さて)進んでいると言う。
こちらでは位相空間の方からやろう。手を替え品を替えで他のチームとも協力して
充実手法化狙い。こちらでも機械工学する。そのために位相から積み上げる。

大学1年のあの解析学が、向学者を弾くような複雑な論理は
抽象化すると位相になるとされる。我々の今年のテーマはニュートン力学であり
それは解析学の誕生をしらしめる分野でもあった。普通の発想として、すると
ニュートン力学や電磁気学を解析学を使って書き下すのではなく、
直結して位相空間扱いするというものがあろう。中間役を除くのである。

一応この思い付きに何か形があるかというのが動機。
物理で何々関数を使うのはあっても開集合によってという論理は少ないからと。
問題から回答付けれれば電磁気の記述が広がり電力になる。
唐突な投入は初等物理に位相の感性を入れてみたいからである。
　
　
そこから雑談で書いて膨らまして話を書いていける。
結構な書く内容がある。最近は中級数学の横のつながりを総合的につかんで
制御や統計や電気や分子生物にというのを目指していて、同ペースな進捗で
日曜はそこから思い付きを適当に文章化というスタイルで、最近やったり思ったことが
と、書けるので、かと言ってつながりのつかみが完成しているわけではなく
そのつかみには3-4か月かかりそうで、わりと適当に書き出す方向でしばし。

というのは、安定→ソリトンの問題→量子群←結び目、その裏側で
カオス→量子乱流のマクロ化で質量の案、
不安定は極が複素右半平面にある場合で、乱流を防止する制御の可能性。
制御では極を複数配置して動かすが、その物理版。
記述に出てくる関数の、位相から押さえる解析学表示。
ソリトンは多時間構成の数値計算アルゴリズムを表現する数学とつながる。

つながりが出来て形になる方が使い道はあるなというのは頷いてもらえると思う。

熱伝導　ut = uxx
バーガース　ut + u ux = uxx
ナビエストークスNS　ut + u ux = uxx - px

KdV　ut + u ux = uxxx
変形KdV　ut + u^2 ux = uxxx

今日は熱伝導、バーガース、ナビエストークス、KdV、変形KdVをする。
ナビエストークスNS以外は非常に容易に解かれその解説。
多ソリトンKdVのベックルント構成、多次元KdVのKP階層構成、
戸田格子に相互作用を入れて各KdVを連続極限導出などはまた。
代表っぽく見えているがこれだけではない。超対称素粒子物理のサイバーグウィッテン解も
可積分系であり、そういう難しい系の方を理論形式にて抑え込んだり
コワレフスカヤ階層、サイバーグウィッテン階層などとして上に発展させる方向を目指す。

プラズマと真空管工学には登場して電磁気の自由系でも重要である。
uxxは拡散現象で、NSは圧力作用の項を拡散と別に足した物。

バーガースを解析的に解いた上で、圧力項を少しずつ入れて
乱流遷移の場所がどこにあるか、何を介入できるか、などは
機械工、航空工、プラント工の研究テーマ。
素人でも手持ちのパソコンでその研究ができるはず。
　
　
uやφを場の量としてtとxを変数とする。微分は右下に単に書いて表す。
φが多いが個人的にφより簡便にuvを使う。タイプも楽だし。
uxxは多次元では△uと書かれる。KdV系3階多次元は2階ほど明らかではない。

熱伝導方程式 ut = uxxの意味。拡散一般であり最初が熱力学だからの名。
放物線の底などでuxxは正の量。式は微小時間後への変化がuxxと主張する。
つまり周囲より低濃度の所は微小時間後に値が増大するような働きを受ける。
ゆえに拡散でありなだらか化な緩和過程のニュアンスが見れた。

シュレーディンガー方程式は i ut = uxx。
導出は、E = (1/2m) p^2 というエネルギーの式に、E u = p^2 u と
作用素と作用空間の含みを入れた式にして、置換 E=-i∂t、p=i∂xをする。
置換の心は u = exp{i (E t - p x)} と想定しその場合なら等値。
u波動型の仮定と作用という思想から出す式だが形は熱伝導の虚数版。

バーガースとKdVも良く似ている。右辺が2階と3階の違いだけ。
本当はそこにexp(∂/∂x) uがあり、基本部分と最初の非自明部分と解釈する。
戸田物理の発想の起源。

KdVの解析解はソリトンだが衝撃波とは関係するか。衝撃波の第一近似はソリトン。
しかしその先は今だにわからない。さらに精密化できれば航空工学に使えて
技術的進展を期待できるため鋭意研究されている。みんなも読んで
自分なりの先のことを考えて貢献してくれれば航空宇宙も進むだろう。

またこの分野の微分方程式の式変形をする。それらをもう一度集めて分析すると
新しい特殊関数が取り出され何種類かのパンルヴェとなる。
　
　
熱伝導方程式の解→。u(t,x) = t^-1/2 exp(a t^-1 (x-x0)^2)。
uに関して線形なので初期値についての重ね合わせ。
uとして基本解のデルタ関数の場合だけを説明すれば十分。
上記はu(0,x) = δ(x-x0)からの時間発展解。検算して確認しよう。

適宜x-x0=cやt^-1=τと置いて見易く。微分に使う変数だけ間違えなければok。
aがカチッと定まれば実際に解だったと言ってよいだろう。
ux = 2 a τ c u　(せめてここまでは素人でも把握。指数関数の中を微分して前に出してる)
uxx = 2 a τ (u + c 2 a τ c u) = a τ (2 + 4 a τ c^2) u

ut = -1/2 τ u - a τ^2 c^2 u = -τ (1/2 + a τ c^2) u
a = -1/4とすればいいはずで実際その時τ (-1/2 + 1/4 τ c^2) u
解になっていることはこれで確認出来た。

バーガース方程式の解→。ut = uxxの解uを使い v(t,x) = a (log u)x = a ux/u。
同じく検算だけしてaが納得レベルで綺麗に決まることを見る。
vt = {a(log u)x}t = a (log u)tx = a (ut/u)x

- v vx + vxx = (-1/2 v^2 + vx)x = [-a^2/2 ux^2/u^2 + a {(uxx u - ux^2)/u^2}]x
a=-2とするとux^2の項が消える。 = -2 (uxx /u)x
uの仮定から3行上と一致していて検算出来た。
　
　
変形KdV方程式の解→。u(t,x) = 1/cosh(x-t)
双曲線関数はsinhとcosh。sとcと略す。s'=c、c'=s、c^2-s^2=1。
ux = - s/c^2
u^2 ux = - s/c^4 = - s(cc-ss)/c^4
ut = s/c^2

uxx = (2ss-cc)/c^3
uxxx = -3s(2ss-cc)/c^4 + (4sc-2cs)c/c^4 = (-6sss+5scc)/c^4
これより
ut + 6u^2 ux + uxxx = {scc - 6s(cc-ss) + (-6sss+5scc)}/c^4 = 0
係数と符号の問題はuの事後(？)修整で。
　
　
KdV方程式の解→。変形KdVの解をuとして v = u^2 + ux。
u^2/2 → u ux → u uxx + ux^2 → ux uxx + u uxxx + 2ux uxx = u uxxx + 3ux uxx
繰り返しx微分を先に見た。さらに単に実直計算。

vt - 6v vx + vxxx = 2u ut + uxt - 6(u^2 + ux) (2u ux + uxx) + (2u uxxx + 6ux uxx + uxxxx)
= 2u ut + uxt - 12u^3 ux - 6u^2 uxx - 12u ux^2 + 2u uxxx + uxxxx

(∂x + 2 u) (ut - 6u^2 ux + uxxx) = utx + 2u ut - 12(u ux) ux - 6u^2 uxx - 12u^3 ux + uxxxx + 2u uxxx
両者は一致している。uが変形KdVの解ならvは(元祖)KdVの解である。

Lax形式という概念。
観測作用素L、時間発展作用素B、がありどの波にも使えるとする。
固有値λが時間保存する時どんな式が成り立つか。
以前のようにtによる微分を右下にtを書くだけで表す。
L φ = λ φ
B φ = φt
λt = 0

Lt φ = λ φt - L φt = λ (B φ) - L (B φ) = B L φ - L (B φ)
最右でλはスカラーなのだからBの内側に入りφとの関係でLに変われる。
φは任意なので Lt = B L - L B という式が状況を特徴付けている。
　
　
KdV方程式 ut + 6u ux + uxxx = 0。
後で解釈や工夫してくれればいいからここではuとφを使う。
L φ = φxx + u φ
B φ = 4 φxxx + 6 u φx + 3 ux φ

B (L φ) = 4 (φxx + u φ)xxx + 6 u (φxx + u φ)x + 3 ux (φxx + u φ)
L (B φ) = (4 φxxx + 6 u φx + 3 ux φ)xx + u (4 φxxx + 6 u φx + 3 ux φ)
引き算なので見え透いている所は始めから落として
上-下 = {4 (u φ)xxx + 6 u (φxx + u φ)x + 3 ux φxx} - {(6 u φx + 3 ux φ)xx + u (4 φxxx + 6 u φx)}

= {4 uxxx φ + 12 uxx φx + 12 ux φxx + 4 u φxxx + 6 u φxxx + 6 u ux φ + 6 u^2 φx + 3 ux φxx}
- {6 uxx φx + 12 ux φxx + 6 u φxxx + 3 uxxx φ + 6 uxx φx + 3 ux φxx + 4 u φxxx + 6 u^2 φx}

= uxxx φ + 6 u ux φ
最後の等号では消えて行って整理されることの確認を。
φの意味はここでは作用素が掛かっていく任意関数のはずだった。
全部φを外してしまい B L - L B = uxxx + 6 u ux。
戻りL = ∂x^2 + u。微分作用素は時不変。Lt = ut。よって上記に定義したLとBでLax形式はKdV方程式。
LとBの中を符号変えたりはあるが463はどの文献も同じでKdV系の時間。

変形KdVと元祖KdVがあり変形KdVの方は1/cosh(x-t)の微分方程式化で作られていた。
この解かれ方を見るとソリトンの代表であるKdVの解かれ方はAIが見つけて来れよう。
coshは中央が凹んで両脇で無限大になる関数で逆数は釣り鐘型の1つ。
その方程式はこうこう、その記述変数を使っての2次関数や何かの微分を足したものは
試行錯誤サーチと数値計算を内部で確認的に使えて元祖KdVの解析解にまで辿り着ける。

変形は関数志向、元祖は式形志向、二段階あったことも啓発的。
そういうシステムを作り、四階でさらに難しくなってるKP方程式を解かせたり、
解がウェーブレットになる方程式を作り、などでAI数学の一つの分野が出来る。
そのシステムは非線形の波よりも問題がさらに難しい一般相対論の解探索に使える。
また統計力学の様々な問題がある。解ける突入点がいくつもあるだろう。

ちなみにut = uxは波動方程式。utt = uxx + sin(u)はサインゴルドン方程式。
単振り子が連なっている時にサインゴルドンが現れるとあった。
物質において場の理論では「各時空点」でこの力学が現れ、調和振動子として
素粒子の生成演算子を導くのが場の量子論な所が、サインゴルドンやKdVがそこであると
より複雑になり、通常の素粒子描像ではない理論が構築できる。
すると深い基礎理論はそういう形になっている場合あるかもしれないしまた
陽子や中間子など複合粒子の生成消滅場の量子論という分野で使える可能性が多いにある。
基礎でも応用でも可能性を持つのが数理物理のようである。
　
　
話題トピックとして基礎と応用における乱流の可能性を先に述べる。
プラズマにおいて乱流のために電気伝導係数と熱伝導係数がかなり増大する。
乱流が無いときは位置が定まって気体型をした固体のような状態だが、乱流では
本質的に動き回り、担体輸送する。これは応用の話題。
定性的に語られるだけの、この現象に数理で正確な数値を当てて、乱流度合いを調整
したりの工学が作れるといいだろう。

基礎物理において質量の起源が同じ可能性があるとされる。

KdVと流体を比べる時、KdVは次々と解析化しているが、ナビエストークスは難しいまま。
ここで啓発的に戦略が誘導され、ナビエストークス現象の傍まで行けるバーガースは
解かれるのだから、解けるような変形ナビエストークスを作ったらという案。

Lax形式により何らかの時間性の方向で、L作用素の固有値は保存し保存量である。
項や数式の形を投入して試行錯誤するシステムを作り、Lax形式に則るという制約を守りつつ
バーガースに項を増やして行く。いざその路線を決めれば、何項も入れてナビエストークス
の相当の近似が出来ると思う。それは流体をよく表しかつ可積分だろう。
個人的にも勉強系が一段落終わったら考察してみようと思う。ロケットとか超流動のために。
　
　
方程式の解と言ってもエネルギーで見ると色々な場合がある。
この視点で複合化していくことも、どの場合にも有用だろう。
方程式の解が求まった。最も高級な言葉では測地線軌道が求まった。それがエネルギーの谷
を辿って行き、少しの揺らぎに戻って来るものだという保証はどこにある？
求解時にそれを考慮に入れていないから何も保証されない。
エネルギーでは谷底が望ましいが、一方向に落ちる踊り場の時もあり
尾根道化して両側に落ちるようになっている場合もある。

そもそも乱流は、ナビエストークスのパラメータ値により、エネルギーでは
谷底型が尾根道型に変わるために起きるとされている。
変形ナビエストークスでそれを示せるだけで結構いいだろう。
いつかオリジナルと主客入れ替わる理論につながるかもしれない。
もともと宇宙の発生も乱流の尾根道型化と同じである。その条件置き。

さてこの戦略では乱流も可積分性を保ったまま、解のエネルギー構造の変化で表現されている。
ナビエストークスの乱流が相対化される。どの方程式の乱流かが雰囲気にも表れる。
そこまで出来れば実用もあろう。くだんの賞金問題も、本物は解けなくても
かなり近い近似で解の定性的なことは全部証明されましたとなるかもしれない。
過冷却や過飽和のようなレイノルズ数を動かす工学の手がかり。

基礎の話をしよう。基礎化学(これはプロセスビッグデータ)もするが基礎物理。
乱流が起きると抵抗が発生する。量子やさらには虚数空間で乱流が起きていること
になっているという設定は、素粒子の質量を与えるに十分な過程である。
結論としてはそういうことで、まず実在乱流には音がする。
乱流を全く発生させずに飛ぶ航空機があれば音がしないはずである。

その音の発生は抵抗の発生と読める。多分。枠組みが定まったら証明していく問題。
量子性により実際に起きている物が丸まる。
流れは虚数方向にあるのだ、というのも界隈ではありがちな話。
そんな状況でも統計量は、実数世界の理論を流用した値であるのもあるともないとも
言えないが、理論の論理の筋としては想定される。
　
　
基礎物理にその過程があり読み取ると、ラグランジアン質量になる。
これならそれが、ファインマングラフの計算などで求まらなかったのが理解できる。
素粒子の質量がこれまでわからなかったのは、どこかに何かの乱流があって
マクロ値化する様子として見なかったからだと思う。

諸乱流がマクロ量化する様子を、数理としてうまくつかんで求めればよく
それを想定乱流に適用すればよいとなる。
確かに乱流には規則があるはずなのである。コルモゴロフ-5/3乗則というような話が
教科書には出て来るが、自動車でも航空でも発電所の各プロフェッショナルでも
計測からこうだと決めつけて実際に信頼運用をしている。

質量もこうして乱流性から発生するマクロ値抵抗と思う。
　
　
我々の設定では乱流も可積分系で近づいて行って表わすし、離散版である戸田理論
では自由なリーゲージ群を乗せる余地がある。すると質量がSU(2)やE系のゲージ群の
表現粒子の物であっても発生させられる。
ヒッグスとの乱流的な反応があるのだろうか。それとももう少し展開した先にあるものか。
ひもから現象論を出す指向性が一つ増えたと言うべきだろう。

これはまだ出来ていないし教科書でそういう視点で書いても無いが
KdVではラックス形式のLとBが、ネーターの保存量演算子なのではと
つかみのそのアイデアは自然と思わせる。現実に保存量をソリトンで具体化してない。

コワレフスカヤには保存量があった。KdVには方程式と解析解があった。
アーノルドの定理は、解析解があること⇔保存量があること。
ネーターの定理は、保存量があること⇔連続変形で不変という対称性があること。
　
　
見えてきてる？言いたいこと(読者に話し掛けてる件)。

系記述用パラメータの連続変化は微分作用素の言葉で書かれる。
即ち位置、温度、電磁荷、それらの微分作用素、さらにそれらの多項式や
関数の引数に任意に取らせたもので書き得る全体、が連続変化を一つずつ書く。

全ての連続変形には、微分作用素の或る多項式関数による表現が一つ対応する。
　
　
初等的にx ∂y - y ∂xは、z軸回りの回転に対応する微分作用素含み多項式。
KdVのLとBはまさに微分作用素の或る多項式の形をしていた。
場合によってはcos(x^2 ∂x)のように微分作用素が関数の中に入っている時も。

以上の思考材料をまとめると、KdVにはある連続変形で不変な対称性が
その解に対してみられるだろう。それは冒頭段落からはLとBの何かの可能性が高い。
KdVでは時間が変なのだったのであり定理をつないで求めなければいけない
ものになるのかもしれない。

一方時間が変であることを許容するのでアルゴリズムという数理が
物理のような顔をして微分作用素による分析を受けることが可能になる。
コワレフスカヤのラックス形式は、定理のつなぎを整理したら簡便導出が見えるだろう。
すると晴れて可積分系の一員と見れる。その方法で多くを扱うことで教程水準の話でも広がる。
数か月前にやった逆二乗重力系のレンツ、ラグランジュ点、冨松佐藤、イジング類のもか。

@dabosuke__

非常に残念なことに福島県白河市表郷でまた猫ちゃんの虐待がありました。

今年の春で推定18歳になったばかりの穏やかなおばあちゃん猫だそうです。
地域猫として長らくお世話していた方たちがここ数日姿が見えないことを心配し探したところ、亡くなっているのを発見しました。

草刈機や鎌のような鋭利なもので背骨も内臓も見えるくらい切られ草むらに遺棄されていました。

前々から猫ちゃんのご飯が投げ捨てられたりお水をひっくり返されたりしていたようですが、犯人の目星はついていたものの、近隣トラブルをたびたび起こすタイプの家族のようで報復が怖くて注意したり話し合ったりはできなかったようです。

今回犠牲になった猫ちゃんはもともとは飼い猫でしたが飼い主に置いていかれ、そのまま地域猫として見守られていたそうです。

最初から最後まで、人間によって苦しめられました。非常に不快。

警察に相談はしたそうですが反応的にとくに動いてはくれなさそうとのこと。
当方からも各関係機関に連絡を入れようと思います。
深夜の迷惑行為で該当宅とトラブルになったことがある方が一生懸命お世話をしていた子らしく、糞尿被害等その猫ちゃんから何かをされて恨みがあってということではなく完全にお世話していた方への嫌がらせとしてやったのではないかなと感じる状況でした。

https://x.com/dabosuke__/status/2070044500349919549

ホモロジー代数ってのをやってみようと思う。意義は途中どこかに。
微分方程式の解が作っている構造をトップダウン側の視点から見るということかな。
界隈で超幾何・パンルヴェ・保型の正則領域で成り立つ連接層について
導来圏をスペクトル系列で解体してどうのこうのと言うが、これを脱掛け声・形にして
ニュートン力学化を目論む。読者を置いていかない解説をする。

我々の考えるのは関数の為す集合である。これをぼんやりと粗視化して
内部構造を見なくしてしまって、其処に土団子程度のカタマリがある程度に読む。
共通して使う性質。関数を集めた物はアーベル群の構造を持つ(8行下)。
重力場の微分方程式の解が宇宙なのだし、時空でも何でもこの文脈に掛かる。
　
　
入門の典型は、点を複数指定して其処で特異性を持たない有理関数全体の集合。
有理関数は因数分解されるが、該当点aに対する(x-a)^-1因子を持たなければいい。
有限領域以外では値が0の関数とか、逆に全領域で微分可能な関数とか
無限遠の両側から同じ値に来るとか。幾つもの変形趣向がある。

そのような通常の集め方（特異性の少ない側という）をした関数の集合について、
fとgがそうならf+gもそうだろうし、0も入ってる、-fも入ってる。
f+gはg+f。自明なようなこの和の可換性をアーベル性と呼ぶ。

和で書いてるから連想ゲームしてしまって自明な気がするだけのもので
抽象公理と適用対象の想定で抽象だけを先に整備する立場だから、
積や何だかわからない後から指定するつもりの演算と書くと連想ゲームをしない。
　
　
関数の集合をぼんやり粗視化する。次に準同型写像を入れる。
ここだけ注記するが我々の群の要素実体は関数なのだけれどその思いも捨てる。そして。

準同型は関数の集合から関数の集合への写像という、メタ性を帯びているもの
なのだけれど、取り立てて言うことのないメタから通常に落とした関数と思う。
構造物を点と思考や形式の中で格下げに置き換えながら先に進むということだね。

設定は、x+yなどで閉じている群Gなどと、その間の準同型写像fなど。
記号の抽象性・関数性をこのように１ランク落としてから開始する。

準同型写像とは線形写像のこと。演算が一個だけなのでスカラーも無く
f(x+y) = f(x)+f(y)
+と書いている演算は抽象物で他で置き換えられる。和の行き先は行き先の和が標語。
言っている通りxとyは関数でfはメタ関数(準同型写像)である。

関数を素朴に集めて来たらこのようなことは成り立っているだろう。
関数の和は普通は、任意の引数sについてx(s)+y(s)=z(s)となるような新関数z。
普通に集めて来たらそのzでもそういう性質が成り立っている内容を明記している。
関数の集合から関数の集合への写像は、演算一つを前後どちらでしても同じ値に
なるような写像だけを考察する。
この辺で心に混乱が起きているだろうが構造を格下げして点化して、また
余計な疑問を探そうとしないで性質の間の約束を動かしていくのだなと受け入れることがこつ。
（混乱とは準同型写像とは足し算等の演算ではないのか…だろう。そうどこか別の空間等へのマップで演算ではない）

(関数の為すアーベル)群と準同型写像を大量に集めて動かすのがホモロジー代数である。
　
　
もう一つの重要な例は多面体である。
四面体をイメージして点・直線・面・体胞の４段階次元の物それぞれが個別にあると見る。
各対象から整数Ｚへの写像がある。通常の意味では１か土俵から消えた０か向きが逆な−１
だけだが任意整数を値にとる写像。
するとそのような写像同士は足し算ができる。

各対象を引数とし引数毎に整数を返す写像。これが（考察している群の）要素である。

境界準同型写像というのを定義する。∂(でる)と名付け、△ABCに対しAB+BC+CA。
その意味は、∂は写像に写像を対応させるメタ写像だが、∂(f)=gと書く場合、
f(△ABC)=aなら、g(AB)=++a g(BC)=++a g(CA)=++a。
g(AB)の値はこれだけではなく別の所から来てそれらの総和になっている。

特異点や関数の有界性の話と、多面体の話。以上で２系統方面を出した。
ここから抽象化が何段階も進むが、先へ進まないでも使えそうを感じ取れればいいな。
そして此処までの抽象世界への入り方の部分は了解して貰えたと思う。

想像するだろう通り解領域の多面体近似やそこからの構造判定をして解解析をする。
楕円関数が解なら複素数に埋込準同型するとトーラスと位相的同型の計算。
楕円関数はコマ解である。単振子・球面振子というのもある。コワレフスカヤとラグランジュ点は超楕円。
近似であるからサンプル点を取って乱数からその「形」をシミュレーションで数値計算同定が可能。
有限要素法の中級理論にも出て来る。

これが微分方程式の話をしていてこんな代数が出てくる理由。
より物理に近い分野では超幾何関数のこの辺の手法による記法を近日中に。
すると摂動や鞍部点法などでもちょうどその場合に落ちて来る題材例が増え力学用に。
解析方面と幾何方面という感じの２系統の出し方である。
　
　
ところで丸で話変わるが重力場の方程式は非線形である。バーガースやKdVも非線形だった。
uが未知関数で、u uxはuについて2次で、或るuが方程式の解でも、その定数倍が解とは言えない。
それは或る項についてはa倍、他の項にはa^2倍等となって等式には辻褄が合わなくなるから。

しかし非線形方程式にも保存量か、もっと一般に時間とは限らない連続変形で
不変な量が存在する時に解がある。そんなニュアンスの総合的には未完成の
定理と理論体系があるのだった。先日までやったことである。

我々は線形をやって其処からの変形で、解ける非線形の数式集一式を得る方向を想定している。

即ち熱伝導からバーガースや変形KdVからKdVの変数変換。
其処にあった線形から或いは別の非線形からの変換による解明法。
整備されたホモロジー代数からの変換は、変換先に解の新性質を出せる可能性が高いだろう。

これが今回ホモロジー代数なわけだが、（抽象化と果実採りの）先は長いから進めよう。
文脈的に必要感でちょっとこんな話を始めてるわけなのだ。

中身の実際に入って行ってみる。個人的には調べ直さなければいけないが
少しでも書いて置く方が復習の回数が増えいいだろう。
様々な設定で(関数の)可換群が大量に作られることがわかった。そこから始める。

f:X→Yを可換群Xから可換群Yへの準同型写像とする。
Yの部分集合であるim(f)を、Xの元の像の全部の集まり。
Xの部分集合であるker(f)を、fでYの0元に行く全体の集まり。

直ぐ上でkerとimを意味で定義したが、意味を捨てて矢印のこういう関係が
あるという性質で特定されるようなXなりYなりの部分集合だ、
という言い方がある。もちろん其処までする。もっと先までする。
随伴だとかいうんである。段階を踏んでするけれど。抽象の階層がある。
　
　
整数で添字付けられる可換群の系列で隣接間には準同型写像が存在しているとする。
X(n-1)→X(n)→X(n+1)→X(n+2)→
f(n): X(n)→X(n+1)
fの添字は写像元の方の対象に合わせるのが流儀である。

すると何回もfを辿り番号を進めて行くことが出来る。
X(n)の要素はf(n)さらにf(n+1)を課すると全てX(n+2)の0元になるような性質。
この性質があるとき系列を複体と言う。

多面体の∂で△ABC→AB+BC+CA→(A-B)+(B-C)+(C-A)。
基底についてこういう写し方をしそれに関数的整数値を取らせた物。
ここでは基底が本質的に消えてしまい∂で構成される構造は複体である。

複体のときim(f(n-1))⊂ker(f(n))⊂X(n)。
なぜなら左からf(n-1)で来た物は次のf(n)で0に行く。この論理を集合で表記
したものが⊂という包含表記だから。

このときim(f(n-1))はker(f(n))の部分可換群であり、そのとき商可換群というのを
作ることが出来る。H(n,X,Z) = ker(f(n))/im(f(n-1)) と書いて
複体XのZ係数n次コホモロジー群と言う。Zは整数の意味。
多面体の時∂は添字が降りる方で添字の方向性によりホモロジー群と言う。

商可換群とは。部分可換群B⊂Aに対して、
Aを x+B = {x+b| b∈B}という形(x∈A)の部分集合の重なりの無い合併と書くことが出来る。
∵)
x∈Aの全部を取ればいいからAの一部が埋まらないことは無い。
y,z∈(x+B)について、z=x+b' = x+b+(b'-b)∈(y+B)。bとb'はy,zに対応して取ったとして。
少しでも重なりがあれば同じグループになることは、上の行から重なりの中にyを取り
他のはy+Bの中に入る。またy+Bの任意の他の元も構成の中であって抜け落ちる状況は無い。□

集合x+Bと集合y+Bに和x+y+Bと定義する。これはきちんと定義されている。
∵)
x+Bの代わりにx'+Bのような書かれ方をするとxとx'の差はBの要素。
すると和の形式にしたときおまけで包括的にするBに吸われてx+y+Bとx'+y+Bは同じ物になっている。□
以上で集合を元、集合と集合についての和演算。これで商可換群と呼ぶべきものが存在している。
　
　
短完全系列。0→X→Y→Z→0。
XYZ0どれも可換群とする。0は0だけの物。
左からのimと右へのkerが何処にても一致しているとする。
im⊂kerは複体だが、im=kerを完全と言う。

次のような設定があり標準的である。本当は今日これの元追いかけ証明を書こうと思って用意
したが説明までになった。係数(多面体の場合関数の行先に相当した物)が整数である記載は省略する。
各可換群X等に標準的な分解展開列があり(多面体の場合を思えば何かありそうはわかる)
それは上の意味のコホモロジー群またはホモロジー群の形状を取り、
一方任意の3つの可換群で短完全列を満たすものに、H(n,X)→H(n,Y)→H(n,Z)→H(n+1,X)→という
無限に続く3×∞の形の写像図が作られる。という内容。調べ直しも加えて来週。