血样分组检验的数学模型_精品文档Word文档下载推荐.doc
- 文档编号:14439848
- 上传时间:2022-10-22
- 格式:DOC
- 页数:6
- 大小:156KB
血样分组检验的数学模型_精品文档Word文档下载推荐.doc
《血样分组检验的数学模型_精品文档Word文档下载推荐.doc》由会员分享,可在线阅读,更多相关《血样分组检验的数学模型_精品文档Word文档下载推荐.doc(6页珍藏版)》请在冰豆网上搜索。
n人群总数
p先验概率
血样阴性的概率q=1-p
血样检验为阳性(患有某种疾病)的人数为:
z=np
发生概率:
检查次数:
平均总检验次数:
3问题的分析
根据题意,由已知的先验概率是一个很小的数值,我们大可不必要一个一个地检验,为减少检验次数,我们通过一次分组,从而可使检验次数大大减少;
然而通过再一次分组,可使结果进一步优化,从而达到一个更佳的结果.
4模型建立与求解
设总人数为n,已知每人血样阳性的先验概率为p,记血样阴性的概率q=1-p
4.1模型一
设分x组,每组k人(n很大,x能整除n,k=n/x),混合血样检验x次.阳性组的概率为,分组时是随机的,而且每个组的血样为阳性的机率是均等的,阳性组数的平均值为,这些组的成员需逐一检验,平均次数为,所以平均检验次数,一个人的平均检验次数为N/n,记作:
(1)
问题是给定p求k使E(k)最小.
p很小时利用可得
(2)
显然时E(k)最小.因为K需为整数,所以应取和,比较E(K),得到K的最优值,见表1.
P
0.01%
0.1%
1%
2%
5%
K
100
32
10
8
5
E(k)
0.020
0.063
0.196
0.274
0.426
表1一次分组检验结果
图一
当p=0.01%时,可用Maple模拟出的图像如图一,曲线是关于k的图像.
图二
同上法,当p=0.1%时,可用Maple模拟出的图像如图二,曲线是关于k的图像.其它情况我们一样可用其所长Maple模拟出类似的图像.
随着p的增加k减小,E(k)变大.只要E(k)<
1,就应分组.
注:
若不取近似,利用导数可得时E(k)最小.
4.2模型二
当E(k)>
1时,不应分组,即:
用数学软件求解得检查k=2,3,可知当p>
0.307不应分组.
4.3模型三
将第1次检验的每个阳性组再分y小组,每小组m人(y整除k,m=k/y).因为第1次阳性组的平均值为,所以第2次需分小组平均检验次,而阳性小组的概率为(为计算简单起见,将第1次所有阳性组合在一起分小组),阳性小组总数的平均值为,这些小组需每人检验,平均检验次数为,所以平均总检验次数,一个人的平均检验次数为N/n,记作(注意:
n=kx=myx)
(3)
问题是给定p求k,m使E(k,m)最小.
P很小时(3)式可简化为
(4)
对(4)分别对k,m求导并令其等于零,得方程组:
舍去负数解可得:
(5)
且要求k,m,k/m均为整数.经在(5)的结果附近计算,比较E(k,m),得到k,m的最优值,见表2.
700
125
22
14
M
25
11
7
4
E(k,m)
0.0028
0.0161
0.0897
0.131
0.305
表2二次分组检验结果
与表1比较可知,二次分组的效果E(k,m)比一次分组的效果E(k)更好.
4.4模型四(平均概率模型)
患病人数:
组的基数:
每组需要检验的人数
平均检验次数:
阳性血样的分组模型:
可分为x组,每组k人
分组要满足的条件:
其中y为患病人数.
4.1分组人数=患病人数(即:
血样呈阳性的人数)时,通过这样的分组模型可以使检验次数达到最优.
4.2当z>
k(k=n/x)时,一组人不能包括所有的病人数,第一次检验的基数较大.
4.3当z<
k时,检验多一组时组的基数会很大,而且每一组的概率相差无几十年来.
具体例子见附录二
5模型推广
本数学模型也可适用于某人民医院要对某地区的居民是否患有某种病(如乙肝)的检验,并对该地区的病情作一定的预测,从而达到预防和及早治疗的效果.乙肝的血样检验只有阴性、阳性两种情况,我们可用本数学模型切实地解决这个问题.
6模型评价
由于血样的先检概率通常很小,为减少检验次数,我们通过先对检验的人群进行分组,引入阳性组的概率,通过阳性组数的平均值作为桥梁,由于阳性组的人需要全部重新检验,最后可得平均总检验次数,进而得到一个人的平均检验次数的一元函数.
然而我们通过对阳性组人群进行再次分组(即对检验人群进行二次分组),从而得到一个关于两次分组人数二元函数,进而得到更为优化的数学模型.
最后,我们引入平均概率模型,再把血样检验中出现的可能性细化,得到当血样检验为阳性的人数等于分组后每一组的人数时,通过这样的分组模型可以使检验次数达到最优,但是我们未能给出确实的理论证明.
【参考文献】
[1]姜启源,谢金星,叶俊 数学模型(第三版).高等教育出版社.2003.2
[2]姜启源等 数学模型(第三版)习题参考解答.高等教育出版社.2003.2
[3]王沫然MATLAB6.0与科学计算.电子工业出版社.2001.9
[4]魏宗舒概率论与数理统计教程. 高等教育出版社.1982.3
附录【1】
假定阳性血样的人群有6个小组时的Matlab的程序如下:
clear;
clc;
counter=0;
z=input('
请输入病人数'
)
forr1=1:
z
forr2=r1:
z-r1
forr3=r2:
z-r1-r2
forr4=r3:
z-r1-r2-r3
forr5=r4:
z-r1-r2-r3-r4
ifr1+r2+r3+r4+r5==z
[r1,r2,r3,r4,r5]
counter=counter+1;
#计数器
end
end
end
end
end
end
counter#输出计数的结果
输入z的值为10,输出计算结果:
couter=7
附录【2】
1.n=1000,p=1%,分100组
阳性组
阴性组
分组可能情况
概率
检验次数
平均检验次数
1
99
P1=1/42
110
2.619
2
98
P2=4/42
120
11.429
3
97
P3=8/42
130
24.762
96
9
P4=9/42
140
30
95
P5=7/42
150
6
94
P6=5/42
160
19.048
93
P7=3/42
170
12.143
92
P8=2/42
180
8.571
91
P9=1/42
190
4.524
90
P10=1/42
200
4.762
=142.9
个人平均检验次数:
E=N/1000=0.1429
2.n=1000,p=1%,分125组,每组8人
124
123
P1=4/40
141
14.100
122
P2=8/40
149
29.800
121
P3=9/40
157
35.325
P4=7/40
165
28.875
119
P5=5/40
173
21.625
118
P6=3/40
181
13.575
117
P7=2/40
189
9.450
116
P8=1/40
197
4.925
115
P9=1/40
205
5.125
=162.8
E=N/1000=0.1628
3.n=1000,p=1%,分50组,每组20人
P1=1/530
70
0.1321
P2=10/530
1.6981
33
P3=33/530
6.8491
64
P4=64/530
15.6981
84
P5=84/530
23.7736
P6=90/530
28.8679
82
P7=82/530
29.3962
P8=70/530
210
27.7358
54
P9=54/530
230
23.4340
42
P10=42/530
250
19.8113
=177.40
E=N/1000=0.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 血样 分组 检验 数学模型 精品 文档
![提示](https://static.bdocx.com/images/bang_tan.gif)