2012年6月8日金曜日

積分の複素数近似

電子の静止エネルギーは電子の電荷の元になる微小電荷を寄せ集めたものになっているであろうとのことです。

微小電荷を寄せ集めるということのイメージは図のようにある力で距離Rが無限遠から0まで押し込んだときのエネルギーの総量ということのようです



実際には荷電粒子間の反発力は距離の逆二乗なので単純には距離0までは押し込めそうもないということですが神様の力で仮に押し込めたとします。

説明の都合上まだ一体にはならないとして話しますが
次に力を緩めると電荷は反発して離れていきます。
力学的には神様の手は反発に伴ってエネルギーを回収しているので無限遠方まで離れたということで全体としてはエネルギーは+-0となるはずと予想できます。

この神様の力の入れ方、抜き方(蛇足ですが神様は手を離すのではなく反発がそうなるように加減しているのかもしれません、離してはエネルギーを回収できないですから・・・)はどうなるのが一番自然かということですが電子は質量があるので、押込むのには+とーの電荷が引き合うときの運動と同じになるのがこれは反発と対称ですし自然だろうと仮定してみます。

また微小電荷には電荷に相当する微小質量があるとします。
そうすることで自然にエネルギーと速度が対応付けできます。

押し込むエネルギーを計算するときに反発力*押込み距離で計算していますが
押込み力=反発力とした場合は力が平衡するので微小電荷は動きません。
反発力*押込み距離で計算するということはつまり
押込み力=反発力 +(+とーの電荷が引き合うときの加速度を生じさせる力)ということです。
押込んだ後は次に反発です、神様が何とかして押込んだ後手を緩めると微小電荷はそのまま反発して無限のかなたに飛んでいくものとします。神様はこのときエネルギーを回収するはずです。
この反発はもともと集まった電荷の作用と質量なのでそのスピードやエネルギーが計算できます。
なので押込みに要するエネルギーは蓄積された電荷と質量に数値として現れていることになります。  つまり、押込みに要するエネルギーはほとんど既に現れていて特になにか余分なものとして
溜まることはないのではないかと思います。

もし余分に蓄積されるエネルギーがあるとすればこの押込みと反発のエネルギー差になるはずです。
ところがエネルギーは逆一乗なので距離0でどうやっても発散は避けられません。
分割して積分するにしても距離0では ∞ー∞(∞から∞を引く)を計算することになりますがこれを単純に0と断定することはできません。

そこで発想を変えて実数から虚数に拡張してルートは実数軸から微小 -iδ(lim ->0とします)だけ虚数側を通します。

これで極を避けることができて連続したルートを通すことができます。
したがって積分範囲は -∞ から ∞までとることができます。(δ=0 でない限り連続になります)
さらに積分範囲が -∞ から ∞ ( うろ覚えですが片側有限の範囲でも極を内に入れていれば
OK・・・?)だと留数定理により円積分できて計算が楽になります。

これで計算できるのはエネルギーの距離積分の押込みと反発分の差なので蓄積される
エネルギーそのものではないですがなにか参考になるかもしれません。

結果、エネルギー積分差は Ke^2/2 となります。

蓄積されるエネルギーそのものは出せませんでした。しかし、もし押込んで一体になっていくとするとエネルギー積分差にかかわる何らかの量が余分に蓄えられているのではないか?ということが伺えて面白いと思います。

それにしてもエネルギーの積分されたもの、これは一体なんでしょ? 力を積分したものがエネルギーポテンシャルならこれはエネルギーを積分した何らかのものでしょうがスカラだしよくわかりません。 ここを見られた方をガッカリさせるとは思いますが、多分エネルギー自身を距離で積分するということ自体あまり意味がなく、まあ・・なんといいますか逆一乗の留数マジックという気がしてます。 その上私自身勉強中なので・・・・
しかしこの留数で極を避けるテクニックは使用制限はないのでしょうか?例えば積分範囲がエネルギーとなる式があっても複素数に拡張して扱えばOK? いくらなんでもいい加減すぎるような気がします。多分この話はどこかに根本的なエラーがあるとは思いますが・・・・。

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
2012/6/12 追記 
(新しくページをつくるのが面倒くさいので追記で対応することにしました。読む方は■■■■■■■ で囲んだここは飛ばして最後に読んだほうが経緯が良くわかると思います)

