作業七

B94611003

一.
第一秒


















第二秒:



















第三秒:



















第四秒:



















第五秒:




















二.

第一桿的
(1)速度 : http://homepage.ntu.edu.tw/~b94611003/7/1v.jpg
(2)加速度 : http://homepage.ntu.edu.tw/~b94611003/7/1a.jpg

第二桿的
(1)速度 : http://homepage.ntu.edu.tw/~b94611003/7/2v.jpg
(2)加速度 : http://homepage.ntu.edu.tw/~b94611003/7/2a.jpg

第三桿的
(1)速度 : http://homepage.ntu.edu.tw/~b94611003/7/3v.jpg
(2)加速度 : http://homepage.ntu.edu.tw/~b94611003/7/3a.jpg


三.

動畫: http://homepage.ntu.edu.tw/~b94611003/7/flash.html


本次作業的函式:

function dyad_run(rho,theta,td,tdd,time)
set(gca,'xlim',[0 100],'ylim',[0 100],...
'NextPlot','replace','Visible','off')
tdT=td;
thetaT=theta;
thetaTEMP=theta;
mov = avifile('walk.avi')
for n=1:length(time)
for p=1:length(rho)
tdT(p)=td(p)+tdd(p)*(n-1);
thetaT(p)=theta(p)+td(p)*(n-1)+0.5*tdd(p)*(n-1)^2;
end
[v,a]=dyad_draw(rho,thetaT,tdT,tdd)
if n > 1
for q=1:length(rho)
vv(q,n-1)=v(q);
aa(q,n-1)=a(q);
end
end
pause(1)
F = getframe(gca);
mov = addframe(mov,F);
end
mov = close(mov);