Phương pháp lặp giải phương trình vi phân cấp cao với hệ số phụ thuộc các phiếm hàm tích phân

ISSN: 1859-2171 TNU Journal of Science and Technology 200(07): 41 - 47 
 Email: jst@tnu.edu.vn 41 
PHƯƠNG PHÁP LẶP GIẢI PHƯƠNG TRÌNH VI PHÂN CẤP CAO VỚI HỆ SỐ 
PHỤ THUỘC CÁC PHIẾM HÀM TÍCH PHÂN 
Vũ Vinh Quang*, Lại Văn Trung 
Trường Đại học Công nghệ Thông tin và Truyền thông – ĐH Thái Nguyên 
TÓM TẮT 
Trong các dạng phương trình vi phân bậc cao, lớp các phương trình có chứa các hệ số phụ thuộc 
các phiếm hàm tích phân mô tả các dạng bài toán đặc biệt trong cơ học đang được các nhà toán 
học quan tâm trong thời gian gần đây. Các kết quả đã nghiên cứu thường đề cập đến tính chất định 
tính của các bài toán, tuy nhiên việc tìm nghiệm số của các dạng bài toán này chưa được quan tâm 
do dạng đặc biệt của bài toán. Nội dung của bài báo gồm hai phần: phần thứ nhất trình bày kết quả 
xây dựng lược đồ sai phân bậc cao đối với phương trình vi phân cấp hai, phần thứ hai trình bày các 
kết quả xây dựng lược đồ lặp giải các phương trình vi phân bậc cao có các hệ số phụ thuộc các 
phiếm hàm dạng tích phân. Các kết quả thực nghiệm số trên các lược đồ sai phân bậc cao khẳng 
định tính hữu hiệu của các sơ đồ lặp đã đề xuất. 
Từ khóa: Phương trình vi phân bậc cao, lược đồ sai phân, phiếm hàm tích phân, sơ đồ lặp. 
Ngày nhận bài: 01/3/2019; Ngày hoàn thiện: 25/3/2019; Ngày duyệt đăng: 07/5/2019 
ITERATIVE METHOD FOR SOLVING HIGHER 
ODER DIFFERENTIAL EQUATIONS WITH COEFFICIENTS 
DEPENDENT ON INTEGRAL FUNCTIONS 
Vu Vinh Quang
*
, Lai Van Trung 
University of Information and Communication Technology - TNU 
ABSTRACT 
In higher order differential equations, classes of equations contain coefficients that depend on 
integral functions describing special types of problems in mechanics that are of interest to 
scientists recently. The results studied often refer to the qualitative nature of the problems. 
However, finding the numerical solution of these types of problems has not been considered due to 
the special form of the problem. The content of the article consists of two parts: The first part 
presents the results of developing a higher oder differential diagram for the second order 
differential equation, the second part presents the results of building a loop iteration of higher oder 
differential equations with coefficients that depend on integral functions. The numerical 
experimental results on higher oder differential diagram confirm the effectiveness of the proposed 
iteration diagram. 
Keywords: Higher order differential equations, differential diagram, integral function, iteration 
diagram. 
Received: 01/3/2019; Revised: 25/3/2019; Approved: 07/5/2019 
* Corresponding author: Email: vvquang@ictu.edu.vn 
Vũ Vinh Quang và Đtg Tạp chí KHOA HỌC & CÔNG NGHỆ ĐHTN 200(07): 41 - 47 
 Email: jst@tnu.edu.vn 42 
1. Giới thiệu 
Khi nghiên cứu về phương trình vi phân phi 
tuyến tính bậc cao, một trong những kết quả 
đã đạt được trong những năm qua có thể kế 
đến các kết quả của nhóm tác giả: Đặng 
Quang Á, Nguyễn Thanh Hường, Ngô Thị 
Kim Quy. Dạng phương trình phi tuyến tổng 
quát nhất được xét là dạng: 
4 ' '' '''
' '
0 1 0
'' ''' '' '''
0 1 0 1
, , , , , ,
, ,
, .
u x f x u u u u a x b
a u a a u a A b u b bu b B
c u a c u a C d u b d u b D
 (1) 