どうにもこうにもエネルギーが出せないのであきらめていたんですが

ひらめきました。  キーワードは対数発散

ワイスコップ( Weisskopf) の計算したもので量子電磁力学を用いて電子の自己エネルギーを計算すると、{m^2C^2e^2/(πhC)}log{h/(mCx)}となるそうです。
この自己エネルギーも 距離 X をゼロに近づけると発散します。
これはくりこみ理論(Renormalization theory)と関連性があるものらしいです。

私が参考にしたサイトは  TOSHIの宇宙 電子の自己エネルギーとDiracの海 です。
<http://maldoror-ducasse.cocolog-nifty.com/blog/2006/12/post_5d9e.html>

この式は極の量子電磁力学的近傍でないと成立しませんが
この近傍では 距離 X で微分すると 力は 1/X と逆一乗になります。
つまり留数でエネルギーを出せます。
ここで問題は今まで円積分はe^iθ と半径1mで計算していましたが
これでは量子電磁力学的近傍じゃないのでだめなんじゃないかということですが、じつは式中でキャンセルされるため計算を簡単にするために半径を省いていました。
なので半径αをかけて線積分路は αe^iθ となります。
そしてこのαは lim -> 0 として絞り上げても良いのです。
したがって、留数は極点近傍だけの式から出すことができます。
なので、極から離れると力は逆二乗分がでてきますが、これは考慮する必要はありません。
あくまでも極近傍の式が重要です。
ということでやっとすっきりした感じです。トリビア1ゲットだぜ。
これでタイトルも「電子に蓄積するものについて」 から 「電子の自己エネルギーについて」にかえました。

追記 終わり
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
追記    2012/6/13,22 ~        
発散しないエネルギーが出たということで喜んでいたんですが
留数と実数軸(距離)との関係を確かめたところ

wiki  の Application of Jordan's lemma
http://en.wikipedia.org/wiki/Jordan%27s_lemma
にあるように使い方を間違えていたようです。

integrate(f(x),-∞,+∞) = i2π*留数(極の位置) となるのでエネルギーが結局虚数になってしまいました。複素数を使った副作用が出た感じです。

この値は、i2π{m^2C^2e^2/(πhC)} で ^-60  オーダと非常に小さく、電子に余分に蓄積されるものは実質0のようです。 力学的なものなのでやはりということです。
蓄積されるものとは結局電荷eと質量m(とそのエネルギーmC^2)そのもの以外にはないということでしょうか?
まあ虚数がでても発散して計算できなくなるよりはましと思うほかなさそうです。
タイトルも大げさなので変えました。

定義域が実数の関数の積分(積分範囲が 実数 -∞ から ∞ まで)は留数をつかって複素数関数に近似することができます。今回の場合は上の追記にあるように 実数0+i 10^-60  オーダの虚数(ω)がでてきました。級数が小さい場合は無視して実質0とできます。

積分定数(Ω)も考慮に入れると Ω + iω となります。(円積分なのでこれは何らかのバイアス値というイメージ?)
積分定数Ωは物理なので計測できるものはそれを使ってもいいのでしょうがこのシナリオではどうでしょう仮に質量エネルギーのmC^2とすると結局 mC^2+ iω = 約mC^2 としてもいいのでしょうか?

留数は  wiki  の Application of Jordan's lemma  http://en.wikipedia.org/wiki/Jordan%27s_lemma
にあるように極が半円の中にあれば適用できます。

ところで積分範囲が 0 から ∞ まで の関数も同じように留数をつかって複素数関数に近似することができるのでは? と想像できます。 積分範囲が 0 から ∞ まで の関数は対数発散の微分ですからψ/x と逆一乗形のものです。 これはまともに積分するとlog(x)  なので ∞ - (-∞)  と発散しますす。
ところが、留数で近似すると半円の中に入っていればいいのなら -∞ から ∞  で積分したものとも同じ結果になっても良いのでは?とも思えます。

なので、まるでパラドックスのようなことになってしまいます。

