D ẠY VÀ HỌC PHƯƠNG TRÌNH VI PHÂN TEACHING AND LEARNING LINEAR EQUATIONS OF THE SECOND ORDER CONSTANT FACTOR WITH THE HELP OF MAPLE MATHEMATICAL SOFTWARE TRẦN QUỐC CHIẾN Trường Đại học
Trang 1D ẠY VÀ HỌC PHƯƠNG TRÌNH VI PHÂN
TEACHING AND LEARNING LINEAR EQUATIONS OF THE SECOND ORDER CONSTANT FACTOR WITH THE HELP OF MAPLE
MATHEMATICAL SOFTWARE
TRẦN QUỐC CHIẾN
Trường Đại học Sư phạm, Đại học Đà Nẵng
TRẦN NGỌC VIỆT
HV Cao h ọc khoá 2005 – 2008
TÓM T ẮT
Mục tiêu của bài báo này là viết chương trình toán học bằng phần mềm MAPLE để phân tích quá trình gi ải phương trình vi phân tuyến tính cấp 2 hệ số hằng với điều kiện ban đầu (Bài toán Cauchy) Từ đó áp dụng giải một số bài toán dạng phương trình vi phân c ấp 2
ABSTRACT
This paper presents a new approach to solve differential equations of the second order constant factor with initial conditions (Cauchy problem) with THE help of program written with the Maple software This program is then applied to some linear equations
of the second order
1 Phương trình vi phân tuyến tính cấp 2
Định nghĩa 1.1 Phương trình vi phân tuyến tính cấp 2 có dạng
'' ( ) ' ( ) ( )
y + p x y +q x y= f x , (1.1) trong đó p x q x( ), ( ), ( )f x là các hàm số liên tục
+ Nếu f x( )=0, thì (1.1) gọi là phương trình vi phân tuyến tính thuần nhất có dạng
'' ( ) ' ( ) 0
y + p x y +q x y= (1.2) + Nếu f x( )≠0, thì (1.1) gọi là phương trình vi phân tuyến tính không thuần nhất có
+ Nếu p x q x( ), ( ) là các hằng số thì (1.1) gọi là phương trình vi phân tuyến tính cấp 2
hệ số hằng
Định lý 1.2 Nếu y x y x là 2 nghi1( ), 2( ) ệm độc lập tuyến tính của phương trình vi phân thuần nhất (1.2) thì y=C y x1 1( )+C y x2 2( ) là nghiệm tổng quát của (1.2), trong đó
Trang 2T ẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 5(28).2008
1, 2
C C là 2 hằng số tùy ý
Chứng minh Xem tài liệu tham khảo [4], định lý 5.3, tr.220
Định lý 1.3 Nghiệm tổng quát của phương trình không thuần nhất bằng tổng của
nghiệm tổng quát của phương trình thuần nhất ( )y và một nghiệm riêng nào đó ( )Y của
phương trình không thuần nhất là y y Y= +
Chứng minh Xem tài liệu tham khảo [4], định lý 5.9, tr.227
2 Phương pháp giải phương trình vi phân tuyến tính cấp 2 hệ số hằng
2.1 Phương trình thuần nhất Cho phương trình thuần nhất có dạng
'' ' 0
trong đó p q, là 2 hằng số
Để giải phương trình thuần nhất (2.1) ta thực hiện như sau:
- Giải phương trình đặc trưng: 2
- Căn cứ vào các trường hợp sau để kết luận nghiệm:
+ Nếu phương trình đặc trưng có 2 nghiệm phân biệt k k thì nghi1, 2 ệm tổng quát của phương trình thuần nhất là 1 2
1 k x 2 k x
y=C e +C e , với C C là 2 h1, 2 ằng số tùy ý
+ Nếu phương trình đặc trưng có nghiệm kép k1 = thì nghik2 ệm tổng quát của phương trình thuần nhất là 1
( ) k x
y= C +C x e , với C C là 2 h1, 2 ằng số tùy ý
+ Nếu phương trình đặc trưng có 2 nghiệm phức k1 = +α βi , k2 = −α βi thì nghiệm
tổng quát của phương trình thuần nhất là .
x
y=eα C βx C+ βx , với C C là 1, 2
2 hằng số tùy ý
2.2 Phương trình không thuần nhất
Cho phương trình không thuần nhất có dạng y''+p y '+q y = f x( ), (2.3) trong đó p q, là 2 hằng số và f x( )≠0
Để giải phương trình không thuần nhất (2.3) ta thực hiện như sau:
Bước 1: Tìm nghiệm tổng quát của phương trình thuần nhất ( )y thực hiện như nêu trên
Bước 2: Tìm nghiệm riêng của phương trình không thuần nhất ( )Y thực hiện như sau:
*Trường hợp 1 Vế phải có dạng .
( ) n( ) x
f x =P x eα , P x n( ) là đa thức bậc n
+ Nếu α (α ≠k k1, 2) không phải là nghiệm của phương trình đặc trưng (2.2) thì tìm nghiệm riêng dưới dạng .
( ) x
n
Y =Q x eα , Q x cùng b n( ) ậc với ( )P x n
+ Nếu α (α ≡ hoặc k1 α ≡k2) là nghiệm đơn của phương trình đặc trưng (2.2) thì tìm
Trang 3nghiệm riêng dưới dạng .
n( ) x
Y =x Q x eα , Q x cùng b n( ) ậc với ( )P x n
+ Nếu α (α ≡ =k1 k2) là nghiệm kép của phương trình đặc trưng (2.2) thì tìm nghiệm
n( ) x
Y =x Q x eα , Q x cùng b n( ) ậc với ( )P x n
* Trường hợp 2 Vế phải có dạng .
( ) x( m( ) co s n( ) sin )
f x =eα P x βx+P x βx , trong đó ( ), ( )P x P x là nh m n ững đa thức bậc m n, và β là hằng số
+ Nếu α ± β không phải là nghiệm của phương trình đặc trưng (2.2) thì tìm nghiệm riêng dưới dạng .
( cos sin ), max( , )
x
Y =eα P βx Q+ βx k = n m + Nếu α± β là nghiệm của phương trình đặc trưng (2.2) thì tìm nghiệm riêng dưới
x( cosk k.sin ), max( , )
Y =x eα P βx Q+ βx k = n m
Bước 3: Vậy, nghiệm tổng quát của phương trình đã cho là y y Y= +
3 Chương trình toán học giải phương trình vi phân tuyến tính cấp 2 hệ số hằng
3.1 Lệnh nhập xuất dữ liệu
+ Hàm readstat("<prompt>"): hiện dấu nhắc <prompt>
+ Hàm
trả về dữ liệu nhập
từ bàn phím
print(data1, data2,… ): hiển thị dữ liệu ra màn hình Lưu ý: xâu ký
tự đặt trong dấu ` `
3.2 Xây d ựng thủ tục trong Maple
+ Maple là một ngôn ngữ lập trình hướng thủ tục (procedure) Chế độ thủ tục được thực
hiện bằng cách đóng gói một dãy các lệnh xử lí cùng một công việc vào một thủ tục duy
nhất, sau đó chỉ cần gọi thủ tục này và Maple tự động thực hiện các lệnh có trong chu trình đó một cách tuần tự và sau đó trả lại kết quả cuối cùng
+ Khai báo thủ tục như sau:
Procedure_name:=proc(parameter_sequence)
[local local_sequence]
[global global_sequence]
[options options_sequence]
statements_sequence
end;
Chương trình giải và vẽ đồ thị nghiệm của phương trình vi phân dạng:
'' ' ( )
y + p y +q y= f x , điều kiện đầu y x( 0)= y0, y x'( 0)= y0'
>
>
restart;
proc_eq2:=proc() # Thu tuc giai pt vi phan
Trang 4T ẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 5(28).2008
local
p,q,f,init_con,pt,ptthuannhat,ptdactrung,dactrung,thuannhat ,phuongtrinh,nktt:
global global_pt:
p:=readstat("p="):
q:=readstat("q="):
f:=readstat("f(x)="):
init_con:=readstat("DKdau="):
global_pt:=pt:
print(` -Bài Giai -`);
pt:=diff(y(x),x$2)+p*diff(y(x),x)+q*y(x)=f:pt:
print(`Phuong trinh vi phan tuyen tinh cap 2 la:`);
print(pt);
ptthuannhat:=diff(y(x),x$2)+p*diff(y(x),x)+q*y(x)=0:ptthuan nhat:
print(`Phuong trinh thuan nhat la:`);
print(ptthuannhat);
ptdactrung:=k^2+p*k+q=0:ptdactrung:
print(`Phuong trinh dac trung tuong ung la:`);
print(ptdactrung);
dactrung:=solve({ptdactrung},{k}):dactrung:
print(`Nghiem pt dac trung la:`);
print(dactrung);
thuannhat:=dsolve(ptthuannhat):thuannhat:
print(`Nghiem tong quat cua phuong trinh thuan nhat la:`); print(thuannhat );
nktt:=dsolve(pt):nktt:
print(`Nghiem tong quat cua phuong trinh khong thuan nhat la:`);
print(nktt);
print(` Dieu kien ban dau la:`);
print(init_con);
phuongtrinh:=dsolve({pt,init_con},{y(x)}):phuongtrinh:
print(`Vay, nghiem cua phuong trinh da cho la:`);
print(phuongtrinh);
with(DEtools):
Trang 5DEplot(pt,y(x),-8 10,{[init_con]},y=-8 10,arrows=NONE,stepsize=0.3,linecolor=blue,title=`VE DO THI NGHIEM CUA PT VI PHAN`);
end;
3.3 Lưu và nạp chương trình
#Ket thuc chuong trinh
+ Maple cho phép lưu chương trình để sử dụng một cách thuận tiện như sau:
>
+ Trong file có sử dụng chương trình thực hiện lệnh read fileproc, trong đó fileproc là tên file (có cả đường dẫn) chứa chương trình biên dịch, được tạo bởi lệnh save:
save proc_eq2 "\\thu_tuc.m":
>
>
restart;
>
read "\\thu_tuc.m";
3.4 S ử dụng chương trình
proc_eq2(-2,2,4*exp(x),y(0)=1,D(y)(0)=2);
Ví dụ:Giải và vẽ đồ thị nghiệm của phương trình vi phân sau
'' 2 ' 2 4 x
y − y + y= e , điều kiện đầu y(0)=1, '(0)y =2
>
p=
proc_eq2(); # khai bao chu trinh
q=
-2; # khai bao he so hang p
f(x)=
2; # khai bao he so hang q
DKdau=
4*exp(x); # khai bao ve phai ham f(x)
y(0)=1,D(y)(0)=2; # nhap gia tri DK ban dau
-Bài Giai -Phuong trinh vi phan tuyen tinh cap 2 la:
=
y'' 2 y' 2 y 4 ex
Phuong trinh thuan nhat la:
=
y'' 2 y' 2 y 0
Phuong trinh dac trung tuong ung la:
=
Nghiem pt dac trung la:
, {k = 1 + I} {k = 1 − I}
Nghiem tong quat cua phuong trinh thuan nhat la:
=
y _C1 e x sin x( ) + _C2 e x cos x( )
Trang 6T ẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 5(28).2008
Nghiem tong quat cua phuong trinh khong thuan nhat la:
=
y ex sin x _C2( ) + ex cos x _C1( ) + 4 ex
Dieu kien ban dau la:
, = ( )
Vay, nghiem cua phuong trinh da cho la:
=
y ex sin x( ) − 3 ex cos x( ) + 4 ex
4 Kết luận
Bài báo này đã trình bày các bước cơ bản phương trình vi phân tuyến tính cấp 2
hệ số hằng bằng chương trình toán học Kinh nghiệm cho thấy là nếu viết một chương trình phức tạp thì trước hết ta viết và chạy từng lệnh một để xem kết quả thế nào, khi
thấy kết quả tốt ta nhóm các lệnh trên thành một chương trình hoàn chỉnh
TÀI LI ỆU THAM KHẢO
[1] Trần Quốc Chiến (2008), Giáo trình Phần mềm toán học, ĐH ĐN
[2] Phạm Huy Điển, Đinh Thế Lục, Tạ Duy Phượng (1998), Hướng dẫn thực hành tính toán trên chương trình MAPLE V, NXB Giáo dục, Hà Nội
Trang 7[3] Lê Ngọc Lăng, Tống Đình Quỳ, Nguyễn Đăng Tuấn, Mai Văn Được(2001), Giúp
ôn t ập tốt môn Toán cao cấp, NXB ĐHQG HN, Hà Nội
[4] Nguyễn Đình Trí, Tạ Văn Đĩnh, Nguyễn Hồ Quỳnh (2001), Toán học cao cấp,
NXB Giáo dục, Hà Nội
[5] K Von Bulow supervised by E.S Cheb-Terrab "Equivalence Methods for Second Order Linear Differential Equations." Master's thesis, Faculty of Mathematics, University of Waterloo (2000)
[6] E.S Cheb-Terrab and T Kolokolnikov "First-order ordinary differential equations, symmetries and linear transformations." European Journal of Applied Mathematics
(2003), 14:231-246
[7] J.A Weil, "Recent Algorithms for Solving Second-Order Differential Equations" (2002), available at http://pauillac.inria.fr/algo/seminars/sem01-02/weil.pdf