OpenGL関数でアニメーション

clear all;

AssertOpenGL;

ListenChar(2);

screenNumber=max(Screen('Screens'));

InitializeMatlabOpenGL;

% 頂点の情報

vertex = [

    0, 0, 0 % A (1)

    1, 0, 0 % B (2)

    1, 1, 0 % C (3)

    0, 1, 0 % D (4)

    0, 0, 1 % E (5)

    1, 0, 1 % F (6)

    1, 1, 1 % G (7)

    0, 1, 1 % H (8)

    ];

% ２つの頂点を結ぶ線分の情報

% MATLABでは行列の参照が０からではなく、１からであることに注意。

edge = [

    1, 2 % ア（A-B）

    2, 3 % イ（B-C）

    3, 4 % ウ（C-D）

    4, 1 % エ（D-A）

    5, 6 % オ（E-F）

    6, 7 % カ（F-G）

    7, 8 % キ（G-H）

    8, 5 % ク（H-E）

    1, 5 % ケ（A-E）

    2, 6 % コ（B-F）

    3, 7 % サ（C-G）

    4, 8 % シ（D-H）

    ];

try

    [windowPtr, windowRect] = Screen('OpenWindow', screenNumber);

    % スクリーンの中心座標

    [centerPos(1), centerPos(2)] = RectCenter(windowRect);

    % OpenGL関数を使い始めるという宣言

    Screen('BeginOpenGL', windowPtr);

    glMatrixMode(GL.PROJECTION); % アンダーバーがドットになっていることに注意。

    glLoadIdentity;

    gluPerspective(30.0, windowRect(3)/windowRect(4), 1.0, 100.0);

    glMatrixMode(GL.MODELVIEW);

    r = 0; % 回転角の初期値

    flag = 1; % flagが1のときrが自動的に加算される。（アニメーション）

    while KbCheck; end; % いずれのキーも押していないことを確認

    while 1 % while 文の中をぐるぐる回ります。

        [ keyIsDown, keyTime, keyCode ] = KbCheck; % キーが押されたか、そのときの時間、どのキーか、の情報を取得する

        if keyIsDown

            if keyCode(KbName('ESCAPE')) % エスケープキーでプログラムを終了

                break;

            end

            if keyCode(KbName('LeftArrow'))

                flag = 1-flag; % アニメーションを有効にしたり、無効にしたり。

            end;

            if keyCode(KbName('RightArrow'))

                r = r + 1; % 手動で回転角度を増やす。

                if r >= 360

                    r = 0;

                end;

            end;

            while KbCheck; end;                    

        end;

        glClearColor(1,1,1,1); % 背景色

        glClear;

        glLoadIdentity;

        gluLookAt(3.0, 4.0, 5.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0); % 視点位置と視線方向

        glRotated(r, 0, 1, 0); % 図形の回転

        % 図形の描画

        glColor3d(0.0, 0.0, 0.0);

        glBegin(GL.LINES); % アンダーバーがドットになっていることに注意。

        for i = 1:12 % 12本の線分を描く

            glVertex3dv(vertex(edge(i,1),:));

            glVertex3dv(vertex(edge(i,2),:));

        end;

        glEnd();

        % OpenGL関数の使用を終わります。

        Screen('EndOpenGL', windowPtr);

        % この段階で、PTBの機能を使って図形を描画することも可能です。

        Screen('FillRect', windowPtr, [255 0 0], [200 300 400 400]);

        Screen('Flip', windowPtr);

        % 次の繰り返しのために

        Screen('BeginOpenGL', windowPtr);

        if flag % flagが1のときはアニメーション

            r = r + 1;

            if r >= 360

                r = 0;

            end;

        end;

    end;

    % OpenGL関数の使用を終わります。

    Screen('EndOpenGL', windowPtr);

    Screen('CloseAll');

    ShowCursor;

    ListenChar(0);

catch

    Screen('CloseAll');

    ShowCursor;

    ListenChar(0);

    psychrethrow(psychlasterror);

end