可能性は
1) そもそも本来留数の適用できないケースである(この可能性が一番高いような気がします)。または使い方が適切でない(近似に使用できる積分範囲は-∞ から ∞ までのみ)
2) もし1)がOKであればこれは対数発散の消去をしたことになる。 複素数の世界ではそれでも案外OK? もしそうなら対数発散消去の手法ゲットだぜ

ということになります。 なので今後これの確認をしようと思います。

以上のことは次のようなイメージ図になります。
※ クリックして拡大











■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
追記    2012/7/4 追記

色々なサイトで教えてもらおうと投稿したのですが進展なしということで、結局自分で考えることにしました。 考えてみれば簡単な線積分なので最初から自己解決を目指せばよかったのです。

結果は
1) 積分範囲が 0~ ∞の場合は発散します。 なので発散しないためには積分は -∞ ~ +∞ に限ります。

次に直線部分を線積分したところビックリする結果がでました。
積分結果は以下ですが

積分範囲(-R~R)が微小でも lim id->0 で iπ となります。 つまりこれはデルタ関数です。
wikiで調べると 佐藤超函数 http://ja.wikipedia.org/wiki/%E4%BD%90%E8%97%A4%E8%B6%85%E5%87%BD%E6%95%B0  というものらしいです。 (英語wikiはhttp://en.wikipedia.org/wiki/Hyperfunction ) つまりこれは 留数が出せます。

翻って元のシナリオに戻って考えて見ますと エネルギーは 力 1/x^2 を積分して出しました。
次にエネルギー(1/x の項になります)の通過前後で何か差が出るかな?ということで発散部分での積分を避けるために複素数で近似して -∞  ~ +∞  でエネルギーを距離積分できるようにした。 それで留数を使おうとしたのです。(これ自体は思い違いでラプラス逆変換にしか適用できないことがわかりましたが・・・)
いかん、頭がどこかにトリップしそうです。 ・・・ところがですね、 直線部分だけでも積分するとデルタ関数で原点でのエネルギーを取り出したということらしいのです。
それで、iπ の意味ですが、これは例えると地質調査のボーリングってありますよね、これはこのパイプのようなもので取り出したエネルギー値は iπ を除いたものになります。

結局、「エネルギーの距離積分の差(押込みと反発でのエネルギー距離積分の差)を出そうとした」 と「エネルギーの原点での値を得る」ということが等価なのではないか?
結果発散しないエネルギー値が出てきた・・・。
出てくるのはいいのだけど本当にそれでOK?ということが問題です。

なんとも不思議な話ですが原点での エネルギー(1/x項)の発散がなくなってしまいました。

しかもこの方法を採った一つ一つの動機はごく自然なものです。

・エネルギー(1/xの項)を -∞  ~ +∞  で連続して距離積分するために複素数に拡張して(つまり連続関数にして) lim id->0  として実数に戻す。そうすれば 1/x は奇関数なので押込みと反発の差に関する何らかの値が出るだろう。それは多分原点での蓄積される何らかの値だろう。

なのにいつの間にかそれはデルタ関数になりエネルギーの原点での値を取り出すという意味になっています。

そうすると、原点でのエネルギー値を出すのに力をわざわざ 1/x の形(エネルギーはlogで対数)の形にする必要はありません。素直に力は 1/x^2でエネルギーは 1/x としても良いようです。

本当にこれでもOKでしょか? 全体を観るとパラドックスのようだけれども、一つ一つ見ると妙につじつまがあうような、なにか狐か狸に化かされているような気がしてきました。
これは、私のレベルでは本当に手に負えません、なのでまた方々に相談してみようと思います。

ところでなぜこんなテーマをやっているのかその理由があるのでついでに記しておきます。
電気力線ってありますよね。 +と-で 力線がつながっている良く見る図ですが、よくよく見ると不思議なところがひとつあります。

それは、中心の外側で無限の先に伸びている力線、これはどこに繋がっているのでしょうか?
中心以外はなんとなく有限の範囲で繋がっているのが想像できるのですが中心だけは想像できません。まああるとしても重ねあわせだし、距離の差分値なのでかなり小さいでしょうけど、無限のかなたに消えていくだけのように見えます。 だけど多分、おそらくこれも繋がっていると考えてもいいはずです。 ところがこれをあらわす方法がわからない。 多分実数の範囲ではむりだろうな・・・ということは予想できますが・・・。

