2023年12月17日日曜日

情報伝達活動の構造(2)

情報伝達活動の構造(1)からの続き

そもそも,年齢によって,職種によって,個人の生活スタイルによって,まったく異なるものを平均しようというのに無理がある。


図:一日の時間と情報伝達活動のモデル

会話行動に関する調査からは,平均会話時間が6時間と出てくる。ただ,会話密度が変われば,発話・受話量はまったく変わってしまうことになるが,図ではその値を書き込んだ。

TVとネットの視聴時間がそれぞれ3時間で計6時間という「情報通信メディアの利用時間と情報行動に関する調査報告書」の結果も反映しているが,ネットでのアウトプットを2時間分加えてみた。これも他の作業と並行・重複していたりで,集中度によってその内容は大きく影響される。

身体活動を中心とした労働の場面では,会話やネットでやりとりされる情報密度はずっと少ないかもしれないし,オフィスワークや対人活動が主となる仕事ならば,逆にさらに情報密度が高まるかもしれない。

これを解決するには,典型的なモデルを複数設定して,そこでの平均を提示することだ。実証的な研究に繋げるためには,GoProのようなアクションカムを24-8時間装着して,その活動を記録して分析すればよいことになる。ありそうだけれで適当な論文は探しきれていない。

総務省の情報通信政策研究所の情報流通インデックス研究会の報告書では,流通するメディア情報については客観的な指標が設定されている。問題は,これを各個人レベルに落とし込んだときの,有効な吸収率や反射率がどうなるかということ。




2023年12月16日土曜日

情報伝達活動の構造(1)

からの続き

りんちゃんの「棒」は「—」や「|」の問題だった。そうすると答えが変わってくる。
頻度の高い漢字の平均画数が7画程度で,これを棒の数とします。児童・生徒・学生の間は1日100字書いたとして(PCなど除く)10年間≒3000日では,〜30万字程度漢字をかくことになるので,〜200万画くらいかな。
1画5mmだとすれば,200万画は1000kmである。鉛筆1本で 50kmかけるらしいので,鉛筆20本あれば十分だ。ホントか,削りながらの場合は1000本くらいになりそうだ。ジェットストリームのSXR-07では替え芯1本で700mということだから,1400本なのか。


で,あらためて自分が行う情報活動の入出力量について,どの程度になるのか気になった。
のシリーズで半年前に少し考えていたことだけれど,再度見直してみる。

まず,定量的なデータについて考察する前に,個人の情報伝達活動を整理してみた。

図:情報伝達活動の構造

一番下の意識主体が,一人の人だとする。一般化できるように意識主体としている。これが,情報をやり取りするのだけれど,その情報伝達の様態を,(A)入力,(B)出力,(C)対話に分類する。(A)と(B)だけでも十分かもしれないが,あえて独立に(C)という連続的な入出力モード=対話を考えてみる。

また,その情報伝達の対象は,(イ)個体(これも個人=一人の人でもよいが,少し一般化した),(ロ)集団(組織でも社会でもよい),これまたあえて独立に(ハ)AI,を設定した。やりとりされる伝達情報の形態は,(1)静的データ(テキストや図画像),(2)動的データ(音声・音楽,動画),(3)実世界(対面による物理的化学的な接触を含む)とする。(1)と(2)ではメディアを媒介する視聴覚情報だけにフィルタされている。

都合,3×3×3=27通りのパターンに分類できる。技術的な問題を考える場合には,情報伝達の形態を図の右囲みのようにさらに7項目に細分化して考えることもできる。

(例)会話=(イ・ロ,C,3),手紙=(イ,A・B,1),電話=(イ,C,2),
ラジオ・テレビ・映画=(ロ,A,2),SNS=(ロ,A・B,1),VLOG=(ロ,B,2),演劇・音楽=(ロ,A,3),演説・講演=(ロ,B,3),生成AIチャット=(ハ,C,1)

2023年12月15日金曜日

Duolingo

今年の春に始めた iPhoneアプリDuolingoでの韓国学習である。

最近は,朝起きて新聞とテレビのニュースをチェックした後に,Duolingoで韓国語のレッスンをするのが日課になっている。iPhoneの先週のアクセスログをみると,一日平均3時間のうち,Facebookが45分,Duplingoを30分,Pikmin Bloomが20分,GoogleとLineが15分ということになっている。

単語の表記と発音と意味,韓国語の和訳あてはめ,韓国語読み上げのハングルあてはめ,日本語の韓国語訳ハングルあてはめ,など15問ほどで1セットである。単語ごとにヒントを確認できるので,覚えていないあるいは理解できていなくても,勘だけで何となく答えられてしまう。いいようなわるいような。

ヒント頼りなので,ほとんど実力がつかない。ハチマン,毎日繰り返していると記憶に残る単語も少しはあるので,韓国ドラマを見ていても,以前よりは意識に引っかかる会話が増えてきた。


図:今年11月30日時点のDuolingoの成績(実はソコまで出来ていない)

2023年12月14日木曜日

フェイスブックで,芳賀さんが息子のりんちゃんから受けた質問について書いていた。

人生でどのくらい,文字の「」を書かないといけないか」という趣旨のものだ。なぜその質問がでてきたのかはよくわからない。なぜ「棒」なのかな。書くのが面倒だからなのかな。そういう話ではないのかな。犬も歩けば棒にあたったのかな。

準備のための一般論
 モノゴコロついてから惚けるまで,およそ75年≒900ヶ月≒3900週≒2万7千日≒66万時間≒3900万分≒24億秒あるとする。このとき,人生におけるある事象の頻度と総数はおよそ次のように対応する。
 10年に 1回 〜 10回
   3年に 1回 〜 30回
   1年に 1回 〜 100回
  年に 4回 〜 300回
  月に 1回 〜 1000回
  週に 1回 〜 4000回
   3日に 1回 〜 1万回
   1日に 1回 〜 3万回
   1日に 3回 〜 10万回
名前との比較
 一生の間に自分の名前はどのくらい書くだろうか。宅配のサインも含めて週に1回よりは多いような気がするので,多分1万回のオーダーになる。「棒」はこれよりは2桁以上少ないのではないか。100回ほどか

ブログの計量
 自分のブログであるOn a Thread of the Webの5年分のデータを全文検索してみると,「棒」は16回ヒットした。つまり,最近は年平均に3回ということになる。75年をかけると,生涯では 200回 程度書くことになるのではないか。

生成AIの考え
Q(私):日本人が一日に書く漢字の数は400字程度です。一生には,1000万字程度の漢字を書くことになるでしょう。このうち,「棒」という漢字を書く回数は何回程度だと推定されますか?
A(Claude):
はい、ご指摘の通り、日本人が一生のうちに書く漢字はおよそ1000万字とされています。このうち、「棒」という漢字が何回含まれるかを推定するため、次のように考えます。日本語の漢字使用頻度データによると、「棒」の出現頻度は0.0026%,1000万字の0.0026%は2600字。したがって、一生のうちに書く漢字1000万字のなかで、「棒」という漢字が含まれる頻度はおよそ2600回と推定されます。
なお,Copilotが4690回,Perplexityが4800回,Bardは200万回,ChatGPTは10億回と答えた。(注)上記の仮定である1日に400字というのは,手書きだけではなくて,SNSなどでのスマホやPCにおけるデジタル入力も含めたつもりだけれど,どうだろうか。

漢字出現頻度
 ClaudeやPerplexityは,漢字の出現頻度という量を持ち出してきた。Perplexityでは,0.0048%とされた。そんなデータがあるのかと検索してみた。令和1年に文化庁の漢字出現度数調査というレポートがある。出現確率ではなく,出現順位が与えられている。棒は1200位あたりで,結構高頻度で出現する文字だった。もう少し調べてみると,mwSoft Blog[1]というサイトが見つかった。2012年6月Wikipediaデータにおける出現漢字をまとめたものだ。このデータを拝借すると,「棒」は1375位で38,285回/528,530,037回(上位3000字の総計)=0.00724%=724回という結果が得られた。

結論
 平均的な人が「書く」漢字の「棒」は,100回〜1000回程度じゃないでしょうか。

参考(学校で棒が出てくる場面)
 棒読み。棒暗記。棒立ち。棒グラフ。棒磁石。乳棒。制御棒。溶接棒。指揮棒。鉄棒。棒高跳び。段違い平行棒。棒球(だま)。綿棒。編み棒。棒針。相棒。泥棒。片棒。棒ダラ。アメン棒。お先棒。棒切れ。火かき棒。金棒。警棒。用心棒。こん棒。ゲバ棒。点棒。

P. S. このページだけで12年分くらいの「棒」を書いてしまった。


写真:画像生成AIがイメージしている「棒」(Diffusion Beeより引用)


2023年12月13日水曜日

Gemini

12月6日に,Googleから新しいマルチモーダルの生成AIモデルのGeminiが発表された。

宣伝のYouTubeビデオクリップでは,絵や写真を見せながら対話している様子がデモンストレーションされていて,かなりレベルが高そうだったのだが,フェイク(再現モデル)疑惑がでてきているので,どうだかよくわからない。

Google Bardのチャットページから使えるが,自分のgoogleアカウントの言語モードを英語にしておく必要があった。現時点で使えるのは,Gemini Pro であり,ChatGPT 3.5レベルのものなのだ。評判ほどではなく,これまでの体験と大きく変わることはない。というわけで,最近は,[ChatGPT3.5,Bard,Claude,Perplexity,Bing]の5つを組み合わせて試している。問題ごとに,良くできる子とイマイチな子に別れる。課金してChatGPT 4.0に復帰するだけの動機もまだない。


