作業七

b94611021 張淳皓
本人4/19有上課 部落格

7.1.1
　在1~5秒內　各桿方向如圖　（藍色線即為桿）
第一秒
第二秒
第三秒
第四秒
第五秒

7.1.2
桿一之速度、加速度
桿二之速度、加速度
桿三之速度、加速度

7.1.3

將過程繪製成影片
影片

本次作業之函式：

function dyad_run(rho,theta,td,tdd,time)
%本函式引用dyad_draw function
%多出的輸入項time是時間矩陣
%ex. dyad_run([31 36 26],[0 0 0],[0.2 0.5 0.3],[0 0.1 0.2],[0 1 2 3 4 5])

thetaTEMP=theta;
tdTEMP=td;
thetaTEMP=theta;
for n=1:length(time)
for m=1:length(rho)
tdTEMP(m)=td(m)+tdd(m)*(n-1);
thetaTEMP(m)=theta(m)+td(m)*(n-1)+0.5*tdd(m)*(n-1)^2;
end
[v,a]=dyad_draw(rho,thetaTEMP,tdTEMP,tdd)
if n > 1
for k=1:length(rho)
vv(k,n-1)=v(k);
aa(k,n-1)=a(k);
end
end
title('b94611021')
pause(1)
end
for j=1:length(rho)
figure
title('桿速度')
plot(vv(j,:),'r+-')
figure
title('桿加速度')
plot(aa(j,:),'cx--')
end