だけど、複素数であればどうでしょうか? 複素数の中ではこの線が繋がっているのではないか?つまり複素数の中では円(またはループ)で実数では -∞  ~ +∞  となるのでは?
物理でこの種の対応が存在する場合、実数の -∞  ~ +∞ というのはその裏に 複素数の円が潜んでいるのでは?という気がしています。 つまり ∞ は複素数世界への接続点・・・(標語としてきれいに決まった ^^)ということなのでしょうか?



追記 終わり
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■

とにもかくにも電子は微小電荷が集まってできたものではなくはじめから電荷eと質量mとして存在しているものです、なので理論の組み立てとして集まるという方法をとると何か話がおかしくなるのではないか?(無限の力で神様が押込めるなど・・・)そのために神様の手とかいろいろなトリック的な方法をとらざるを得ないのでは・・・という気がします。   それよりも逆に電子というものは微小電荷を神様が切り離していったら(切り離したら反発力で飛んでいきますから)エネルギーの積分されたものが最後に残る・・・そういうものと定義してもよいのだが、しかしそれを人間が知りたいがため集めるというストーリーを取らざるをえなかったというのが真相のような気がします。 なので微小電荷が集まったらエネルギーの積分されたものも集まっていたというストーリでも、まあなんといいますか、とにかくも計算できさえすれば良いのではないでしょうか。

と言っても数学的に表さなければなりません。 ここで気になるのは留数の積分範囲です。 蓄積を出すのに -∞ から ∞ まで積分しましたが積分範囲が 0 から ∞ まででも良いとすればこれは押込みだけを計算するのと押込から反発までを積分するのと等価ということになります。
これは確認しなければなりませんが片側有限の範囲でも極を範囲内に入れていればOKであれば数学的にもOKということになります。


更に面白いのが
1.  相対論のエネルギー式γmc^2に適用したときでも、元のエネルギー式KQdQ/R
に適用しても結果は同じということです。

KQdQ/R はRの符号で+-があるので感覚的に納得できるのですが γmc^2 は符号がありません。 考えれば複素数だから当たり前ですがこれも感覚的に面白いです。

2. さらに 級数などで 1/R^2 以上の逆べき項や R^n 等の正則なものは0になります。
つまり 1/R 以外は無いものとして扱えます。不思議なものです。

3. 式15(%o15)をみると円積分の結果  e^2*K  と  m*C^2  部 ときれいに分離されているのがわかります。

4. 式8(%o8) を見ると 微小電荷の速度V=C となるのは R=0 です。つまり原点ぎりぎりまで加速は可能です。      また、微小電荷dQとその微小電荷対応分の質量dmは関係なしです。

Vのグラフは  (見やすいように距離Rは対数、 電荷Qが 1.609*10^-37 まで集まった時点のもの)


--------  式  ------------------------------------

電子の電荷 e:1.61*10^-19;        クーロン定数  K:8.976*10^9;     光速 C:3*10^8;
電子の素の質量  m :(注 電子の質量は9*10^-31ですがこのシナリオでの値は蓄積エネルギー分が +αされるため素の質量は多分少し異なるはずです);
Q:微小電荷が集まった電荷  dQ: 微小電荷量  dm:微小電荷に対応する質量  R:距離






2012年5月27日日曜日

引力とエーテル その2

相対的には重心の位置にいる観測者が重心点から発する場と運動する質量の関係を観ているということですが

重心は質量の逆比の内分点となります。質量比ですから重心点は固定してます。
引力によりMとNが引き合いそれぞれ運動していきます。
MとNがそれぞれ無限遠から引き合いある間隔Rとなったときの様子を観てみます
積分範囲が換算質量となっていますので同様に空間縮小変化率を導き出すことができます。


2012年5月26日土曜日

引力とエーテル

引力による空間の縮小について思ったこと。

