カテゴリー: 日記

速度に比例する抵抗を受けながらの落下について、Pythonを使って、Euler法、odeint、Runge-Kutta法の3種類の計算をやってみた。グラフは比較用の解析解を加えた4つをプロットしているけど、一目で分かる違いなんてない。ちなみにグラフは上が変位、下が速度である。まぁ、終端速度がありますね、というくらいの話である。

しばらくプログラミングは何もしていなかったので、タイピングがどこかぎこちなくて、打ち間違いが滅茶苦茶多かったが、ようやくマシになってきた。

仕事で使うためにのPythonコードを書いていて、大きな勘違いをしていることに気がついた。csvファイルを読み込んで、各行ごとに文字列中の全角かなから半角英数に切り替わる箇所を検出しようとしていたが、対象となるファイルを確認したら、その場所にはタブが入っていたのである。これならタブをセパレータとしてsplitを使えば良いだけの話である。すげー間抜けだ…

なかなか恥ずかしい勘違い（見落としというべきか）だが、これが分かれば後にやる作業ははっきりしている。メインの業務ではないので、割ける時間は多くないが、早ければ今週中、遅くとも来週中には結果を出せるだろう。

随分前に購入した、「計算物理I」（夏目雄平・小川建吾、朝倉書店）を参考に、4次のRunge-Kutta公式の導出を行った。卒論で一度導出はしているはずだけど、もはや記憶にない（苦笑）

上記書籍をメインに計算を追い掛けて、1階の常微分方程式の簡単な問題について、解析解との比較までは上手く行ったが、理解の程度はまだまだなので、それ以上は繰り返し読み込むしかないだろう。Simpson公式と一致するところまでは、しっかり理解したいところである。

その後は、2階の微分方程式への適用である。「数値計算による量子力学」（桜井捷海、裳華房）を参考に、もう少し考える必要がある。ブランク空きすぎだ（苦笑）

偉い人が、自社の技術者のレベルが下がっていることを嘆いていた。一つは過去の経験に固執して新たな挑戦をしないこと、もう一つは、数学力がないことである。

数学力についての話を聞いていると、本質は、数学に限らず自分の頭で考えていないことのように感じられた。過去の経験に基づいて仕事をすることも、楽をしたいという気持ちの現れだと思う。

帰宅後、速度に比例する抵抗を受けながらの落下について、紙と鉛筆で計算をした。符号の向きで考え込んだりと、初歩的なところで躓いたりしたが、一通り計算して、matplotlibで解析解をプロットするところまではできた。明日以降は複数の方法で数値計算を行って、それぞれの方法を比較しようと思う。

Pythonでプログラミングをしていると、取り敢えず本に書いてあるコードをEmacsで打ち込んでも意味はなく、真面目に理論を追い掛けたり、odeintなどの使い方を理解しないと何も残らないのだと、改めて思う。こんなことは敢えて書くまでもないかも知れないが。

小さなことからコツコツと、である。

ここ数日、平面内の荷電粒子の運動を計算をするためにPythonでコードを書いていたのだが、明かに挙動がおかしい。参考にしている書籍と比較すると、一箇所符号が逆になっているようなのだが、その理由が分からなくて困っていた。

何か符号に関する勘違いをしているのだろうと思い、上記書籍に繰り返して目を通していたら、最初のところに「質点の持っている電荷は−1であるとします」って書いてある…電荷は+1だと勘違いしてた…俺あるあるだ。

その一方で、リストやタプルを使ったfor文の扱いは参考になった。最初は愚直にインデックスをつけたりしていたのだが、そうしなくてもできることなのね。自分で考えて書いたコードよりも格段にすっきりしている。

まぁ、こうした方が知識が身に付くと考えよう。

“Don’t think. Feel!”は「燃えよドラゴン」に出てくる有名なセリフだが、”Don’t think. Move!”というのはその派生なのだろうか。いずれにしても、「ストレスフリー超大全」（樺沢紫苑）の「不安を行動で取り除く」というのはもっともな指摘である。

