HW7

B94611023許惠善
本人(4/19)有上課

本次作業其實可以利用一個簡短的程式同時完成，所以我只做了一個程式。將三小題的答案同時呈現在同一張圖表中。紅色為各桿之加速度，藍色為速度，綠色為合速度。

另外也可將這隻程式碼做動畫以及圖片的切換。 將12行clf以及23行pause刪掉就變成圖片(最下方依序為起始點到五秒過後的對應圖片)。


以下為程式碼：

%輸入a長度，我的學號為b94611023，故取a=3(最末碼)+10
a=3+10;

%輸入長度，角度，角速度，角加速度之初始值
rho=[a, a+5, a-5];
theta=[0, 0, 0];
td=[0.2, 0.5, 0 .3];
tdd=[0, 0.1, 0.2];

%本圖形的角度以及角速度會受角加速度影響，亦即受時間改變
for t=0:1:5
    clf;
    for i=1:1:3  
        %在此我們必須先做角度的改變，才能作角速度的改變
        %因為角度的公式是theta2=theta1+v1t+ωt^2，所以我們要用角度在前一秒的值以和速度前一秒的值來做
        %又t間隔都為一秒，所以在此t不代1進去計算
    theta(i)=theta(i)+td(i)+0.5*tdd(i);
        %角度則為td2=td1+ωt，為前一秒的角度加上角加速度*t，所以同樣t間隔1
    td(i)=td(i)+tdd(i);
    end;
        %繪圖
    dyad_draw([a, a+5, a-5],[theta(1),theta(2),theta(3)],[td(1),td(2),td(3)],[tdd(1),tdd(2),tdd(3)])
    pause(2);

end;
axis equal;