プログラム＆数理
こんばんは皆様、三頌亭です。今日は以前、朝日選書で出ていたものがちくま学芸文庫に入ったのでご紹介いたしましょう。ノーマン・マクレイの「フォン・ノイマンの生涯」です。フォン・ノイマンの伝記の決定版といっていいのじゃないかと思ってます。この伝…
こんばんは、皆様、三頌亭です。数値的逆ラプラス変換のGaver-Stehfest 法についての作例をあげておきます。係数はn=14に固定してあります。逆変換の例はもっとも簡単な指数関数をやってみました。細野法と2つを今までの記事で説明してきましたが、算法の解…
こんばんは、皆様、三頌亭です。数値的逆ラプラス変換のvbaによるユーザー定義関数のプログラムです。これは以前掲載しました有理関数の作例ですが、複素数計算の関数も併せて貼り付けるだけの形にしました。 Public Type Complex x As Double y As DoubleEn…
こんばんは、皆様、三頌亭です。さらに小ネタ集の蛇足ですｗ。今回使った為替レートの生データ（1分足）ですが、どこからどのようにして持ってきたかを具体的に申し上げておきます。Forexiteというサイトでロシアのサイトです。Meta stockの部分でダウンロー…
こんばんは、皆様、三頌亭です。さて小ネタ集の続きですｗ。というか補足でございます。為替レートの変動（1分足）6カ月分をストレージしていきなり指数分布に当てはめました。こういった時系列変動がどのような分布を取るかについては様々な研究がありまし…
こんばんは、皆様、三頌亭です。これも小ネタ集のひとつです。今回はVBAを使わなくてもできるでしょうｗ。といいますか1種のフィールドワークのようなものです。まず表題通り、為替の分足データをゲットしてきます。分足データはいろいろあるのですが下記の…
こんばんは、皆様、三頌亭です。今日は小ネタ集の一つですｗ。以下の関数は縦に並んだデータを一定量シフトさせて足し合わせていく関数です。何に使うのかよくわからないという方が大半でしょうが、必要な人には自分でいうのもなんですが「便利」ですｗ。 Re…
こんばんは、皆様、三頌亭です。前回はあまり得意でない矩形波を逆変換した例をお見せいたしましたので、今回はベッセル関数を例にとって細野法で数値的逆ラプラス変換行ってエクセルの関数値と比較してみましょう。ところで最近ではフリーの数式処理システ…
こんばんは、皆様、三頌亭です。細野法でオイラー変換の項数と未変換の項数を自由に設定できるようにしたVBAのコードです。これでおそらくラプラス変換した関数部分の記述を変えるだけでオッケーでしょうＷ。これらの数値的逆ラプラス変換法のアルゴリズムが…
こんばんは、皆様、三頌亭です。このまえから解説してます逆ラプラス変換方法ですが、ご参考までに細野法ではない、Gaver-Stehfest 法について簡単にご説明いたします。まあ大体、利用者の立場としては「能書きはいいから早くプログラムをよこせ！」というの…
こんばんは、皆様、三頌亭です。今日は数値的な逆ラプラス変換をエクセルに関数として組み込むプログラムを紹介いたします。もともとこのプログラムは実験データを拡散方程式の解に当てはめるために作成したものです。もうだいぶ前に作成したものですが、最…
Public Type Complex x As Double y As DoubleEnd Type Function ToComplex(x As Double, y As Double) As Complex ToComplex.x = x ToComplex.y = yEnd Function Function Cadd(z1 As Complex, z2 As Complex) As Complex Cadd.x = z1.x + z2.x Cadd.y = z1…
Rem 数値的逆ラプラス変換filt06(t) Hosonoのアルゴリズム----------------------Rem atstp:分母 btstp:分子---ラプラス関数値の算出に必要なパラメータをセル参照Rem 作例 F(s) = ht(8次多項式)/gt(8次多項式) : 有理関数Public Function filt06(time As Do…
こんばんは、皆様、三頌亭です。これまたプログラム関係のお話。ちょっと画像解析をやる必要があったりするのでPythonを勉強してみたいとよく思います。大した解析ではないので今はほかのプログラムにて代替してますｗ。ところでPython関連の数理モデルのペ…
こんばんは、皆様、三頌亭です。珍しくリアルお仕事関連のプログラムのお話です。多群（3群以上）の平均値を比較する統計手法に多重比較法という手法があります。もう３０年以上前のお話になりますが、当時あまり使われていなくて適当なパッケージがすぐ手に…
