108
D=[0.025, 0.050, 0.100, 0.200]; % cac he so can Lehr D khac nhau
ii=1;
t=linspace(0, 4*pi/Omega, 100);
A=q/om_0^2/sqrt((1.0-eta_0^2)^2+4*D(ii)^2*eta_0^2);
% bien do
Psi=atan2(2.0*D(ii)*eta_0, 1-eta_0^2);
% goc pha
x=A*cos(Omega*t-Psi); % ham x(t)
xp=-Omega*A*sin(Omega*t-Psi); % ham xp(t) hay x_dot(t)
figure('name','3D')
plot3(t, x, xp, 'k','Linewidth',2), hold on
% chuyen dong khong gian
plot3(t, x, 1.2*min(xp)*ones(size(t)), 'k','Linewidth',1)
% hinh chieu trong mat phang (t, xp)
plot3(t, 1.2*max(x)*ones(size(t)), xp, 'k','Linewidth',1)
% hinh chieu trong mat phang (t, x)
plot3(t*0, x,xp, 'b','Linewidth',2)
% qui dao pha (x, xp)
axis([t(1) t(end), 1.2*[min(x) max(x)], 1.2*[min(xp) max(xp)]]);
box on
xlabel('t'); ylabel('x'); zlabel('xp');
view(28,42);
Chạy chương trình sẽ cho ta đồ thị như hình 3-46.
0
0.2
0.4
0.6
0.8
1
1.2
−0.2
−0.1
0
0.1
0.2
−2
−1
0
1
2
x
t
xp
Hình 3-46. Đồ thị dao động trong không gian chuyển động, cùng các hình chiếu
