映画”ドリーム”(Hidden Figures)を観てきた。数値計算好きな私にはめっちゃ面白かった!フリーデンの計算機をいつ使うか?と思ったら、、、あそこで使うか!そしてIBM7090でFORTRANを使い、オイラー法で計算する。 #ドリーム
| Determination of Azimuth Angle at Burnout for Placing a Satellite Over a Selected Earth Position |
« 2017年8月 | トップページ | 2017年10月 »
| Determination of Azimuth Angle at Burnout for Placing a Satellite Over a Selected Earth Position |
こんなコンデンサ見たことない!
本日、加速器の新型電源用のコンデンサーの1號機が搬入されました
このコンテナの中にコンデンサーが入っていて、その容量は1F!
1mFじゃないですよ、1F!
メイカーの工場内でコンテナにコンデンサーを組み込んでもらう、というのは僕の提案ですので、それがそのまま実現して嬉しいです pic.twitter.com/xOcK0tbPnF— 多田 将 (@sho_tada) 2017年9月26日
さらになんという大電流。
そしてこれがその電磁ホーンにその大電流を流すバスバー(銀色の板)とトランス(左の緑色のもの)
電源の出力が32,000Aであるため、このトランスで10倍の320,000Aにします pic.twitter.com/nK0IB9jFtU— 多田 将 (@sho_tada) 2017年9月28日
私も電子部品業界にいて長いですが、こんな容量と電流の話を聞いたことすらない、、、 加速器の世界はすごいなあ、、、、 まあ小さな容量は小さな容量で2001サイズとか、10005サイズとかそれはそれで大変ですのだ。
※電気二重層コンデンサなら1Fとかは普通にあるが、スペックが全然違うんでしょうね。

| Determination of Azimuth Angle at Burnout for Placing a Satellite Over a Selected Earth Position |
そこでビジネスの専門家であるジョンさんは会社を辞め、財団を立ち上げ、研究者に連絡を取り、ウイリアムズ・キャンフィールド博士と知り合う。

iPhone8/8plus, Apple Watch Series3に続いて今度はApple TV 4KをiFixitが分解していた。
ルンゲクッタ8次のDOP853で三体問題の周期解を1223個見つけた話をこの前見て検証してみたりしてましたが、
14次の公式があるのを知ってますか?
これだ!
http://sce.uhcl.edu/rungekutta/
35段14次!
これは倍精度くらいではそもそも使う意味がなさそうなので、任意精度の計算ができるPARI/GPに実装してみる。
まずはいつものようにローレンツ方程式をやってみよう。式はこちら。
dx/dt = -σ*(x-y)
dy/dt=-y-x*z+r*x
dz/dt=x*y-b*z
でプログラムのソースコードはこちら。60桁で計算してますよ。
使い方は、PARI/GPを立ち上げて、
\r Lorenz
でロード。
rk14("ファイル名")
で実行されて、ファイル名のところに保存される。散布図がどうも描きにくそうなので、ファイルに落としてExcelで描いてみた。それがこれ。
なるほど、計算できてそうだ。
で、このブログではルンゲクッタ8次のDormand & Princeルーチン(DOP853)をExcel VBAに移植したりまた普通のルンゲクッタ4段4次をカシオの高精度計算サイトkeisan.casio.jpにUPしたり、結構いろいろやったので全部比較してみよう。
お題は、やはりローレンツ方程式でσ=10,r=26, b=8/3, 時間刻み幅0.01, そして初期値はx=y=z=0.1としてやってみる。
まずは4次、8次、14次の比較。
t=25くらいですでに4次がやばくなってきた。次行こう。
もうt=40過ぎると8次もやばい。
t=100くらいまでいくともうぐだぐだ。
ということで、常微分方程式なら全部4段4次のRunge kuttaでいいや、とか思っているといろいろと落とし穴がでたり。注意しよう。
え?オイラー法?t=1くらいでもう終わってるんだけど、、、
この前、三体問題の周期解が1223個も見つかったという論文を見て、、、
ExcelのVBAにルンゲクッタ8次のDOP853 (Dormand&Prince)を移植したもので計算したりしていたが!
論文はこれですね。
監修の白石さんが第四話はポピュラーな式、と書かれていた通りに最もポピュラーな式でした。波動方程式を不変にする変換としてローレンツ変換を導出。
関連リンク:
だいぶ前に、
なんじゃこれ!力技で1223個も見つけたのか!
1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにDOP853移植ルーチンを使って計算、GIFアニメにしてみた。今回はI.A68(0.5)というもの。えぐい動きをしてるな、、、計算の詳細は続くツイートにて。 pic.twitter.com/BjCu7uJFIQ
— tomo (@tonagai) 2017年9月23日
うーん、なんちゅう動きをしてるんだ(ちょっと表示の刻みが粗すぎるか、、、)
次はIID2(2)ってやつ。
1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAへのDOP853移植ルーチンを使って計算、GIFアニメにしてみた。今回はII.D2(2)というもの。 pic.twitter.com/npOPHjfBSe
— tomo (@tonagai) 2017年9月23日
今年6月に発売されたモッチッチ。生めんのような食感で人気とのこと。
毎度毎度オーストラリア(時差の関係で最も早く買える)に行って分解しているiFixitのスタッフ。
72人の統計学者が連名でp値の閾値を0.05から0.005にしようと言っている。
'One-size-fits-all’ threshold for P values under fire
~ Scientists hit back at a proposal to make it tougher to call findings statistically significant. ~
Justify Your Alpha: A Response to “Redefine Statistical Significance”


