黒木玄 Gen Kuroki
- いいね数 389,756/311,170
- フォロー 995 フォロワー 14,556 ツイート 293,980
- 現在地 (^-^)/
- Web https://genkuroki.github.io/documents/
- 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
2019年05月12日(日)
@sekibunnteisuu 重力加速度gが高さに依存する場合のモデルは
∂²Φ/∂t² = ∂²Φ/∂x² - g(x)
ではなく、
∂²Φ/∂t² = ∂²Φ/∂x² - g(Φ)
です。xはバネ上の座標で、xに対応するバネ上の点の時刻tでの高さがΦ=Φ(t,x)なので、g=g(x)ではなく、g=g(Φ)です。
ここは結構混乱し易い所だと思う。非線形偏微分方程式。
タグ：
posted at 01:07:31
@sekibunnteisuu φ_{tt} = φ_{xx} - m²φ は例外的に線形でKlein-Gordon方程式。
φ_{tt} = φ_{xx} - sin φ はシャレでSine-Gordon方程式と呼ばれている。Sine-Gordon方程式はソリトン方程式なので色々特殊解の作り方が知られていますが、私は十分にいじったことがないです。今回の件に参考になる特殊解があるかも。
タグ：
posted at 01:12:48
あまりにも可笑しかったので記録に残しておく。
ミュート済みなので発見が遅れた。
すでに笑っている人達がいて存在に気付いた。
【理解のレベルの差ですねー】発言の返信先に注目(笑)。
また一つ伝説が生まれた！！！
この人はどこの何者？ pic.twitter.com/1Ip6H8c87i
タグ：
posted at 01:42:45
非公開
タグ：
posted at xx:xx:xx
小池一夫【ベビーカーがしばしば問題提起されますが、江戸時代から拝一刀は大五郎をベビーカーに乗せてましたよッ！】
クソリプ【「子連れ狼」はフィクションなので、現実のベビーカー云々の根拠に使うのはやめたほうがいいと思いますよ】
この件もクッソ笑えた。
twitter.com/koikekazuo/sta... pic.twitter.com/Xx4cCT8KKQ
タグ：
posted at 02:04:29
非公開
タグ：
posted at xx:xx:xx
@genkuroki >∂²Φ/∂t² = ∂²Φ/∂x² - g(x)
ではなく、
そこはわかっているつもりです。
＞∂²Φ/∂t² = ∂²Φ/∂x² - g(Φ)
です。
Φ(x,t)の取り方にもよると思うけど、xからのずれをΦ(x,t)としたので、
x+Φ(x,t)の位置にあることになるから、
∂^2Φ/∂t^2＝v^2*∂^2Φ/∂x^2-g(x+Φ(x,t)) としました。
タグ：
posted at 04:34:26
@genkuroki でも、どうせ2階微分だから、xがある位置、xからずれ、どちらをΦ(x,t)としてもよさそうですね。それなら、多少でも形が単純なg(Φ(x,t))の方がよさそう。といっても、やはり難しそう。
タグ：
posted at 04:40:05
@genkuroki 結論としては、
この場合も、遠くで起こったことの影響が波及する速度は、∂^2Φ/∂t^2＝v^2*∂^2Φ/∂x^2-g(Φ(x,t)) の場合、速度vを超えない
0<x、t=0で、v^2*∂^2Φ/∂x^2-g(Φ(x,0))が成り立っているなら、
vt<x、t=tで、Φ(x,t)＝Φ(x,0)が成り立つ
となるものなのでしょうか？
タグ：
posted at 04:44:01
@genkuroki バネが不均一に帯電していて、電界がある場合まで考えると
∂^2Φ/∂t^2＝v^2*∂^2Φ/∂x^2-g(Φ(x,t))σ(x) まで一般化できそう。
とても手に負えそうもない。
タグ：
posted at 05:04:39
@genkuroki 物体をバネで数珠つなぎにした離散モデルの方、係数が0でない最低次数の項の次数が2ずつ増えていくというのも、ばね定数と質量が不均一の場合まではやったけど、重力が一定でない場合は、手を付けていない。
そっちを先に考えた方がいいのかな？
タグ：
posted at 05:06:15
@genkuroki 離散モデルと連続モデルの関係もまだ十分に理解していません。
連続モデルでの、影響が速度vを超えて波及することない、というのを、離散モデルの次数が2ずつ上がるという結果の方から求めることは出来るのか？とかいろいろ気になります。
タグ：
posted at 05:08:48
@sekibunnteisuu なるほど、バネの自然長がゼロでない場合を考えているのか。確かにその場合の運動方程式は
u_{tt} = u_{xx} - g(λx+u)
の形になります。ここでλはバネ上の座標と自然長を関係付けるパラメーターです。バネの自然長が0ならλ=0になる。
そうか、この場合はまだ数値計算もしていません。続く #数楽
タグ： 数楽
posted at 05:27:02
@sekibunnteisuu #数楽 続き。離散版は両端を除けば
d²y_i/dt² = (k/m)(y_{i+1}-y_i-L) - (k/m)(y_i-y_{i-1}-L) - g(y_i)
の形になります。Lはバネの自然長。これを
y_i= Li + u_i
で変数変換すると
d²u_i/dt² = (k/m)(u_{i+1}-u_i) - (k/m)(u_i-u_{i-1}) - g(Li+u_i).
続く
タグ： 数楽
posted at 05:34:19
@sekibunnteisuu #数楽 Δx=1/N, (k/m)Δx²=1, x=iΔx, L=λΔx, u(t,x)=u_i(t) とおいて、N→∞とすると
u_{tt} = u_{xx} - g(λx+u)
が得られます。y=λx+u とおくと、
y_{tt} = y_{xx} - g(y).
もとの離散版でも d²y_i/dt² の式内のLはキャンセルして消える。
結局 y_{tt} = y_{xx} - g(y) の基本形に戻ります。
タグ： 数楽
posted at 05:39:11
@genkuroki 自然長の設定が本質的じゃないのはわかりつつも、考えるときは一応自然長を考慮しないと不安になる。何度もやるうちに「自然長は0としても構わない」と確信できるという感じです。
タグ：
posted at 05:52:23
@genkuroki 今現在の私が理解した範囲
重力なし、一様重力で
離散モデル → 次数が2ずつ上がる
↓
極限を取る
↓
連続モデル → 波及速度はvを超えない
タグ：
posted at 05:56:12
@genkuroki 未解決および手を付けていない部分
次数が2ずつ上がる
↓
極限を取る
↓
波及速度はvを超えない
と出来ないか？
重力が一様ではない場合の、離散、連続に関しての考察
タグ：
posted at 05:58:15
@genkuroki 物理をそれほど理解しているわけじゃないので、この人が言っていることが妥当かどうかさっぱりわからなくてスルーしていたのですが、ここから続くツイートを見て、まあ、いろいろ思いました。
twitter.com/sekibunnteisuu...
タグ：
posted at 06:07:03
中央公論の呉座先生の記事読みました。
・百田→文筆家として知的体力が無い。口喧嘩には強いが論理的な文章は苦手
・I沢氏→あくまでも小説家
・H郷先生→学者というよりもはやタレント
三者のうち、H郷先生への批判が「ここ十年学界で言及されてない」「一般書の室も疑わしい」と一番厳しい
タグ：
posted at 06:09:26
@sekibunnteisuu 物理を分かっている人はもっとクリアにかつ地道に説明してくれます。私はすぐにミュートしました。
私は「久々の大型新人登場」だと思っていたのですが、案の定やらかしてくれました(笑)。
twitter.com/genkuroki/stat...
タグ：
posted at 10:39:23
「史料による裏付けのない思いつきを『東大の日本史の教授』というかたがきで世間にばらまくことは、見方によっては百田氏や井沢氏よりも罪深い」
「学界で支持されていない自説を押し売りする前に、学界の主流説を世間に紹介することに努めるべき」
H郷先生は同じ立場なので容赦ない。
タグ：
posted at 10:53:51
テープのり便利ですよね〜、でも案外細かいところはスティックのりの方が使いやすかったりしませんか……ヽ(´o｀;。そんなこんなでまだ完全にテープのりに移行できずにいるのです……。
タグ：
posted at 11:18:08
大学の学費免除が家族4人で270万未満ということは、其れ位の収入の家庭の子供が大学に進学すらほぼ考えられないであろうことを知ってての設定だと思われ。つまりはほとんど国は払う気ないんだよ。
今の授業料免除を受けてる人たちの切り捨て。
タグ：
posted at 11:51:07
きくちさんの所にはまだ来てないのか。
そのツイート1個だけじゃない点もすごいです。
小池一夫さんに【「子連れ狼」はフィクションなので、現実のベビーカー云々の根拠に使うのはやめたほうがいいと思いますよ】と言った人を見かけて以来の快挙。
伝説級の大型新人。
twitter.com/kikumaco/statu...
タグ：
posted at 12:20:20
Ninja DAO | CryptoNi @CryptoNlnjaNFT
うーん、だったら何でその後次官辞任まで何もしなかったんだ。 twitter.com/Kyoko_quetta/s...
タグ：
posted at 14:09:47
Ninja DAO | CryptoNi @CryptoNlnjaNFT
前川さんの一連の話は、いろいろな事について「政権に問題があった」と言いながら、その後、影響力があったはずの次官退任までほとんど何もしていないので（退任はご存知天下り問題）、個人的にはどうなのかなぁ、と思っている。
タグ：
posted at 14:15:15
非公開
タグ：
posted at xx:xx:xx
How to use redirect_stdout(f, [stream])? It's not easy to guess from the document...😑 #julialang
タグ： julialang
posted at 14:28:56
#数楽 #Julia言語
twitter.com/woody_egg/stat...
と同じようなことをやってみました。
nbviewer.jupyter.org/gist/genkuroki...
Lorenz方程式の数値解のアニメーション
DifferentialEquations.jl と Plots.jl を使った簡単な演習問題なのですが、数値解を求める部分よりもプロットや動画作成の方が難しい。 pic.twitter.com/HbsN5bTnwf
posted at 14:36:05
言ってることは特に異論ないけど、そう言ってる人がどんな結果を仕事で出してるのかさっぱり判らんのではどうしようもないな。 twitter.com/paruparu100/st...
タグ：
posted at 14:45:54
Hmm, I could use another overload:
output = nothing
backup = stdout
rd, wr = redirect_stdout()
try
# Do something using stdout
finally
redirect_stdout(backup)
close(wr)
output = String(read(rd))
close(rd)
end
# Use `output` here
タグ：
posted at 14:48:01
For unit testing, it is often useful to redirect stdout and/or stderr to a memory buffer so that we can examine whether the program writes expected message or not. Although I found a solution as above, there should be an easier way... Hmm... #julialang
タグ： julialang
posted at 14:53:41
地毛の「黒染め指導」は廃止を 大学生らが署名活動 www.kyobun.co.jp/news/20190508_... @kyoiku_shimbunさんから
タグ：
posted at 14:55:22
OokuboTact 大久保中二病中年 @OokuboTact
ツイッターの世界では、こういう無敵人がたまに現れるけど、実際に会うと普通の人だと思う。 pic.twitter.com/ydiyQM5dkm
タグ：
posted at 15:01:50
OokuboTact 大久保中二病中年 @OokuboTact
辛い環境で頑張っているのに報われない日々が長く続くと、SNS世界で匿名を使って暴れたくなる人はいるからなあ
タグ：
posted at 15:06:27
#Julia言語 初期条件を変えてもう1つ動画を作りました。N=256
nbviewer.jupyter.org/gist/genkuroki...
Lorenz方程式の数値解のアニメーション
こちらの方は最初から退屈せずにすむ。こういう動画を作る作業は微調節の試行錯誤の時間を結構取られる。
ぶわーっとなるところが結構快感。 pic.twitter.com/duftQo5nHe
タグ： Julia言語
posted at 15:09:59
Wrote a utility (overload) function of `redirect_stdout(f, stream) which can be used for an IOBuffer:
gist.github.com/sgryjp/088f57b...
#julialang pic.twitter.com/j6LL80dtcP
タグ： julialang
posted at 15:14:24
元呟きの二つの動画より少し動きが早いようでより素敵感が。ただ色目にもう少し差がある方が特に別れてからの挙動が分かりやすそうに思う。 #自分でやれ #やっとこないだ環境だけ再構築しただけやし twitter.com/genkuroki/stat...
タグ： やっとこないだ環境だけ再構築しただけやし 自分でやれ
posted at 15:22:14
#Julia言語 についてググるとちょっとアレな状態になって困ると思っている人はGoogleの設定を
地域の設定：アメリカ合衆国
にするとよいです。結構おすすめ。
あと、Julia関係の検索時にキーワードとして
jl
を入れておくという技もあります。 pic.twitter.com/46gPzuIQk2
タグ： Julia言語
posted at 15:25:24
面従腹背とか言いながら、現役時代は上に言われた通り従って、後で自分は上の命令に従っただけと言い訳するのなら、気概も矜持も感じない。本当に面従腹背を座右の銘だというなら表向き上からの命令をハイハイと聞いておいて密かに方向を是正するとかそれくらいの芸当を見せないとただの責任逃れ。
タグ：
posted at 15:29:06
都槻郁稀Totsuki, Yuki @totsukiyuki
@rouninsyurachi 不安…… pic.twitter.com/D6t7euoH3a
タグ：
posted at 15:38:37
非公開
タグ：
posted at xx:xx:xx
OokuboTact 大久保中二病中年 @OokuboTact
@genkuroki ＞ 仲正さんって「いいひと」ですよね。
続きがアップされました。
twitter.com/meigetsus/stat...
タグ：
posted at 15:48:08
ロシアの路上ライブでワガノワの生徒が飛び入りダンスするの巻。
ミュージシャン、こんなシットインされたらテンション上がりまくりだろうな〜（そしてバレエじゃなくても膝から下の動きが完璧ワガノワっ子） pic.twitter.com/mN8PHzpM6i
タグ：
posted at 15:50:03
OokuboTact 大久保中二病中年 @OokuboTact
仲正氏の隔月ネット連載が更新された（以前は毎月更新してたけど、多忙のため隔月に変更）。
「山川ブラザーズ」という言葉は見当たらず、「山川」だけになってしまっているのは悲しい。
雑魚には触れないという姿勢に変更してしまったんだろうか？ pic.twitter.com/CTP90oMCme
タグ：
posted at 15:57:36
OokuboTact 大久保中二病中年 @OokuboTact
仲正氏は無名＆匿名の批判にもいちいち反応する大人気なさ（褒め言葉）が良かったのに、多忙のせいかそういうことをしなくなったのは悲しい。
タグ：
posted at 16:09:30
注）を付けると嫌がられる、売れないという版元が多いというツイがTLに流れているが、そんなことはないぞよ。少なくとも私は、注があって、章末・巻末に参考文献や参照URLがちゃんとあると、なんか得した気分でうれしいぞ。
タグ：
posted at 16:12:44
OokuboTact 大久保中二病中年 @OokuboTact
私は仲正氏ウォッチャーを何年もしているけど、専門よりも副業的な方が無茶苦茶忙しいように思える。
現代思想の啓蒙活動とか俳優業とか。
それはそれで人類への偉大な貢献だけど、やはり以前のように無名＆匿名からの仲正批判への応答に労力と時間を割いて欲しい。
タグ：
posted at 16:14:56
非公開
タグ：
posted at xx:xx:xx
OokuboTact 大久保中二病中年 @OokuboTact
注：「山川ブラザーズ」と言っても知らない人が多いと思うので、簡単に説明すると、サブカルチャー研究者の山川賢一氏が仲正氏を批判して賛同者がかなり出たので、仲正氏がまとめて「山川ブラザーズ」と命名した。
タグ：
posted at 16:18:55
OokuboTact 大久保中二病中年 @OokuboTact
私も山川賢一氏の仲正氏批判には賛同したのだけど、山川賢一の著書は一冊も読んだことはなく、
それに対して仲正氏の本はかなり購入しているので、どちらかと言えば、仲正ブラザーズではある。
仲正ブラザーズというよりは、仲正ウォッチャーあるいは仲正ストーカーの方が正確
タグ：
posted at 16:22:26
#Julia言語 まさにそれ。そういう細かいことをやるのが難しい。やり方を調べ上げること自体が結構面倒。だから、「デフォルトの設定を上書きしてのカスタムプロット！」がもっと公開されてくれないと不便なのだ。
やってみました。続く
twitter.com/tsatie/status/...
タグ： Julia言語
posted at 16:29:07
#Julia言語 まず pyplot() バックエンドから gr() バックエンドに変えただけ。色の設定はデフォルトのまま。gr()バックエンドだと背景の表示が弱すぎ。しかし、動画作成時間は1/3に減る。 pic.twitter.com/g9Wp0IqKKK
タグ： Julia言語
posted at 16:29:07
#Julia言語 gr()バックエンドで色の設定を自前で行ったもの。かなり派手になっている。確かに見易くなったような気がする。 pic.twitter.com/NmYx1rku3m
タグ： Julia言語
posted at 16:29:08
#Julia言語 pyplot()バックエンドで色の設定を自前で行ったもの。plot()の中で
color_palette=(色の巡回リストの配列)
とするだけ。
2つの傾いた円盤のあいだの「相互作用」の様子が見ていて楽しい。
ソースファイル↓
nbviewer.jupyter.org/gist/genkuroki...
Lorenz方程式の数値解のアニメーション pic.twitter.com/VGPOKfLO5z
タグ： Julia言語
posted at 16:35:26
#Julia言語 ううむ、pyplot()版は線が細過ぎるのかも。切りがないのでやめます。その辺をいじって自分のために完璧な動画を作りたい人は
nbviewer.jupyter.org/gist/genkuroki...
の lw=0.2 の部分を lw=0.3 とか lw=0.4 に変えて実行し直してみて下さい。変更後のファイルを私に無断で公開配布して構いません。
タグ： Julia言語
posted at 16:39:50
玉井克哉（Katsuya TAMAI） @tamai1961
よくここまで無責任なことが言えるな。初等中等教育の責任者だったわけだろ。大臣が思い込みで変なこと言ったら、止めるのが仕事。それしないで辞めてから批判するとは。
大臣を止め、止まらなければ職を賭す。辞職して大臣を諌める。それをしない者が、何で局長の高給を食み、次官にまでなるのか。
タグ：
posted at 17:04:34
「「江戸しぐさ」をまだやってるんですか！ あれはやめた方がいいですよてんあれはねえ、大失敗。僕が初等中等局長のとき、下村さん（下村博文・元文科相）に言われて作った。あんなインチキなもの…」他人事感が… / “親を悩ます「PTA問題」 前川喜平さんに聞いた - 前川…” htn.to/2KxrYP6P4g
タグ：
posted at 17:16:14
#Julia言語 の Plots.jl パッケージを使い倒したい人は特に
docs.juliaplots.org/latest/attribu...
Attributes
の熟読が必須。私はブラウザでブックマークしてある。
タグ： Julia言語
posted at 17:20:55
#Julia言語 の Plots.jl ではlegendの表示を半透明にするには、
default(:bglegend, plot_color(default(:bg), 0.5))
default(:fglegend, plot_color(ifelse(isdark(plot_color(default(:bg))), :white, :black), 0.6))
nbviewer.jupyter.org/gist/genkuroki...
タグ： Julia言語
posted at 17:24:30
とんでもニセ歴史の「江戸しぐさ」を文科省が作った道徳教材に入れる仕事をしたのは前川喜平氏(元・文部科学官僚)だった！
なんちゅう迷惑なことを！
webronza.asahi.com/national/artic... pic.twitter.com/RTRzcU4TTe
タグ：
posted at 17:35:04
@genkuroki 読みたくないでしょうが^^、資料として
twitter.com/sekibunnteisuu...
に続くツイートです pic.twitter.com/SVWPCG6SHg
タグ：
posted at 17:44:08
別の記録来た！
後で検索で見つけ易いようにキーワードを並べておく
パルジャミーノ
理解のレベルの差ですねー
小池一夫 子連れ狼 フィクション
移項
みんなに笑われている
大型新人
twitter.com/sekibunnteisuu...
タグ：
posted at 17:52:27
子連れ狼の「ベビーカー」がチート性能を持っていることについては
lonewolfcub.jimdo.com/%E5%AD%90%E9%8...
などを参照せよ(笑)。
沢山の強力な武器が収納されており、拝一刀自身の手で完成させられたガトリングガンまで装備。
底には鉄板。防御にも使える。
船としても利用可能で激流も降れる。
タグ：
posted at 18:00:24
物理クソ苦手マンだからバネの話頑張って読んだが
レッジャーノさんはスリンキーは理想ばね（ばね質量を無視する）とは違うから完全な証明とは限らないだろ 現実ばねと理想ばねのギャップを無視して理論適用できるのかよ証明しろって延々絡んでるんだなこれ
タグ：
posted at 18:06:13
そのレベルの話してないやん、ぶっちゃけ
もちろん理想ばねとは違うが性質をおおまかに捉えた現象をスリンキーで見られるよ！って本文に書いてあんのに
その「おおまか」部分が異常に許せないらしい
タグ：
posted at 18:07:35
で、理想ばね想定理論を現実ばねに当てはめていいのか、に納得行く説明もらえなかったからってそこが面白い部分なのに理解してない時代遅れ大したことない無名研究者だの言いたい放題
性格に難有りだけどこの偏執的な完璧主義は学者向きな気もしますね
タグ：
posted at 18:15:36
非公開
タグ：
posted at xx:xx:xx
高 史明(TAKA, Fumiaki) @Fumiaki_Taka
元文科の前川氏って頭はいいのかもしれんけど本当にしょーもないな…。自分が責任ある地位にいたときにバカな政治家にほいほい従っといて後になって「面従腹背だった」もクソもない。心の中で反感持つだけなら誰でもできるし「面従」したことこそが自分の遺産だろ…。
タグ：
posted at 19:20:42
前川喜平さんが「江戸しぐさ」導入の担当者だったのが意外という人がおられるようだけど『私たちの道徳』作成当時の初等中等局長が前川さんだったことはすでにわかっていたのだからむしろ当然だったわけで、この度のインタビューではそれが本人の証言で裏付けられただけ
タグ：
posted at 19:23:28
高 史明(TAKA, Fumiaki) @Fumiaki_Taka
「面従できなくて1話のBパートあたりで殺されそうな人」という定評のある僕はああいう散々長い物に巻かれておいて後になって「あれは間違っていると思ってました」とか言い出す種類の人間が嫌いなんだ…。
タグ：
posted at 19:40:54
T先生が釈迦やニュートンやアインシュタインではない、は文字通りの意味では良いとして（相対的な差の比喩としては近いと思いますが）、だから自分のほうが物理を理解している、というロジックはすごい
タグ：
posted at 19:45:37
ごまふあざらし(GomahuAzaras @MathSorcerer
こういう問題N回見ていてもう疲れたので berrycondaでいいじゃないかって思ってまキュ。 twitter.com/termoshtt/stat...
タグ：
posted at 20:12:38
権力側にいるときは天下りの斡旋もやってて、追放されると今度は左翼的に受けのいいことを言ってるわけで、要するに保身がうまい口先だけ上手な人で全く信頼できないと思うな。 twitter.com/tamai1961/stat...
タグ：
posted at 20:56:38
blogos.com/article/376481/
なんで役に立つかどうかとか数学の話をしているふうをよそおって小学校の算数みたいな例を出すんだろう？
なんで、高校数学で教えられている（or教えられていた）整数とか三角関数とか微積分とか行列とかベクトルとか、スゲー社会に役立ってるって話をしないんだろう？ pic.twitter.com/4HiDWLfKL6
タグ：
posted at 21:13:09
面の皮をぶ厚くして「面従腹背」でのし上がっていかないと、ああいう人は出来上がらないってことなのかも。 twitter.com/tamai1961/stat...
タグ：
posted at 21:29:41
Kuri TsukaΦ (くり つかお) @20_0855369231
本日の演習授業にて
誰か一人くらいは「こんなん10/49に決まってますよね」と即答してほしかったのですが残念ながらそれはおらず…（条件付き確率の定義に従ってゴリゴリ計算ばかり笑）
いくつか布石は打ったつもりだったんですが難しい pic.twitter.com/vJFqeIrgxG
タグ：
posted at 21:43:38
前川の悪事は天下り斡旋に関わったことだけではなかったということです。この人物を「識者」扱いするのは愚かな行為ですよ twitter.com/genkuroki/stat...
タグ：
posted at 21:56:28
@kankichi57301 @sekibunnteisuu もう面倒なので私は期待の新人さんをブロックしました。
件の新人さんとやりとりする場合には返信先から私を外して下さい。
タグ：
posted at 22:45:16
ほぼ半々に分かれている点が何とも。
アンケート問題を読まずに無作為に答えたのと区別がつかない(笑)。
twitter.com/musicisthebest...
タグ：
posted at 22:48:14
噂のパルミジャーノなる人物のツイートをぱらぱらと眺めてみたのだけど、ようわからん。フォロワーを増やしたくて炎上を狙ってる高校生くらいにしか思えん。炎上より、プロフィールを書くほうがいいんじゃないか？
タグ：
posted at 23:25:47
#数楽 #Julia言語
nbviewer.jupyter.org/gist/genkuroki...
非一様重力下で落下するバネの偏微分方程式のFFTを使った数値解法
u_{tt} = u_{xx} - g(u) 型の偏微分方程式(一般には非線形)を数値的に解いてみました。普通に離散化すると質点をバネで繋いだ場合と同じになりそうなので、FFTを使った別の方法を採用。 pic.twitter.com/kQAiVKXTLF
posted at 23:28:24
#数楽 #Julia言語 2つ前のツイートの動画は u_{tt} = u_{xx} - 1 の解です。これは一様重力(重力加速度1)の場合で、初期条件の x=0 での高さ(初期状態の下端)はしばらく全く動かない。
以下は非一様重力下の場合。次の動画は
u_{tt} = u_{xx} - u
の解です。数値計算が破綻するまで計算している。 pic.twitter.com/axryC0iONy
posted at 23:35:24
#数楽 #Julia言語
u_{tt} = u_{xx} - 1/u
数値計算が破綻する直前まで。 pic.twitter.com/MdZ7owHcU8
posted at 23:39:11
#数楽 #Julia言語
u_{tt} = u_{xx} - 1/u²
これも数値計算が破綻する直前まで pic.twitter.com/ZtfEXyXwn7
posted at 23:39:11
#数楽 #Julia言語
u_{tt} = u_{xx} - sin(π u)
所謂 SIne-Gordon 方程式の場合です。
数値計算が破綻するまでプロット。 pic.twitter.com/AwxHUqkSqp
posted at 23:39:12
#数楽 #Julia言語 これも
u_{tt} = u_{xx} - sin(π u)
所謂 SIne-Gordon 方程式の場合。初期状態が 2 sin(π x) です。
段々 sine カーブが歪んで行く。
x軸をもっと長くして別の解を求めることもできるのですが、まだやっていません。
ソースファイル↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/tFbW4zFv6W
posted at 23:41:57