引力に限らずある遠隔力が働いて物が運動するという場合その運動の様子が遠隔力の働く空間の定義のひとつになりえます。
引力そのものは F=Gmn/r^2 という有名な式です。
物体が運動すると相対性理論によるとその進行方向に空間が縮みます。
これでひとつ遠隔力(場?)を定義してみます。
ある二物体 MとNがありNは無限のかなたから距離 R までMの引力に引かれてやってきたとします。

無限のかなたでの速度は 0 、空間の縮小率は 1(縮小なし)で速度Vに応じて進行方向の
空間が縮小します。 縮小率はローレンツ収縮なので 1/γ = (1-V^2/C^2)となります。

距離Rでの運動エネルギーは力と距離の積分で
En : integrate(-fn(X)*Qn*N, X,R,∞)   fn(x)は単位質量あたりの引力で加速度、Qn は換算質量
なのでこれは重心点を起点にしているということになります。

速度Vでの運動エネルギーは N*V^2/2 (計算の簡便のため近似領域でのシナリオとします)
これで速度が出せます。 速度が出ると空間の縮小率がでます。

一方引力は加速度だから距離Rに応じて空間の縮小率が異なるはずです。
なので空間縮小変化率ψは 1/γ をRで微分すれば良いということになります。

この空間縮小変化率ψが加速度=引力に対応していると仮定すると。
次のことが推測できます。

1、相対的には重心の位置にいる観測者が重心点から発する場とNとの関係を観ている
     ということになります。(こういう関係は対称的なのでたぶんMに対しても同じ・・・のはずです)

2、これは重心点の及ぼす空間の定義で  M>>Nでは Mが重心点となる。

3、加速度にするには C^2*(1/γ)で補正する必要がある

4、C^2が静止質量 M+N に係るとするとそれは静止エネルギーである、場は質量ではなくエネルギーが作っているようだ

5、γは ≧1なので空間縮小変化率は実際の引力の加速度より大きいが
    これは何かしらイメージとしてはNは単位質量あたりにスリップや抵抗を受けているとも思える。