http://www.cc.kyoto-su.ac.jp/project/MISC/slide/seminar-s/2011/110604Tsutsui01.pdf
仮面ライダービルド第三話の話数を表す数式 (2^n + 1)/n^2 が整数であるときのnを求めよ、は
IMO(国際数学オリンピック)1990の第三問に出てた。 https://t.co/u6U2d2yQqq— tomo (@tonagai) 2017年9月17日
追記:合ってたみたいです。
黒板の式は、指摘の通りでLeggett-Garg不等式です。巨視的な量子重ね合わせ状態(シュレディンガーの猫状態)があったときに、それが本当に巨視的な量子重ね合わせ状態なのかを検証するための式として提案されたものです。式の構造そのものはベルの不等式と同じです #仮面ライダービルド
— シータ (@Perfect_Insider) 2017年9月17日
関連リンク:
今年のイグノーベル賞が発表になった。
いつものように普通にアップデートしたら、、、
えー! また変な仕様変更を、、、
ちょっと前にAppleの領収書を騙ったフィッシングメールが流行って、私のところにも来た。
奈良に行ったときにたまにたちよるイトーヨーカドーさんですが、
によると、観光型複合商業施設「奈良平城プラザ」とのこと。
仮面ライダービルドの第二話を見て、
仮面ライダービルドの今週の黒板の式は「Arthurs-Kelly-Goodmanの不確定性関係(測定の誤差と擾乱の不確定性関係)」です。日本語の教科書だと沙川貴大、上田正仁『量子測定と量子制御』には出てるけど、多分だいぶマニアック。
— シータ (@Perfect_Insider) 2017年9月10日
これは知らなかった、、、
ここら辺に説明があります。この次の章は小澤の不等式だった!
109シネマズのIMAX次世代レーザーでダンケルク観てきました。70㎜の作品をちゃんと見ようと、expocityのあの超でかいスクリーンで。
トミーズ雅さんの奥さんも使っていてめちゃくちゃ楽だとのこと。
今週から始まった仮面ライダービルド。数式の解説などをしてましたが↓
ゲストブレーンは西野 精治さん(スタンフォード大学教授)。
モリタさんは私の実家の近くにある会社。消防車の売り上げ世界第三位だとか。
http://www.morita119.jp/index.html



遺伝の「優性」「劣性」使うのやめます 学会が用語改訂 というニュースを見た。
| ||||
| ||||
とでている。
斜方晶というのがある。こんなの。
| 斜方晶 (直方晶) |
単純 | 底心 | 体心 | 面心 |
いやどうみても直方体だろう、、、ということでこれ。
http://www.crsj.jp/journal/Vol56/index-56-3-j.html
日本結晶学会誌のこの号の会報を見てみると、
216ページに、
6.晶系名称変更について
坂田会長より, 現在用いられているOrthorhombic の
日本語訳「斜方晶系」が, 適切な訳語でないとの指摘が
あり, 世界結晶年を機会に適切な訳語を日本結晶学会か
ら提案することが提案された. 議論の結果, より適切で
あると思われる「直方晶系(斜方晶系)」の学会決議を諮
り, その使用を地道に提案していくことが承認された.
とある。なので、直方晶系に今後はなっていくんだろう、、、ってなんで最初に斜めの字を使ったんだろう。。。
というニュースを見た。バードストライクではない、ということでしたが、ではバードストライクがあるとどうなるかのシミュレーションは結構なされている。
例えば、これとか、、、
ミズノが面白いボールの発表をしていた。

立川裕二さんのTweetで知った話です。
EEtimesにこんな記事が出てました。
おお、これ私も行ったことがあるんですよね。
モンティホール問題のモンテカルロシミュレーションをExcel VBAと高速メルセンヌツイスタで。
日本では2006年、千葉県船橋に開店したのち、今では8店舗のイケア。日本だけで767億円、世界では4兆円以上の売り上げだそうだ。
ゲストはイケア・ジャパン株式会社 代表取締役社長の ヘレン・フォン・ライスさん。
アグナリッド村のエルムタリッド農場で生まれた。
ゲストブレーンはSF作家の山本弘さんでした。以下メモ。
« 2017年8月 | トップページ | 2017年10月 »
最近のコメント