「Pythonによる数値計算とシミュレーション」と「なんでもPythonプログラミング」を取っ替え引っ替え、たまに「独習Python」を読みながら、放物運動のグラフを描くところまでやってみた。敢えて本に書いてあることから逸れてみることで、理解が深まったところもある。とはいえ、独立して「趣味れーしょん」に載せるレベルには程遠いので、ここで計算結果を上げるに留めておく。45度方向に投げるのが一番遠くまで飛ぶよねという、ありきたりな結果でしかない。出射角は、青が35度、赤が45度、緑が55度である。

「趣味れーしょん」の方は、10年間未更新になることを回避できれば良いと思っているが、あと2年ほどあるので、それは大丈夫だろう…多分（笑）

ついでと言っては何だが、これまでMacのメモを使っていたここの下書きを、nDiaryを使っていた頃のようにEmacsに戻した。プログラミングはEmacsですることになるので、プログラミングの時間が増えるのであれば、エディタはEmacsをメインにした方が都合が良い。

ここ数日、Pythonで文字列の処理をしようとあれこれやっている。文字列中に、全角かなから半角英数に切り替わる箇所が複数あって、そこで文字列を分割したいのだが、今のところ、切り替わる箇所を検出して、そこにセパレータになる文字を挿入するところまでは、できている。その後にsplitを使えば良いだろう。

次の問題は、今は上記の場所が1つだけある場合にしか対応していないので、それを複数箇所対応にすることだが、もう数日かかりそうだ。その後は、csvファイルを行毎に読み込んで、そうして得られるリストの特定の要素に、上記の処理を行うことが最終目標である。

10年くらい真面目にプログラミングをしていないので、この程度の処理でも結構大変だし、自宅ではEmacs、職場ではxyzzyの操作が少々ぎこちない。しかしながら、色々と心に火がついた（ような気がする）し、仕事で使うためのものなので、もう一踏ん張りである。

HDD2個とEye TV HDを、小型家電のリサイクルBOXに出してきた。片道2kmの散歩であった。

HDDの1つはTime Machineで使っていたのだが、先日動作中に異音を発したので、新しいHDDと入れ替えて、お役御免である。多分、5年以上使っていたと思う。もう一つはLinux Boxで使っていたものなので、Time Machineで使っていたものよりも古い。

Eye TV HDは、数年前にDTCP-IP対応のソフトを導入したことで不要になったが、取り敢えずSTBとMacに接続した状態で放置されていた。Time Machine用のHDDの隣に置いてあったので、丁度良い機会だと思い、処分した次第である。

また、やんごとなき事情でモニターのHDMIポートが一杯になりそうなので、Mac miniをDisplayportで接続するため、USB-C/Displayport変換ケーブルを注文した。

一本ケーブルがあったと思ったて確認したら、mini-Displayport/Displayport変換ケーブルだった。Mac miniとモニター間のHDMI接続が不調の頃に買ったのだが、変換ケーブルを買ったら不具合が収まったので放置したため、間違った買い物をしていたことに今頃気がついたという、間抜けな顛末である。

実際に使ったかどうかは記憶が定かではないのだが、多重ループを一気に抜ける場合の方法は、gotoを使うことしか知らなかった。

ところが、である。Pythonの本を読んでいたら、elseとcontinueを使う方法と、フラグ変数を使う方法の2つが出ていて、軽く衝撃を受けた。前者はelseだけを使うという意味で、後者はC言語でもgotoより一般的な方法らしいのに、今まで知らなかったことに…

年末年始で仕事関係の本を一冊読み切った。頭の整理に役立ったと思う。

Python、物理、深層学習、データ分析、ゲノム編集といった本にも手を伸ばしたが、いずれもつまみ食い程度でしかない。こういうことにもっと時間を使えたらなぁと思うのだが、現実を見なければいけないだろう（苦笑）