D ẠY VÀ HỌC CÔNG THỨC XÁC SUẤT BAYES VỚI SỰ TEACHING AND LEARNING BAYES’ FORMULA WITH THE HELP OF THE MAPLE MATHEMATICS SOFTWARE Tr ần Ngọc Việt Trường Cao đẳng GTVT II Ph ạm Văn Tiến
Trang 1D ẠY VÀ HỌC CÔNG THỨC XÁC SUẤT BAYES VỚI SỰ
TEACHING AND LEARNING BAYES’ FORMULA WITH THE HELP
OF THE MAPLE MATHEMATICS SOFTWARE
Tr ần Ngọc Việt
Trường Cao đẳng GTVT II
Ph ạm Văn Tiến
Trường Cao đẳng Công nghệ - Kinh tế và Thủy lợi miền Trung
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 áp d ụng công thức xác suất toàn phần - công thức xác suất Bayes và điều quan tr ọng hơn hết là phải biết vận dụng các kiến thức đó để giải quyết các bài toán, các tình huống cụ thể Từ đó áp dụng giải một số bài toán xác suất dạng nâng cao hơn Như vậy mới
g ọi là nắm vững và hiểu thấu đáo môn học, đồng thời đưa những thành tựu nổi bật của công ngh ệ thông tin để hỗ trợ việc đổi mới phương pháp dạy và học theo chủ trương của Bộ Giáo
d ục & Đào tạo
ABSTRACT
This paper presents a new approach to the computing of probabilities – Bayes’ formula with the help of a program written in the Maple math software It is important that we should know how to apply such knowledge to solving mathematical problems and dealing with practical situations Accordingly, this can be applied to the solution of various probability problems In this way, the subjects are thoroughly mastered so that great IT achievements in support of the teaching and learning innovation issued by the Ministry of Education and Training can be recorded
1 Đặt vấn đề
Từ mấy năm lại đây, môn học Xác suất thống kê đã được giảng dạy rộng rãi ở các trường đại học, cao đẳng khối tự nhiên, kỹ thuật, kinh tế, …Để đáp ứng nhu cầu dạy
và học đó, qua kinh nghiệm giảng dạy nhiều năm môn học này cho các đối tượng khác nhau, chúng tôi nghiên cứu viết chương trình xác suất Bayes và sử dụng nó như một công cụ làm việc thì thường thu được kết quả bất ngờ Thực tế cho thấy cán bộ kỹ thuật
có thể tự mình giải quyết một bài toán nẩy sinh từ thực tiễn công việc
Định nghĩa (xác suất cổ điển)
Cho Ω = { w , , w1 n} là không gian các sự kiện cơ bản Giả sử các sự kiện cơ bản { } w , , w1 { }n
có cùng khả năng xuất hiện, tức có cùng xác suất, nghĩa là
1 ( w ) ( w )n
n
Khi đó với A = { w , , wi1 ik}, ta nói w , , wi1 ik là
ự kiện cơ bản thuận lợi cho
Trang 2( ) k
P A
n
= (với k là các sự kiện ngẫu nhiên cơ bản thuận lợi cho A; n là các
sự kiện ngẫu nhiên cơ bản của không gian mẫu)
Định nghĩa (xác suất có điều kiện)
Giả sử P B ( ) ≠ 0 Xác suất có điều kiện của A khi điều kiện B được định nghĩa bởi:
( / )
( )
P AB
P A B
P B
= (1.1)
M ệnh đề (xác suất tích)
(a)P AB ( ) = P A B P B ( / ) ( ) = P B A P A ( / ) ( ) (1.2) (b) P A A ( 1 2 An) = P A ( 1 / A2 A P An) ( 2 / A3 An) ( P An−1 / A P An) ( n)
Chứng minh Xem tài liệu tham khảo [3], mệnh đề 2.5.1 , tr.34
Định nghĩa (sự kiện độc lập)
a) Hai sự kiện A và B được gọi là độc lập với nhau nếu
P AB ( ) = P A P B ( ) ( )
b) Các s ự kiện A1, , An ( n > 2) được gọi là độc lập từng đôi nếu
P A A ( i k) = P A P A ( ) (i k), 1 ≤ ≠ ≤ i k n
Giả sử A1, , An là xung khắc từng đôi và đầy đủ Khi đó với bất kỳ sự kiện A,
ta có:
1 ( ) ( / ) ( )
n
i
P A P A A P A
=
= ∑ (1.3)
Chứng minh Xem tài liệu tham khảo [3], mệnh đề 2.6.1 , tr.36
Mệnh đề (công thức xác suất Bayes)
Giả sử A1, , An là xung khắc từng đôi và đầy đủ, và A là một sự kiện ngẫu nhiên Khi đó:
( / ) ( ) ( / )
( )
i
P A A P A
P A A
P A
= (1.4)
Chứng minh
Áp dụng công thức (1.1), (1.2) và (1.3) ta có:
Trang 3( ) ( / ) ( ) ( / ) ( )
( / ) ( )
k
P A A P A A P A P A A P A
P A A P A
=
∑
2.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
ký tự đặt trong dấu ` `
2.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 trong
thủ tục đó và 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;
>
>
restart;
proc_xs2:=proc() # Thu tuc xs Toan Phan-Bayes
local
P1,P11,P2,P21,P3,P31,xstoanphan,toanphan,tph,xsbayes,bayes1 ,bayes11,bayes2,bayes21,bayes3,bayes31,bien:global
global_toanphan:
P1:=readstat("P(A[1])="):
P11:=readstat("P(A/A[1])="):
P2:=readstat("P(A[2])="):
P21:=readstat("P(A/A[2])="):
P3:=readstat("P(A[3])="):
P31:=readstat("P(A/A[3])="):
global_toanphan:=toanphan:
print(` -Bài Giai -`);
Trang 4print(`Goi A[i]la Su kien lay san pham cua phan xuong thu i:`);
print(`P(A[1])=`);
print(P1);
print(`P(A[2])=`);
print(P2);
print(`P(A[3])=`);
print(P3);
print(`Goi A la Su kien san pham lay ra la Phe Pham:`); print(`P(A/A[1])=`);
print(P11);
print(`P(A/A[2])=`);
print(P21);
print(`P(A/A[3])=`);
print(P31);
xstoanphan:=P(A)=Sum(P(A[i])*P(A/A[i]),i=1 n):xstoanpha: print(`+Cong thuc xac suat TOAN PHAN:`);
print(xstoanphan);
toanphan:=P(A)=P(A[1])*P(A/A[1])+P(A[2])*P(A/A[2])+P(A[3] )*P(A/A[3]):toanphan:
print(toanphan);
tph:=P(A)=P1*P11+P2*P21+P3*P31:tph:
bien:=P1*P11+P2*P21+P3*P31:bien:
print(`Vay, xac suat cua Bien Co A la:`);
print(tph);
xsbayes:=(P(A[i]/A)=(P(A[i])*P(A/A[i]))/P(A),i=1 n):xsba yes:
print(`+Cong thuc xac suat BAYES:`);
print(xsbayes);
bayes1:=(P(A[1]/A)=(P(A[1])*P(A/A[1]))/P(A)):bayes1:
print(bayes1);
bayes11:=(P(A[1]/A)=(P1*P11)/bien):bayes11:
print(bayes11);
bayes2:=(P(A[2]/A)=(P(A[2])*P(A/A[2]))/P(A)):bayes2:
print(bayes2);
Trang 5bayes21:=(P(A[2]/A)=(P2*P21)/bien):bayes21:
print(bayes21);
bayes3:=(P(A[3]/A)=(P(A[3])*P(A/A[3]))/P(A)):bayes3: print(bayes3);
bayes31:=(P(A[3]/A)=(P3*P31)/bien):bayes31:
print(bayes31);
end;
2.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_xs2 "\\xac_suat.m":
>
>
restart;
>
read "\\xac_suat.m";
2.4 S ử dụng chương trình
proc_xs2(0.36, 0.12, 0.34, 0.1, 0.3, 0.08);
Ví d ụ: Một Công ty có 3 phân xưởng sản xuất cùng một loại sản phẩm Trong
đó phân xưởng I chiếm 36%, phân xưởng II chiếm 34%, phân xưởng III chiếm 30% sản lượng toàn Công ty Biết tỷ lệ phế phẩm tương ứng của các phân xưởng là 12%, 10%, 8%
a)Tính tỷ lệ phế phẩm chung của toàn Công ty ?
b)Lấy ngẫu nhiên từ kho một sản phẩm, biết sản phẩm lấy ra là phế phẩm Tìm xác suất để sản phẩm ấy thuộc phân xưởng I, phân xưởng II và phân xưởng III ?
>
P(A[1])=
proc_xs2();
P(A/A[1])=
0.36;
P(A[2])=
0.12;
P(A/A[2])=
0.34;
P(A[3])=
0.1;
P(A/A[3])=
0.3;
-Bài
Giai -0.08;
Goi A[i]la Su kien lay san pham cua phan xuong thu i:
Trang 60.36
P(A[2])=
0.34
P(A[3])=
0.3
Goi A la Su kien san pham lay ra la Phe Pham:
P(A/A[1])=
0.12
P(A/A[2])=
0.1
P(A/A[3])=
0.08
+Cong thuc xac suat TOAN PHAN:
=
( )
P A ∑
=
i 1
n
( )
P A i
P A
A i
=
( )
P A P A( 1) + +
P A
A1 P A( 2)
P A
A2 P A( 3)
P A
A3 Vay, xac suat cua Bien Co A la:
=
( )
P A 0.1012
+Cong thuc xac suat BAYES:
, =
P
A i
A
( )
P A i
P A
A i
( )
P A i = 1 n
=
P A1
A
( )
P A
1
P A
A1
( )
P A
=
P
A1
A 0.4268774704
=
P
A
2
A
( )
P A
2
P A
A
2
( )
P A
Trang 7=
P
A2
A 0.3359683794
=
P
A3
A
( )
P A3
P A
A3
( )
P A
=
P
A3
A 0.2371541502
3 K ết luận
Bài báo đã trình bày các bước cơ bản bài toán xác suất toàn phần - xác suất Bayes 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 Hy vọng qua chương trình xác suất Bayes sẽ giúp cho người học có thêm phương pháp và tư duy mới về xác suất
[1] Trần Quốc Chiến, Giáo trình Phần mềm toán học, ĐH ĐN, 2008
[2] Phạm Huy Điển, Đinh Thế Lục, Tạ Duy Phượng , 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, 1998
[3] Nguyễn Văn Hộ, Xác suất thống kê, NXB Giáo dục, Hà Nội, 2006
[4] Đặng Hùng Thắng, Bài tập Xác suất thống kê, NXB Giáo dục, Hà Nội, 2006
[5] Trần Quốc Chiến, Trần Ngọc Việt, Tạp chí khoa học & công nghệ số 28-năm 2008, ĐHĐN
[6] Trần Quốc Chiến, Phạm Văn Tiến, Tạp chí khoa học & công nghệ số 31-năm 2009,
ĐHĐN