6、γは例えば 地表ではほぼ1 (M(地球)>> N として)
     G*M/(C^2*R) は (7*10^-11 )*(6*10^24)/((9*10^16)*(6*10^6) = 0.8 * 10^-9 <<<< 1

特に4番目は何かしらが空間縮小変化率に抵抗している感じでこれがあのエーテルの流れかもしれないと一瞬思いましたが今までのことを思い出すとやはりちがうのではないかと思います。

どういうことかというと、つまり空間の縮小で場を定義すると言うことで、縮小は近くでは密で
無限遠の先で縮小率1(縮小なし)ということになります。これはエーテルの流れではなく
イメージとしては力線にフィットするという感じがします。第一エーテルの流れだとエーテルが
どこにいくかもイメージできません。
力線は広大無辺の宇宙に薄まり消えていく、悪く言えばゴミ捨て場にすてて後は知ったこっちゃないと無責任でいられます。

一方エーテルの流れというのは思うに実際の運動ベクトルの方向が力線の方向(空間の縮小では密から疎)と反対であるというイメージであることから発想されているような気もします。
どちらを実態と思うかということであり、結局哲学的な話になってしまいます。
あと、速度を出すのに 近似の N*V^2/2 ではなくちゃんと γNC^2-NC^2 で出すと結果は以下のようになります。


結局最後は空間縮小変化率を引力の加速度に合わせこむということで相当怪しい感じで全体的にもそこはかとなくイタイ感が漂ってます。

まず、距離Rでの運動エネルギーを出すのに En : integrate(-fn(X)*Qn*N, X,R,∞)   としていますが
Nの速度による運動質量の増加を入れなくてもよいのか? これはR点より逆に引き離すことを考えるとこの引き離しの速度でエネルギーが場合により変わってくるということであり指標になるエネルギーがでないし、不思議なことにエネルギーは速度と時間が無くても計算できてしまいます。(もし考慮するとしても速度が無限小、変わりに時間が無限大 なので質量増加は考慮しなくてもよいのでは・・・とも思えます)、また一般相対性理論では引力で光が曲がるということで運動質量の無いものにも適用できるもの・・・なのでひょっとしたらポテンシャルエネルギーというのはそういう性質のものではないのかもしれません。
等等いろいろ突っ込みどころ満載で結局トリビアにもならなかったのですが、それでもなんとなく捨てがたしという感じではあります。


2012年3月18日日曜日

4元数と回転

ベクトルの回転公式はコンピュータグラフィックスをする方はご存知の方も多いと思いますが、
ロドリグの公式では (A(軸),B(回転させるベクトル)は単位ベクトル)
C=Bcosθ+A<A , B>(1-cosθ) - (B×A)sinθ  または
C= A<A,B>+(B- A<A , B>)cosθ - (B×A)sinθ
にあるようにクロス積で直交を出して軸回りにcos(θ)とsin(θ)の組み合わせで回転させています。

ロドリグの公式から
式1 : Bcosθ - (B×A)sinθを抜き出すとどうなるでしょうか?

これは下図に示すように楕円回転のベクトル式になり、演算規則はクロス積ですから対応は  ( i , j , k ) = ( jk - kj , ki – ik , ij – ji )  で
ij = –ji = k  ,   jk = –kj = i  ki = –ik = j 

代数的に総積でき、結合則があると拡張すると
k = ij = (jk –kj)( ki – ik) = jkki – jkik – kjki + kjik = jkki – jjk – kii – kkk
なので kk = jj = ii = –1 となります。これは内積(符号は逆)になります。

したがって式1は 代数的に出すと
式2: B * A = B * [cosθ – Asinθ]  =    Bcosθ – (– < B , A >  + (B×A )sinθ)
=  < B , A > + [Bcosθ – (B×A)sinθ] と 内積分(スカラ)が増えて計算されます。

式2は反時計回り回転ですが、その先にベクトルD(Bが2θ回転したもの)がありこれが逆回転してベクトルCになるとします。
式3 : C =  A * D =  [cosθ – Asinθ] *  D = Dcosθ  – (–< A , D > + (A×D)sinθ)
 = < A , D > + [Dcosθ – (A×D)sinθ]

式2と3をみてみると 軸との内積は同心円上にあるので同じです。ベクトルは逆回転になっています。したがって
C =  A * D = B * A  となるので Aの逆演算 を{A}とすると {A} = 1 / A と なるものがあればDが求まるわけです。 これは式2と同じ要領で計算すると
式4 :  {A} = 1 / A = 1 / (cosθ – A sinθ) = cosθ + A sinθ となります。

なのでDは
式5: D = {A} * B * A     となります、回転角は2θ。

しかし、ほんと4元数って機械みたいですね。

2012年3月6日火曜日

三つ子のパラドックス

宇宙船の中の時間は 航路図 と 光路図  http://akimpotos.blogspot.com/?v=0p
にある外から観た玉の飛行時間α と宇宙船中の玉の飛行時間の比で
遅れます。(1/γとなる)
ところがこれは対象的な関係なのでいわゆる双子のパラドックスとなります。
http://ja.wikipedia.org/wiki/双子のパラドックス

で、双子でダメなら三つ子にしてしまえということで、
三つ子の1人を地球に、あとの二人をそれぞれ宇宙船で旅行させます。
宇宙船が帰ってきた時に二人の兄弟にそれなりの年齢差が出るように
それなりの速度差で飛ばしまた同時に帰ってこれるように…など条件を整える
ものとします。

そうすると、もし地球の時間が遅れるのならば帰ってきた時宇宙船は
速度差があるのだから地球の時間の遅れにも差がでます、
つまり帰ってきた兄弟は同時に年齢差のある地球の兄弟を見ることになり
これではオカルトになってしまいます。


なので、時間が遅れるのは宇宙船の方ということになります。



宇宙船の時間(T)の遅れは地球の時間をtとすると
T=t/γ ですがこの分母分子にm*c^2を乗じると
T=t*静止エネルギー/運動エネルギー とエネルギー比の形にできます。
なので時間が遅れるのは実際にエネルギーを与えられた方(つまりwikiにあるように
加速度を持った方)という解釈でもよさそうです。


-----------------------------------------------------------------------------------------------
続きがあります。
There is a continuance.
<http://akimpotos.blogspot.jp/2012/12/2-triplets-paradox-2-deterministic.html>

2012年2月28日火曜日

航路図 と 光路図

「速度合成 と MAXIMAへの式の貼付け」では、玉の航路図を用いましたがこの際なので光路図との関係をまとめておきます。
シナリオは速度Uで飛行する宇宙船の中を速度V=C/3で飛行する玉と光の様子を観るものです。
光は船の中で反射して往復しています。
光が逆行する場合は図中の左にある場合の座標変換関係を流用して表しています。
したがって、本来は玉の飛行中に黄色の矢印で示す宇宙船の図を置くことは適切ではありません。
本来の意味内容を考えると宇宙船の図を置くということは、その場所で座標変換の様子を観たい・・・・・ということであり、単に位置関係を知りたいならそれは斜交座標中のピンクの矢印で示す宇宙船の距離軸(上の方の)への射影でわかります。
また座標関係が反転しており、黄色とピンクで示す関係が繋ぎ目のところで突然入れ替わるという物理的にオカルトチックなものになっています。 地図でいえばメルカトル図法のように角度を保存するために面積が犠牲になっているという感じでしょうか。(この図では事象の一致と座標変換の保存)


2012年2月25日土曜日

速度合成 と MAXIMAへの式の貼付け

「(物理)運動量保存はエネルギー保存則? その12  」で特殊相対性理論での速度合成を逆行列で求めましたが、式に冗長というか不要なものがくっついていたので、結果には影響は出ないのですが、まああまりにひどいので整理して再掲しました。
くっついていた不要なものとは逆変換で求める際に付加していたγで、下図でベクトルA,B を C,Dから求めるのに使っていたものですがまったく必要ないものでした。
Tsuを変換行列とすると
[C*t, t*V]=[A , B]=Tsu*[α , β]  したがって         INVERT(Tsu)*[A,B]=[α,β] で 合成速度は  C*β/α となる。 [C*t, t*V]は更に時間 t  が共通なので γ と同じくこれも省略できますがあまり省略しすぎると物理的意味がわからなくなります。

更に玉にも座標系を与えると [C*t, 0]=[A , B]=Tsv*Tsu*[α , β]  したがって  
INVERT(Tsu)* INVERT(Tsv)*[A,B]=[α,β] で 合成速度は  C*β/α となる。(図はありません)
この C*t も片側が 0 なので実際は0以外の何でもよいです。



ところでいままでMAXIMAの式を添付していたのだけど、MAXIMAへのペーストのテストをウインドウズのメモ帳で確認して問題ないと思っていたのだけど、不思議なことに少しでもメモ帳を編集すると例えば

/* [wxMaxima: input   start ] */
beu:u/c;
/* [wxMaxima: input   end   ] */

のようなブロック単位ではではペーストできなくなってしまうことがわかった。 「/* [wxMaxima: input   start ] */」の部分を除いて 式「beu:u/c;」だけコピーは可能だがこれだと大変だ。どうも文字コードか何かの関係のような気もする。これは掲載した式でも同じ。

仕方ないので下記のように式だけにすると全てひとつのセルに入ってしまう。
動作はちゃんとします、がコメントも入れられなくなってしまうので  簡単なものならいいけど長いものだと使いにくくてどうしようも無い。
要所要所分割して入れると少しは見やすいかもと試してみたのだけどこれだと動作がどうもおかしい。
一行毎に徹底的に分けるか1セルに全て入れるしかないようです。
もうね、ガッカリです、お手上げ。

速度合成の式(MAXIMA)  
     ・全てひとつのセルに入れるようにしてください。(全行一括でコピーしてペースト)
---------------------------------------------
kill(all);
beu:u/c;
gmu:1/sqrt(1-beu^2);
tsu:gmu*matrix([1,-beu],[-beu,1]);
bev:v/c;
gmv:1/sqrt(1-bev^2);
tsv:gmv*matrix([1,-bev],[-bev,1]);
outvct:matrix([c*t],[0]);
invert(tsv).invert(tsu).outvct;
%[2]/(%[1]/c);
ratsimp(%);
outvct2:matrix([c*t],[t*u]);
invert(tsv).outvct2;
%[2]/(%[1]/c);
ratsimp(%);
---------------------------------------------