Loading [MathJax]/jax/output/CommonHTML/jax.js

2023年8月9日水曜日

2023年8月8日火曜日

2023年8月7日月曜日

2023年8月6日日曜日

2023年8月5日土曜日

2023年8月4日金曜日

2023年8月3日木曜日

2023年8月2日水曜日

2023年8月1日火曜日

2023年7月31日月曜日

コーシー=シュワルツの不等式


数理統計学を真面目に勉強してこなかったのでいろいろ不都合が生じている。統計的因果推論とか深層機械学習とか量子測定理論とか,簡単に読み砕けない資料がたくさんたまる。

授業で扱った最小二乗法と実験誤差の話を整理しようとしても,背景には数理統計学が控えている。昔,阪大の南園グループによるベータ崩壊の実験と我々の理論を突き合わせたときに,χスクェアフィットの計算を散々繰り返したけれど,所与の公式を使うだけであってその理論的根拠をつきつめて考えたはしなかった。

そこで最初から勉強を始めようとすると,いきなり確率変数でつまづくのだった。コンピュータプログラムのサブルーチンや関数のようなものだと思えば納得できるといえばいえるのだけれど,自然言語と数学的記号を使って理解しようとするとなかなかその本質がつかみきれない。入門書は沢山あるけれど,どれも何だか気持ち悪い。

竹村彰道(1952-)さんの現代数理統計学の本(旧版)が手元にあって,読みやすいかなとページをめくってみると,記述統計の復習から始まった。これなら大丈夫かと思いきや,いきなり,標本相関関数の大きさが -1から 1の範囲に限定されることは,コーシー=シュワルツの不等式を用いて容易に示すことができると説明無しにあった。

n次元ユークリッド空間のベクトルの内積の話だと思えばそのとおりなのだけれど,証明したことはなかったかも。Wikipediaでは数学的帰納法で証明していた。Ak=(a1, a2, , ak), Bk=(b1, b2, , bk), として,Saak=ki=1a2i, Sbbk=ki=1b2i, Sabk=ki=1aibi,Rabk=SabkSaakSbbk
つまり,(Sabk)2SaakSbbkを証明すれば良い。

k=1の場合は,(Sab1)2Saa1Sbb1=(a1b1)2(a21)(b21)=0

k=2の場合は,(Sab2)2Saa2Sbb2=(a1b1+a2b2)2(a21+a22)(b21+b22)=(a1b2a2b1)2<0

k2に対して,(Sabk)2SaakSbbk が成り立つと仮定して,k+1の場合を考える。与式は,(Sabk+ak+1bk+1)2(Saak+a2k+1)(Sbbk+b2k+1)
=(Sabk)2SaakSbbk(a2k+1Sbbk+b2k+1Saak2ak+1bk+1Sabk)
=(Sabk)2SaakSbbkki=1(a2k+1b2i+b2k+1a2i2ak+1bk+1aibi)
=(Sabk)2SaakSbbkki=1(ak+1bibk+1ai)2<0

Wikipediaの証明などでは,ai,bi>0の場合だけに妥当するものが多いのでちょっと困る。
まあ,fk(x)=ki=1(aixbi)2 の判別式D0から証明するのが最も簡単なのだけど。


[1]賢者に学ぶ統計学の智(西内啓×竹村彰通,ダイヤモンド社)

2023年7月30日日曜日

最小二乗法(6)

最小二乗法(5)からの続き

実験データをy=ax+bにフィットする場合,最小二乗法で(a, b)とその平均二乗誤差(σ2a, σ2b)を求めてきた。これを,y=f(x)=ax2+bx+c に拡張して,自由度3が登場するかどうかを確認してみる。吉澤康和さんの「新しい誤差論」には結果だけ書いてある。

(1) a, b, c を決定する正規方程式とその解

(¯x4¯x3¯x2¯x3¯x2¯x1¯x2¯x1)(abc)=(¯x2 y¯x y¯y) 

