複素力学系で有名な図にマンデルブロ集合がありますが、同じ漸化式を異なる条件で計算するとジュリア集合を計算する事が出来ます。ここでは、Pythonによるジュリア集合を描画し、さらに変数を連続的に変化させた場合で動画にする方法を紹介します。
続きを読むPythonでジュリア集合を描画！条件変更で動画にしてみた
工学
Pythonで描くマンデルブロ集合!フラクタルの旅を体感してみる
複素力学系の範囲に大変不思議な図形を描画する「マンデルブロ集合」というものがあります。ここではPythonを使って複素数の漸化式をシミュレートし、マンデルブロ集合を実際にプロットしてみました。
続きを読むPython/sympyとnumpyで書くニュートン-ラフソン法
方程式系を数値解析的に解く方法としてニュートン-ラフソン法があります。有名な手法で記事も沢山ありますが、ここではPythonのsympyを使って記号的に本手法を書いてみます。
続きを読むsympyで定義した2変数関数をマクローリン展開してプロットする
テイラー展開やマクローリン展開は1変数関数では有名ですが、2変数関数の場合でもよく使います。ここではsympyを使って記号ベースで2変数関数をマクローリン展開し、matplotlibによるプロットをする所までをやってみます。
続きを読むPython/sympyでテイラー展開した結果をグラフ化する方法
テイラー展開は任意点周りの近似計算をすることが出来るため、様々な工学的・物理的場面で活躍します。ここではsympyを使って面倒なテイラー展開をPythonにやらせ、グラフプロットまで自動で行うことを目標とします。
続きを読むPython/pydubでmp3をwavファイルに変換する方法
mp3ファイルは様々なアプリで作成することができるファイルフォーマット形式ですが、信号処理ではwavファイルの方が汎用性があります。ここではPythonのpydubでmp3をwavに変換する方法を紹介します。
続きを読むNumpyだけで書いたガウスの消去法で連立1次方程式を解いてみた
ガウスの消去法は掃き出し法とも呼ばれ、機械的に連立1次方程式を解く方法で有名です。Pythonには様々な便利ライブラリがありますが、ここではNumpyだけを使ってガウスの消去法を書き、アルゴリズムの理解をすることを目標とします。
続きを読むPythonの4次ルンゲ・クッタ法で多自由度連成振動を解く方法
多数の質点系から成る多自由度系は各点の振動が影響し合う連成振動をします。連成振動を解く方法は色々ありますが、ここでは有名な4次のルンゲ・クッタ法をPythonで作成して解いてみます。
続きを読むPythonの過渡応答解析で力を滑らかにかけて応答の違いを見る
振動の時刻歴応答シミュレーションを行う場合、加振の信号は多大に影響します。今回は滑らかに振幅増加する信号を使った場合とそうでない場合で、応答波形がどう変化するのかをPythonで確認していきます。
続きを読むPythonで滑らかに振幅増加する正弦波を作る1つの簡単な例
実験やシミュレーションには理想的な信号を用いることがよくあります。ここでは1つの例として、滑らかに振幅が増加する正弦波をPythonで生成する方法を紹介します。
続きを読む12