分数階微分(Fractional Calculus)の計算をPythonのdifferint(Riemann-Liouvilleの定義で計算するもの)でやってみる。
たまたま非整数階微分のことを調べていて、
を見ていた。昔、非整数階の拡散方程式(∂u/∂t = D∂^α u / ∂u^α )を計算したことがありますが、
その時とは違うRiemann-Liouvilleの定義で計算しようと思った。
こういうのは絶対Pythonでやっている人がいるだろう、と思ったら、やっぱりあった。こちら。differint。
https://pypi.org/project/differint/
では早速やってみる。
f(x)=x
f(x)=sin(x)
f(x)=exp(-x^2)
なるほど直感的にもそんな感じだ。なぜこれを調べていたかというと(続く)。
« 京都の猿ラインへ行く(その3) 猿ヶ辻(京都御所) 何度も見ていた光景なのに、なんでこんな形か全然疑問に思ったことがなかった。 | トップページ | ココイチで牛カレー(ぎゅうかれー)+やさい、400gの2辛を食す。かなりコクがあって美味しい。 »
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース2025年5月23日 HUBER+SUHNERが76-81GHzのミリ波レーダ向け3D waveguide antenna発表、Silicon LabsがIoT向けシリーズ3 SoC発表、GSMAがM360ユーラシアでAIと5Gのイノベーションと協業について発表、ロームがAIサーバー向けMOSFET発表(2025.05.23)
- 高周波・RFニュース 2025年5月22日 三星電機(SEMCO)が165℃対応の車載インダクタ発表、KYOCERA AVXがリップル電流についての技術文書発行、QualcommとXiaomiの契約15年目、OmdiaがNokiaをPrivate 5Gの2025の王者と決定(2025.05.22)
- Google ColabのJulia言語でFPUT問題(Fermi–Pasta–Ulam–Tsingou、非線形結合した振動子が最初に与えたモードに戻る再帰現象)をDifferentialEquations.jlの2階の常微分方程式ソルバーDynamicalODEProblemでシンプレクティック8次のKahanLi8で計算、振動子の動きも動画にしてみる。(2025.05.22)
- 高周波・RFニュース 2025年5月21日 TDKが0201のRFインダクタ発表、InfineonがUWBのFiraコンソーシアムの理事会に、ubloxがロボット用GNSSモジュール発表、FibocomがMediaTekのT930を使った5Gモジュール発表、Motolora Edge 60 Pro分解動画(2025.05.21)
« 京都の猿ラインへ行く(その3) 猿ヶ辻(京都御所) 何度も見ていた光景なのに、なんでこんな形か全然疑問に思ったことがなかった。 | トップページ | ココイチで牛カレー(ぎゅうかれー)+やさい、400gの2辛を食す。かなりコクがあって美味しい。 »
コメント