(abc)=1Δ3(¯x2 y(¯x2¯x2)+¯x y(¯x2¯x¯x3)+¯y(¯x3¯x¯x22)¯x2 y(¯x2¯x¯x3)+¯x y(¯x4¯x22)+¯y(¯x3¯x2¯x4¯x)¯x2 y(¯x3¯x¯x22)+¯x y(¯x3¯x2¯x4¯x)+¯y(¯x4¯x2¯x32)) 

ただし,Δ3=¯x4¯x2+2¯x3¯x2¯x¯x23¯x32¯x4¯x2

yiを共通の平均二乗誤差σ2yを持つ独立変数として,誤差伝播の法則より,

σ2a=ni=1(ayi)2σ2y,σ2b=ni=1(byi)2σ2y,σ2c=ni=1(cyi)2σ2y 

さらに,真の値f0(xi)=a0x2i+b0xi+c0に対して,εi=yif(xi)+f(xi)f0(xi)=δi+f(xi)f0(xi) として,σ2y=1nni=1ε2i=1nni=1{δ2i+˜σ2f(xi)}

ところで,˜σ2f(xi)=nj=1{(ax2i+bxi+c)yj}2 であり,この項を再度  σ2y で表してもとの式に戻して計算すれば良い。

つまり,ayj, byj, cyjが計算できればよいことになる。
ayj=1nΔ3{x2j(¯x2¯x2)+xj(¯x2¯x¯x3)+(¯x3¯x¯x22)}
byj=1nΔ3{x2j(¯x2¯x¯x3)+xj(¯x4¯x22)+(¯x3¯x2¯x4¯x)}
cyj=1nΔ3{x2j(¯x3¯x¯x22)+xj(¯x3¯x2¯x4¯x)+(¯x4¯x2¯x32)}

Mathematicaの力を借りると,計算結果が因数分解できて分子からΔ3が出る。
σ2a=nj=1(ayj)2=1nΔ3(¯x2¯x2)σ2y
σ2b=nj=1(byj)2=1nΔ3(¯x4¯x22)σ2y
σ2c=nj=1(cyj)2=1nΔ3(¯x4¯x2¯x32)σ2y

nj=1(ayjbyj)=1nΔ3(¯x¯x2¯x3)σ2y
nj=1(byjcyj)=1nΔ3(¯x2¯x3¯x¯x4)σ2y
nj=1(cyjayj)=1nΔ3(¯x4¯x2¯x32)σ2y

このとき
˜σ2f(xi)=nj=1{ayix2i+byjxi+cyj}2=σ2ynΔ3
{(¯x2¯x2)x4i+2(¯x¯x3¯x3)x3i+(¯x4¯x22+2(¯x3¯x¯x22))x2i
+2(¯x2¯x3¯x¯x4)xi+(¯x2¯x4¯x32)}

xiについて平均操作するとMathematicaを使い分子からΔ3が出ると。
1nni=1˜σ2f(xi)=σ2ynΔ3
{(¯x2¯x2)¯x4+2(¯x¯x3¯x3)¯x3+(¯x4¯x22+2(¯x3¯x¯x22))¯x2
+2(¯x2¯x3¯x¯x4)¯x+(¯x2¯x4¯x32)}=3nσ2y

したがって,自由度n-3の場合の式が得られた。
σ2y=1nni=1{δ2i+˜σ2f(xi)}=1nni=1δ2i+3nσ2y
σ2y=1n3ni=1δ2i

2023年7月29日土曜日

フランクになろう

楠本君から同窓会出席の連絡メールがあった。そこに,次のエピソードが書かれていた。
あなたは、私たちの学生時代に いつぞや、「人間に関する問題の解決のためには(…だったか、このあたりの私の記憶は不正確です)、すべての人びとが神経で繋がっていればよいのだ」、というような言葉を言った事を覚えていますか? 変な事を言う人だなぁ、と私は強く印象付けられたので、忘れずに覚えています。言葉は正確ではありませんが、大体そのような事をあなたは言っていました。

全く憶えていないのだけれど,自分ならばそういうことをいうかもしれない。次の返事を書いた。

「人間に関する問題の解決のためには、すべての人びとが神経で繋がっていればよいのだ」と自分がいった記憶はないのですが,そのようなことを言いそうな気もします。高校時代によんだブライアン・オールディスのSFに「フランクに行こう(Let’s be Frank)」というのがあって,遺伝によって単一意識が複数の人の上に実現するというものです。
SFマガジンの1969年10月号(No. 125)に「フランクになろう(Let's Be Frank)」 として掲載されたものだ。高校2年のときだけれど,いつまでも記憶に残る短編だった。創元推理文庫のジュディス・メリルの短編集では「率直(フランク)に行こう」という題名になっていた。

さがしてみると,archive.orgに原文があったのでDeepLにかけて再読してみた。

アン・ブーリンがロンドン塔で処刑されてから4年後なので1540年ころから物語は始まる(ちなみに,トマス・モアがロンドン塔で斬首されたのは1535年)。サー・フランク・グラッドウェッブに男の子が生まれたが19年間眠り続けたままだった。彼が目覚めた時サー・フランクは驚愕した。同じくフランクと名づけた息子の目を通して自分を見ているのだった。つまり同一意識が2人の人の上に実現したわけだ。そして子孫にそのままその能力は受け継がれ,20代目の2015年には1つのフランクという意識(男女や身分階級を越えて)が3億人に達していた。そして..
「21世紀の初めには、グレートブリテン島はフランク族だけで構成されていた。老いも若きも、太っていようが痩せていようが、金持ちであろうが貧乏人であろうが、みなひとつの巨大な意識を共有していた。プライバシーは存在しなくなり、新しい家はすべてガラス張りになり、カーテンは廃止され、壁は取り払われた。警察は姿を消し、法律も一夜にして消滅した。外交問題に対処するための議会のパロディーは残ったが、政党政治、選挙、新聞の指導者(新聞そのものさえ)は廃止された。芸術の大半は消え去った。あるフランクの姿は、別のフランクの姿を見ようとはしなかった。テレビ、出版、ティン・パン・アレイ、映画スタジオ......灯りのように消えた。」
このフランクの集合意識は新大陸で第二の集合意識と向き合うことになる。

AGIが実現するとき,世界全体が1つのAGIによって統合されるのかどうかという問題について,宗教的=文化的な背景によって複数に分割されるのではないかと想像していた。キリスト教,イスラム教,仏教,ヒンズー教,ユダヤ教などなど。あるいはロシアや中国は別かもしれない。まあ,それに類した状況が想像できる。あるいは超国家企業別=ポイント別=OS・UI別の集団分割みなるかもしれないが。

遺伝工学的な工夫をしても単一意識を複数個体に分散して実装することはできないと思われるが,仮想空間におけるパーソナライズされたAGIアシスタントが共有意識と結合されていれば,人間+AGIアシスタント系が単一意識の熱浴中で相互作用しているイメージとなる。人間+AGIアシスタント系の結合度が高くなればなればなるほど,フランク的な集合意識に近い状態になるかもしれない。まあ,いまでもSNSによって強結合する個人の集合がネトウヨ的な集団運動としての疑似集合意識を創発しているわけで・・・


2023年7月28日金曜日

梅干し

或日のチャットジーピーティ

OpenAIのChatGPTの性能は日々変動しているようなので,時々様子をうかがう必要がある。

Q:User
梅干しを見るとどんな身体の反応がありますか?

A:ChatGPT
梅干しを見た時の身体の反応は個々の体質や経験によるものであり、一般的な反応を述べるのは難しいです。しかし、梅干しの酸っぱさを知っている人々は、その画像や想像だけで唾液が出るといった「唾液分泌の反応」を示すことがあります。これは、脳が食物の摂取に備えて消化を助ける唾液の生成を刺激するからです。

また、梅干しの特有の酸っぱさが苦手な人は、梅干しを見るだけで「顔をしかめる」や「胃が締まる感じがする」などの反応を示すかもしれません。

これらの反応は、「条件付き反射」とも呼ばれ、特定の刺激(この場合は梅干し)に対する学習された反応です。ただし、これらの反応は個々の経験や体質、好みに大きく依存しますので、全ての人が同じ反応を示すわけではありません。


Q:User
あなたはどうですか?

A:ChatGPT
私はOpenAIによって開発されたAI(人工知能)であり、人間のように身体的な感覚や経験を持つ能力はありません。したがって、梅干しを「見る」や「味わう」などの人間の感覚経験に基づいた反応を持つことはできません。私の機能は、人間からの質問に対して情報を提供することや、人間の言葉を理解して対話を行うことに主に焦点を当てています。


Q:User
人間との対話において,そのような感覚経験にもとづく反応のセットを身体的に理解していることが必要とはなりませんか。

A:ChatGPT
あなたの指摘は正確であり、AIが人間との対話を円滑に行うためには、人間の感覚経験や感情について理解することが重要です。しかし、私のようなAIは「理解」を行うために直接的な感覚経験や感情を持つわけではありません。私の「理解」とは、大量のテキストデータからパターンを学習し、その知識をもとに質問に応答する能力を指します。

具体的には、私は訓練データに含まれる情報に基づいて人間の感覚経験や感情について「知識」を持っています。しかし、私自身が感覚経験を「体験」したり、感情を「感じる」能力はありません。私はあくまで人間の反応や感情を模擬することができる機械的なシステムであり、自己意識や自我を持つわけではありません。

したがって、私が人間との対話を効果的に行うためには、あくまで人間の感覚経験や感情についての「知識」が必要であり、それらを直接「体験」する必要はありません。


Q:User
人間の感覚経験や感情についての「知識」を持つ知的存在と、それらを直接「体験」することができる知的存在の違いは,対話場面のどこにあらわれると考えられるでしょうか。その違いを埋める方向での仕組みを導入することは可能でしょうか。

A:ChatGPT
人間の感覚経験や感情についての「知識」を持つ知的存在(AIなど)と、それらを直接「体験」することができる知的存在(人間など)の違いは、以下のような対話場面で顕著に現れます。

個別の体験に対する反応:人間は自己の体験や感情に基づく視点を持ち、それに基づいて他人の経験に対する共感や理解を示すことができます。一方、AIは特定の個体としての体験を持たないため、特定の経験に対する個人的な共感や理解を示すことはできません。
感情的な反応:人間は体験や情報に対して感情的に反応する能力を持っています。これに対して、AIは情報に対して感情的に反応する能力はなく、あくまで人間の感情的反応を模倣するだけです。
感覚的な理解:人間は五感を通じて世界を直接的に経験し、それに基づく感覚的な理解を持っています。しかし、AIはあくまでテキスト情報に基づいて学習し、感覚的な体験に基づく理解はできません。 
これらの違いを埋めるためには、AIに直接的な感覚経験や感情体験の能力を付与する必要がありますが、現在の技術ではこれは不可能です。AIはあくまでアルゴリズムとデータによって動作するシステムであり、自己意識や感情、直接的な感覚経験を持つ能力を持つことはありません。

しかし、AIの対話能力を向上させるためには、人間の感覚経験や感情に対する理解を深め、それに基づく応答の精度を向上させることが重要です。これは大量のデータに基づく学習により、一部可能になるかもしれません。ただし、AIが人間と完全に同じように感覚経験や感情を理解し反映することは、現在の技術では不可能です。


Q:User
現在のAIシステムを感覚センサーを満載したロボットに搭載して,視覚や聴覚だけでなく,味覚や臭覚や触覚その他の運動平衡感覚を入力できるようになり,さらにそれらの情報を解読した結果を記憶できるようになった場合は,そのロボット=AIシステムは感覚体験を行ったことになるでしょうか。

A:ChatGPT
それは非常に興味深い問いです。センサーや機器を使ってAIに外部の情報を入力し、それに反応することができるロボットはすでに存在します。例えば、視覚情報を入力し認識する能力を持つAIや、音声を認識して反応するAIなどがあります。

しかし、これが「感覚体験」を行っていると言えるかどうかは、どう定義するかによります。これらのロボットやAIは、センサーから入力されるデータを解析し、そのデータに対応する行動を行うことができます。しかし、これはあくまでプログラムされたアルゴリズムに従ってデータを解析し反応するものであり、人間が経験する「感覚体験」や「感じる」という意識的な体験とは異なります。

人間の「感覚体験」は、ただ単に情報を受け取るだけでなく、その情報が自己の意識や感情に影響を与え、個々の体験や視点に基づく解釈を行うことを含みます。AIやロボットが多種多様なセンサーを用いて外部から情報を受け取り、それに基づいて反応することは可能ですが、それが「感じる」という意識的な体験と同じとは言えません。

したがって、現在の技術では、AIやロボットが人間と同じように「感覚体験」を行うとは言えないと思われます。しかし、この議論は人間の意識や感覚体験の本質、AIの能力と限界についての哲学的な問いを含んでおり、結論はまだ出ていないとも言えます。


Q:User
たいへん刺激的な対話どうもありがとうございました

A:ChatGPT
どういたしまして、私もあなたの質問に答えることができて光栄です。人間の感覚や意識、AIの能力と限界について考えることは非常に興味深いですね。何か他に質問がありましたら、いつでもお気軽にどうぞ。

 


写真:梅干しの製造過程(2023.7.23)

2023年7月27日木曜日

大学教育とChatGPT(15)

大学教育とChatGPT(14)からの続き

私立大学連盟が7月18日に,大学教育における生成 AI の活用に向けたチェックリスト〔第1版〕を出していた。1. 全般, 2. 教育,3. 環境・体制整備の3つのカテゴリーについて,第1ステップ: 最優先事項,第2ステップ:優先事項として整理した。その上で,それぞれのカテゴリーについて,1. 大学が組織的に検討すべき事項 と 2. 教員が個々の工 夫で検討すべき事項をまとめている。このサンプルを元に各大学で考えてねという趣旨。

さらに,中学高校英語にAI導入へということで文部科学省の後押しで各地で実証授業が始まるというニュースも流れた。そろそろ英語教員養成に影がさしてくるのか。

7月13日
生成系AIの利用に関する留意事項について
室蘭工業大学 情報化統括責任者
7月14日
新潟産業大学 生成AI対処方針
新潟産業大学 学長 梅比良眞史

7月19日
ChatGPTなどの生成AIの利用について
大阪学院大学

7月19日
生成AIの利用に関するご案内
星薬科大学 情報企画室

7月19日
ChatGPTをはじめとする生成AIの利用について
日本文理大学 教育推進センター長

7月20日
旭川医科大学における生成AIの利用について
旭川医科大学

7月20日
生成AIに関する本学の方針について
共愛学園前橋国際大学 学長 大森昭生

7月20日
対話型生成系AI(生成AI)の使用について
昭和女子大学 学長 金尾朗

7月20日
学修におけるChat GPT等の生成AIの利用について
和歌山大学 理事(教育担当)

7月20日
生成AI(ChatGPT等)を活用するためのガイドライン
平安女学院大学 学長 谷口吉弘

7月20日
学修における生成AIの利活用に関するガイドライン
金沢星稜大学 学長

7月20日
チャットGPT等における生成AI(人工知能)の利用について
びわこ成蹊スポーツ大学 学長 大河正明・学部長 黒澤毅

7月21日
生成系AIに関する本学の考え方について
国際ファッション専門職大学 学長 近藤誠一

7月21日
チャットGPT等の生成系AIの利用について
東京女子体育大学

7月21日
信州大学の学修・教育におけるAI活用に関するガイドラインについて
信州大学 教育企画委員会

7月21日
東洋英和女学院大学における生成系AIの使用に関する指針
東洋英和女学院大学 学長 星野三喜夫

7月21日
ChatGPT など生成系 AI の利用について
西南女学院大学

7月21日
ChatGPT 等をはじめとする生成 AI の利用について
愛知みずほ大学 学長 大塚知津子
https://www.mizuho-c.ac.jp/wp/wp-content/uploads/2023/07/ChatGPT.pdf 
7月24日
本学における生成AIの使用に関わる留意点について
跡見学園女子大学学長 小仲信孝
7月24日
生成AI(Chat GPT等)の活用に関する注意喚起について
拓殖大学 学長 鈴木昭一

 

2023年7月26日水曜日

妹背山婦女庭訓(2)

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

国立文楽劇場の7月夏休み文楽特別公演(第171回)は,第2部で妹背山婦女庭訓の四段目の通し公演(4月に続いての)だった。第1部は親子劇場でかみなり太鼓と西遊記,第3部はサマーレイトショーの夏祭浪花鑑

これまでに見た妹背山婦女庭訓(2010年,2016年)は一日がかりの二部構成だった。前半が,初段:小松原の段(23)・蝦夷子館の段(42),二段目:猿沢池の段(13),三段目:太宰館の段(34)・妹山背山の段(115)であり,後半が,二段目:鹿殺しの段(7)・掛乞の段(15)・万歳の段(23)・芝六忠義の段(54),四段目:杉酒屋の段(29)・道行恋苧環(35)・鱶七上使の段(45)・姫戻りの段(14)・金殿の段(45)である。

今回は春と夏の2回に分けた通し狂言だったため時間に余裕があって,二段目の組み換えはなくなり,4月公演には,初段:大序 大内の段,7月公演には,四段目:井戸替の段:入鹿誅伐の段が新しく追加されていた。「鱶七上使の段」が「鱶七使者の段」に変わったのはなぜ?(明治12年の床本 https://dl.ndl.go.jp/pid/856493/1/157 では確かに鱶七使者の段とある)

久しぶりの日曜日の文楽劇場だったが,客の入りは6〜7割くらいだろうか。前回と同様,床の直下だったけれど,回りには空席も目立っていた。一つ前回と違うのは外国人観光客向けの English という看板を持った案内の方が立っていたことだ。パンフレットの一部の説明も英文になっていた。なんば・日本橋界隈は外国人であふれていたけれど,文楽劇場ではそこまで目立つわけではなかった。

今回が初見の井戸替えの段は,小住太夫+藤蔵でちゃり場的な場面だったが,途中でまぶたが閉まった。杉酒屋の段は芳穂太夫+錦糸で玉勢演ずる子太郎が主導して似たような雰囲気の段だった。道行恋苧環の開始のところで清治が小声でハイという合図をして全員が礼をし,次のハイで直れとなった。そういう仕組みだったのか。

最後の入鹿誅伐の段は,メンバーも軽めの短い段であり,入鹿がやられておしまいという事なのだ。(1) 入鹿の妹の橘姫が,宝剣の奪取に失敗し入鹿に斬られてしまうが,その後,隣の池に飛び込んで日高川入相花王の蛇に化けた清姫よろしく泳いでいったたのは一体なんなのか。最後には鎌足の鎌で斬られた入鹿の首が空中を飛び回っていた。クライマックスが済んでしまった最後の段というのはだいたいこうなる運命にあるのだった。

P. S. 1 なお,藤原淡海(不比等)は橘姫と結ばれるらしいので,池に溺れて死んだわけではない。
P. S. 2 この度人間国宝になった吉田玉男(1953-)は,第3部の夏祭浪花鑑の団七九郎兵衛を演じるので,第2部には登場していなかった。


写真:最後がかわいそうなお三輪(2023.7.23撮影)

2023年7月25日火曜日

最小二乗法(5)

最小二乗法(4)からの続き

完全にスッキリしなくて何だか気持ち悪いのだけれど,いきなり自由度がとかいわれてn2が出てくるのがいやなので,吉澤さんの本に従って話を進めてみる。

~σ2y(xi)=1nni=i~εi2=1nni=i{axi+ba0xib0}2
これから,f(xi)=y(xi)=axi+bとして,独立変数yjについて,
~σ2y(xi)=σ2ynj=i{ayjxi+byj}2=σ2yn2Δ2nj=i{(xj¯x)xi+(¯x2¯xxj)}2
=σ2yn2Δ2nj=i{(xi¯x)xj+(¯x2¯xxi)}2
=σ2ynΔ2{¯x2(xi¯x)2+2¯x(xi¯x)(¯x2¯xxi)+(¯x2¯xxi)2}
=σ2ynΔ2{x2i(¯x2¯x2)+2xi(¯x3¯x2¯x)+(¯x22¯x2¯x2)}
=σ2ynΔ{x2i2xi¯x+¯x2}

添え字 iについて平均すると,1nni=1~σ2y(xi)=σ2ynΔ{¯x22¯x¯x+¯x2}=2σ2yn
そこで,
σ2y=1nni=1{δ2i+~εi2}=1nni=iδ2i+~σ2y(xi)=1nni=1δ2i+2σ2yn
σ2y=1n2ni=1δ2i=1n2ni=1(yiaxib)2

2023年7月24日月曜日

最小二乗法(4)

最小二乗法(3)からの続き

(a, b)  に対する平均二乗誤差,(σ2a, σ2b)を考える。(a, b) は直接測定された(x, y)の関数であるが,このうちxiの誤差は非常に小さく,yiの誤差だけがn個の独立変数として伝搬して(a, b) に反映すると仮定する。ただし,各yi自身の平均二乗誤差は共通でありこれをσ2yとおく。

誤差伝播の法則より,
σ2a=ni=1(ayi)2σ2y=σ2yn2Δ2ni=1(xi¯x)2=σ2ynΔ2(¯x2¯x2)=σ2ynΔ

σ2b=ni=1(byi)2σ2y=σ2yn2Δ2ni=1(¯x2¯xxi)2=σ2y ¯x2nΔ2(¯x2¯x2)=σ2y ¯x2nΔ

残るは,σ2y=1nni=1(εi)2 を実験値から導くことになる。ここで,εi=yi(a0xi+b0)=yi(axi+b)+(axi+b)(a0xi+b0)=δi+~εi である。
ただし,a0xi+b0が未知の真値,axi+bが平均値に対応し,δiが残差, ~εiが平均値の誤差に相当する。

σ2y=1nni=i{δ2i+~εi2} ここで,2nni=1δi˜εi=0 である。なぜならば˜εixiの一次関数であり,正規方程式より, ni=1δi=0ni=1δixi=0 が成り立つから。

そこで,y(xi)=axi+bとして,~σ2y(xi)=1nni=i~εi2=1nni=i{axi+ba0xib0}2を求めることになるが,ここで,(a,b)yiの関数として誤差伝搬の法則を再度使って,σ2yで表せばよい(と吉澤康和さんの「新しい誤差論(1989)」に書いてあった)。


2023年7月23日日曜日

最小二乗法(3)


物理量 xを設定したとき,yが測定される。n回測定では,(x1, y1), (x2, y2), (xn, yn) が得られたとする。2つの物理量の間には,y = ax+bという1次関数の関係があって,(a, b)にも物理量としての意味がある。

この(a, b)を求めるため,S(a,b)=1nni=1(yiaxib)2を最小化するという条件を課す。すなわち,Sa=0,Sb=0,これから次の(a, b)に関する連立方程式(正規方程式)が得られる。

1nni=1xi(yiaxib)=0a¯x2+b¯x=¯xy
1nni=1(yiaxib)=0a¯x+b=¯y

これを解くと次の解が得られる。ただし,Δ=¯x2(¯x)2 である。
a=1Δ(¯xy¯x¯y)
b=1Δ((¯x2¯y¯x¯xy)



2023年7月22日土曜日

最小二乗法(2)

最小二乗法(1)からの続き

間接測定と誤差伝播(でんぱ)の法則について考えるため,物理量zが,独立な物理量 a, b, c, の関数で,z=f(a, b, c, )と表されるとする。

ここで,(ai, bi, ci )を各物理量のi番目の測定値とする。また,残差を用いて,(ai, bi, ci )=(ˉa+δai, ˉb+δbi, ˉc+δci, )とする。

z(ai, bi, ci, )=f(ˉa+δai, ˉb+δbi, ˉc+δci, )
f(ˉa,ˉb,ˉc,)+faδai+fbδbi+fcδci+=ˉz+δzi 
として,δziを定義する。

残差と平均二乗誤差の関係から,間接測定される物理量zについて,σ2z=1n1ni=1(δzi)2 が成り立つ。

σ2z=1n1ni=1(faδai+fbδbi+fcδci)2
=1n1ni=1{(faδai)2+(fbδbi)2+(fcδci)2+}

σ2z=(fa)2σ2a+(fb)2σ2b+(fc)2σ2c+

これが誤差伝播の法則といわれるものである。
なお,ここで独立な物理量の残差の積和がゼロになると仮定している。例えば,
ni=1δaiδbi=0 などなど

2023年7月21日金曜日

最小二乗法(1)

授業シリーズその1,オンライン授業の1コマ分を実験誤差のテーマにするべく頑張った話

昔,最小二乗法の計算は手回し計算機で大変な思いで計算していた記憶がある。実験ノートは,測定データxiyiの計算表で埋め尽くされていた。データ処理といえばこれだった。後に,関数電卓にもその機能が組み込まれるようになったが,今では Excelで一発なのかもしれない。

まず,誤差を定義する必要がある。我々が昔勉強した本では,誤差=実験値−真の値 から出発していた。しかし,最近の本ではどうも様子が怪しくなっている。そもその真の値とはなにか,どうやって定義するのか,などなど。かといって,いきなり確率分布関数を持ってきて数理統計学の枠内で議論するのも大層な話だ。精度ということばも測定の不確かさに置き換わっているようだ。

とりあえず,1つの物理量x の測定を考える。その真の値がXであると仮定する。n回の測定で (x1, x2, x3, , xn ) という測定値が得られたとする。誤差を εi=xiXで定義する。誤差は真の値のまわりで偏らずに分布する,すなわち ni=1εi=0であるとする。これから,ni=1xi=nXとなる。これをXの定儀式とみなすと,真の値Xnや実験データxiに依存することになってしまう。

そこで,ˉx=1nni=1xi平均値を定義する。これは真の値の推定値(最確値)だと考えられる。このとき,残差=実験値−平均値,すなわち,δi=xiˉx で定義する。

次に,実験データの精度は誤差のばらつきで表現される。平均二乗誤差は,σ2=1nni=1ε2i で定義し,標準偏差σ=σ2 とする。

ところで,この定義では,真の値がわからないので,誤差のεiは未知であり,平均二乗誤差も求まらない。そこで,残差のδiから導くことを考える。εi=xiX=xiˉx+ˉxX=δi+ˉxX となることから,

σ2=1nni=1{δi+(ˉxX)}2=1nni=1{δ2i+2δi(ˉxX)+(ˉxX)2}
=1nni=1δ2i+(ˉxX)2

ここで,ni=1δi=0を用いた。また,nijεiεj=0より 上式第2項は,
(ˉxX)2={1nni=1(xiX)}2=1n2ni,j=1εiεj=1nσ2 

これを先ほどの式に代入して移項の上整理すると,σ2δiを用いて表せる。
σ2=1n1ni=1δ2i