MATLAB LECTURE - Trang 38

16/09/2015

1

Nguyen Quang Hoang

Department of Applied Mechanics

Giải phương trình vi phân trong Matlab

Nguyen Q.Hoang

Department of Applied Mechanics

Hanoi University of Science and Technology

Lecture 6

1

Nguyen Quang Hoang

Department of Applied Mechanics

Chương 6.

Lập trình trong Matlab

2

6.1 Giải phương trình vi phân bậc nhất với ode23 và ode45

6.2 Giải các phương trình vi phân bậc hai

6.3 Bài tập thực hành

Nguyen Quang Hoang

Department of Applied Mechanics

Ví dụ giải phương trình vi phân thường

3

Tìm hàm y(t) thỏa mãn phương trình vi phân

Với phương trình đơn giản này ta có thể viết ra được nghiệm chính xác
(nghiệm giải tích) ,

cos( )

dy

t

dt



( )

sin( )

y t

t

C





(0) 2

y



với điều kiện đầu

với C là hằng số tích phân phụ thuộc vào điều kiện đầu. Từ điều kiện
đầu nhận được hằng số C:

(0)

sin(0)

2

2

y

C

C











( )

sin( ) 2

y t

t





Nguyen Quang Hoang

Department of Applied Mechanics

Giải số bằng lệnh ode23

4

Cú pháp lệnh ode23

[t, y] = ode23(

‘funct_name’, [t0, t_end], y0)

Giá trị điều kiện đầu

Khoảng thời gian cần giải

t từ t0 đến t_end

Tên m-file chứa phương trình vi
phân cần giải

Kết quả tính được ghi vào hai véctơ
t và y chứa các cặp [t(i), y(ti)]