0.3 Symbolic và simulink  (Page 2/2)

Vẽ mặt 3D bằng hàm ezsurf: Ví dụ vẽ mặt $S=f(x,y)=\frac{y}{1+x^2+y^2}$ trên miền xác định: -5<x<5; -2<y<2.

>>syms x y

>>S=y/(1+x^2+y^2)

>>ezsurf(S, [-5 5 –2*pi 2*pi])

Hình 3.2 – Vẽ đồ thị 3D cho hàm số symbolic

Simulink

Simulink 5.0 (Simulation and Link - R13) được MatWorks giới thiệu vào tháng 6 năm 2003. Nó cho phép phân tích, mô hình hóa và mô phỏng các hệ thống động tuyến tính và phi tuyến, liên tục và rời rạc một cách trực quan trong môi trường giao tiếp đồ họa, bằng các thao tác chuột đơn giản. Có thể nói, không tận dụng được Simulink là một thiệt thòi lớn cho người làm công tác mô phỏng!

Khởi động Simulink bằng một trong các cách sau:

nhập:>>simulink

hoặc nhấp chuột vàotrên menubar của Matlab

Thư viện simulink hiện ra như hình 3.3:

Trước tiên, sinh viên hãy nhấp chuột vào các thanh cuộn của thư viện để có cái nhìn thân thiện về simulink.

Từ đây, để có thể tạo mô hình bằng simulink, hãy:

• nhấp chuột vào biểu tượngcủa thư viện simulink
• chọn: File – New – Model trong Menu của thư viện Simulink
• chọn: File – New – Model trong cửa sổ lệnh của Matlab

Hình 3.3 – Cửa sổ chính của thư viện Simulink

Hình 3.4 – Môi trường soạn thảo của Simulink

Cửa sổ này (hình 3.4) cho phép ta ‘nhấp - kéo - thả’ vào từng khối chức năng trong thư viện simulink. Ví dụ, đặt vào đây khối ‘Sine Wave’ trong thư việncủa(hình 3.5):

Hình 3.5 – Lấy một khối từ thư viện

Sau khi đặt tất cả các khối cần thiết của mô hình, ta nối chúng lại bằng cách ‘nhấp - giữ’ và kéo một đường từ ngõ ra của khối này đến ngõ vào của khối khác rồi thả phím trái chuột, một kết nối sẽ được thiết lập.

Xây dựng mô hình hệ thống xe tải:

Hình 3.6 – Mô hình xe tải

cho bởi phương trình: $m\frac{\text{dv}}{\text{dt}}=u-\text{bv}$ hay $\frac{\text{dv}}{\text{dt}}=\frac{1}{m}(u-\text{bv})$ .

Trong đó m là khối lượng xe, u là lực tác động của động cơ (ngõ vào của mô hình), b là hệ số ma sát và v là vận tốc đạt được (ngõ ra của mô hình).

Xuất phát điểm của việc xây dựng các mô hình hệ thống từ các phương trình vi phân tương ứng là các khối tích phân (Integrator). Nếu trong phương trình mô tả hệ thống có vi phân bậc n thì ta sẽ đặt vào mô hình n khối tích phân, do quan hệ $\int \frac{\text{dv}}{\text{dt}}=v$ .

• Mở một cửa sổ mô hình mới.
• Đặt vào mô hình khối ‘Integrator’ từ thư viện ‘Continuous’ và kẻ các đường thẳng nối đến ngõ vào và ngõ ra của khối này.
• Đặt nhãn ‘vdot’ (dv/dt) cho cho đường nối đến ngõ vào và ‘v’ cho đường nối đến ngõ ra bằng cách nhấp đúp chuột ngay phía trên các đường này.

Từ phương trình hệ thống ta thấy dv/dt (vdot) bằng tích của thành phần (1/m) và thành phần tổng (u-bv), nên ta thêm khối 1/m ngay trước khối tích phân:

• Đặt vào khối ‘Gain’ trong thư viện .
• Nhấp đúp chuột vào khối này để thay đổi độ lợi thành 1/m.
• Đặt nhãn ‘inertia’ cho khối này để tượng trưng cho quán tính của xe (nhấp đúp vào nhãn ‘Gain’ bên dưới khối).

Bây giờ ta đặt khối tổng với 2 ngõ vào ‘+-‘, ngõ vào ‘+’ sẽ được nối với u, ngõ vào ‘-’ sẽ được nối với thành phần bv để được (u-bv)

• Đặt vào khối ‘Sum’ trong thư viện
• Nhấp đúp vào khối này để đổi ngõ vào từ ‘++’ sang ‘+-’

***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.***

Để được thành phần bv ta chỉ cần đặt thêm khối ‘Gain’ với độ lợi b:

