数学建模 对汽车保险问题的研究文档格式.docx

1、6.支出不变。7.后五年的收支平衡。8.注销人数为自然退保人数和死亡人数之和。9.每一类别的人没有索赔时补贴比例不变。10.注销时每人平均的偿还退回金额不变。11.下一年平均修理费不变。12.近几年物价不变。四、符号说明符号意义本年度的人数实施安全带法规之后第一年的人数各部分费用金额医疗费下降系数基本保险费用取0，1，2，3，4，5，代表年份死亡人数下降比例取0，1，2，3，代表类别总保费说明：其中，下标“总”表示总投保人数或总收入，“续”表示续保人数，“新”表示新投保人数，“注”表示注销人数，“索”表示索赔人数或索赔费用，“死”表示死亡人数或死亡赔偿费用，“修”表示修理费，“医”表示医疗费。

2、五、模型建立与求解（1）首先我们需理清各投保人及其之间的关系，我们想到利用Venn图。图示如下： 图一 Venn分析图（2）根据题意，本年度的统计报表如下：表1 本年度发放的保险单数基本保险费：775元类别没有索赔时补贴比例（%）续保人数新投保人数注销人数总投保人数112807083846201826416653282517648972824017648982401154461138573508760058324114总收入：6182百万元，偿还退回：70百万元，净收入：6112百万元；支出：149百万元，索赔支出：6093百万元，超支：130百万元。表2 本年度的索赔款索赔人数死亡司机人数平

3、均修理费（元）平均医疗费（元）平均赔偿费（元）58275611652102015263195582463233151223123138861158572292947823294170087270138058142321总修理费：1981（百万元），总医疗费：2218（百万元），总死亡倍偿费：1894（百万元），总索赔费：6093（百万元）分析表1，表2可得：下一年0类续保人数=0类索赔人数-0类死亡人数+1类降为0类的人数（1类索赔人数-1类死亡人数）+2类降为0类的人数（2类索赔人数-2类死亡人数），即： （1）下一年1类续保人数=0类升为1类的人数（0类总投保人数-0类索赔人数-0类注销人

4、数+0类死亡人数）+3类降为1类的人数（3类索赔人数-3类死亡人数），即： （2）下一年2类续保人数=1类升为2类的人数（1类续保人数-1类索赔人数-1类注销人数+1类死亡人数），即： （3）下一年3类续保人数=2类升为3类人数（2类续保人数-2类索赔人数-2类注销人数+2类死亡人数）+3类续保人数-3类索赔人数-3类注销人数+3类死亡人数），即： （4）现在分析新投保人数。由人口阻滞模型：可令新投保人数为，其中为时间；为新投保人数的增长率。设为的线性函数，这里为固有增长率。为了确定系数的意义，引入投保人数最大值。因此，新投保人数增长率可表示为： （5）在上式假设下建立阻滞增长模型如下： （6

5、）用分离变量法求解，结果为： （7）带入本年度新投保人数并取，得出，所以新投保人数的阻滞增长模型即为： （8）由上式即可得今后五年的新投保人数如下表：表3 五年的新投保人数第一年第二年第三年第四年第五年396549404144408895411831413632然后分析索赔人数。查阅相关资料可知，在总投保人数中每一个人的索赔次数服从泊松分布，所以他索赔次的概率为： （9）所以，他至少索赔一次的概率 （10）所以总人数中个人索赔的概率为： （11）所以： （12）由此可得到索赔人数与总投保人数之间的关系为： （13）带入数据计算可得到各类索赔人数与总投保人数的关系，见下表所示。 表4 索赔人数与

6、总投保人数之间的关系函数第0类0.4307第1类0.4005第2类0.1057第3类0.0834最后，我们认为死亡司机人数以及注销人数和总投保人数无关，只因为政府实施安全带法规而发生改变。（3）运用上述公式（1），（2），（3），（4），（8），（13），以及本年度此公司的统计报表，我们使用C语言编程运算得到数据，记录至下列表格中： 表5 实施安全带法规后第一年发放的保险单数124381713603164036617698191891411775091294011775108769124321309表6 实施安全带法规后第一年的索赔款减少20%减少40%57404069911221916584

7、0671398998573911811113756584947016784208651488由表5，表6由于本年度0，1，2，3类总注销人数为，所以每人偿还退回的金额为根据假设11，总修理费不变，为1981百万元。附加保费等于总保费减去总收入，且保持不变。由本年度可知： （14） （15） （16） （17） （18） （19） （20） （21） （22）由此，我们可以得到具体的模型： （23） （24）在医疗费下降20%40%的情况下，我们计算出基本保险费的金额为677、623。（4）同以上方法相同，我们利用循环语句可计算出第二，三，四，五年度的预计统计报表，表格如下： 表7 第二年度发放

8、保单数没有索赔补贴比例（%）12538631658008175718411808278798178 表8 第三年度发放保单数12561981665093177097711723628827893 表9 第四年度发放保单数12623811674212177796111816038847614 表10 第五年度发放保单数12688041682436178546711862828874071在医疗费下降20%的情况下，此四年的基本保险费依次为：676、675、 675、 674。在医疗费下降40%的情况下，此四年的基本保险费依次为：622、622、621、621。我们可利用MATLAB分析此结果。

9、图二 本模型两种结果对比图六、模型误差分析由于汽车保险涉及到很多社会因素和人为因素。在模型中不可能完全考虑在内，所以此模型存在误差。1.社会因素引起的误差汽车保险公司收取的保险费受很多因素影响。诸如每年的经济发展情况、车价、交通部门发布的一些政策、各个汽车保险公司之间的竞争等等。这些因素都会影响汽车保险费用。这些因素是我们不能精确考虑在内的，所以此模型与实际有很大的误差。2.建模过程中的误差 此题建立的模型简化了诸多波动因素，如死亡司机人数、注销人数，所得结果可能与现实还有很大的误差，同时我们还做了很多的假设，这些假设也会给模型带来误差。还有在数据优化方面也有所限制，可能略微影响到结果地准确性

10、。且运算本身存在误差。以上为模型无法逃避的误差。七、模型评价及推广 本模型是对基本保险费的计算和估测。模型中，我们利用通过Venn图分析个模型量之间的关系，建立出较为合理的数学关系，同时本模型结合部分实际情况，以实际问题为背景，考虑到了人口阻滞增长的问题和索赔人数与总人数的关系，具有一定得科学性。还有此模型用到了大量的数学模型，具有一定的理论基础，可信度高。但此模型也存在了一些缺点，此模型较为复杂，不能对其进行全方面的考虑，使得结果与实际有一些偏差。还有我们的模型有需改进之处，本题认为保险公司收取的费用对不同的人是一样的。但现实生活中，保险公司应对不同人采取不同的策略。由于每个人驾驶经验，年龄

11、，家庭条件不一。所以他们需要交的汽车保险费应不同。同时不同人的车也是不同的，车的价位应影响到保险的费用，因为一旦出现事故，汽车保险公司就必须索赔，但由于车的价位不同，修理的费用自然就不同，所以保险公司应该根据车的名贵程度来衡量汽车保险费用。本模型可以应用于很多其他多因素决定单变量的问题,如：人口增长问题，其他类型保险业务等等，本模型尤其对费用的估算与预测有很大帮助。参考文献1 吴建国，数学建模案例精编，中国水利水电出版社，20052 李玉泉，责任保险与索赔理赔，人民法院出版社，20023 楼顺天，MATLAB7.x程序设计语言（第二版），西安电子科技大学出版社，20074 谭浩强，C语言程序设

12、计，清华大学出版社，2000附录1C语言程序#include math.hstdio.hstdlib.hvoid year1（）;void year2（）;void year3（）;void year4（）;void year5（）;void money（）;main（） year1（）; year2（）; year3（）; year4（）; year5（）; money（）; printf（The End!n）; getch（）;void year1（）float BR46,SP4047=0,a; int i,j; FILE *outfile,*sp,*br,*s40; float A46=

13、0; float SP46=0; outfile=fopen（BR0.txt,r sp=fopen（SP0.txt br=fopen（BR1.txtw s40=fopen（SP1.txt for（i=0;i4;i+） for（j=0;j6;j+） Aij=0; fscanf（outfile,%f,&a）; Aij=a; BRij=a; SPij=0; fscanf（sp, SPij=a; if（j=5） SP40ij+1=a;else SP40ij=a; BR02=SP01-SP02+SP11-SP12+SP21-SP22; BR12=A03+A02-SP01-A04+SP02+SP31-SP

14、32; BR22=A12-SP11-A14+SP12; BR32=A22-SP21-A24+SP22+A32-SP31-A34+SP32; BR13=0; SP40i2=0.6*SP40i2; SP40i5=0.6*SP40i4; SP40i4=0.8*SP40i4; BR03=416440/（1+0.0827*exp（-0.5）; BRi5=BRi2+BRi3; BRi4=BRi4-0.4*SPi2; SP4001=BR05*（1-exp（-0.4307）; SP4011=BR15*（1-exp（-0.4005）; SP4021=BR25*（1-exp（-0.1057）; SP4031=BR

15、35*（1-exp（-0.0834）; for（j=0; fprintf（br,%8.0f,BRij）;7; fprintf（s40,SP40ij）; fclose（outfile）; fclose（sp）; fclose（br）; fclose（s40）;void year2（）float BR46,a; FILE *outfile,*br,*s20,*s2; float SP47=0;BR2.txt s20=fopen（ s2=fopen（SP2.txt fscanf（s20, BR03=416440/（1+0.0827*exp（-1）; SP01=BR05*（1-exp（-0.4307）; SP11=BR15*（1-exp（-0.4005）; SP21=BR25*（1-exp（-0.1057）; SP31=BR35*（1-exp（-0.0834）; fprintf（s2,SPij）; fclose（s20）; fclose（s2）; void year3（） FILE *outfile,*br,*s20,*s3;BR3.txt s3=fopen（SP3.txt BR22=A12