エクセルマクロVBAを使ってパラメータの変化をリアルタイムで確認する

資料請求番号:SH43 SH44 TS41

微分方程式パラメータが解曲線に与える影響を考察するために

大学や会社で研究開発活動をしていると、微分方程式における様々なパラメータが変化すると解曲線はどのように変化するのかをストレスなく確認したいという場面があると思います。
また、パラメータを求める実験を行う際には、ターゲットとしている微分方程式がパラメータによってどのように変動するのかを予め知っておけば、実験計画をスムーズに立てられると思います。

今回は、微分方程式のパラメータの変化が解曲線にどのような影響を与えるのかを理解するためのアニメーションの作成方法をまとめたいと思います。

今回作成したいアニメーションは以下の通りです。こちらは振動現象の微分方程式

をルンゲ・クッタ法で解いた曲線で、「γ=0.1の粘性によって減衰してしまうところ、どのくらいの強制振動F0を与えると振動を安定に続けることができるのか」を考察しています。

こちらのページは振動に関する上記の微分方程式をルンゲクッタ法で解くプログラムがあることを前提にしています。ルンゲクッタ法のプログラムは以下の記事で解説しています。

※VBAソースコード 振動

はじめに

ストーク
・・・ああでもねぇな・・・こうでもねぇな・・・。いくつが適切なんだ・・・・それとも、このモデルじゃ、実験結果は再現できねぇのか・・・・?
シママ
どうしたの?なんか、しんどそうね。
ストーク
微分方程式の解って、パラメータが変われば変わるやろ?どのくらいが適切なのか、試行錯誤して何度もルンゲクッタ回してて、途中でなんか疲れてきて、気づいたらさっき試したパラメータと同じパラメータ試したりしてて、アタマん中がしっちゃかめっちゃかで、どうもスッキリせんたい。
シママ
あ~そっか・・・。そうなのね・・・。
ストーク
てかオマエ、研究か仕事で微分方程式解いて現象を表現したり、現象から微分方程式立てたりしたこと、ねぇやろ?
シママ

なんでそんなケンカ腰なのよ!?イライラしてるからってワタシに当たらないでちょうだい!!

シママ

まぁ、確かにないけど。ワタシはプログラム書くことが仕事だら?一応微分方程式を解くプログラムは書けるけど、その微分方程式の中身はよく分かんない~なんてよくあることだし?
ただね、例えば・・・この微分方程式のF0ってパラメータあるでしょ?

ストーク

ああ。

シママ

これを0.1ずつ変動させながら、繰り返しグラフを描画させられる方法なら知ってるわよ。

ストーク

おおっ、ってことは、いちいち0.1ずつ変化させてルンゲクッタループしなくてもええってこと?

シママ

まぁね。1回アニメーション見て、気になったパラメータで刻み幅を小さくするなり、他のパラメータ考察するなりすれば、効率よさそうじゃない?

ストーク

せやな!教えてくれよ!

シママ

お腹すいたな~。今日、テンイチ行きたいな~

ストーク

・・・しょうがねぇな。ええよ。行こう。

シママ

やった!

説明

シママ

それで、今回やりたいことっていうのは、例えばF0を0.5から4まで動かしたときの解曲線の挙動が見たいっていう目的なら、
0.5でルンゲクッタループ回して、グラフ見て、0.6でルンゲクッタループ回して・・・を続けて4までたどり着くようにすればいいのよね。

ストーク

せやな。

シママ

なら、基本形として、C7セルにあるF0の値を代入しては計算させてプロットさせての繰り返しになるわね。

/* –vibration constant simulation——– */

For force = 0.5 To 1 Step 0.1
Cells(7, 3) = force
Application.Run “rkvib”
Next force

/* –Rose curve Realtime plot ——End– */

ストーク

このApplication.Run “rkvib”ってのは・・・ルンゲクッタループか?

※Sub rkvib()
シママ

そうそう!0.5でルンゲクッタループ回して、Step0.1だから次は0.6で回して・・ってこと。

ストーク

ほう。繰り返しの制御変数に整数型以外も使えるのか。便利やな。

シママ

まぁ、一応、テストとして、とりあえず1までにしてやってみたら?

ストーク

わかった。

ストーク

・・・・。

ストーク

数字が変化するだけで、リアルタイムでグラフが変化せんのやけど。
もうテンイチ行かんぞ!!

シママ

まぁまぁ待って待って。
あのね、VBAのプログラム回している間はExcelの制御を完全にVBAに渡している状態になるの。
だから、時間のかかるマクロを回している間はExcelがフリーズしたみたいな感じになるのよね。

ストーク

で?

シママ

それで一旦制御をOSに返しますっていう命令を書かなきゃいけないの。それがDo Events関数。
これを入れてやってみなよ。制御がOSに返ったらグラフを描くことができるから。

完成形(サンプルプログラム)

/* –vibration constant simulation——– */

For force = 0.5 To 1 Step 0.1
Cells(7, 3) = force
Application.Run “rkvib”
DoEvents: DoEvents: DoEvents
Next force

/* –Rose curve Realtime plot ——End– */
※Sub rkvib()

ストーク

おおっ!グラフがうごきよった!!最初からこれ教えろよ~(笑)

シママ

でしょ?スゴイでしょ?そしたら、F0 = 4までまわしてごらん?

ストーク

これでF0=3がちょうどいいってことがわかるし、F0 = 2.1~2.2で重ね合わせ特有の不安定な振動がなくなるのもわかるたいね!

ストーク

いやぁ~。よかよか!ありがとう!

シママ

じゃ、テンイチ行きましょ!

ストーク

はいはい・・・。

2 Comments

フロスト

こんにちは、

貴殿のサイト読ましてもらっています。

とても詳しくわかりやすく説明されていて、VBA初心者の私にとってはとても学びやすいです。

ところで、以下のコードなんですが、これをどこに書けば最後の動画のようになるのでしょうか?

rkvib内のどこかなのかモジュール(Module1 やModule2)内に導入すればいいのでしょうか?

For force = 0.5 To 1 Step 0.1
Cells(7, 3) = force
Application.Run “rkvib”
DoEvents: DoEvents: DoEvents
Next force

返信する
shimakei8364

こんにちは。
質問の回答ですが、同一モジュール内にrkvibとは別の新しいSubプロシージャを作成して
For force = 0.5 To 1 Step 0.1
Cells(7, 3) = force
Application.Run “rkvib”
DoEvents: DoEvents: DoEvents
Next force
の構文を書けばよろしいかと思います。

返信する

コメントを残す

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です