手慰み

マンデルブロ集合を計算してみる3

前々回:マンデルブロ集合を計算してみる - 超ウィザード級ハッカーのたのしみ 前回:マンデルブロ集合を計算してみる2 - 超ウィザード級ハッカーのたのしみ マンデルブロ集合は以下示される複素数列の挙動について表したものである。 そこで、この複素数列…

マンデルブロ集合を計算してみる2

前回:マンデルブロ集合を計算してみる - 超ウィザード級ハッカーのたのしみ 通常のマンデルブロ集合は、 でm=2として、n→∞のときにz_nが発散しない複素数cの集合です。 mを変えてみたらどうなるのでしょうか?試してみました。 m=2の場合: m=3の場合: m=4…

マンデルブロ集合を計算してみる

マンデルブロ集合を計算してみた。 マンデルブロ集合とは で表される漸化式でn→∞としたときにz_nが発散しない複素数cの集合である。 コード マンデルブロ集合を計算してPPM形式の画像を出力します。 結果 画像は複素平面を表しており、横が実数軸、縦が虚数…

ターミナル上でセル・オートマトン

ウルフラムの一次元セル・オートマトンをターミナル上で計算するスクリプトを手慰みに作った。定期的にセル・オートマトンとかで遊びたくなる。 ターミナルなので幅が少ないから、あまり見た目的に面白くはない。作業に飽きた時に和むためのスクリーンセーバ…

水面波の数値解析における安定性について

どこかのだれか(たぶん世界に2〜3人ぐらい)にはきっと理解されて役に立つであろうこと。闇に葬られるのも耐え難いので、英語にして公開する。 Stability of numerical analysis about free surface flow

ニューラルネットワークで囲碁の終局を判定2

前回 ニューラルネットワークで囲碁の終局を判定 - 超ウィザード級ハッカーのたのしみ 何パターンか試してみた。 入力値は、終局図から黒石を取り除いたものと白石を取り除いたものを並べたものである。石がある点を1.0、石がない点を0.0としている。黒石が1…

ニューラルネットワークで囲碁の終局を判定

ニューラルネットワークで囲碁の終局を判定する.終局したかどうかではなくて,終局状態でどこが自分の地かどうかを判定する. これだけなら,わざわざニューラルネットワークを使わずに,手続き的にプログラミングしたほうがいい.しかし,これが出来ないと…

数学定数でブラウン運動 2 - 1億桁で計算

前回 数学定数でブラウン運動 - 超ウィザード級ハッカーのたのしみ 1億ステップでやってみた. Mathematical Constants - Millions of Digits から数学定数のデータはもらった. メルセンヌツイスター 2の平方根 3の平方根 自然対数の底 円周率 に加えて,…

数学定数でブラウン運動

擬似乱数を可視化するためにブラウン運動を計算している. 今回は,数学定数の各桁を外力としてみる. 数学定数のn桁目をdとしたときに,時間n-1からnまでに粒子に与える外力を (F_x, F_y) = (cos(2 * pi * d / 10), sin(2 * pi * d / 10) とする. 試した数…

ブラウン運動をシミュレーションしてみる

題名の通り.擬似乱数を可視化するために,ランジュバン方程式に基づくブラウン運動を計算してみた.*1 結果の絵. コードは以下. import math import random class Particle(object): ''' A particle moving in 2D space. ''' position_x = 0.0 position_y…