凄い速さで驚きです。

数値ラプラス逆変換

　　数値ラプラス逆変換 島村1243 2008/07/22 07:33:55  (修正2回)
　　├（２）について 山中和義 2008/07/23 21:22:50 
　　│└約3分の１になりました 島村1243 2008/07/24 19:50:08 
　　├できるだけ演算回数を減らしてみました。 荒田浩二 2008/07/23 23:03:12  (修正1回)
　　│└更に速い 島村1243 2008/07/24 19:54:40 
　　└（１）について 山中和義 2008/07/24 19:50:54 
　　　├凄い速さで驚きです。 島村1243 2008/07/24 20:14:16

Re: （１）について	返事を書くノートメニュー
島村1243 <bjllmpcujp> 2008/07/24 20:14:16
凄い速さで驚きです。山中さん、手間のかかるコードを作成していただき有り難うございます。私のPCで約0.９秒でした。やはりFFTを使う（だけでは無いようですが）と速いのですね。私のプログラムは、元々は、ラプラス変換の物理的イメージを捉える方法を考えている中で気づいたものです。具体的には、電気的歪み波の複素フーリェ係数(An-jBn)から、AnとBnを抜き出して時間関数を復元出来る手法をフーリェ変換に拡張利用し、更にラプラス変換に拡張利用出来れば、ラプラス変換のAw、Bwは、複素フーリェ係数のAn、Bnと同じような物理的イメージとして捉えることが出来ると考え付いたものです。十進BASICを知る前から、もしかしたら複素フーリェ係数(An-jBn)と同じ様な考え方でフーリェ変換F(jw)やラプラス変換F(s)も扱えるのかな？というぼんやりした物理イメージは有ったのですが、テストする手段が無くて具体的な式を考えるには至りませんでした。しかし十進BASICを知って、複素数を扱う計算が簡単なので、今回作成を試したものです。結果的には、数値ラプラス逆変換を行っていると言う事をウェブで知りました。私のプログラムは物理的なイメージを分かり易く示していると思うのですが、処理精度と速度を真の数値逆変換プログラムを作る事が出来ないので比べる事が出来ずにいました。ご教示頂いた、真の数値逆変換プログラムを、これから勉強したいと思います。お忙しい中、本当に有難うございました。

　　　└Ganmaと誤差 島村1243 2008/07/25 09:14:17 
　　　　└精度を上げるのは難しいようです。 山中和義 2008/07/25 14:53:36 
　　　　　└良く分かりました 島村1243 2008/07/25 21:02:36

インデックスへ EXIT
新規発言を反映させるにはブラウザの更新ボタンを押してください。