• Đặt khối ‘Gain’ có độ lợi b
• Đặt nhãn là ‘damping’ tượng trưng cho thành phần lực cản của xe.

***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.***

Đến đây việc xây dựng mô hình xe tải với ngõ vào u và ngõ ra v coi như hoàn thành. Tuy nhiên, để mô phỏng mô hình này, ta cần đặt thêm khối ‘Step’ vào u và hiển thị v trên khối ‘Scope’

• Đặt khối ‘Step’ trong thư việnbiên độ u ngay ngõ vào.
• Đặt khối ‘Scope’ trong thư viện ngay ngõ ra v.

***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.***

Nhớ rằng m, b và u là các biến cần được gán trị trước khi mô phỏng.

>>m=1000

>>b=50

>>u=500

Thời gian mô phỏng hệ thống tùy thuộc vào thông số ‘Stop time’ trong menu ‘SimulationSimulation paramrters’, giả sử đặt 120 (Hình 3.7).

Để chạy mô phỏng, ta có thể thực hiện bằng 1 trong các cách:

• nhấp chuột vào biểu tượngtrên menubar của mô hình
• chọn: Simulation  Start
• Ctrl-T

Nhấp đúp vào ‘Scope’ để xem kết quả mô phỏng.

Hình 3.7 – Thay đổi thông số mô phỏng

Xây dựng mô hình hệ thống điều khiển vị trí motor DC cho bởi phương trình vi phân sau:

$\frac{d^2\theta }{{\text{dt}}^2}=\frac{1}{J}K_{t}i-b\frac{\mathrm{d\theta }}{\text{dt}}$

$\frac{\text{di}}{\text{dt}}=\frac{1}{L}-\text{Ri}+V-K_e\frac{\mathrm{d\theta }}{\text{dt}}$

Trong đó:

J = 0.01 Kgm2/s2 là moment quán tín của rotor

b = 0.1 Mms là hệ số ma sát của các bộ phận cơ khí

K = Ke = Kt = 0.01 Nm/A là hằng số sức điện động

R = 10 ohm là điện trở dây quấn

L = 0.5 H là hệ số tự cảm

V là điện áp đặt lên cuộn dây của motor

 là vị trí trục quay (ngõ ra của mô hình)

i là dòng điện chạy trong cuộn dây của motor.

Hình 3.8 – Mô hình toán hệ điều khiển vị trí motor DC

Quan sát từng phương trình mô tả hệ thống ta thấy cấu trúc của chúng cũng tương tự như phương trình của câu . Sinh viên lần lượt thực hiện trên từng phương trình để được (hình 3.9):

(Hình 3.9)

Kết hợp 2 phương trình:

(Hình 3.10)

Đặt vào mô hình khối ‘Step’ để làm tín hiệu tham khảo, khối ‘Scope’ để quan sát đáp ứng. Sinh viên hãy gán trị cho tất cả các thông số của mô hình, thực hiện mô phỏng và quan sát đáp ứng (Hình 3.11).

Hãy thay đổi ngõ vào mô hình bằng khối tạo xung vuông. Mô phỏng, quan sát kết quả và nhận xét.

Tự chọn

Sinh viên hãy tính đạo hàm cấp 2 của hàm $y={\text{xe}}^{1-x^2}$ bằng tay và kiểm chứng kết quả bằng symbolic.

Tính tích phân sau và kiểm chứng kết quả bằng symbolic: $\underset{0}{\overset{\infty }{\int }}{e}^{-x^2}$

Hình 3.11 - Mô hình Simulink hoàn chỉnh của hệ điều khiển vị trí motor DC

Xây dựng mô hình hệ thống xe lửa cho bởi phương trình:

$M_1\frac{d^2{x}_1}{{\text{dt}}^2}=F-k(x_1-x_2)-{\mathrm{\mu M}}_{1}g\frac{{\text{dx}}_1}{\text{dt}}$

$M_2\frac{d^2{x}_2}{{\text{dt}}^2}=k(x_1-x_2)-{\mathrm{\mu M}}_{2}g\frac{{\text{dx}}_2}{\text{dt}}$

(Hình 3.12 Photo courtesy: Dr. Howard Blackburn)

Trong đó các thông số tượng trưng như sau:

M1=1 kg là khối lượng toa kéo;

M2=0.5 kg là khối lượng toa khách;

k=1 N/sec là độ cứng lò xo kết nối giữa 2 toa;

F=1 N là lực tác động của đầu máy (ngõ vào mô hình);

=0.002 sec/m là hệ số ma sát lăn;

g = 9.8 m/s^2 là gia tốc trọng trường

x1, x2 vị trí 2 toa (ngõ ra).

Hình 3.13 – Mô hình toán của hệ thống xe lửa