2007年3月30日 星期五

機動學作業四

4.1
程式碼
AXIS([-80 80 -80 80])
L=34+10
x=[L/sqrt(3)*cosd(30) 0 -L/sqrt(3)*cosd(30) L/sqrt(3)*cosd(30)]
y=[-L/sqrt(3)*sind(30) L/sqrt(3) -L/sqrt(3)*sind(30) -L/sqrt(3)*sind(30)]
TRI=line(x',y')
xlabel('X');
ylabel('Y');
grid on
for n=1:360
rotate(TRI,[0 0 1],-1,[x(3) y(3) 0])
pause(0.0034)
end;

for n=1:360
rotate(TRI,[0 0 1],-1,[x(2) y(2) 0])
pause(0.0034)
end;

for n=1:360
rotate(TRI,[0 0 1],-1,[x(1) y(1) 0])
pause(0.0034)
end;
三角形連桿
4.2
程式碼
x=[0 0]
y=[10 0]
for n=1:12
m=[cosd(30*n) sind(30*n); sind(-30*n) cosd(30*n)];
Y=y*m
linkshape([Y(1) Y(2)],[x(1) x(2)],4)
axis equal
line([Y(1) 15],[Y(2) 0],'Color','g','LineWidth',1.5)
end
連桿羽彈簧
藍色線為彈簧的變化

4.3
4.3.1
程式碼
A=[0 0]
B=[3 4]
C=[13 4]
D=[10 0]
linkshape(A,D,2)
linkshape(A,B,3)
linkshape(B,C,1.5)
linkshape(C,D,2)
四連桿

4.3.2
程式碼
A=[0 0]
B=[3 4]
C=[13 4]
D=[10 0]
linkshape(A,D,2)
for n=[1:12]
m=[cosd(30*n) sind(30*n); sind(-30*n) cosd(30*n)];
b=B*m
m1=[cosd(30*n) sind(30*n); sind(-30*n) cosd(30*n)];
c(1)=b(1)+10
c(2)=b(2)
linkshape([c(1) c(2)],D,2)
linkshape([b(1) b(2)],[c(1) c(2)],1.5)
linkshape(A,[b(1) b(2)],3)
axis equal
pause(2.0)
end
四連桿軌跡

沒有留言: