207
Trong phần sau ta sẽ khảo sát chuyển ñộng của con lắc có khối lượng
0.25
m
=
kg, chiều dài dây
0.5
l
=
m, với các góc lệch ban ñầu khác nhau từ 10 ñến 110 ñộ.
Các kết quả mô phỏng ñược ñưa ra như trên hình 9-19, 9-20 và 9-21.
global m l k
m = 0.25; l = 0.5;
k = 0.0; % khong can
% k = 0.1; % co can
t_start = 0; t_end = 6.0;
phi = [10 : 20 : 110];
for i = 1:length(phi)
y1 = phi(i)*pi/180; y2 = 0; y0 = [y1 y2]';
[t,y] = ode45('conlac_don',[0: 0.02: t_end],y0);
plot(t, y(:,1),'k-', 'linewidth',1), grid on
xlabel('t [s]'), ylabel('\phi [rad]'),
hold on
end
0
1
2
3
4
5
6
- 2
- 1
0
1
2
t [s]
φ
[r
a
d
]
Hình 9-19. ðồ thị theo thời gian của góc lắc khi không có cản, k = 0
0
1
2
3
4
5
6
- 2
- 1
0
1
2
t [s]
φ
[r
a
d
]
Hình 9-20. ðồ thị theo thời gian của góc lắc khi có cản, k = 0.2
