初めて訪れる方は「ごあいさつ」からご覧ください。
------------------------------------------------------------
目次
------------------------------------------------------------
ごあいさつ
第１回 ｚ＝ｆ(ｘ，ｙ)のグラフ／準備編
第２回 ｚ＝ｆ(ｘ，ｙ)のグラフ／プログラミング前編
第３回 ｚ＝ｆ(ｘ，ｙ)のグラフ／プログラミング後編
ごあいさつ(2021/04/18掲載)(*1)
このページは「超小型記号言語」を使った数値計算とゲームについて書いてあります。
非力な整数型プログラミング言語でどこまでできるか、おおよそ４０年前のテーマの再挑戦です。
このページは、他のページが忙しくなると放置されるということを繰り返してきました。(他のページも似たようなものだけど)
気長に見守ってやってください。
なお「超小型記号言語」の開発段階の記事は無かったことに(削除)してありますので、悪しからず。
--------------------
(*1)元々タイトルは無かったのですが「こあいさつ」とタイトルを付けました。
掲載日はタイトルを付けた日付にしてあります。
■第１回 ｚ＝ｆ(ｘ，ｙ)のグラフ／準備編(2018/12/09掲載,2020/04/12最終変更)
最初のテーマは、z=f(x,y)のグラフです。
このテーマは、パソコン黎明期にグラフィック機能をアピールするにはうってつけでした。
この記事も、当時のパソコン雑誌を参考にしています。(*1)
ただ、超小型記号言語にはグラフィック機能がありません。
そこで、コマンドプロンプト画面に"*"を表示することでグラフィック機能を代用することとします。
幸いなことに、現在のパソコンではコマンドプロンプト上で400×150以上の文字表示が可能で、当時のパソコンのグラフィック機能と比べ遜色ありません。
さて、今回は３次元空間に、水平方向には直交するx軸、ｙ軸、高さ方向にはｚ軸があるものとします。
この場合、関数z=f(x,y)は、任意のx,yに対してz軸(高さ)方向にf(x,y)だけ上方向にプロットされます。
このプロットされた点の集まりが、関数z=f(x,y)で表された曲面ということになります。
非常に単純な例ですが、f(x,y)=1であれば、高さ1のところにある水平な曲面(平面)を表します。
次に、２次元であるパソコンのディスプレイ上では、縦方向(下向き)にX軸、水平方向(右向き)にY軸があることとします。
このパソコンのディスプレイ(２次元)に立体的なz=f(x,y)の曲面を表示させるには、曲面をある方向から見た場合に限定してディスプレイに表示 するしかありませ ん。
ここで、任意の方向から曲面を見た時のディスプレイ表示を考えると、たぶん行列とか三角関数が登場してきます。
そっち方面(ＣＧ)は私の専門外ですし、どうせ超小型記号言語にはそんな関数(行列とか三角関数)はありません。
ということで、見る方向はある特定の斜め方向に限定することとします。
今回の話は、単純化すると３次元空間の任意の点(x,y,z)を、ある関係式で２次元空間(＝パソコンのディスプレイ)の点(X,Y)に変換することです。
では、その変換式は次の図1-1ではどういう式になるでしょうか。
図1-1
まだ、プログラミグしていないので、１００％の確証はないのですが、たぶん
X = +x/2 -z
Y = +y -x/2
でしょう。(*2)
図をみれば判ります。(数学は不要です)
それなら、３次元空間の曲面f(x,y)は２次元座標表示では
X = +x/2 -f(x,y)
Y = +y -x/2
となるはずです。(*3)
唐突ですが、今回はここまで。
次回は、実はそう簡単にはいかない諸問題について説明します。
--------------------
「エンサイクロペディア・アスキーＶｏｌｕｍｅ３、Ｈ６８ＴＲ ３次元グラフィックス(アスキー出版)」
(*2)(*3)最初はxに係数(1/2)を付けていなかったですが、x軸が真下ではなく斜め左下方向に延びているためx軸の増減は他の軸の半分の効果しか無いと考え、ｘに係数(1/2)を付けました。
■第２回 ｚ＝ｆ(ｘ，ｙ)のグラフ／プログラミング前編(2018/12/24掲載,2019/01/03最終変更)(*2)
前回(第１回)は、関数z=f(x,y)のグラフを表示する仕組みについて説明しました。
理屈としてはこれですべてなのですが、グラフィックライブラリのない超小型記号言語でこれを実現しようとすると、実はそう簡単にはいきません。
今回は、具体的な仕様と設計過程を説明しながら、実はそう簡単にはいかない諸問題の解決方法を示してゆくことにします。
【手順１】画面設計
曲面z=f(x,y)は、縦(x軸)×横(y軸)×高さ(z軸)=100×100×100の立方体の範囲で変化させることにします。
パソコン画面には、この曲面を左斜め上(じゃなくて、右斜め上でした)から見下ろした時の見え方を、縦(S軸)×横(T軸)=50×100の長方形内に表示させます。
なお、前回は画面上の座標を大文字で(X,Y)としていましたが、３次元空間内での座標(x,y,z)と混同しないよう、今回から(S,T)とすることにします。
以上を図2-1に示します。
図2-1
【手順２】３Ｄ⇒２Ｄ変換式
図2-1より、３次元空間内の任意の点(x,y,z)は、パソコン画面上の点(S,T)としては
S = +0.7x -f(x,y)
T = +y -0.7x
となります。
前回は、この式のx係数は0.5だったはずですが、実際に0.5でプログラミングすると理屈通りのグラフとならず、最終的に0.7に変更しました。
0.7の根拠は図2-1でx軸とy軸との関係が直角二等辺三角形の長辺と短辺の関係(1.41:1)にあり、そのおおまかな補正(1.41×0.7=0.987)ということで納得ください。
正直なところ私自身半信半疑で、xの係数が1未満なら0.5でも0.7でも何だって良いような気がするのですが、0.5だとうまくいかないので仕方ありません。
ところで先程の式は原理的な式であって２次元座標と３次元座標の原点位置の差と、各座標の変動範囲を考慮してスケーリングを行うと
S = (+0.7x -f(x,y))/4 +25
T = (+y -0.7x)/2 +50
となります。(*1)
これで、３次元空間内でx,y,f(x,y)がそれぞれ0～99の範囲で変化しても、画面(S=0～49、T=0～99)からはみ出すことは無いはずです。
【手順３】配列の設計とアクセス方法
次に、今回のプログラムで必要な配列を考えます。
まず、画面表示バッファとして要素数50×100=5000個の配列が必要となります。これを
画面表示バッファ Z(0,0)～Z(49,99)
とします。
例えば、配列Z(34,87)に数値42が書き込んである状態で表示ルーチン(まだ説明していないけど)を実行すると、パソコン画面上の座標(34,87)にはASCIIコードで42である"*"が表示されます。
それから、立体を画面表示させる場合「障害物に隠れて見えない面は表示しない」という処理(隠面処理)が必要となります。
次回説明しますが、この隠面処理に横軸(T軸)の各点毎の最大高さと最小高さを記憶する配列を用意します。
最大高さ記憶配列 MAXH(0)～MAXH(99)
最小高さ記憶配列 MINH(0)～MINH(99)
とします。
なお、現状の超小型記号言語(Ver1.4)の配列の要素数は約1000個ですから、これらの配列エリアを準備するのは到底無理です。
ここはあっさり諦めて、超小型記号言語のバージョンアップ(Ver1.4→Ver1.5)で配列エリアを拡張することにしました。
次に、超小型記号言語での配列へのアクセス方法について説明します。 単純な変数Zへのアクセスであれば
#42 Z! (変数Zに数値42を代入)
Z@ ￥ (変数Zの内容を文字表示)
で代入、文字表示できます。
この例では、変数Zの内容は文字コードと解釈され、ディスプレイ上には"*"が表示されます。
１次元配列Z(X)の場合、その実体は「変数Zの番地と変数Xの値を番地と見做して加算した番地」のメモリです。
したがって、配列Z(55)なら変数Zの代わりに「Z #55 +」を使い(逆ポーランド記法ね！)
#42 Z #55 +! (配列Z(55)に数値42を代入)
Z #55 +@ ￥ (配列Z(55)の内容を文字表示)
とします。
超小型記号言語ではアルファベット大文字(A～Z)が単独で現れた場合、数値(0～25)と解釈されることさえ知っていれば問題ないと思います。
では、２次元配列Z(X,Y)はどうするかというと、これを配列Z(X+50*Y)と考えます。(ここではX,Yの範囲が、X=0～49,Y=0～99という前提です)
１次元配列になってしまえば理屈は簡単で、例えば配列Z(13,77)に数値42を代入するには
#42 Z #13 #50 #77 *++! (配列Z(13+50*77)に数値42を代入)
Z #13 #50 #77 *++@ ￥ (配列Z(13+50*77)の内容を文字表示)
とします。
２次元配列Z(13,77)を１次元配列Z(13+50*77)に置き換えているのが判るかと思います。
今回の説明はここまでとします。
次回、後編で表示方法、隠面処理、エッジの強調などの説明をしたいと思います。
理屈はともかく、すでに完成してる３Ｄグラフィクスプログラム(3dg.vts)を動かしてみましょう。
Windowsパソコン上で動作させるとして、必要なファイルは
超小型記号言語 vts15.zip
サンプルプログラム vts-prg1.zip
です。両方とも「超小型記号言語」のページからダウンロードできます。
２つのファイルを解凍後、同一フォルダに入れ、超小型記号言語の実行ファイルvtsxx15.exeをダブルクリックします。
コマンドプロンプトが立ち上がって"?"を表示したら、"3dg.vts"Enterをキーインすると動作します。
なお、コマンドプロンプトの文字が大きいと正しく曲面が表示されません。
この場合は、１回だけは次の操作を行ってください。 コマンドプロンプトの上側のタイトルバーを右クリックし「プロパティ」を選択
ここで
フォントサイズ=10
フォント=MSゴシック
画面バッファのサイズ=幅120以上、高さ60以上
ウィンドウのサイズ=幅120以上、高さ60以上
であることを確認(必要に応じて変更)
3dg.vts内の関数FNZを書き換えることで、色々な関数の曲面を表示することができます。
サンプルプログラムに曲面定義集が同梱されていますので参考にしてください。(三角関数とかは無理なので悪しからず)
最後に、実行画面を図2-2に示します。(ディスプレイの解像度は1280×1024)
図2-2
では、次回まで。
--------------------
(*1)実際のプログラミングでは、この式を変形して
S = (+7x -10f(x,y))/40 +25
T = (+10y -7x)/20 +50
で計算します。この方が整数演算で誤差が小さくなります。超小型記号言語のコードは逆ポーランド記法なので
X@ #7 * X@ Y@ :FNZ #10 * - #40 / #25 + S!
Y@ #10 * X@ #7 * - #20 / #50 + T!
となります。
(*2)2019/01/02に全面改訂しました。
■第３回 ｚ＝ｆ(ｘ，ｙ)のグラフ／プログラミング後編(2019/01/27掲載,2019/02/04最終変更)
前回(第２回)は、具体的な仕様と設計過程(途中まで)を説明しながら、最後に完成したプログラムを実行してみました。
すでにプログラムは完成しているので「実はそう簡単にはいかない諸問題」はすべて解決済みなのですが、一応簡単に概念だけでも解説しておきます。
すでに勝敗が決まった消化試合みたいな気分ですが、頑張りましょう。f(^_^;
【手順４】画面表示と隠面処理(*1)、エッジの強調
まずは画面表示バッファ(配列)に表示データを書き込み、このデータを画面に表示することから始めます。
(1)画面表示
前回、曲面z=f(x,y)とパソコン画面上の点(S,T)の関係は、整数演算で誤差を小さくすることを考慮すると(脚注で)
S = (+7x -10f(x,y))/40 +25
T = (+10y -7x)/20 +50
と説明しました。
そこで、この式に従ってx,yをそれぞれ0から99(*2)まで変化させて、この時の画面表示バッファ(S,T)に"."(ドット)をプロット(文字コードを書込む)ことにします。
画面表示バッファ(S,T)の初期値を" "(空白)の文字コードにしてあるので、全部プロットした時点では、" "または"."のいずれかの文字コードが書き込まれていることになります。
この配列を画面にプリントすれば、" "でも"."でもカーソルは右に進むので画面上に"."がプロットされた曲面が表示されます。
これが画面表示ルーチンで簡単ですが100文字ごとに改行を挟み込みます。この理由は説明不要ですよね。
(2)隠面処理
隠面処理とは、目の前に物体があればその後の物体は隠れて見えないようにする処理です。
現実世界ではこの事実をあたりまえに受け入れていますが、今回の曲面z=f(x,y)も自分で自分の曲面の他の部分を隠す場合があります。
これをきちんと「見えない」ように処理してあげないと立体には見えません。(動いていればまだマシかもしれませんが…)
具体的には(1)で説明した画面表示バッファに書き込むとき
・書き込む順番は、視点から見て手前から奥の順(x=99,98,97,…,2,1,0)とする(*2)
・このとき手前に曲面があれば、隠れて見えないはずなので"."は書き込まない
とします。
手前の曲面とは、現時点までに画面に書き込んだSの最小値～最大値の範囲です。
前回ちょっとだけ触れた
最大高さ記憶配列 MAXH(0)～MAXH(99)
最小高さ記憶配列 MINH(0)～MINH(99)
は、この判定を行うためのもので、もし点(S,T)が
MAXH(T)≦S≦MINH(T)
なら書き込みを行わないようにします。
反対に点(S,T)がこの範囲を外れていれば視点から隠れていないことになるので"."を書き込み、配列MAXH(T)またはMINH(T)更新を更新します。
図3-1は、手前(x=99)から奥(x=0)へ短冊上の面が山形に折れ曲がっている曲面です。
縦位置Tに注目すると、おもて面は通常通り書き込みますが、裏面(図で白丸の部分)は、MAXH(T)≦S≦MINH(T)が成立するため書き込みません。
図3-1
(3)エッジの強調
これはおまけみたいな機能ですが、隠面処理でMAXH(T)より大きい、すなわち手前のプロットのときだけ"."ではなく"*"を使用すると、手前の曲面の端が強調され曲面の境界が見やすくなります。
試しに追加した処理ですが、あまりにも見やすいのでそのまま残してあります。
前回の実行画面(図2-2)でその効果が判ると思います。
さて、以上で３回に渡った「ｚ＝ｆ(ｘ，ｙ)のグラフ」は今回で終了です。
次回のテーマは決まっていませんが、また近いうちにお会いしましょう。
--------------------
(*1)参考文献
「エンサイクロペディア・アスキーＶｏｌｕｍｅ３、Ｈ６８ＴＲ ３次元グラフィックス(アスキー出版)」
「エンサイクロペディア・アスキーＶｏｌｕｍｅ３、ａｐｐｌｅⅡ ３Ｄグラフィックス(アスキー出版)」
(*2)実際のコードでは、すべてのx,yで計算をするのではなく、xは4つ飛びに、yは2つ飛びで計算しています。
画面の解像度が低い(粗い)のでこれで十分ですし、計算時間も短くなります。
それからxは99ではなく100から始めています。
すなわち、xは100,96,92,…の順に、yは0,2,4,…の順に計算しています。
このソフトウェアに関するご意見・ご質問等をお待ちしています。
ｔｏ＠ｔａｋａｍｉ.ｃｏｍ （アドレスは半角で入力ください）
[トップページに戻る]