lapinv_fx
例:ラプラス関数 (π/(4*s^3))^0.5*e^(-1/(4*s))(原関数 = sin(sqrt(x)))の逆変換
- 機能:
- 任意式F(s)のラプラス逆変換
- 書式:
- lapinv_fx(F(s),[ot],[結果],t始,t終,t増)
- 解説:
- FILT法に基づき、ラプラス関数F(s)の逆変換を行い、指定区間[t始,t終]と刻み値での逆変換値(原関数値)を数値組[ot],[結果]に出力します。
原関数が不連続である場所付近(例えば、ステップ関数の立ち上がるところは不連続)を除けば、本関数の逆変換結果の精度は充分に高いです。
F(s) ・・・・・・・・・・・・・・・ ラプラス関数の複素数型式表現
[ot] ・・・・・・・・・・・・・・・ 結果のt値の出力数値組
[結果] ・・・・・・・・・・・・・・・ 逆変換結果の出力数値組
t始,t終 ・・・・・・・・・・・・・・・ 逆変換の出力区間(t始>0)
t増 ・・・・・・・・・・・・・・・ 正:刻み値、負整数:均分点数
lapinv_fx(sqrt(PI/(4*s^3))*exp(-1/(4*s)),[1],[2],0.2,10,0.2); → 最大計算誤差 = 0.000000298
return |