Advertisement
tdttvd

BTL GT

Dec 31st, 2021
74
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 3.05 KB | None | 0 0
  1. function B9
  2. clc;
  3. clear;
  4. close all;
  5.  
  6. %Phần đầu (Khai báo biến, hàm...)------------------------------------------
  7. %Khai báo biến symbolic t
  8. syms t;
  9.  
  10. %Khai báo các hàm x và y theo t
  11. x = 3 * t^2 - 4/3 * t^3;
  12. y = 8 * t;
  13.  
  14. %Tìm vận tốc theo phương x và y = đạo hàm của x, y theo t
  15. vx = diff(x);
  16. vy = diff(y);
  17.  
  18. %Tìm gia tốc theo phương x và y = đạo hàm của vx, vy theo t
  19. ax = diff(vx);
  20. ay = diff(vy);
  21.  
  22. %Phần sau (Thực hiện yêu cầu các câu hỏi)----------------------------------
  23.  
  24. %Nhập khoảng thời gian cần vẽ đồ thị:
  25. disp('Nhap khoang thoi gian can ve do thi');
  26. t1 = input('Moc thoi gian dau tien: t1 = ');
  27. t2 = input('Moc thoi gian thu hai: t2 = ');
  28.  
  29. %----------------Hết phần nhập dữ liệu vẽ đồ thị câu a--------------------%
  30.  
  31. %Câu a: Vẽ đồ thị x, y theo t với giá trị t từ 0 tới 5
  32.  
  33. %Đồ thị 1: Đồ thị x, y ----------------------------------------------------
  34. subplot(3,1,1)
  35. fplot(x, y, [t1 t2]);
  36.  
  37. %Ký hiệu các trục x, y và tiêu đề
  38. title('Quy dao cua vat trong khoang thoi gian tu t1 den t2');
  39. xlabel('x');
  40. ylabel('y');
  41.  
  42. %Hiển thị lưới
  43. grid on;
  44.  
  45. %Vẽ 2 trục tọa độ
  46. xline(0);
  47. yline(0);
  48.  
  49. %Đồ thị 2: Đồ thị x, t ----------------------------------------------------
  50. subplot(3,1,2);
  51. fplot(t, x, [t1 t2]);
  52.  
  53. %Ký hiệu các trục và tiêu đề
  54. title('Hoanh do x cua vat trong khoang thoi gian tu t1 den t2');
  55. xlabel('t');
  56. ylabel('x');
  57.  
  58. %Hiển thị lưới
  59. grid on;
  60.  
  61. %Vẽ 2 trục tọa độ
  62. xline(0);
  63. yline(0);
  64.  
  65. %Đồ thị 3: Đồ thị y, t ----------------------------------------------------
  66.  
  67. subplot(3,1,3);
  68. fplot(t, y, [t1 t2]);
  69.  
  70. %Ký hiệu các trục và tiêu đề
  71. title('Tung do y cua vat trong khoang thoi gian tu t1 den t2');
  72. xlabel('t');
  73. ylabel('y');
  74.  
  75. %Hiển thị lưới
  76. grid on;
  77.  
  78. %Vẽ 2 trục tọa độ
  79. xline(0);
  80. yline(0);
  81.  
  82. %----------------------------Hết câu a------------------------------------%
  83.  
  84. %Nhập thời điểm cần xét cho các câu b, c, d
  85. time = input('Nhap thoi diem can xet, t = ');
  86.  
  87. %----------------Hết phần nhập dữ liệu cho câu b, c, d--------------------%
  88.  
  89. %Câu b: Tính độ lớn vận tốc tại t = 1s
  90. disp('Cau b: Do lon van toc chat diem tai t la:');
  91. v = sqrt(subs(vx, t, time)^2 + subs(vy, t, time)^2);
  92. disp(v);
  93.  
  94. %----------------------------Hết câu b------------------------------------%
  95.  
  96. %Câu c: Tính độ lớn gia tốc tại t = 1s
  97. disp('Cau c: Do lon gia toc chat diem tai t la:');
  98. a = sqrt(subs(ax, t, time)^2 + subs(ay, t, time)^2);
  99. disp(a);
  100.  
  101. %----------------------------Hết câu c------------------------------------%
  102.  
  103. %Câu d:
  104. %Tính bán kính cong của quỹ đạo tại t = 1s
  105. disp('Cau d: Ban kinh cong quy dao tai t la:');
  106.  
  107. %Vận tốc tại thời điểm t
  108. v = sqrt(vx^2+vy^2);
  109.  
  110. %Tính gia tốc tiếp tuyến
  111. at = diff(v);
  112.  
  113. %Tính gia tốc hướng tâm
  114. an = sqrt(a^2 - at^2);
  115.  
  116. %Tính bán kính quỹ đạo
  117. R = subs(abs(v^2/an), t, time);
  118. disp(R);
  119. end
  120.  
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement