彩票中的数学建学建模文章.doc

《彩票中的数学建学建模文章.doc》由会员分享，可在线阅读，更多相关《彩票中的数学建学建模文章.doc（18页珍藏版）》请在三一文库上搜索。

1、高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的

2、话）： 所属学校（请填写完整的全名）： 重 庆 大 学 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：彩票中的数学摘要本模型讨论的是如何评判传统型彩票和乐透型彩票的一般评奖方案的合理性问题。本文首先根据彩票中奖规则，利用古典概率求出了这两种类型彩票的各种奖项出现

3、的可能性。把每注彩票中奖与否看成贝努利试验，得到当期销售的n注彩票内中第i项奖的注数服从二项分布。在假设每期彩票的销售量足够多的前提下，由贝努利大数定律归结为正态分布，从而求出了每注彩票的平均收益率。在此基础上，结合公平尺度，利用彩民的博彩心理变化构造了评判方案合理性的判别函数，利用MATLAB6.1软件编程计算，判别出题目所给方案的奖金设置的优劣。并且利用这个判别函数，我们建立了求解最优方案的非线形规划模型。通过求解所建立的模型，找到在给定的彩票销售注数下最优方案，极其奖项和奖金额的设置。本模型可操作性强，它使彩票运作有章可循，在今后类似活动中有科学的指导作用。关键词 彩票 二项分布 期望收

4、益率 博彩心理 判别函数一、 问题的重述近年来“彩票飓风”席卷中华大地，巨额诱惑使越来越多的人加入到“彩民”的行列，目前流行的彩票主要有“传统型”和“乐透型”两种类型。“传统型”采用“10选6+1”方案：先从6组09号球中摇出6个基本号码，每组摇出一个，然后从04号球中摇出一个特别号码，构成中奖号码。投注者从09十个号码中任选6个基本号码（可重复），从04中选一个特别号码，构成一注，根据单注号码与中奖号码相符的个数多少及顺序确定中奖等级。以中奖号码“abcdef+g”为例说明中奖等级，如表一（X表示未选中的号码）。 表一中 奖等 级10 选 6+1（6+1/10） 基 本 号 码 特别号码说

5、明一等奖abcdef g选7中（6+1）二等奖abcdef 选7中（6）三等奖abcdeX Xbcdef 选7中（5）四等奖abcdXX XbcdeX Xxcdef选7中（4）五等奖abcXXX XbcdXX XXcdeX XXXdef 选7中（3）六等奖abXXXX XbcXXX XXcdXX XXXdeX XXXXef 选7中（2）“乐透型”有多种不同的形式，比如“33选7”的方案：先从0133个号码球中一个一个地摇出7个基本号，再从剩余的26个号码球中摇出一个特别号码。投注者从0133个号码中任选7个组成一注（不可重复），根据单注号码与中奖号码相符的个数多少确定相应的中奖等级，不考虑号码

6、顺序。又如“36选6+1”的方案，先从0136个号码球中一个一个地摇出6个基本号，再从剩下的30个号码球中摇出一个特别号码。从0136个号码中任选7个组成一注（不可重复），根据单注号码与中奖号码相符的个数多少确定相应的中奖等级，不考虑号码顺序。这两种方案的中奖等级如表二。 表二中 奖等 级33 选 7（7/33）36 选 6+1（6+1/36）基 本 号 码 特别号码说 明基 本 号 码 特别号码说 明一等奖 选7中（7） 选7中（6+1）二等奖 选7中（6+1） 选7中（6）三等奖选7中（6） 选7中（5+1）四等奖 选7中（5+1） 选7中（5）五等奖选7中（5） 选7中（4+1）六等奖

7、选7中（4+1） 选7中（4）七等奖选7中（4） 选7中（3+1）注：为选中的基本号码； 为选中的特别号码； 为未选中的号码。 以上两种类型的总奖金比例一般为销售总额的50%，投注者单注金额为2元，单注若已得到高级别的奖就不再兼得低级别的奖。现在常见的销售规则及相应的奖金设置方案如表三，其中一、二、三等奖为高项奖，后面的为低项奖。低项奖数额固定，高项奖按比例分配，但一等奖单注保底金额60万元，封顶金额500万元，各高项奖额的计算方法为：(当期销售总额 总奖金比例) -低项奖总额 单项奖比例 （1）根据这些方案的具体情况，综合分析各种奖项出现的可能性、奖项和奖金额的设置以及对彩民的吸引力等因素评

8、价各方案的合理性。（2）设计一种“更好”的方案及相应的算法，并据此给彩票管理部门提出建议。（3）给报纸写一篇短文，供彩民参考。 二、模型的假设及符号说明21 模型的假设 （1） 设彩票定期开奖，本文只考虑一期的奖项;（2） 每注彩票只兑付最高奖级奖金，不可兼得。（3） 假设彩票的规则是以公正公平为原则；（4） 假设彩票的发行费用不计，彩票总奖金比例一般为销售总金额的50%；（5） 假设彩票的宣传工作做得很到位，在某一区域，一段时期内潜在的购买力是固定的；（6） 假设高项奖按给定的百分比分配，且按当期各奖项实际中奖注数平均分配该奖项奖金；（7） 假设彩民大都具有博彩心理。22 符号说明(1) 表

9、示等奖是否被取走，表示等奖没被取走，表示等奖被取走；(2) 表示第等奖中奖注数；(3) 表示第等奖的奖金额；(4) ：表示当期已售出的彩票注数；(5) ：表示取走的奖金额； (6) ：表示一注中奖彩票被取走的奖金比率，即单注彩票的平均收益率；(7) ：表示第等奖的中奖概率；(8) 表示第等奖奖金分配的百分比；(9) 表示每注彩票中第项奖的奖金额；(10) :表示博彩心理函数；(11) ：表示低项奖奖额在总奖金中所占的比例；(12) ：判别方案合理性的判别函数。三、问题的分析3.1 求两种类型彩票每注中各奖项的概率首先，对某一方案而言，每注彩票中各等奖的概率是可求的，分别用来表示中一等奖到七等奖

10、的概率。由古典概率问题求得传统型和乐透型概率如下(共可分为4类)：（1）传统型 （10选6+1）= = = = = = （2）乐透型 (N选M )= = = = = = =（3）乐透型 （N选M+1）= = = = = = =（4）无特殊号型 (N选M)= = = = = 所给的29种方案可归纳为15类，利用上述公式对不同的和，具体概率值计算见附录中的程序（一）3.2 设当期销售注彩票，研究每注彩票的收益率由于当期彩票的总奖金与售出彩票注数有关，我们取不同的n值来研究方案的合理性。因为单注彩票中第等奖服从0-1分布，因此注彩票中获得第等奖的中奖注数随机变量服从二项分布，即由贝努利大数定律可知，

11、当足够大时，近似服从正态分布，其期望与方差分别为记 由的独立性高项奖的第等奖的奖金额由题中给出的计算方法得 表示第等奖是否被取走，服从（0-1）分布，表示第等奖没被取走，表示第等奖被取走，其分布率为 若表示取走的奖金总额，则它是一个随机变量，则可表示单注彩票的收益率。基于彩票的总奖金为销售总额的50%，单注彩票金额为2元。首先分析的数学期望是否接近于1来评价彩民抽奖的收益率，这是因为表示一注彩票的平均收益金额。（1） 求每注彩票的平均收益率假设表示高项奖的项数，表示总奖项数，则 取，利用 （1）求出题目表三中29种方案，可求出不同下每注彩票的平均收益率（见表0），从运算结果看，随着的增大，平均

12、收益率趋于稳定值1。表0收益分析表序号100万150万200万300万12345678910111213141516171819202122232425262728290.48520.60530.57980.54980.50830.61790.67760.54570.48090.67530.34890.42350.41200.37380.39090.41230.33320.37860.26970.45380.48290.30010.03090.47310.40430.35900.37700.17130.33570.55130.65010.62140.59240.63140.69320.6920

13、0.62490.57010.73000.44150.50300.48240.44750.45120.46940.39790.43320.32700.49120.51740.34580.07320.44830.39620.41250.2530. 0.3989 0.39620.59600.68470.65840.63220.74310.77740.76840.71370.67140.78800.55110.59350.56840.53830.52810.54280.47840.50370.40230.54190.56470.40880.13960.55100.49000.44690.45990.3

14、6020.48310.66920.74190.72040.69890.86560.88270.85990.82610.80030.85830.69680.7.740.68560.66330.63620.64680.5873 0.60310.50940.61810.63630.50520.25830.6.350.54950.52570.53300.51380.6004根据收益分析表1，我们得到以下排序方式：10 ，24， 6， 7， 2， 3， 4， 8， 16， 15， 5 ，12， 21， 1， 9， 20， 12， 13， 18， 21， 11， 17， 26， 29， 22， 19， 2

15、8， 23。该顺序反映了这些方案的平均收益，但是，单纯考虑单注奖金的平均收益，并没有综合考虑博彩心理和公平因素，因此这只是一种不全面的排序方式，还需要综合考虑其他因素。（2）求风险概率若销售注彩票，则等奖被取走的概率为，平均有注彩票获一等奖，因此一等奖能取走的金额为，又一等奖单注保底金额60万元，封顶金额500万元，因此,， 记，对任意的风险率的计算结果不超过0.2（见程序一），说明彩票发行有风险的可能性非常小。3.3从公平因素出发进行分析：利用公平原则确定奖金在各个奖项中分配情况，中奖的概率与第项奖单注中奖金额应该具有相关性，即中奖可能性越小，奖金额越高，例如，在打麻将中，和边与和夹的几率相

16、差近2倍，获得报酬也差2倍。再者中奖几率和成本有直接关系，假设一等奖概率是0.001,二等奖概率是0.01,则一等奖平均抽奖1000注（需2000元）才能中一次，而二等奖需要100注（需200元），可见，中奖收益和中奖概率应成反比，最理想的情况为中奖概率和奖金金额乘积等于常数。考虑到收益率（50%）和博彩因素，我们首先将总奖金额分成低项、高项两部分，在每项内部考虑公平因素，若考虑奖项为1到7等奖的方案：， 我们称该公式为公平尺度。若低项奖的百分比为，则高项奖的百分比为，下面仅对L=3，J=7讨论，由低项奖内部公平尺度有以下公式： （2）由高项奖内部公平尺度有以下公式： 则 = （当n是定值时，

17、其值是常数）。下面研究如何刻画公平程度的数学表达式：因为当，同理也反映了低项奖的公平程度。 我们可取可作为衡量公平程度的数学表达式，但由于远小于1, 很小，我们可设，记为总公平因子。3.4博彩心理的因素分析：彩民的兴趣大小与单注收益率、公平性相关。对一个方案，彩民的看法主要是由收益决定的，我们用这个量的某个函数来衡量该方案对彩民的吸引力，而彩民的看法是一个心理因素，是一个很难准确衡量的量。根据心理学的知识，人的心理情况的变化可用博彩心理函数来近似刻画,博彩心理与单注的收益率有关，因此，取博彩心理函数为其中 ，表示销售注数。四、模型的建立与求解综合以上的分析，影响方案合理性的因素有三个：收益期望

18、，公平因子，搏彩因素。由此来构造给定方案合理性的判别函数，并且利用各个因素处于理想状态时，得出一个理想状态函数，通过求解这个函数的最大值可以构造最优方案。（1）评价函数的构造（第一问的解答）方案的吸引力函数=博彩心理函数*公平因子，即： 根据彩票目前每期开奖的有关资料，以及题目中对单注一等奖奖金的约束，取彩票销售注数=200万注进行计算。 将29个方案的有关数据代入判别函数中，求出的值（见表1），通过比较 值的大小。，我们可以确定给定的各种方案的优劣（见排序表，即表2）。 表1 判别函数的值表 表2排序表排序序号12345678910111213141516171819202122232425

19、262728299511871920162410622182512261314231517212922732841序号方案值（）12345678910111213141516171819202122232425262728296+1/106+1/106+1/106+1/107/296+1/297/307/307/307/317/317/327/327/327/337/337/347/347/357/357/357/357/356+1/366+1/367/367/376/405/603.23203.95113.64483.36546.12385.07915.90235.99256.30265.7

20、2356.09145.45775.37355.30074.24475.79724.15645.6235 5.89685.82114.10785.67504.27315.74495.51385.44433.70993.36734.0083 由此得到对29个方案的排序为： 9，5，11，8，7，19，20，16，24，10，6，22，18，25，12，26，13，14，23，15，17，21，29，2，27，3，28，4，1。 经过计算得出了给定方案的固定部分奖金占总额的比例（）：0.0132 0.5181 0.4979 0.5081 0.2393 0.4205 0.5378 0.3074 0.2

21、964 0.6839 0.3577 0.4341 0.4391 0.4368 0.4850 0.4404 0.3130 0.4036 0.1980 0.3941 0.5659 0.2260 0.4013 0.4887 0.3977 0.3009 0.3167 0.1517 0.121429个方案的公平因子：（） 下面的数据反映了对应的方案的公平程度 0.1982 0.4078 0.4015 0.3886 0.5021 0.4223 0.4709 0.4844 0.5049 0.4521 05003 0.4369 0.4314 0.4307 0.4684 0.0022 0.4002 0.4769

22、 0.4751 0.0075 0.4494 0.4513 0.4198 0.4527 0.4384 0.4348 0.4012 0.3984 0.4103（2）第二问的模型及解答由问题一的求解结果可以看出，传统型和无特别号（23号）的方案已不可能列入最优方案中，因此，问题二中只需在乐透型彩票的两种方案（/和/）进行讨论，我们需要求出当，取何值，取几个奖项，高项奖的百分比及低项奖的奖金额为多少时方案最优。 最优方案的目标函数仍然为问题一的判别函数，其中的变量，根据题目的表3给出的数据和问题一的讨论，满足如下的条件：，。，。由公平原则，因此，而最大为60，因此，可约束 。综合以上的分析可得到问题二

23、的求解模型为： 把乐透型的两种情况的概率代入上式，利用MATLAB6.1编程计算得到此非线形规划问题的最优解为：，。 上述结果是当期彩票的注数时，应选择的方案。利用此模型，我们求得了当期彩票的销售注数注时，应选择的最优方案：计算结果注数最优方案()100万6/326.3210.800.110.091403082150万7/316.8090.780.120.102101850250万6+1/325.9340.770.100.132041920300万6/367.0690.720.120.161213560350万7/336.9640.740.150.111902550400万7/355.8750

24、.750.170.082051620500万7/364.8960.800.100.101432810五、模型优缺点的讨论1 优点：本模型全面考虑了彩票的收益、公平性及博彩心理，建立了一个对彩票方案进行评价的模型。结合实例，用本模型可以寻找给定摸彩方法的最优方案，并能对给定方案进行评价。它利用评价函数，运用公平因子，结合收益原则，建立了优化模型，这种方法是本模型最突出的特色。它使奖金的设置有章可循，在今后类似活动中有较科学的指导作用。2 缺点：本模型是对固定的销售注数的基础上进行研究的，分析了给定摸彩方法的分析评价并求其最优方案，没有能够脱离销售注数考虑问题。而且我们没有解决设计不同的方案还会影

25、响到销售注数。通过本模型的适当调节，当然可以处理这种情况。在本模型的两个函数中，我们把概率也看作变量的情况下，运用条件极值的有关知识，即可解决这个问题。六、给报纸写的一篇短文把握机会 理智博彩在高节奏的现代社会，人们的生活节拍加快，紧张的精神状态渴求得到放松。彩票悄然走进人们的生活，正在默默地发挥着心理按摩医生的作用。彩票是建立在机会均等的基础上，公平竞争的娱乐性游戏。它把穷人和富翁的距离变得不再遥远。它将成为社会保障基金多元化来源的一种，可以集中利用社会的闲散资金，“取之于民，用之于民”，用大家乐于支出的钱，办大家希望办成的事。不过值得注意的是,当前人们在对彩票业的认识和相关宣传中存在着诸多

26、似是而非的误区。为了促进这项具有重大意义的事业发展,我们通过对“传统型”和“乐透型”两种彩票的各种奖项出现的可能性、奖项和奖金额的设置的数学分析，给广大彩民如下几点建议：一、 购买彩票，奉献爱心彩票是筹集资金的有效手段，在市场经济条件下，我国要大力发展教育、体育事业，搞好福利事业，完善社会保障体系，这些都需要大量资金，根据民政部门有关规定，彩票销售总额的30%要作为社会福利基金专项专用，福利基金的三分之二用于本地区的社会福利事业，三分之一上交中央或本地政府调剂使用。所以购买一注彩票，既为自己增加了一次发财的机会，也为社会奉献了一片爱心。二、 把握尺度，合理购彩较高的预期收益总和较大的预期风险相

27、联系。投资彩票游戏，并不是获取收益的正常途径。把彩票作为一种消遣，一种娱乐，把握尺度,合理购彩。忌当成职业，透支赊帐博彩。彩民应谨慎理智购彩，每期适当投入几元，既能买个希望梦想回家，也能享受开奖时刻的心动，又不至于因盲目的投入引发“家庭经济危机”。三、 讲究投注的科学性，树立彩民良好形象彩号的出现随意性强，本身具有不可预测，是一种科学刺激、往往能带给人惊喜的游戏。因此，彩票本身就是玩的心动与潇洒。做一个文明彩民，将科学购彩进行到底，把买彩当作一项娱乐，重在参与。至今为止，还没有发现哪一期特等奖的号码是人为地被预先“研究”出来的。事实上，买彩者“心想”未必一定“事成”。购买者一定要理智对待。 从

28、我们对“传统型”和“乐透型”彩票中各种奖项出现的可能性、对彩民的吸引力等因素的综合分析中，可以看出，每注“乐透型”彩票比“传统型”彩票的平均收益率大，而且乐透型彩票的趣味性也很高，它正逐渐成为世界彩票业的主流。周而复始的梦想与希望，带给人们的是一种轻松愉悦的心理体验过程，它能使人们在平时工作和生活中长期绷紧的神经得以放松，使人们在不知不觉中拥有一种快乐的心情。这就是彩票游戏的乐趣。 七、参考文献1 高强 从经济学角度审视博彩现象 陕西经贸学院学报 2001.102 吴珊娜 彩票的运行机理 渭南师范学院学报 2001.63 何文章 数学建模与实验 哈尔滨工程大学出版社 2002.34 何文章 大

29、学数学实验 哈尔滨工程大学出版社 2000.85 母丽华 数学实验 黑龙江科学技术出版社 2002.86 吕盛鸽 概率统计在彩票选号中的应用 统计与决策 2001.87 许乘 概率论与树理统计 哈尔滨工业大学出版社 2002.3 附录： 程序一：为了求各方案各奖项的概率，编写了jiecheng.m , c.m和gl.m文件：%求阶层 jiecheng.m文件function y=jiecheng(x)jiecheng=1;if x=0 jiecheng=1;else for i=1:x jiecheng=jiecheng*i; endendy=jiecheng;%求解组合数 c.m文件func

30、tion y=c(n,m)y=jiecheng(m)/(jiecheng(n)*jiecheng(m-n);%求解概率 gl.mp=zeros(15,7);p(1,:)=1/(5*106) 8/(106) 18/(106) 261/(106) 342/(105) 419.95/(104) 0;p(2,:)=1/c(7,29) c(7-1,7)/c(7,29) c(7-1,7)*c(1,29-7-1)/c(7,29) c(7-2,7)*c(1,29-7-1)/c(7,29) c(7-2,7)*c(2,29-7-1)/c(7,29) c(7-3,7)*c(2,29-7-1)/c(7,29) c(7

31、-3,7)*c(3,29-7-1)/c(7,29);p(3,:)=1/c(6+1,29) c(7,7)*c(1,29-6-1)/c(6+1,29) c(6-1,6)*c(1,29-6-1)/c(6+1,29) c(6-1,6)*c(2,29-6-1)/c(6+1,29) c(6-2,6)*c(2,29-6-1)/c(6+1,29) c(6-2,6)*c(3,29-6-1)/c(6+1,29) c(6-3,6)*c(3,29-6-1)/c(6+1,29);p(4,:)=1/c(7,30) c(7-1,7)/c(7,30) c(7-1,7)*c(1,30-7-1)/c(7,30) c(7-2,7)

32、*c(1,30-7-1)/c(7,30) c(7-2,7)*c(2,30-7-1)/c(7,30) c(7-3,7)*c(2,30-7-1)/c(7,30) c(7-3,7)*c(3,30-7-1)/c(7,30);p(5,:)=1/c(7,31) c(7-1,7)/c(7,31) c(7-1,7)*c(1,31-7-1)/c(7,31) c(7-2,7)*c(1,31-7-1)/c(7,31) c(7-2,7)*c(2,31-7-1)/c(7,31) c(7-3,7)*c(2,31-7-1)/c(7,31) c(7-3,7)*c(3,31-7-1)/c(7,31);p(6,:)=1/c(7,

33、32) c(7-1,7)/c(7,32) c(7-1,7)*c(1,32-7-1)/c(7,32) c(7-2,7)*c(1,32-7-1)/c(7,32) c(7-2,7)*c(2,32-7-1)/c(7,32) c(7-3,7)*c(2,32-7-1)/c(7,32) c(7-3,7)*c(3,32-7-1)/c(7,32);p(7,:)=1/c(7,33) c(7-1,7)/c(7,33) c(7-1,7)*c(1,33-7-1)/c(7,33) c(7-2,7)*c(1,33-7-1)/c(7,33) c(7-2,7)*c(2,33-7-1)/c(7,33) c(7-3,7)*c(2,

34、33-7-1)/c(7,33) c(7-3,7)*c(3,33-7-1)/c(7,33);p(8,:)=1/c(7,34) c(7-1,7)/c(7,34) c(7-1,7)*c(1,34-7-1)/c(7,34) c(7-2,7)*c(1,34-7-1)/c(7,34) c(7-2,7)*c(2,34-7-1)/c(7,34) c(7-3,7)*c(2,34-7-1)/c(7,34) c(7-3,7)*c(3,34-7-1)/c(7,34);p(9,:)=1/c(7,35) c(7-1,7)/c(7,35) c(7-1,7)*c(1,35-7-1)/c(7,35) c(7-2,7)*c(1,

35、35-7-1)/c(7,35) c(7-2,7)*c(2,35-7-1)/c(7,35) c(7-3,7)*c(2,35-7-1)/c(7,35) c(7-3,7)*c(3,35-7-1)/c(7,35);p(10,:)=1/c(7,35) c(7-1,7)*c(1,35-7)/c(7,35) c(7-2,7)*c(2,29-7)/c(7,35) c(7-3,7)*c(3,29-7)/c(7,35) c(7-4,7)*c(4,29-7)/c(7,35) 0 0;p(11,:)=1/c(6+1,36) c(7,7)*c(1,36-6-1)/c(6+1,36) c(6-1,6)*c(1,36-6-

36、1)/c(6+1,36) c(6-1,6)*c(2,36-6-1)/c(6+1,36) c(6-2,6)*c(2,36-6-1)/c(6+1,36) c(6-2,6)*c(3,36-6-1)/c(6+1,36) c(6-3,6)*c(3,36-6-1)/c(6+1,36);p(12,:)=1/c(7,36) c(7-1,7)/c(7,36) c(7-1,7)*c(1,36-7-1)/c(7,36) c(7-2,7)*c(1,36-7-1)/c(7,36) c(7-2,7)*c(2,36-7-1)/c(7,36) c(7-3,7)*c(2,36-7-1)/c(7,36) c(7-3,7)*c(3

37、,36-7-1)/c(7,36);p(13,:)=1/c(7,37) c(7-1,7)/c(7,37) c(7-1,7)*c(1,37-7-1)/c(7,37) c(7-2,7)*c(1,37-7-1)/c(7,37) c(7-2,7)*c(2,37-7-1)/c(7,37) c(7-3,7)*c(2,37-7-1)/c(7,37) c(7-3,7)*c(3,37-7-1)/c(7,37);p(14,:)=1/c(6,40) c(6-1,6)/c(6,40) c(6-1,6)*c(1,40-6-1)/c(6,40) c(6-2,6)*c(1,40-6-1)/c(6,40) c(6-2,6)*c

38、(2,40-6-1)/c(6,40) c(6-3,6)*c(2,40-6-1)/c(6,40) c(6-3,6)*c(3,40-6-1)/c(6,40);p(15,:)=1/c(5,60) c(5-1,5)/c(5,60) c(5-1,5)*c(1,60-5-1)/c(5,60) c(5-2,5)*c(1,60-5-1)/c(5,60) c(5-2,5)*c(2,60-5-1)/c(5,60) c(5-3,5)*c(2,60-5-1)/c(5,60) c(5-3,5)*c(3,60-5-1)/c(5,60);p程序二：%求解出风险概率， fengxian.m文件pa=zeros(4,10);p

39、a1=zeros(20,10);pa2=zeros(3,10);pa3=zeros(1,10);n=10000 20000 100000 200000 500000 1000000 2000000 5000000 10000000 20000000;p=0.0000002 0.000009 0.000018 0.000261 0.00342 0.0423 0; 0.0000002 0.000009 0.000018 0.000261 0.00342 0.0423 0; 0.0000002 0.000009 0.000018 0.000261 0.00342 0.0423 0; 0.000000

40、2 0.000009 0.000018 0.000261 0.00342 0.0423 0;v=0.50 0.20 0.30 50 0 0 0; 0.60 0.20 0.20 300 20 5 0; 0.65 0.15 0.20 300 20 5 0; 0.70 0.15 0.15 300 20 5 0; (由于篇幅原因，此处的5个矩阵及输出部分略去：)p1=; v1=; p2=; v2=; p3=; v3=1.00 2000 20 4 2 0 0;for i=1:10 for j=1:4 t=0; r=0; for k=4:7 t=t+v(j,k)*n(i)*p(j,k); r=r+v(j,k)2*n(i)*(n(i)-1)*