分数階微分(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辛を食す。かなりコクがあって美味しい。 »
「日記・コラム・つぶやき」カテゴリの記事
- Qwen3.6-35B-A3Bが発表され、Ollamaでも使える。そこで電子レンジの動作原理(2.45GHzは水分子の共振周波数でない)と隕石が大気圏突入で燃える原理(摩擦熱ではない)を聞くと、誘電緩和と断熱圧縮について正しく答えられた。今までのローカルLLMで一番賢い回答と思う。(2026.04.17)
- 高周波・RFニュース 2026年4月17日 atisの3GPP Rel.20ウェビナー動画公開、MWCバルセロナ2026でのGSMA Device Enablement Summit資料公開、ハリファ大学が無線周波数AI言語モデルRF-GPT発表、レドームの解説など(2026.04.17)
- ExcelのOfficeスクリプト(TypeScript)で数値計算ライブラリmath.jsを使う(1) Officeスクリプトは外部API呼び出せるし、math.jsは RESTful APIで呼び出せることがわかった。まずは選択したセルのデータを読み、行列演算。LU分解で一次方程式を解き、逆行列と行列式を求める。(2026.04.17)
- 高周波・RFニュース 2026年4月16日 AmazonがGlobalstarを買収、GSMAが日本のデジタル化をレポート、Mini-Circuitsがケーブルアセンブリを動画で解説、Kymetaが米国海軍研究局と衛星通信で契約、PerasoがドローンIFF向け60GHzモジュール出荷、SEMCOが1500V耐圧MLCC発表(2026.04.16)
- 高周波・RFニュース 2026年4月15日 Microwave Journalはアンプと発振器特集、Signal Integrity Journalは100GHz越えのインターコネクトのAIを使うHFSSモデル化、ローデ・シュワルツが潜水艦通信をUDT2026で発表、Xiaomi Poco X8 Pro分解動画、atisの5Gポリシーレポート(2026.04.15)
« 京都の猿ラインへ行く(その3) 猿ヶ辻(京都御所) 何度も見ていた光景なのに、なんでこんな形か全然疑問に思ったことがなかった。 | トップページ | ココイチで牛カレー(ぎゅうかれー)+やさい、400gの2辛を食す。かなりコクがあって美味しい。 »





コメント