Về mặt lý thuyết, sự hội tụ của sơ đồ lặp trên 
đã được chứng minh bằng lý thuyết của 
phương trình toán tử. Việc giải số các bài toán 
cấp hai đã thực hiện bằng việc xây dựng các 
lược đồ sai phân với độ chính xác cấp 4 cho 
bài toán cấp hai. 
''
'
0 1
'
0 1
( ) ( ), ,
( ) ( ) ,
( ) ( ) .
u x f x a x b
u a u a A
u b u b B (2) 
Các kết quả đã được công bố trong [1]. 
Trong các công trình gần đây, một số tác giả 
trên thế giới đã đề cập tới các mô hình bài 
toán cơ học có sự phụ thuộc tích phân. Cụ 
thể, trong [2], các tác giả: N. Kachakhidze, N. 
Khomeriki, J. Peradze, Z. Tsiklauri đã nghiên 
cứu mô hình bài toán được mô hình hóa bởi 
phương trình phi tuyến cấp hai 
1
'' 2 ''
0
(w ) w ( ), 0 1,
w(0) w(1) 0.
dx f x x
 (3) 
Trong [3], các tác giả Q.A. Dang, T.L.Vu đã 
nghiên cứu mô hình bài toán cấp bốn phi 
tuyến dạng 
4 '' ' 2 ''
0
'' ''
2
( ) ( ), 0 1,
(0) ( ) 0, (0) ( ) 0.
y y y dx y p x x
y y y y
 (4) 
Trong [4], các tác giả T.F. Ma, A.L.M. 
Martinez đã xét dạng bài toán biên phi tuyến 
cấp 4 dạng 
24 ' '' '
0
'' '' '
( ( ) s) ( ) ( , , u ),
(0) , ( ) ,
(0) , ( ) ( ( )).
L
u M u s d u x f x u
u A u L B
u C u L g u L
 (5) 
Có thể thấy rằng điểm chung của các mô hình 
bài toán (3), (4), (5) mà các tác giả đã nghiên 
cứu đều có các hệ số của phương trình phụ 
thuộc tích phân của hàm cần tìm. Để nghiên 
cứu các bài toán này, chúng ta có thể sử dụng 
các phương pháp biến đổi để chuyển các 
thành phần phụ thuộc tích phân sang vế phải 
của phương trình chuyển các bài toán đang 
xét về dạng bài toán (1). Khi đó có thể nghiên 
cứu sự hội tụ của phương pháp lặp bằng các 
phương trình toán tử và việc tính toán số sử 
dụng các lược đồ sai phân với độ chính xác 
bậc bốn của dạng bài toán (2). Hiển nhiên khi 
thực hiện phép biến đổi, thành phần vế phải 
sẽ trở nên phức tạp hơn vì sẽ chứa thêm thành 
phần hàm ( )f x . Một cách tự nhiên, chúng ta 
cần nghiên cứu phương pháp giải các bài toán 
trên mà không cần chuyển các thành phần 
chứa tích phân sang vế phải. 
Để giải quyết vấn đề trên, trong bài báo này 
chúng tôi trình bày việc xây dựng phương 
pháp lặp tổng quát cho mô hình các bài toán 
có hệ số phương trình chứa thành phần tích 
phân của đạo hàm cần tìm và lược đồ sai phân 
với độ chính xác cấp bốn cho bài toán cấp hai 
dạng tổng quát hơn. 
Bài báo gồm 5 phần, sau Phần giới thiệu là 
Phần 2, trình bày lược đồ sai phân với độ 
chính xác bậc cao; Phần 3, trình bày mô hình 
bài toán biên với hệ số phụ thuộc phiếm hàm 
tích phân và sơ đồ lặp để giải bài toán; Phần 
4, trình bày cá kết quả thực nghiệm và Phần 5 
là phần Kết luận. 
2. Lược đồ sai phân với độ chính xác bậc cao 
Vũ Vinh Quang và Đtg Tạp chí KHOA HỌC & CÔNG NGHỆ ĐH