for文で得た値をplotしたい

for文で0~3000回姿勢計算を行い、その一回ごとの結果をplotしてグラフにしたいのですが、どのように配列すればいいのでしょうか。

for i = 0 : 0.001 : 3

R = R + r *0.001

end

a=~;b=~;c=~;

上のa,b,cを3000ずつplotしたいです。

[Hernia Baby]

rはiでよろしいでしょうか？ またa,b,cはRに格納されており、i-1は保持されますか？ それとも更新し、過去の情報は消しますか？

[Jun Ozaki]

すみません、説明不足でした。

r ＝ R（単位行列：3×3行列）×W（角速度:3×3行列）の行列です。

MATLABに関して知識が薄いためi-1を保持する場合としない場合の違いが分からないのですが

iが０から0.001づつ変化するごとにプロットしたい場合、更新して過去の情報を消すと、iが3になった結果のみがグラフになるのでしょうか。

[Hernia Baby]

返信ありがとうございます。

自分の勘違いで更新はwhile文にする必要があると思い込んでいました。

確認しましたが、for文でも可能でした。

プロットはhold on状態でなければ、for文内に組み込むことでプロット可能です。

hold onの場合は事前のプロットが保持されます。

遷移をゆっくり見たい場合は、回答に入れました通り、pause関数でしばらく次の計算を待つことができます。

Answer 1

[Hernia Baby]

for文というよりwhile文のほうが便利かなと思い、サンプルコードを書きました。

刻み幅が0.001だと、plotがゆっくりなので0.01刻みにしています。

clear,clc,close all;

R = zeros(3,1);

cnt = 0;
while cnt <= 3
    r = 0.001*[cnt 10^3*sin(2*pi*cnt) 20^3*cos(2*pi*cnt)]';
    R = R + r;
    
    for i = 1:3
        subplot(3,1,i)
        scatter(cnt,R(i));
        xlim([0 3])
    end
    pause(0.001);
    cnt = cnt + 0.01;
end

[Hernia Baby]

今確認取ってみましたが、for文でも可能でした。

またplot3で3次元描写する方法も以下に示します。

clear,clc,close all;

R = zeros(3,1);

cnt = 0;
for j = 0:0.01:3
    r = 0.001*[j 10^3*sin(2*pi*j) 20^3*cos(2*pi*j)]';
    R = R + r;

    plot3(R(1),R(2),R(3),'ro')
    
%{
    % a,b,cのそれぞれで描写する場合
    for i = 1:3
        subplot(3,1,i)
        scatter(cnt,R(i));
        xlim([0 3])
    end
%}
    pause(0.001);
end