ON-OFF調整の動作をエクセルでシミュレーションしたけれど、もちろん、適当なプログラム言語を使っても同じ事ができる。
それを示すため、十進basicを使ってプログラムを書いてみた。何で十進basicかというと、知りあいの勤めていた高校の課題研究で生徒さんが使っていたのがこれだったから。無料で入手できるbasicだし、日本語の簡易マニュアルがあるので、Scilabなんかよりは高校生さんあたりだと手が出しやすいだろうと思う。 ON-OFFのプログラムは LET StageT=20 LET RoomT=20 LET TargT=100 LET HPower=160 LET HNKeisu=1 LET NetuYoryo=160 OPEN #1:NAME "C:/TempD.txt" FOR i=1 TO 1000 IF StageT <= TargT THEN LET StageT=StageT+(HPower-(StageT-RoomT)*HNKeisu)/NetuYoryo ELSE LET StageT=StageT-(StageT-RoomT)*HNKeisu/NetuYoryo END IF PRINT #1: StageT NEXT i CLOSE #1 END な感じ。 最初のLETの部分は変数の定義。StageTがステージの温度。RoomTが室温で一定。TargTが目標温度でこれも一定。HPowerがヒーターパワー、HNKeisuが放熱係数、NetuYoryoが熱容量である。 次のOPEN文は、結果を書き込むためのファイルを作る命令。以下、PRINT文とCLOSE文はファイルに書き込むのとファイルを閉じる操作である。 というわけでプログラム本体はFOR i=からNEXT i までの間である。やっていることは IF StageT <= TargT THEN でステージの温度が目標温度以下であるかを調べて、以下なら LET StageT=StageT+(HPower-(StageT-RoomT)*HNKeisu)/NetuYoryo で、新しいステージ温度として、現在のステージ温度に加熱要素と放熱要素を加える。 もし、ステージ温度が目標温度より高いなら LET StageT=StageT-(StageT-RoomT)*HNKeisu/NetuYoryo で、新しいステージ温度として、現在のステージ温度から放熱要素文のみを引く。 そして、これを1000回繰り返すというものだる。 エクセルファイルと比べると、少しばかり取っつきにくいかもしれない。また、結果をグラフにするのには、できたファイルをエクセルなどで読ませてグラフ化する必要がある。でも、その一方で、ファイルサイズは遙かに小さいし、また、エクセルだとゲロゲロになるような、面妖な処理も比較的楽に実装できる。 コンピュータは、ソフトも含めて道具である。ので、目的のためにとりあえずは使えて、少ない手間で結果が得られるものを使えば良いのだろうと思う。ただ、そのためには、ヒマなときに、いろんなことを試して、それぞれのソフトの得意範囲などを実感しておかないといけなかったりするのだけれど。 ところで、エクセルでは、繰り返しの計算はしていないように見えるかもしれないけれど、 E3のセルの式を等号を使って書き表すと E3=E2+$B$5*(F2*$B$15-G2*$B$11)/$B$9 と一つ前の代のセルの値に足す格好になっており、上の LET StageT=StageT+(HPower-(StageT-RoomT)*HNKeisu)/NetuYoryo とまったくパラレルになっていることが分かる。もうすこしパラレル度を上げたかったら、配列を使って十進basicのプログラムを書き直して LET StageT(i)=StageT(i-1)+(HPower-Tdiff(i-1)*HNKeisu)/NetuYoryo と、明示的に一つ前の世代を使って計算する形になる。プログラムの全体は一番下につけてあるけれど、エクセルの場合は、世代をiという文字で表す代わりに LET StageT(2)=StageT(1)+(HPower-Tdiff(1)*HNKeisu)/NetuYoryo LET StageT(3)=StageT(2)+(HPower-Tdiff(2)*HNKeisu)/NetuYoryo という具合に、具体的な数が入って、繰り返しでなく、下へ下へと続いていっている訳である。エクセルは、上から下へと計算を行っていく。この流れが、プログラムのステップを定めているのである(それ故に、上下を逆にするとまったく動かなくなる) エクセルというと、表やグラフを作るためのソフトと思っている人もいるようだけれど、それは、完全な誤解である。 配列版ON-OFFの十進basicプログラム DIM StageT(1001) DIM TDiff(1001) LET StageT(1)=20 LET TDiff(1)=0 LET RoomT=20 LET TargT=100 LET HPower=160 LET HNKeisu=1 LET NetuYoryo=160 OPEN #1:NAME "C:/TempMD.txt" FOR i=2 TO 1001 IF StageT(i-1) <= TargT THEN LET StageT(i)=StageT(i-1)+(HPower-Tdiff(i-1)*HNKeisu)/NetuYoryo ELSE LET StageT(i)=StageT(i-1)-Tdiff(i-1)*HNKeisu/NetuYoryo END IF LET Tdiff(i)=StageT(i)-RoomT PRINT #1: StageT(i) NEXT i CLOSE #1 END ※DIMは配列の定義命令である
by zam20f2
| 2010-05-04 07:38
| 科学系
|
Comments(0)
|
リンク
お気に入りブログ
カテゴリ
タグ
SCI(841)
道具(405) 組織写真(314) プロジェクトT(216) 顕微鏡(180) Ch(130) 散歩(119) 製図(99) N相(70) キラル(65) のほのほ(57) Sm相(48) 表面・界面(42) プロジェクトLC(40) B/W(30) Iso相(30) 新興織物撰集(21) SmA相(20) 結晶相(12) 転傾(11) 最新のコメント
記事ランキング
以前の記事
2024年 03月
2024年 02月 2024年 01月 2023年 12月 2023年 11月 2023年 10月 2023年 09月 2023年 08月 2023年 07月 2023年 06月 2023年 05月 2023年 04月 2023年 03月 2023年 02月 2023年 01月 2022年 12月 2022年 11月 2022年 10月 2022年 09月 2022年 08月 2022年 07月 2022年 06月 2022年 05月 2022年 04月 2022年 03月 2022年 02月 2022年 01月 2021年 12月 2021年 11月 2021年 10月 2021年 09月 2021年 08月 2021年 07月 2021年 06月 2021年 05月 2021年 04月 2021年 03月 2021年 02月 2021年 01月 2020年 12月 2020年 11月 2020年 10月 2020年 09月 2020年 08月 2020年 07月 2020年 06月 2020年 05月 2020年 04月 2020年 03月 2020年 02月 2020年 01月 2019年 12月 2019年 11月 2019年 10月 2019年 09月 2019年 08月 2019年 07月 2019年 06月 2019年 05月 2019年 04月 2019年 03月 2019年 02月 2019年 01月 2018年 12月 2018年 11月 2018年 10月 2018年 09月 2018年 08月 2018年 07月 2018年 06月 2018年 05月 2018年 04月 2018年 03月 2018年 02月 2018年 01月 2017年 12月 2017年 11月 2017年 10月 2017年 09月 2017年 08月 2017年 07月 2017年 06月 2017年 05月 2017年 04月 2017年 03月 2017年 02月 2017年 01月 2016年 12月 2016年 11月 2016年 10月 2016年 09月 2016年 08月 2016年 07月 2016年 06月 2016年 05月 2016年 04月 2016年 03月 2016年 02月 2016年 01月 2015年 12月 2015年 11月 2015年 10月 2015年 09月 2015年 08月 2015年 07月 2015年 06月 2015年 05月 2015年 04月 2015年 03月 2015年 02月 2015年 01月 2014年 12月 2014年 11月 2014年 10月 2014年 09月 2014年 08月 2014年 07月 2014年 06月 2014年 05月 2014年 04月 2014年 03月 2014年 02月 2014年 01月 2013年 12月 2013年 11月 2013年 10月 2013年 09月 2013年 08月 2013年 07月 2013年 06月 2013年 05月 2013年 04月 2013年 03月 2013年 02月 2013年 01月 2012年 12月 2012年 11月 2012年 10月 2012年 09月 2012年 08月 2012年 07月 2012年 06月 2012年 05月 2012年 04月 2012年 03月 2012年 02月 2012年 01月 2011年 12月 2011年 11月 2011年 10月 2011年 09月 2011年 08月 2011年 07月 2011年 06月 2011年 05月 2011年 04月 2011年 03月 2011年 02月 2011年 01月 2010年 12月 2010年 11月 2010年 10月 2010年 09月 2010年 08月 2010年 07月 2010年 06月 2010年 05月 2010年 04月 2010年 03月 2010年 02月 2010年 01月 2009年 12月 2009年 11月 2009年 10月 2009年 09月 2009年 08月 2009年 07月 2009年 06月 2009年 05月 2009年 04月 2009年 03月 2009年 02月 2009年 01月 2008年 12月 2008年 11月 2008年 10月 2008年 09月 2008年 08月 2008年 07月 |
ファン申請 |
||