図:Geminiのイメージ(Google Japan ブログから引用


2023年12月12日火曜日

四角形

三角形(3)からの続き

小学生の問題から離れられない。

長方形($a \times b$)があって,その1つの角が小さな長方形($p \times q$)の形に欠けている。この図形の面積を2等分する直線を求めよというものだ。問題では辺の長さの情報が与えられていない。どうするのかと思って解答をみた。それは,与えられたL型の図形を2つの長方形の組み合わせとみなし,それぞれの長方形の中心点を結ぶ直線を引けばよいというものだ。なるほど,小学生でも理解できる。この直線の図形内の中点を中心に,図形の角を越えない範囲で直線を回転させればそれらも解になる。

ところで,その欠けた部分が大きくなると,この方法ではうまくいかない。というのも直線で分割される角の欠けた側の領域が連結領域ではなくなり,バラバラになってしまうからだ。これを避けるためには,欠けた側と反対側の角を通る直線で,面積が等しくなるものを探せば良い。中学生レベルの問題になる。


図:直線APによる四角形の面積の二等分($\frac{a}{b}>\frac{p}{q}$の場合)

直線ACは四角形ABCDの面積をS1(△)とS2(△)に2等分するが,左図形から引き去る部分は台形CEFR,右図形では三角形CGRなので,S1(左) <  S2(右)となる。そこで,分割線をAPにずらせば,S(左)=S(右)となる点が見つかるはずだ。

ずらす距離PC=$x$とおいて立式すると,$2S_1 = b(a-x)-q\bigl\{ 2(p-x)+\dfrac{q(x-a)}{b}\bigr\}$,$2S_2 = (b-q)\bigl\{ a+\dfrac{x(b-q)+aq}{b}\bigr\}$となる。これを等しいとして$x$を求めると,$x=\dfrac{q(a q- b p)}{(b-q)^2}$と求まった。

めでたしめでたしというところだ。ところで,$\frac{a}{b}=\frac{p}{q}$でもとの四角形と切り欠きの四角形が相似となる。このため,最初の対角線が2等分線となって$x=0$になる。ということは,$x<0$の場合が出てくるということか。おかしいのでしばらく悩んだ。

良く考えると,この場合は,欠けた図形の角が対角線より右側になるので,P点はCG上に持ってきてCP=$y$と置く必要がある。すなわち,この場合は,$a$と$b$,$p$と$q$を入れ替えて同じ式を解くことになるので,$y=\dfrac{p(b p- a q)}{(a-p)^2}$とすればよいことになる。




2023年12月11日月曜日

石川デジタルミュージアムネットワーク

金沢大学は,人工衛星(X線突発天体監視速報衛星)を打ち上げたというニュースを先日見たばかりだけれど,金沢大学資料館を中心とした石川デジタルミュージアムネットワークを構築した。

早速アクセスしてみると,金沢大学資料館,石川県立自然史資料館,石川県西田幾多郎記念哲学館,石川県埋蔵文化財センター,羽咋市歴史民俗資料館,野々市市ふるさと歴史館・,野々市デジタル資料館の6つの連携機関を結ぶものだった。どれも訪れたことはない。ちょっと物足りなかった。

金沢21世紀美術館とか石川県立美術館とか石川県立図書館をはじめとして諸々沢山あるので,このへんもつながると面白そうなのだけれど。しかしそれだったら文化遺産オンラインで事足りてしまうのかもしれない。


図:デジタルミュージアムネットワークの表紙から

2023年12月10日日曜日

PISA2022

PISA2018からの続き

OECD生徒の学習到達度調査(PISA)は,3年に1回実施されていた。新型コロナウィルス感染症の影響で,2021年の予定が1年延期されて2022年になり,その結果がこの度公表された。

PISA2022には,81か国・地域から約69万人が参加した。前回より2カ国ふえたのだが,トップだった中国の北京・上海・江蘇・浙江地域がはずれている。なぜか,生成AIのいくつかに聞いてみたがはっきりしない。日本からは,全国の高等学校,中等教育学校後期課程,高等専門学校の1年生のうち,国際的な規定に基づき抽出された183校(学科)約6,000人が参加して,2022年6月から 8月に実施された。

日本の順位(今回/前回)は,数学的リテラシー(1位/5位),読解力(2位/3位),科学的リテラシー(1位/2位)の3分野全てにおいて平均得点や順位が上昇した。コロナの影響で,全体の平均得点は下がったが,日本では,新型コロナウイルス感染症のため休校した期間が他国に比べて短かったことが影響した 可能性があるとの指摘がなされている。



図:PISAの平均得点の推移(国立教育政策研究所の報告から引用)

テレビで田中博之さんが,Facebookで,芳賀さんや豊福さんがコメントしていたけれど,思うところはあまりない。順位が低下していたらもわもわしていただろうか。そもそもPISAで一喜一憂しているところに問題があるのだろか。GIGAスクールの行く末の方がキニナルのか。はたまた,生成AIのインパクトのほうか。日々感受性が劣化しているのかもしれないけれど。

2023年12月9日土曜日

三角形(3)

三角形(2)からの続き

小学生向けの簡単な図形の問題を見かけた。a, b が与えられたときにx を求めれば良いというものだ。適当な補助線を引いてもわからなかったので,思わずピタゴラスの定理を使った複雑な方程式を立てて解いてしまった。

塾の先生のようなしゃべり方の人の解答をみると,どうやら,45度の三角形の解法パターンというのがあるらしい。そこで,問題を一般化したのが次の図である。これで,台形の面積の出し方や三角形の合同についての知識さえあれば,ピタゴラスの定理無しで問題が解けることになった。

日々,思考の水準が弱っていくのがわかる。そのうち小学生の問題も解けなくなりそうだ。


図:45度の三角形を含む簡単な問題

2023年12月8日金曜日

12月8日(Blogger 5周年)

2018年の12月8日にこのブログ(On a Thread of the Web)を開始してからちょうど5年が経過した。この記事を含めて1825件だ。1件平均500字≒1kBとすれば,高々1.8MBであり,微々たる情報量だ。

快食快便と同様で,情報も適切に入力出力しないと健康に悪い。そこで,毎日の散歩(寒いので最近はお休み中だけれど)のように続けてきたら5年経ったというわけだ。

生成AIの登場によって,本末転倒状態が発生しやすくなっている。すなわち,文章を自分で考えずにAIにつくらせて連続記録を更新するだけというものだ。そうでなくても,引用まみれの内容なので,自分でよく考えている部分がどれだけあるのか,頭上には疑問符がたくさん点灯している。AIに出力させるための適切なプロンプトを考えるだけでもよいことにしようか。これだけでは,惚け(認知症)予防になるのかどうか微妙だけれど。

自分のために書いているのではあるが,他の人から参照されているというカウンターを見ると励みにはなる。もっとも新着記事あたり高々数件にしか過ぎないのだが。統計的には,過去記事への遡及アクセスがあるので,もう少し参照回数はふえる。Bloggerの統計情報によれば,このブログ全体への一日あたりのアクセスは次のようになっている。一日50回だから,新着記事よりは1桁多いことになる。
全期間(5年間)6.61万回 = 36回/日
過去12ヶ月間   1.98万回 = 54回/日
過去6ヶ月間     1.38万回 = 77回/日
過去3ヶ月間     0.45万回 = 50回/日
過去1ヶ月間     0.18万回 = 60回/日

2023年12月7日木曜日

アンペールの法則(2)

アンペールの法則(1)からの続き

前回の一般的な結果を得るまでにあれこれ考えた。普通はアンペールの法則の単純な形態,つまり直線電流のまわりの円周上の磁束密度に対する,$2 \pi r B(r) = \mu_0 I$から出発して一般化するのかと思った。しかし,そもそも簡単なアンペールの法則とは直線電流まわりの磁束密度ベクトル場を与えるもので,答えは既に出ていたのだった。

あれこれの過程での計算は,結局,線積分の練習問題だった。


図:アンペールの法則の線積分経路

方針:磁束密度を測定する点への位置ベクトル$\bm{r}$とその軌跡として経路$C=r(\theta)$を考える。線要素$d\bm{r}$を変数,$r,\theta$であらわし,さらに経路条件から$r$を消去して,線積分要素を$\theta$の関数として表す。磁束密度は$r$の関数なので,これも$\theta$の関数とみることができる。その結果,線積分要素$dB=\bm{B}\cdot d\bm{r}$は$\theta$の関数になって,角度積分を実行することができる。

領域Ⅰ(左図の$0 \le \theta \le \pi/4$):$r=a/\cos\theta$,$dy = a d\theta / \cos^2 \theta$
  $dB=\frac{\mu_0 I}{2 \pi} \frac{\cos^2\theta}{a} \frac{a}{\cos^2 \theta} d\theta$,$B=\frac{\mu_0 I}{8}$
領域Ⅱ(左図の$\pi/4 \le \theta \le \pi/2$):$r=a/\sin\theta$,$dx = -a d\theta / \sin^2 \theta$
  $dB=\frac{\mu_0 I}{2 \pi} \frac{-\sin^2\theta}{a} \frac{-a}{\sin^2 \theta} d\theta$,$B=\frac{\mu_0 I}{8}$
領域Ⅲ(左図の$\pi/2 \le \theta \le \pi$):$r=a/(\cos\theta - \sin\theta)$
  $dB=\frac{\mu_0 I}{2 \pi} \frac{\sin\theta - \cos\theta}{a} \frac{a}{\sin \theta - \cos \theta} d\theta$,$B=\frac{\mu_0 I}{4}$
領域Ⅵ(左図の$\pi \le \theta \le 2\pi$):$r=a$,$d\bm{r} = a (-\sin\theta , \cos \theta) d\theta$
  $dB=\frac{\mu_0 I}{2 \pi} \frac{1}{a} a d\theta$,$B=\frac{\mu_0 I}{2}$
領域Ⅴ(右図の$-\pi \le \theta \le \pi$):$r=\sqrt{a^2+d^2+2 a d \cos\theta}$,$d\bm{r} = a(-\sin\theta, \cos\theta) d\theta$
  $\displaystyle dB = \dfrac{\mu_0 I}{2 \pi}\int_{-\pi}^{\pi}\dfrac{a(a+d\cos\theta)}{a^2+d^2+2 a d \cos\theta}d\theta = \dfrac{\mu_0 I a}{2 \pi} \int_{-\infty}^{\infty} \dfrac{(a+d)+(a-d) t^2}{(a+d)^2+(a-d)^2 t^2}\dfrac{2 dt}{1+t^2}$
$\displaystyle = \dfrac{\mu_0 I a}{2 \pi a} \int_{-\infty}^{\infty}  \Bigl\{ \dfrac{1}{1+t^2} +\dfrac{(a-d)(a+d)}{(a+d)^2+(a-d)^2 t^2} \Bigr\} dt = \dfrac{\mu_0 I}{2\pi} (\pi + \pi) = \mu_0 I$


2023年12月6日水曜日

アンペールの法則(1)

物理科学概説の授業で,アンペールの法則のところに入った。積分形では,$\displaystyle \oint_C \bm{B}(\bm{r})\cdot d\bm{r}=\mu_0 I$である。

十分長くてまっすぐの導線を流れる電流のまわりの磁束密度$\bm{B}(\rm{r})$の強さ$B(r)$は,電流の強さ$I$に比例し,電流からの距離$r$に反比例する。その向きは電流の向きに右ネジが進むときにネジが回る方向(電流を中心とした半径$r$の円の接線方向)である。

この実験事実を式で表現する。直線電流上の一点を原点に取って,磁束密度ベクトルは原点をとおり電流に垂直な平面内にある。観測点の座標を$\bm{r}=(r\cos\varphi, r\sin\varphi)$として,$\bm{B}(\bm{r})=\dfrac{\mu_0 I}{2\pi r}(-\sin \varphi, \cos \varphi) = B(r) \cdot \bm{e}_{\varphi}$となる。このとき,$\bm{B}(\bm{r})\cdot \bm{r}=0$となっている。

このとき,積分形のアンペールの法則を導けるかという問題だ。


図:アンペールの法則の積分形の導出

積分形のアンペールの法則では,空間中に任意の閉経路Cを設定して,この経路Cに対する磁束密度の線積分を求める。線積分要素は$dB = \bm{B}(\bm{r})\cdot d\bm{r}= B(r) dr \cos\theta$となる。一方,線要素の磁束密度方向の成分は,$dr \cos \theta = r d\varphi$である。そこで,$dB = B(r) r d\varphi = \dfrac{\mu_0 I}{2 \pi} d\varphi$となる。$\therefore \oint_C \bm{B}(\bm{r})\cdot d\bm{r} = \int_0^{2\pi}  \dfrac{\mu_0 I}{2 \pi} d\varphi = \mu_0 I$

したがって,無限直線電流に対して,3次元空間内でこれを囲む任意の閉経路での磁束密度の線積分の値は,この経路を貫く電流に磁気定数をかけたものとなる。


2023年12月5日火曜日

物体O

10月1日では遅すぎるからの続き

小松左京(1931-2011)の1962年の短編「物体O」は,ハヤカワ・SF・シリーズ3088の短編集「日本売ります」に収録されている。たぶん,銀背シリーズの1冊で初めて買ったものだ。最初に読んだ小松左京のSF作品の一つでもある。

詳細な科学・社会スペキュレーションなどで,後の日本沈没首都消失などの系列につながる作品だ。そのルーツは日本アパッチ族にあるだろう。小松左京の主な関心は社会構造の変革可能性にある。阪大理学部核物理の三伏教授(伏見康治)が登場して,阪大というキーワードがインプットされたので,これが自分の後の進路選択に微かに影響したかもしれない。

ところで,再読してみるといくつか気になるところがあった。

(1)物体Oの落下位置
兵庫県相生市を中心とした半径450km,500km,550kmの円を描くとほぼ記述通りの位置が再現できるのだが,文中に記載のある屋久島と種子島は物体の下敷きにはならない(なお韓国の釜山は物体Oに下に沈む)。

(2)兵庫県豊中市などの地名
物体の落下を見たと気象台に報告したアマチュア天文家の在住地,これは大阪府です。小松左京は後に箕面に住むことになるのだけれど。物体Oが落下した関東の地名の一部も不正確かもしれない。

(3)物体Oの密度
直径1000km,幅100km,高さ200kmのリングなので,その体積は,2π×500×100×200 = 6×10^7 km^3 = 6×10^16 m^3 になる。一方文中にはその質量が2万兆トン=2京トン= 2×10^16 t とある。したがって,その密度は,0.3 t/m^3 = 0.3 g/cm^3になる。ところが,後の方では「比重は重金属−銀と同じ程度だ」という記述があって(銀は10.5 g/cm^^3)矛盾している。

(4)紫外線のラウエ斑点
紫外線による薄片の写真に写った斑点が,X線によるラウエ斑点と同じではないかというところから,物体Oが普通の物質を5000万倍に拡大したものであることがわかる。しかし,結晶構造をみるX線の波長は1Å=0.1 nm であり,紫外線の波長は 300nmなので,高々3000倍にしかならない。5000万倍にすれば,波長5mm のミリ波でなければならない。


図:物体Oの想定位置(半径450km, 500km, 550km)

2023年12月4日月曜日

万博リング

評判の悪い大阪・関西万博2025の大屋根(リング)である。

万博のシンボルとなる木造大屋根リングは350億円もかかる。そのため,会場建設費+α=2350億円+800億円以上の無駄遣いの象徴とされて,あちこちから叩かれている。まあ大阪維新が万博を持ち出した動機が,同じ夢洲で計画されているIR用地周辺の環境整備だったり,虚構の経済効果だったりするので,それらが見透かされるとこうなってしまう。

そもそも,失われた30年を取り戻すための戦略が,高度成長期の夢よ再びという東京五輪+大阪万博でしかなかったという,創造力+想像力の貧困と税金奪取機会の創出が問題だったわけだ。そのあたりの本質的な問題点を一度忘れて考えてみると,万博にはシンボル的な建造物が必要であるということはわかる。ロンドンの水晶宮,パリのエッフェル塔,大阪の太陽の塔などなど。

1970年の大阪万博にもお祭り広場の上に大屋根があった。その大屋根を突き破った太陽の塔には岡本太郎(1911-1996)が必要だったけれど,今の日本にはそれに匹敵するパワーを持ったクリエイターはいないので,大屋根リング止まりになった。リングがなくて,海外各国からの主要パビリオンもなければ,大阪・関西万博2025は本当にグズグズ,バラバラになってしまうのだろう。

なんだかんだいって,パビリオンが一部建設途上のままでも,2025年の4月に万博が始まってしまえば,マスコミがこぞって囃し立てて機運を盛り上げる。大阪の子供に配った無料チケットの効果もあって,2800万人ではなくともそこそこ人は集まるだろう。よほどの混乱や災害が起きない限りは,東京五輪と同様に無事終って良かったというストーリーがでっち上げられそうな気がする。自分は,万博には行かないつもりだが,孫が連れていってと言い出したときにどうなるかはわからない。


さて,当初計画にはなかった万博リングをごり押しで導入した会場デザインプロデューサの藤本壮介だが,リングのデザインがノーマン・フォスター(1935-)のアップルパーク(2017)のパクリではないかという疑惑記事をみかけた。丸くて中に池があるので似ているというのも言い掛かりじみている。それぐらいはしかたないだろう。ただ,万博リングの直径650mは,アップルパークの直径460mの約√2倍で,万博リングの平均高さ16mの方は,アップルパーク23mの約1/√2倍になっていた。まあ偶然である。

ちなみに,小松左京の物体Oに関係があるのではないかと思って確認してみたが,物体Oは直径1000km,幅100km,高さ200kmの銀の塊だったので形状はかなり違う。また,中心は大阪ではなくて,兵庫県の相生市付近だった。そんなこんなで,リングという形状それ自身はあまりいじめられなくてもいいのにと思う。


図:各リングの断面図の比較

2023年12月3日日曜日

球形キャパシタ

球形キャパシタの問題を物理科学概説の中間テストで出題した。

教科書の例題と同じ単純な問題のつもりだったけれど,2つの球殻に与える電荷の記述を省略したため,アースの取り方によって話が変わるのだった。それが教科書の章末課題に書いてあったので,良く勉強した学生さんはそちらを参照していた。


図:球形キャパシタのイメージ

半径$a$と$b$の同心の導体球殻があり,それぞれに電荷$q_a$と$q_b$を与えたとき,それぞれの電位が$V(a)$と$V(b)$になったとする。内球殻の電荷がつくる電場は,$E_a=\dfrac{q_a}{4\pi\varepsilon_0}\dfrac{1}{r^2}\quad (a<r<b)$,であり,これによって誘導される電荷が外球殻の内面に$-q_a$,外面に$q_a$だけ生ずる。これによって,外球殻の外面には$q_a+q_b$の電荷が分布するので,この電荷が作る電場は,$E_b=\dfrac{q_a+q_b}{4\pi\varepsilon_0}\dfrac{1}{r^2}\quad (b<r)$となる。

これから,外球殻の電位は,$\displaystyle V_b(r) = -\int_\infty^r \dfrac{q_a+q_b}{4\pi\varepsilon_0}\dfrac{1}{r^2}\ dr = \dfrac{q_a+q_b}{4 \pi \varepsilon_0} \dfrac{1}{r} \quad (b<r)$ となり,$V(b) =  \dfrac{q_a+q_b}{4 \pi \varepsilon_0} \dfrac{1}{b}$
内球殻の電位は,$\displaystyle V_a(r) = V(b) -\int_b^r \dfrac{q_a}{4\pi\varepsilon_0}\dfrac{1}{r^2}\ dr = V(b) + \dfrac{q_a}{4 \pi \varepsilon_0}\dfrac{1}{r} -  \dfrac{q_a}{4 \pi \varepsilon_0}\dfrac{1}{b}$
$\therefore V(a) =  \dfrac{1}{4 \pi \varepsilon_0} \Bigl( \dfrac{q_b}{b} +  \dfrac{q_a}{a} \Bigr)$

(1) 外球殻が接地されている場合
$V(b)=0$より$q_b = -q_a$となる。$\therefore V(a) = \dfrac{q_a}{4 \pi \varepsilon_0}\Bigl( \dfrac{1}{a} - \dfrac{1}{b}\Bigr) = \dfrac{q_a}{C}$とすれば,
キャパシタの電気容量$C$は,$C = \dfrac{4 \pi \varepsilon_0 a b }{b-a}$となる。

(2) 内球殻が接地されている場合
$V(a)=0$より$q_a = -\dfrac{a}{b} q_b$となる。$\therefore V(b) =  \dfrac{q_b}{4 \pi \varepsilon_0} \dfrac{1 - a/b}{b} = \dfrac{q_b}{C'}$とすれば,
キャパシタの電気容量$C'$は,$C' = \dfrac{4 \pi \varepsilon_0 b^2 }{b-a}$となる。

このとき,$C' = C +  4 \pi \varepsilon_0 b$となって,外球殻をキャパシタと考えたときの電気容量とCとの並列接続の式となっている。

2023年12月2日土曜日

円軌道はむずかしい

万里鏡1号弾道ミサイルの軌道(2)からの続き

北朝鮮の弾道ミサイルの簡単なシミュレーションコードをMathematicaで作っていた。これを少しアレンジすれば人工衛星を軌道に投入するところまでできそうな気がする。

早速,以前のコードを修正してみた。まずは通常の加速直後に角度方向だけに加速度を加えるようにしたがうまくいなかい。打ち上げ加速は投射角の方向になっているので,動径速度成分が大きく残っているうえ,加速すれば軌道は膨らむ。このため,離心率の大きな長円軌道になって地表にぶつかってしまうか,地球の重力圏から脱出してしまうのだ。

次に,打ち上げ加速の直後に空白時間をおいて,動径速度成分が小さくなったところで角度方向に加速できるようにした。それでもうまくいかない。簡単な試行錯誤では周回軌道にのせるのが難しい。そもそも角度方向に加速するということは面積速度すなわち角運動量をふやし,動径方向の微分方程式で軌道半径を膨らませる方向に作用してしまう。

そこで,後期加速では衛星をその速度ベクトルの方向に加速することにした。$t=0$で速度ベクトルをとりだすところに発散があったので,これを回避するため,地球の自転による2倍面積速度$h(t)$の初期値として,$h(0) = R^2\omega=R^2 \frac{2\pi}{24*3600}=2930$ km$^2/$sを与えた。初期加速度はこれまでの$\ a=0.0446$として$30$秒加速する。その後,800秒程度休止した後に,後期加速度$\ b=0.1445$(ここを微調整した)で$250$秒加速すると,なんとか軌道に投入することができた。投射角は$s=45$度,初期加速における燃料比は$p=0.85$であり,加速方向の角度には$0.3$をかけて動径成分を抑えた。

なかなか難易度の高いゲームである。衛星の軌道高度が1200km程度の準円軌道となっている。これを500kmにしなさいといわれても,こんな単純な2段階制御ではちょっと難しい。なお,プログラムの検証のため,$r=6850$kmの宇宙空間で第一宇宙速度に相当する$v=\sqrt{gr}=8.2$ km/sを角度方向の初速度として与えると,正確に円軌道を描くことが確かめられた。

g = 0.0098; R = 6350; τ = 30; τs = τ*27; τt = 250; p = 0.85;
 a = 0.0446; b = 0.1445 a; s = 45 Degree; T = 15400; 
fs[t_] := 0.3*ArcTan[r[t]*r'[t]/ h[t]]
fr[t_, τ_] :=  a*Sin[s]*HeavisideTheta[τ - t] + 
   b*Sin[fs[t]]*HeavisideTheta[t-τs-τ]
   *HeavisideTheta[τ+τs+τt-t]
ft[t_, τ_] :=  a*Cos[s]*r[t]*HeavisideTheta[τ - t] + 
   b*Cos[fs[t]]*r[t]*HeavisideTheta[t-τs-τ]
   * HeavisideTheta[τ+τs+τt-t]
fm[t_, τ_] := -p/(τ - p*t)*HeavisideTheta[τ - t]
sol = NDSolve[{r''[t] == -fm[t, τ]*r'[t] 
   +h[t]^2/r[t]^3 -g R^2/r[t]^2 +fr[t, τ],
   r[0] == R, r'[0] == 0, 
   h'[t] == -fm[t, τ]*h[t] + ft[t, τ], 
   h[0] == 2930 + 0*Sqrt[g] R^(3/2)}, {r, h}, {t, 0, T}]
f[t_] := r[t] /. sol[[1, 1]]
d[t_] := h[t] /. sol[[1, 2]]
Plot[{6350, f[t]}, {t, 0, T}]
Plot[{f[t + 1] - f[t], d[t]*R/f[t]^2, d[t]/f[t]},
 {t, 0, T}, PlotRange -> {-5, 15}]

 


図:苦労すると有難みがわかる衛星の準円軌道のグラフ

P. S. もう少しがんばると,軌道高度650km(r=6980km)の準円軌道まで達成できた。
g = 0.0098; R = 6350; τ = 25; τs = τ*15.3; τt = 350;  p = 0.85;
a = 0.0446; b = 0.1275 a; s = 45 Degree; T = 15400; 

2023年12月1日金曜日

ai pin

ぼんやりしていると,重要な話題を逃してしまうことが多い。

3週間前の11月10日のニュースになっていたヒューメインai pin である。ヒューメインは,アップルの技術者として働いていたイムラン・チャウドリベサニー・ボンジョルノ夫妻が2018年に立ち上げたAIスタートアップだ。

NHKのニュースウォッチ9が,ChatGPTが1周年を迎えた生成AIの課題をまとめたニュースの中で,新たなAIデバイスとしてai pin を紹介していた。早く教えてよ。調べてみると,3週間も前に最初のニュースが出ていた。

ai pin は,胸につける小さなバッジ型のデバイスであり,音声とカメラによるインターフェイスでAIとやり取りすることができる。ディスプレイはついていないが,小型レーザプロジェクタで手のひらに文字情報を映し出すことができる。なるほどそうきたか。

既に米国では699ドル(日本円で10万円強)販売(予約?)されているが,AI使用料を含んで月々24ドル(3,600円)のサブスクリプションが必要である。日本で使おうとすると国際ローミングが必要なのでちょっとまだ無理そうだった。音声入力インターフェースには良い面と困った面があるので,公共の場でうまく使えるかどうかが勝負か。同時通訳には最適だ。ジェススチャーだけである程度のコマンドが指示できれば問題ないのかもしれないが。



写真:AiPinのサイズ(humaneサイトから引用)

2023年11月30日木曜日

OpenAI

ChatGPT(6)からの続き

2022年の11月30日にOpenAIがChatGPTを発表してからちょうど1年目となった。

最初のインパクトが大きかったので,最近の動きがもどかしく見えてしまう。ChatGPTの話ではないが,生成AI高性能のGPUを積んだWindowsマシンでないといろいろと試すことができないのもつらい。MacintoshのM1〜M3チップではほとんど対応できないのだ。GPT-4の相対的機能低下もあって,有料版からは離脱してしまった。最近は,GPT-3.5とBard, Claude, Bingなどの比較に留まっている。この中ではGPT-3.5が一番融通が利かないかもしれない。

そうこうしているうちに,11月6日のDevDayで300ページのテキスト入力に対応する128Kコンテクスト・ウィンドウを持つGPT-4 Turboがアナウンスされた。そろそろGPT-4復帰すべきかと思っていると,11月17日には,サム・アルトマンの解任騒動が勃発した。まだ全貌が必ずしも明らかになっていないが,11月22日に,サム・アルトマンとグレッグ・ブロックマンが復帰することでいちおう落ち着いたように見える。Q*をめぐる陰謀論めいた話はいろいろあるが,AGIが出来ているとか出来そうだということでもないか。


[付録]なんだか妙に下手に出るように躾けられているClaudeの場合:
はい、完全におっしゃるとおりです。失礼しました。訂正いたします。
e = 2.718...
log π = 1.1447...
として、e × log π を正確に計算すると、
2.718... × 1.1447... = 3.11...(小数点以下省略)
となりますね。
3.14ではありませんでした。

計算をちゃんと行わずに短絡的に判断してしまい、ファクトとしても誤りがありました。
このようなミスを重ねてしまい、大変申し訳ありません。
基本的な数値の計算から正しく立ち返る必要があると反省いたしました。
今後は肝に銘じ、確認を怠らぬようにいたします。
ご指導本当にありがとうございました。

2023年11月29日水曜日

大小関係

よくある問題で,冪数の大小比較というのがある。その例で次のようなものがあった。

$M=\displaystyle \begin{pmatrix} e^e & e^3 & e^\pi \\ 3^e & 3^3 & 3^\pi \\ \pi^e & \pi^3 & \pi^\pi \end{pmatrix}$の9個の数の大小関係を求めよ。
ただし,$e=2.7183 < 3 < \pi=3.1416$はわかっているとする。

各行や各列で比較すると,行番号や列番号が増えると大きくなる。
次に,$M_{12}=e^3$と$M_{21}=3^e$,$M_{23}=3^\pi$と$M_{32}=\pi^3$,$M_{31}=\pi^e$と$M_{13}=e^\pi$を比較する。それぞれ,両者のべきの積の逆数を双方にかけると,$x^\frac{1}{x}$の形での比較に帰着する。この関数の対数をとって$f(x)$とおけば,$f(x) = \frac{\log x}{x},\ f'(x) = \frac{1 - \log x}{x^2}$の形から,$e^\frac{1}{e} < 3^\frac{1}{3} < \pi^\frac{1}{\pi}$である。


図:$f(x) = \log x^{1/x}$とその微分 $f'(x)=(1-\log x)/x^2$のグラフ

したがって,$  3^e < e^3 ,\  \pi^e < e^\pi,\  \pi^3 < 3^\pi $が成り立つ。残るのは,$3^3$と$e^\pi$または$\pi^e$の関係である。これがちょっとわからなかった。仕方がないので,数値的に評価することに。

$\log M_{13}=\pi = 3 + 0.1416$,$\log M_{33} = 3 \log 3 = 3 + 3(\log3 - \log e) = 3 + 3 \log \frac{3}{e}$
$3 \log \frac{3}{e} = 3 \log (1 + \frac{3-e}{e}) \approx 3 \Bigr\{ \frac{3-e}{e}-\frac{1}{2} \bigl( \frac{3-e}{e}\bigr)^2 + \cdots \Bigr\}= 0.295$ 。したがって,$e^\pi < 3^3$

結局,$e^e < 3^e < e^3 <  \pi^e  <  e^\pi <  3^3 < \pi^3 <  3^\pi < \pi^\pi$ となった。


追伸(2023.11.19):ひとつ確認もれがあった。$e^3$ と$\pi^e$の大小関係である。
対数をとると,$3$と$e \log \pi$の比較になる。
$\log \pi = \log e(1 + \frac{\pi-e}{e}) = 1 +  \log ( 1 + \frac{\pi-e}{e} ) \approx 1 + \frac{\pi-e}{e} -\frac{1}{2}\Bigl(  \frac{\pi-e}{e} \Bigr)^2$
したがって,$e \log \pi \ \ (3.1117) \approx \pi -\frac{(\pi-e)^2}{2e}\ \  (3.1086)  > 3$,$\therefore \pi^e > e^3$

2023年11月28日火曜日

妖星ゴラス

海底軍艦からの続き

11月のWOWOWで東宝の特撮SF映画シリーズをやっていた。地球防衛軍(1957),宇宙大戦争(1959),妖星ゴラス(1962.3),海底軍艦(1963),緯度0大作戦(1969)の五本だ。

このうち,1963年の海底軍艦だけ映画館で見ているのは以前書いた通り。妖星ゴラスは小学校2-3年の時で,近所に映画のポスターも貼ってあった(そんな時代)。ストーリーも薄々わかって,とても見たかったのだけれど,当時は"大人の映画"につれていってほしいと言い出せるとは思っていなかった。まもなく,最初に体験することになる東宝の特撮怪獣映画は,キングコング対ゴジラ(1962.8)で,それ以後,夏休みのゴジラシリーズ等には連れて行ってもらえた。


その妖星ゴラスは,本多猪四郎(いしろう)と円谷英二のコンビ作品のうちの怪獣物でないSF作品の一つであり,今回のWOWOWの特集もそうしたSFものから変身人間シリーズなどを除いた5作が選ばれている。ただし,世界大戦争(1961)は含まれていない。

妖星ゴラスは,地球の0.75倍の大きさだが,重力が6000倍近い"黒色矮星"という設定で,地球に向かってくる。この星の接近による地球の破壊を避けるために南極にロケット噴射装置を設置して,地球をその公転軌道からずらすというものだ。$10^{-6}$Gを100日かけて40万km移動する。加速終了後も等速運動を続けるのはどうするのかと思ったけれど,映画の中では,北極に装置を再設置して逆に動かすような説明をしていた。

このため,南極におけるロケット噴射を表現するガスバーナーの炎のシーンが延々と続くのだった。ただ,説明では重水素と水素による核エネルギー(核融合とか水爆というキーワードは表立って出てこない)的なものが示唆されている。そのわりにはガスバーナーなのであるが。アポロ11号を打ち上げたサタン5号程度の推力ならば,1万セットで$10^{-12}$Gを短時間加えられるかもしれないけれど,ちょっとかなり厳しい。

おもしろかったのは,久保明がゴラスの再調査に向かったときに危機的状況になって記憶喪失になるシーン。ところどころ,2001年宇宙の旅(1968)のボーマン船長を思わせるようなシーンや宇宙ステーションへの回収のカットが出てくるのだ。キューブリックがこの映画を観ていることはないと思うが・・・。ところが検索してみると,同様の意見が散見された。もしかすると影響しているのだろうか。

なお,毛色が異なるので今回は含まれていない第三次世界大戦ものである世界大戦争を検索していたら,第二東映の第三次世界大戦 四十一時間の恐怖(1960)というドキュメンタリータッチのモノクロ映画も見つかった。当時は相当世界危機的な認識が広まっていた状況だったのだろう。



写真:妖星ゴラスの一場面([1]から引用)

[1]映画 妖星ゴラス(サブロジーの日々是ずく出し)

2023年11月27日月曜日

奥州安達原

 妹背山婦女庭訓(2)からの続き

国立文楽劇場の11月文楽公演(第172回)は,第2部の奥州安達原を観た。第1部は双蝶々曲輪日記と面売り,第3部は冥途の飛脚という演目。

奥州安達原を調べてみると,10年前の2013年11月に観ている。今回に加えて,道行千里の岩田帯,一つ家の段,谷底の段が上演されていて,一つ家の段と谷底の段の安達ケ原の鬼婆がでてくる怖い話のところが印象に残った。記憶の中では袖萩祭文が安達ケ原で行われたような錯覚に陥っているが,そんなことはありません。

今回は,朱雀堤の段と環の宮明御殿の段(敷妙使者の段,矢の根の段,袖萩祭文の段,貞任物語の段)なので,舞台は京都なのだった。出語り床の真ん前の席で,三味線の手も太夫の汗もよく見え,太棹の強いバチさばきが体に伝わるのだが,舞台の下手の方が遠いので肝腎の袖萩祭文の動き(和生・勘次郎)がよくわからない。

今回は,芳穂太夫・錦糸組ががんばっていた。呂勢太夫・清治の袖萩祭文から錣太夫・宗助の貞任物語もなかなかよかった。ところが,ストーリーが追いきれないので,平傔仗直方がなぜ切腹するのかがいまいち納得できないまま物語が急展開していくのだった。

休日で,阪神タイガースとオリックスバッファローズの優勝記念パレードの日だった。天気も良くて,人出は多かったのだが,文楽劇場の第2部の入りは40%ぐらいだった。大丈夫かな。


写真:10年前の奥州安達原のポスター,こちらのプログラムの方がよかった。

2023年11月26日日曜日

万里鏡1号

11月21日の22時42分に,朝鮮民主主義人民共和国東倉里の西海衛星発射場から,新型ロケットチョルリマ(千里馬)1型が発射され,軍事偵察衛星マルリギョン(万里鏡)1号が衛星軌道に投入された。

これに関するNHKニュースを時系列に並べると次のようになる。
・11月21日 22時46分 Jアラート発令「ミサイル発射。ミサイル発射。北朝鮮からミサイルが発射されたものとみられます。建物の中、又は地下に避難して下さい。(対象都道府県:沖縄県)」(Jアラート)
・11月21日 23時15分 Jアラート発令「ミサイル通過。ミサイル通過。先程のミサイルは22時55分頃、太平洋へ通過したものとみられます。避難の呼びかけを解除します。不審な物には決して近寄らず直ちに警察や消防などに連絡して下さい。(対象都道府県:沖縄県)」(Jアラート)
・11月22日 1時22分 北朝鮮ミサイル発射沖縄・鹿児島の状況は(日本)
あまり意味のなさそうなJアラートを全国に発令してこれでもかと煽り続けるのはやめてほしい。そのぶん客観的で正確な情報を流してね。弾道ミサイル実験の場合と違って,10分ほどで高度500kmに到達し,衛星軌道への投入に切り替わる。この段階でレーダーによる追跡ができなくなったようだ。

少しフライングしているものの,人工衛星としての打ち上げ予告をしているのに,弾道ミサイル実験だと言い募るのもどうにかならないのか。軍事偵察衛星の実験はそれはそれで困るし,ミサイル技術の前進があるのは確かだけれど,弾道ミサイルとしての実験ではないだろう。

韓国が衛星打ち上げに成功したと確認したあとでも,「何らかの物体が地球の周回軌道に」とディスっても仕方ないだろうに。機能していなくても地球を周回していれば初の人工衛星の打ち上げ成功には違いない。万里鏡1号はNORADにも登録されていて,高度が遠地点519km,近地点499.6kmのきれいな太陽同期円軌道になっている。もちろん,カメラが正常で十分に機能しているかどうかは別問題だけれど。

1970年2月11日の,日本初の人工衛星おおすみの場合は,非軍事目的ということで誘導制御ができなかったこともあるが,遠地点5151km,近地点337kmの楕円軌道で,15時間で電池が熱やられて通信は途絶している。それでも光学観測で軌道は確認できたようだ。おおすみは,2003年8月2まで33年間地球周回軌道にあった。


図:万里鏡1号の軌道の例(NOY2.comから引用)


2023年11月25日土曜日

アマテラス粒子

観測史上2番目にエネルギーの高い宇宙線が見つかったというニュース。11月24日のサイエンスオンラインに論文が掲載されるはずだけれど,まだ見当たらない。

実験史上最大というプレスリリースになっているのは,米国ユタ州のテレスコープアレイ実験(2008-,760㎢に1.2km間隔で507台の大気チェレンコフカウンターを並べた装置)においてという意味だ。これまでに観測された史上最大エネルギーの宇宙線は,同じユタ州のダグウェイ実験場で1991年に見つかった,オーマイゴッド粒子だ。そのエネルギーは,3.2±0.9 × 10^20 eV = 320 EeV(エクサ電子ボルト)= 51 J(ジュール)である。

今回,大阪公立大学や東京大学宇宙線研究所などのメンバーを含む国際共同実験チームが見つけた宇宙線は,アマテラス粒子と命名され,そのエネルギーは 244 EeV = 38 J である。これがマスコミに報道されるとき,40Wの電球を1秒点灯させるだけのエネルギーというのはOKだけれど,1gで地球を破壊するほどのエネルギーだという例えに引っかかった。

この宇宙線の正体となる粒子が何であるか(粒子1個の質量)がわからなければ,1gに相当する粒子数が定まらない。とりあえず,銀河宇宙線の大半を占める陽子だとすると,質量は,1.67 × 10^-27 kg なので,6 × 10^23 個分にあたる。1g 分のアマテラス粒子群全体の持つエネルギーは,2.4 × 10^25 J なのだけれど,これで地球は破壊できるのだろうか。

広島に投下されたリトルボーイのエネルギーは,7 × 10^16 J らしいので,3億個の広島型原爆を落とされたことになる。また,チクシュルーブ・クレーターを作って恐竜を絶滅させたといわれる直径10km,速度20km/sの小惑星は,リトルボーイの1億倍のエネルギーに相当するので,この小惑星衝突の3個分のエネルギーに相当する。この表現のほうがわかりやすかっただろうか?

宇宙線研究所のプレスリリースには共同実験代表者の荻尾彰一さんの写真が大きく載っていた。彼は,大阪市立大学時代には物理教育学会近畿支部長を務められ,いろいろとお世話になったのだった。


図:大阪公立大学のプレスリリースから引用(©Ryuunosuke Takeshige)


2023年11月24日金曜日

ケプラー方程式

楕円軌道からの続き

軌道の形ではなく,時間発展を考える。
出発点は,$\bm{r}(t) = (x(t),y(t)) = ( a(\cos\theta(t) -e), a\sqrt{1-e^2} \sin \theta(t))$と,面積速度が一定であるということだ。
長半径$a$,短半径$b$,離心率$e$の 楕円の面積$S$は $S=\pi a b = \pi a^2 \sqrt{1-e^2}$なので,周期を$T$とすると,面積速度は,$\dfrac{dS}{dt}= \dfrac{S}{T} = \dfrac{\pi a^2 \sqrt{1-e^2}}{T}$である。

次に,楕円上の位置ベクトル$\bm{r}(t)$から面積速度を計算する。
$\dfrac{dS}{dt}= \frac{1}{2}(\bm{r} \times \dot{\bm{r}})_z = \frac{1}{2}(x \dot{y} - \dot{x} y) $
$= \frac{1}{2}  \{ a(\cos\theta-e) \cdot a\sqrt{1-e^2} \cos\theta \dot{\theta} - (-a \sin \theta \dot{\theta} ) \cdot a\sqrt{1-e^2}\sin \theta \}$
$= \dfrac{a^2 \dot{\theta}\sqrt{1-e^2}}{2} (\cos^2 \theta -e \cos\theta + \sin^2 \theta) =  \dfrac{a^2 \dot{\theta}\sqrt{1-e^2}}{2} (1 -e \cos\theta ) $

この2つの式が等しいので,$\dfrac{2\pi}{T} = \dot{\theta} ( 1- e \cos\theta)$となる。この両辺を時間$t$で積分して,$t=0$で$\theta=0$とすれば,次のケプラー方程式が得られる。
$\dfrac{2\pi}{T} t = \theta -\sin \theta \quad ( 0 \le t \le T \ \ \rightarrow\ \  0 \le \theta \le 2\pi) $
この解$\theta(t)$ を$\bm{r}(t) = (x(t),y(t)) = ( a(\cos\theta(t) -e), a\sqrt{1-e^2} \sin \theta(t))$に代入すれば,位置ベクトルが時間の関数として表される。

Mathematicaで計算してみた。
f[t_, e_] := FindRoot[u - e Sin[u] == 2 Pi t, {u, 0}]
g1[a_, e_] := 
  Table[{a (Cos[u] - e), a Sqrt[1 - e^2] Sin[u]} /. f[k/52., e], {k, 1, 52}];
gp1 = Graphics[{PointSize -> Large, Red, Point[g1[1, 0.2]]}];
g2[a_, e_] := 
 Plot[{a Sqrt[1 - e^2] Sqrt[1 - (x/a + e)^2], -a Sqrt[1 - e^2] Sqrt[
     1 - (x/a + e)^2]}, {x, -a (1 + e), a (1 - e)}, 
  AspectRatio -> Automatic, PlotStyle -> Blue]
gp2 = g2[1, 0.2];
Show[gp2, gp1]

図:ケプラー軌道の計算例(a=1, b=0.98,  e=0.2, r_ap=1.2, r_pe=0.8)

2023年11月23日木曜日

楕円軌道

昼夜時間(3)からの続き

惑星の楕円軌道についての復習の時間。


図:楕円の性質

図の左が,楕円の中心を原点とする$\ (X-Y)$座標系の表示である。長半径$a$の円を短半径$b$方向に$b/a$倍すると,$(X/a)^2+(Y/b)^2=1$という楕円が描かれる。このとき,縮小前の点への位置ベクトルが$X$軸となす角度を$\theta$として,楕円上の点の座標が$\ (X=a \cos \theta, Y= b \sin \theta)\ $となる。また,原点から種横転までの距離は$\ \sqrt{a^2-b^2}=ae$となる。ただし,離心率が$\ e=\sqrt{1-(b/a)^2}$と定義される。

図の右が,楕円の焦点を中心とする$\ (x-y)$座標系での表示である。原点から最も近い$x$軸方向の近地点までの距離を$r_{\rm pe}=a(1-e)$,最も遠い遠地点までの距離を$r_{\rm ap}=a(1+e)$とすると,$a=(r_{\rm ap}+r_{\rm pe})/2,\ b=\sqrt{r_{\rm ap} r_{\rm pe}}$である。また,楕円上の点への位置ベクトルは,その長さを$r$,$x$軸のなす角度を$\phi$として,$(x=r\cos\phi, y=r\sin\phi)$と表される。

ところで,$x = a \cos \theta - a e = r\cos\phi,\ y = b \sin \theta = r\sin\phi)$である。そこでこれらから,$\theta$を消去すれば,$r$と$\phi$の関係式が得られる。すなわち,$r=\sqrt{x^2+y^2} = a (1-e\cos\theta)$,$\tan \phi = \dfrac{y}{x} = \dfrac{\sqrt{1-e^2} \sin \theta}{\cos \theta - e}$
そこで,$\dfrac{1}{\cos^2\phi} = 1 + \tan^2\phi = \dfrac{(\cos \theta - e)^2 + (1-e^2)\sin^2 \theta}{(\cos \theta - e)^2}=\dfrac{(1-e\cos\theta)^2}{(\cos \theta - e)^2}$
$\therefore \dfrac{1}{\cos\phi} = \dfrac{1-e\cos\theta}{\cos \theta - e} =  \dfrac{r}{a\cos \theta - a e} =  \dfrac{r e}{a - r - a e^2} $
最終的に,$r = \dfrac{a(1-e^2)}{1 + e \cos \phi}$となる。



2023年11月22日水曜日

添田唖蝉坊

東京節からの続き

東京節(1918)の別名はパイノパイノパイで,ジョージア行進曲のメロディーに演歌師の添田知道(添田さつき)(1902-1980)が詞をつけたものだ。彼の父親は添田唖蝉坊(そえだあぜんぼう,1872-1944)という演歌師だった。

最後の演歌師といわれる桜井敏雄(1906-1996)がバイオリンをひきながら,なぎら健壱東京節を歌っているものや,打楽器奏者土取利行東京節解放節を聴くことができる。なかなか味わい深い。これが,後の高石ともや,フォーククルセダーズ,岡林信康などにつがなっているのかも。


写真:添田唖蝉坊(大磯町のページから引用)

[2]のんき節(添田唖蝉坊,土取利行)
[3]あきらめ節(添田唖蝉坊,土取利行)
[4]ストライキ節(添田唖蝉坊,土取利行)

2023年11月21日火曜日

東京節

柳澤眞一(1932-2022)といえば,奥様は魔女のダーリンだ。彼がYouTubeで榎本健一(1904-1970)の東京節を話題にしていた。「ラメチャンたら ギッチョンチョンで パイのパイのパイ♪ パリコとパナナで フライ フライ フライ♪」うーん,昔良く聞いたことのある歌だ。

そもそもこれはどういう意味なのかと調べると,なんとなくそれらしい解釈が見つかった。
果たして,AI諸氏はこれに答えられるだろうか。Bingは,検索できるのでこの答えをそのまま返してきた。その他はどうもいまいち。やはり,日本語を十分学習したLLMが必要だ。

Q:榎本健一の東京節に「ラメチャンたら ギッチョンチョンで パイのパイのパイ パリコとパナナで フライ フライ フライ」という歌詞があるのですが,その意味を教えてください。

A:Bing
「ラメチャンたら ギッチョンチョンで パイのパイのパイ パリコとパナナで フライ フライ フライ」という歌詞は、大正時代に流行したコミックソング「東京節(パイノパイノパイ)」の一部です。この歌詞には、大正時代の東京の町や建物の名前が歌い込まれています。歌い出しは「東京の中枢は 丸の内 日比谷公園 両議院」。コーラス部分の「ラメチャンタラ ギッチョンチョンで パイノパイノパイ」は、日本の民謡の囃子言葉(はやしことば)のように、歌の調子を整える意味のない言葉です。 「ラメチャン」の「ラメ」は、「デタラメ」の「ラメ」に由来しています。 「ギッチョンチョン」は、お座敷歌『ぎっちょんちょん』から転用されたものかもしれません。 「パイ」は、大正時代の日本で当時新しかった洋食の一つです。 後述の「フライ」も同様です。 「パリコト」の「パリコ」とは、1918年に終結した第一次世界大戦の戦後処理として、1919年に締結されたパリ講和会議を意味しています。

A:Bard
榎本健一の東京節の歌詞「ラメチャンたら ギッチョンチョンで パイのパイのパイ パリコとパナナで フライ フライ フライ」の意味は、以下のとおりです。
ラメチャン:デタラメの「ラメ」のこと。デタラメちゃん。深い意味はない。
ギッチョンチョン:擬音語。何かが動く音を表す。
パイのパイのパイ:囃子言葉。意味はない。
パリコ:当時流行していた菓子の名前。
パナナ:当時珍しい果物の名前。
この歌詞は、歌の途中で挟まれる囃子言葉です。深い意味はなく、単にリズムをとったり、雰囲気を盛り上げたりする役割があります。
「ラメチャン」は、当時流行していたアメリカの流行歌「テネシーワルツ」の歌詞に出てくる「ラメイ」という単語を、日本語風に言い換えたものです。この「ラメイ」は、もともとは「デタラメ」という意味の俗語です。
「ギッチョンチョン」は、何かが動く音を表す擬音語です。この場合は、ラメちゃんが何かを動かしている様子を表現していると考えられます。
「パイのパイのパイ」は、囃子言葉としてよく使われる言葉です。特に意味はありません。
「パリコ」と「パナナ」は、当時流行していた菓子と果物の名前です。この歌詞では、ラメちゃんがこれらの食べ物を手に持っている様子を表現していると考えられます。
このように、この歌詞は、当時の流行や風俗を反映した、ユーモラスで軽快な囃子言葉です。

A:Claude
榎本健一さんの「東京節」は、東京の街とそこに住む人々の生活をユーモラスに歌った曲です。
ご指摘の歌詞部分の意味は以下のように解釈できると思います。
「ラメチャンたら」は、ラーメン屋さんの景色を描写しています。
「ギッチョンチョンで」は、客がラーメンを食べている様子。ギッと音を立てて食べるさま。
「パイのパイのパイ」は、ラーメン屋の中で流れているBGM。「パイパイパイ」と重なる音。
「パリコとパナナで」は、ラーメンの具材。チャーシュー(パリコ)とメンマ(パナナ)。
「フライ フライ フライ」は、揚げ物料理をいうスラング。ラーメン屋のメニューを描写。
このように、東京のラーメン屋の雰囲気を特徴的な言葉で表現した歌詞だと思われます。
言葉遊び的な面白さがある反面、一見意味不明にも見える歌詞ですが、東京の街空間を体感できる造形力があると思います。

A:perplexity
「ラメチャンたら ギッチョンチョンで パイのパイのパイ パリコとパナナで フライ フライ フライ」という歌詞は、榎本健一の「東京節(パイのパイのパイ)」の一部であり、東京の観光名所や食べ物を歌ったものです。以下が歌詞の全文です。
(略:全文に近いが間違っているものを引用している)
したがって、「ラメチャンたら ギッチョンチョンで パイのパイのパイ パリコとパナナで フライ フライ フライ」という歌詞は、東京節の中で、東京の観光名所や食べ物を歌ったものであると言えます。

A:ChatGPT
「ラメチャンたら ギッチョンチョンで パイのパイのパイ パリコとパナナで フライ フライ フライ」という歌詞は、榎本健一が作詞した歌「東京節」の一部です。この歌詞は、単語やフレーズが連なっていて、直接的な意味があるようには見えません。
一部の歌詞やフレーズは、単にリズムや響きを重視した言葉遊びや音の連なりである場合もあります。したがって、この歌詞の具体的な意味は作者の意図や解釈者によって異なる場合があります。
もしこの歌詞が特定の文脈や背景を持っている場合は、それに基づいた意味があるかもしれませんが、一般的な日常会話や表現としては認識されていないようです。


そうこうしているうちに,OpenAIのChat-GPTは更なる進化を遂げているようだ。

[1]榎本健一(NHKアーカイブス人物)

2023年11月20日月曜日

チェレンコフ放射(2)

チェレンコフ放射(1)からの続き

チェレンコフ放射の理論的取り扱いは,JacksonのClassical Electrodynamicsをみればよいらしいけれど,それは宿題ということで・・・。

ここでは,原子炉から核燃料から放出される典型的なβ崩壊の高エネルギー電子の速度を求めてみる。キセノン133,ヨウ素131,セシウム137などでは,ベータ線の運動エネルギーは0.5MeVから1MeV程度である。この電子の速度が 媒質中の光速度 $c'=c/n$ を超えるかどうかが知りたい。つまり$T$求めた $v/c$が,$1/n$を超えれば,チェレンコフ光が観測される。ただし,$n$は媒質(水)の屈折率である。

電子の静止エネルギーを$mc^2$,電子の相対論的な全エネルギーを$E=\dfrac{mc^2}{\sqrt{1-(v/c)^2}}$,相対論的な運動エネルギーを$T=E-mc^2$とおく。これから$\ v/c\ $を$T$の関数として表せばよい。
$(T + mc^2)^2=\dfrac{(mc^2)^2}{1-(v/c)^2}$であるから,$1-(v/c)^2 = \dfrac{1}{(T/mc^2 + 1)^2}$
$\therefore (v/c)^2 = 1 - \dfrac{1}{(T/mc^2 + 1)^2} = \dfrac{(T/mc^2)^2 + 2(T/mc^2)}{(T/mc^2 + 1)^2}$

$\therefore v/c = \dfrac{\sqrt{T^2 + 2mc^2 T\ \ }}{T+mc^2} = \dfrac{\sqrt{\tau^2 + 2\tau\ \ }}{\tau + 1}$

ここで,$\tau=T/mc^2$である。例えば,$\tau=\{0.5, 1, 2\}$に対して,$v/c=\{0.74,\ 0.86,\ 094\}$なので,水の屈折率の逆数$1/n= 1/1.33 = 0.75$ をほぼ超えることがわかる。


写真:チェレンコフ放射の例(Wikipediaから引用)


2023年11月19日日曜日

チェレンコフ放射(1)

チェレンコフ放射は知っている。

屈折率が,$n=\sqrt{\varepsilon_r \mu_r} > 1\ $である媒質中の光速度は,$c' = c/n$と真空中より遅くなる。ここで,$\varepsilon_r,  \mu_r\ $は無次元の比誘電率と比透磁率である。この媒質中を進む荷電粒子が媒質中の光速度を超える場合,波面の作る包絡線に垂直な方向に生ずるのがチェレンコフ光である。典型的な例は水に浸かった原子炉中の核燃料が出す放射線から生ずる青白い光である。

ところで,荷電粒子が電磁波を放出するのはそれが加速度運動している場合である。上記の放射線(高エネルギーのベータ線)は媒質の水の中を等速度で運動している。

砂川さんの理論電磁気学によれば,点電荷の座標を$\bm{r}(t_0')$,観測点の座標を$\bm{x}$,粒子の位置から観測点に向かう単位ベクトルを$\bm{n}(t_0')=\dfrac{\bm{x}-\bm{r}(t_0')}{|\bm{x}-\bm{r}(t_0')|} = \dfrac{\bm{x}-\bm{r}(t_0')}{R(t_0')}$とする。
さらに次の量$\ \bm{\beta}(t_0') = \bm{\dot{r}}(t_0')/c\ $と$\ \alpha(t_0')=1-\bm{n}(t_0')\cdot \bm{\beta}(t_0')\ $を定義した。
ただし,$t_0'\ $は$\ t_0'=t -|\bm{x}-\bm{r}(t_0')|/c \ $の解であり,$t_0'$のなかに$\bm{x}$が含まれる。

スカラーポテンシャル$\phi(\bm{x},t)$とベクトルポテンシャル$\bm{A}(\bm{x},t)$は,次式で与えられる。
$\phi(\bm{x},t)=\dfrac{e}{4\pi\varepsilon_0}\dfrac{1}{\alpha(t_0') R(t_0')}$
$\bm{A}(\bm{x},t)=\dfrac{\mu_0 e}{4\pi}\dfrac{\bm{\dot{r}}(t_0')}{\alpha(t_0') R(t_0')}$
また,電場$\bm{E}(\bm{x},t)$と磁場$\bm{B}(\bm{x},t)=\dfrac{1}{c}\bm{n}(t_0')\times \bm{E}(\bm{x},t))$は,
$\bm{E}(\bm{x},t)= \dfrac{e}{4\pi\varepsilon_0}\Biggl[ \dfrac{(\bm{n}-\bm{\beta})(1-\bm{\beta}^2)}{\alpha^3 R^2}+\dfrac{(\bm{n}-\bm{\beta}) (\bm{n}\cdot \bm{\dot{\beta}}) - \alpha \bm{\dot{\beta}}  \} }{c\alpha^3 R} \Biggl]_{t_0'}$
$\bm{B}(\bm{x},t)= \dfrac{\mu_0 e}{4\pi\varepsilon_0}\Biggl[ \dfrac{(\bm{\beta}\times \bm{n})(1-\bm{\beta}^2)}{\alpha^3 R^2}+\dfrac{ (\bm{\beta}\times \bm{n})(\bm{n}\cdot \bm{\dot{\beta}}) + \alpha \bm{\dot{\beta}} \times \bm{n} \} }{c\alpha^3 R} \Biggl]_{t_0'}$

加速運動する荷電粒子から生ずる電磁波は$\bm{\dot{\beta}}$の項からくる。これを含まない項は,遠方で$R^{-2}$で減衰するのでエネルギーの放射には関係しない。一方,媒質中で光速を超える場合は,$\bm{\dot{\beta}}=0$ではあるが,同時に$\alpha=0$になる可能性がある。そこでこの項が消えずに残るというのが,ものの資料[1]の説明だったが,イマイチよくわからない。フーリエスペクトル以降の計算を追えていない。

結局,チェレンコフ放射についても自分はよくわかっていなかった。まあそんなものだ。

図:チェレンコフ放射のイメージ(github-nakashoから引用)





2023年11月18日土曜日

三角形(2)

三角形(1)からの続き

半径$r$の円が内接する直角三角形で,円の接点が斜辺を$a,\ b$に分割するものの面積が,簡単な表式 $S=ab$で与えられるので,ピタゴラスの定理を経由せずに幾何学的に説明できそうな気がする。


図:長方形への図形断片の埋め込み

そこで,一辺が$a,\ b$の長方形を対角線で分割した△AQB=$ab/2$に,前回の図における図形の断片がきれいに埋め込めるのではないかないかと思ってトライしてみる。2種類の三角形は底辺の$a, \ b$と高さ$r$をそのままにして頂点の位置をずらせばきれいにおさまる。すなわち,△BCQ=$ar/2$と△ACQ=$br/2$である。

長方形を対角線で分割した三角形△AQBからこれらの面積を除けば,薄い三角形△ABCが余る。したがって,この面積は,△ABC $= \frac{1}{2}\{ab -(a+b)r\}$である。前回の円を内接する直角三角形の面積条件は,$S=r(a+b+r)$だったので,△ABC $= \frac{1}{2}\{(ab -S) + r^2\}$となる。つまり,△ABC =$r^2/2$を満足する場合に$S=ab$となって,断片が長方形に収まることになる。

何だか回りくどい話になって,図形からすんなりと説明できたとはいいにくかった。上の例では,面積を保ったままA→Eに変形すれば,△ABC=△EBCになっているのだけれど,それは特別な場合である。

2023年11月17日金曜日

三角形(1)

三角形の面積を求めるという小中学生向けの問題があったので,ちょっと一般化してみた。逆行進化している自分のレベルにふさわしいかもしれない。

次のような直角三角形OABと内接円Cがあって,斜辺ABと円Cの接点をQとする。与えられているのは,AQとBQの長さであり,このとき直角三角形ABCの面積を求めるのが課題だ。


図:三角形の面積の問題

円の半径を$r$として三平方の定理から$r$の方程式を作れば簡単に解ける。この方程式は次の形になる。$(a+r)^2+(b+r)^2= (a+b)^2$。これを整理すると,$r*(r+a+b) = ab$となる。これから二次方程式の解の公式を使ってモチャモチャしていたのだけれど,その必要はなかった。

三角形OABの内接円の中心CからA,B,P,Q,R点と結ぶ線を切り離して並べ替えると,高さがrで幅が(r+a+b)の細長い長方形ができる。この長方形の面積$ r*(r+a+b)$がもとの三角形OABの面積と等しい。すなわち先ほどの式の右辺である$ab$が答えの面積になるのだった。

2023年11月16日木曜日

相似形

同志社大学田辺キャンパスの交隣館でおにぎりをほおばりながら,Facebookをみていたら,芳賀さんが,盛大に算数の教科書をディスっていた。

教科書にあるのは,相似の例題だ。校舎の高さを影の長さと相似の関係で求めたり,校舎を挟んで立っている2本の木の距離を,運動場の一点から見た角度を測って求めるといった問題だ。目的がその値を求めることなら,屋根から降ろしたロープの長さをはかるとか,学校の設計図面で確認するほうが早いだろうという主張である。みんなで便乗して,算数教育の在り方を批判していた。

そこまで,いじめなくてもよいのにと思ってしまう。自分がこどものときに読んだ「数のふしぎ・形のなぞ」では,ピラミッドの高さの図り方としてこの方法が説明されていた。ピラミッドの横に棒を立ててその影の長さと棒の長さの比を求めておく。ピラミッド頂上の影の長さに相当する距離を測って,先ほどの比を当てはめればピラミッドの高さが推定できるというものだ。

小学校の算数の時間に運動場で測量を行ったときは,確かに角度を正確に求めるのが難しかった。なので,気持ちがわからないこともない。しかし,この方法を理解しておけば,簡単にいろいろな量の概算ができるので,それはそれでよいのではないでしょうか。最も効率のよい方法だけに固執する必要はない。環境次第では,何かの役に立つこともあるかもしれない。

そんなことを考えながら,キャンパスを歩く学生さんを見ていたら,みんな自分の身長の1.3倍くらいの影をずるずる引きずりながら,スイスイ進んでいた。今ごろの太陽高度は南中時刻でも40度はないのだろう。ところが自分の視線の方向を変えると,短くなった影はお団子のように人の後に張り付いていた。そうだ,影を観察する方向によってその見かけの長さを制御できるのだ。

ということで,身長と影の長さが等しくなる方向を選択すれば,ややこしい計算なしに,目標とする高いものの高さがほぼ水平距離に置き換えることができる。それはそれで学びの成果になるだろうか?

図:高さ2の影が4のとき,方向によって長さ2にできる

[1]太陽高度(一日の変化)(CASIO高度計算サイト)

2023年11月15日水曜日

忘れ物

天理市では火曜日と金曜日が燃えるゴミ(生ゴミ)の回収日だ。

それをすっかり忘れていて,非常勤の授業に向かおうとして玄関を出ると,まだゴミが回収されずに積んであった。さっそく,取って返して台所の生ゴミをまとめて・・・とバタバタして駅に着いたら,iPhoneとボールペンを携帯するのを忘れていることに気がついた。

切符は予備のPiTaPaで対応し,ボールペンは大学に予備が置いてある。ところが,電車の待ち時間や電車の移動時間の間が持たない。普段は,常時iPhoneに没頭して情報入力しているので,それがないと調子が狂うのだ。

かばんの中には文庫本(ハヤカワ文庫JA,AIとSF)があったので,早速これを読んで無事に過ごすことができた。最近,本を読めなくなったことの最も大きな原因が,インターネットにつながった端末からの情報入力に取って代わられていることをあらためて痛感する。

大学のエスカレータを上りながら縷々考えてみると,毎日の食事によって身体の再生産を維持するのと同様に,人間は日々情報をインプットしてアウトプットしないと生きていけないのでは思い至った。そして,人やその食物が地球物質循環の中に位置づくように,人と世界の相互作用の中にある情報循環のイメージが浮かんできた。


2023年11月14日火曜日

宮腰

NHKの午後の列島ニュースでは,地方局からの話題が2局ずつペアで紹介されていく。

先日,青森局と金沢局のニュースがあって,後者は石川県立美術館国立工芸館の共催による「皇居三の丸尚蔵館収蔵品展」の入場者数がこのシリーズで過去最多だったといいうもの。

問題は,前者のほうだ。青森県中泊町の旧家で大正ステンドグラスが残る宮越家の美術品や書状などが展示されたという話題だ。同時に宮越家の離れと庭園の秋公開も行われている。宮越というのは金沢市金石(かないわ)の旧地名なので,そのルーツは金沢なのではないかと思って調べてみた。

自分が勘違いしていたのは,金石の旧称は宮越ではなく宮腰(みやのこし)だったこと。しかしルーツが石川県というのは正しかった。中泊町の宮越家の初代は江戸時代に加賀国江沼郡(石川県加賀市)の宮ノ越から移住して代々庄屋などをつとめていた。

なお,金沢の宮腰は,初代中村歌右衛門(1714-1791)の出身地でもある。


写真:宮越家の大正ステンドグラス(青森県観光情報サイトから引用)

2023年11月13日月曜日

昼夜時間(3)

昼夜時間(2)からの続き

昼夜時間の式はできたものの,現実の天理市(緯度:34.6°,経度:135.8°)の日の出,日の入りとは差がある。日本標準時との違いは,0.8°だから3分程度になるが,それだけでは足りないような気もする。

そもそも,冬至(12月22日)に向けて,日の出が最も遅くなる日(1月8日ごろ)と日の入りが最も早くなる日(12月6日ごろ)が違うのは,単純な円軌道モデルでは説明できないのではないか。

とりあえず,国立天文台計算室のこよみの計算によって,天理市の日の出時刻(7:04〜4:43),日の入り時刻(16:46〜19:14)を求めてみた。±2時間半くらいは変化する。この結果,昼時間は9時間50分から14時間30分までかわる。南中時刻はこんなにフラフラするものだったのか。


図:天理市の日の出時刻,日の入り時刻,昼時間,南中時刻

[1]こよみの計算(国立天文台計算室)

2023年11月12日日曜日

奈良の神社寺院

融通念仏宗からの続き

散歩の犬棒で思い出したが,うちの近所の寺には融通念仏宗が多いような気がする。早速,調べてみよう。日本の神社・寺院検索サイト「八百万の神」という,センスのよい有難い場所がある。誰が運営しているかというと,株式会社 INFO UNITE というあまり聴いたことのない会社だった。とりあえずなんとなくニュートラルな印象で,他にも日本の住所というデータベースを運営している。なんなのだろうか。

その結果を以下に整理した。まず。奈良県の神社・寺院数を市町村別に見たベスト7だ。五條とか宇陀などがかつては相対的に栄えていたのかもしれない。
奈良市   396   12%
五條市   238   7%
宇陀市   229   7%
橿原市   214   7%
大和郡山市 193   6%
桜井市   156   5%
天理市   155   5%
その他   1602  50%
また,神道の系列別のベスト7は次のようになる。八王子神社も散歩でよく見かける。
春日系列  146  20%
八幡系列  139  19%
祇園系列  83  11%
天神系列  81  11%
八王子系列 39  5%
伊勢系列  35  5%
稲荷系列  24  3%
その他   188  26%
問題の仏教の宗派別のベスト7は次のとおりで,融通念仏宗は第3位,奈良県ではけっこう存在感を示していた。
浄土真宗本願寺派 424  24%
浄土宗      318  18%
融通念仏宗    205  11%
高野山真言宗   164  9%
真宗大谷派    100  6%
曹洞宗      74  4%
真宗興正派    67  4%
その他      441  25%
なお,天理市には融通念仏宗の寺が25あった。なるほど近所でよく見かけるわけだ。


写真:融通念仏宗本山は大阪市平野区の大念彿寺(Wikipediaから引用)

2023年11月11日土曜日

融通念仏宗

朝の散歩でたまに珍しいものに出会う。犬も歩けば棒に当たる

11月8日の水曜日,北東に向って天理喜殿町セブンイレブンまで往復の散歩。ひまわり保育園の近くの八幡神社の向いに融通念仏宗大念仏寺の大きめののマイクロバスが停まっていた。中には運転手だけが待っている。こんな朝早くから法事でもあるまいし?そのまま北東に向かった。

帰り道,どこかから鉦の音が響いてくる。やっぱり法事なのだろうか。しかし,あいかわらずマイクロバスは空のままだ。さらに進むと,止んだと思った鉦の音が大きくなってきた。そちらの方をみると遠くに黒衣の坊さんが歩いているようだ。そのあたりには確か墓地があるが,そこに向かって田んぼの中のあぜ道を5人ばかりが鉦を鳴らしながら進んでいた。

うーん,これはいったい何だろう。本当にあのマイクロバスと関係があるのか。あのなんの変哲もない田んぼの中の小さな墓地に一体何があるのだろうか。謎は謎をよぶばかりだ。

かくして,人間は僅かな情報から,これらを結びつける説明や物語を紡ぎ出そうとする。


写真:鉦の音のなる方向に見えたのは(撮影 2023.11.08)

[1]融通念仏縁起絵巻(Wikipedia)
[2]融通念仏縁起絵巻(クリーブランド美術館)
[3]融通念仏縁起絵巻(シカゴ美術館)


2023年11月10日金曜日

昼夜時間(2)

昼夜時間(1)からの続き

早速これを使ってグラフを描いてみる。Mathematicaで次のコードを打ち込んだ。純粋な三角関数との違いを調べるために,ピーク値を採寸して,そのコサインを比較のために描かせてみた。
In[1] = f[τ_] := ArcSin[Cos[τ]*Sin[-23.4/180*Pi]]
In[2] = t[τ_, θ_] :=  12 (1.0 - 2.0/Pi*
ArcTan[(Tan[f[τ]]*Tan[θ])/Sqrt[1 - (Tan[f[τ]]*Tan[θ])^2]])
In[3] =  (12 - t[0, 34.58/180*Pi])/2
Out[3] = -1.15705
In[4] = g1 = Plot[t[τ, 34.58/180*Pi], {τ, 0, 2 Pi}, 
  PlotStyle -> {Blue}]
In[5] = g2 = Plot[{12 + 2*1.15705 Cos[t]}, {t, 0, 2 Pi}, PlotStyle -> {Red}]
In[6] = g3 = Plot[
  10*(t[τ, 34.58/180*Pi] - (12 + 2*1.15705 Cos[τ])) + 
   12, {τ, 0, 2 Pi}, PlotStyle -> {Green}]
In[7] =Show[{g1, g2, g3}]
図:地球の夜時間(冬至〜夏至〜冬至)

グラフの縦軸は夜時間で,中心値が12時間になっている。横軸は1年を位相2πに換算したもので,冬至−春分−夏至−秋分−冬至に対応する。違いがわかりにくいが,青が今回の結果,赤が三角関数,緑が(青−赤)の10倍である。つまり,純粋な三角関数からの ズレは,たかだか10分程度ではあるが,少しモジュレーションがかかったような関数になっていた。


2023年11月9日木曜日

昼夜時間(1)

11月8日は立冬。朝の散歩で,日の出がしだいに遅くなっていく。

日の出から日没までの昼時間の一年における変化は,三角関数になっているはずだと思い込んでいた。秋分から冬至にかけての位相はπ//2で立冬がπ/4だ。そうならば今ごろの日の出時間は午前6時(秋分)と午前7時(当時)の間の1/√2の午前6時40分ごろになるはずだが,実際には午前6時20分ごろ。もしかして,三角関数ではない?

小学校のころ緯度と太陽高度の関係ははさんざん勉強したはずだ。小学校5年生の理科では,春夏秋冬の良く晴れた一日が太陽高度の観察日にあてられた。授業中でも1時間毎に運動場に出て,画用紙をおいた画板に立てた棒の影の長さを記録するのだ。それでも昼の時間が1年でどう変わるかの定量的議論には至らなかった。小学生には計算できません。


図:昼夜時間計算のための座標

地球の自転軸を$z$軸とする$xyz$座標系を考える。$x$軸$\phi$方向すなわち$\bm{s} = (\cos\phi, 0, \sin \phi)_{xyz}$から太陽光線がくるものとする。ただし,公転面に対する自転軸の傾き$\phi_0=23.4^\circ$として,$-\phi_0 \le \phi \le \phi_0$の範囲で振れることになる。

この角度$\phi$は,地球の公転面を$XY$平面とした座標系において,自転軸の方向ベクトル$\bm{a}=(\sin \phi_0, 0, \cos \phi_0)_{XYZ}$と,太陽から見た公転軌道上の地球への方向ベクトル$\bm{t}=(\cos \tau, \sin \tau, 0)_{XYZ}$の内積の角度を$\pi/2$から引いたものになる。つまり,$\sin \phi = \sin \phi_0 \cos \tau$となる。


昼夜時間を求めるために,半径を1とした地球の$xyz$座標系で考える。

太陽入射光線の方向ベクトル$\bm{s}=(\cos\phi,0,\sin\phi)$に垂直な平面と地表面が交わる大円を考えると,大円上の点$(x,y,z)$は$x\cos\phi + z \sin \phi = 0$を満たす。なお,この大円より$x$軸負方向側が夜である。

緯度$\theta$の観測点は,$z=\sin\theta$なので,$x^2+y^2=\cos\theta^2$の小円上にある。さきほどの大円との交点が昼夜分界点となるので,これらの連立方程式を解いて,$(x_b,y_b)$を求めればよい。その結果,$x_b = -\tan\phi \sin\theta,\ \  y_b=\pm \sqrt{\cos^2 \theta-\tan^2 \phi \sin^2 \theta\ \ }$となる。

昼夜分境界点までの角度$\alpha$は,$\tan\alpha = \dfrac{x_b}{y_b} = \dfrac{\tan\phi \tan \theta}{\sqrt{1-\tan^2 \phi \tan^2 \theta\ \ }}$となる。
$\alpha$はラジアン単位なので,$\alpha \cdot \frac{180}{\pi}$で度になおし,さらに $\frac{24}{360}$をかけて,$\alpha \cdot \frac{12}{\pi}$が時間単位の値だ。この2倍が12時間からの夜時間の余剰部分に相当する。
これから夜時間の長さは,$T= 12\Biggl\{1 - \dfrac{2}{\pi} \tan^{-1}\Bigl( \dfrac{\tan\phi \tan \theta}{\sqrt{1-\tan^2 \phi \tan^2 \theta\ }} \Bigr) \Biggr\}$で与えられる。


2023年11月8日水曜日

レジカート

近所のマーケットに行ったら,レジカートなるものが導入されていた。

新しもの好きなので,早速トライアルしてみた。セルフレジは使ったことがある。重さを検知していて待たされたり,やり直しになったりと,一見便利そうなのだけど実際には店員さんが入力したほうが早く終る。

レジカートは普段のカートより大きめで,タブレットサイズの入力チェック画面とがっちりしたバーコードスキャナーがついている。使い方は特に教わる必要もなくスイスイと入力できる。なかなかいいぞ。キーウィ6個パックをカゴに入れたところで,隣に半値の4個パックがあった。早速前者を取り消したが,これも問題なくできた。スキャンしてカゴに入れた商品はタブレット画面で確認できるし,総額も常にわかって便利だ。

問題は,プリペイドカードで払う必要があることだけれど,これもアプリに切り替わりつつある。いいような悪いような。ICOCAもiPhoneにしているので,何かのときに困るような気もする。とりあえず,チャージしてレジカート専用窓口に並んだところ,対応のお姉さんが一瞬でチェックして通過することができた。

設備投資にはコストがかかりそうだし,レジカートがかさ張るのが問題かもしれないが,労働人材不足の折り,さらにこの手のシステムの普及は進むのだろう。


写真:はじめてのレジカート(撮影 2023.11.6)