2013深圳杯数学建模D题

合集下载

2013年“深圳杯”数学建模夏令营

2013年“深圳杯”数学建模夏令营

如果事件A、B互斥,那么()()()P A B P A P B +=+如果事件A、B相互独立,那么(.)().()P A B P A P B =如果事件A在一次试验中发生的概率是P,那么n 次独立重复试验中恰好发生k 次的概率是()(1)k k n kn n P k C P P -=-第Ⅰ卷(选择题 共60分)一:选择题(本大题共12小题,每小题5分,共60分。

在每小题给出的四个选项中,只有一项是符合题目要求的。

)1、已知集合{}b P ,0=,{}Z x x x x Q ∈<-=,032,若φ≠⋂Q P ,则b 等于( )A 、1B 、2C 、1或2D 、82、 3、4、在一个口袋中装有5个白球和3个黑球,这些球除颜色外完全相同。

从中摸出3个球,至少摸到2个黑球的概率等于 (A )27 (B )38 (C )37 (D )9285、.已知二面角βα--l 的大小为030,m n 、为异面直线,m n αβ⊥⊥且,,m n 则、所成的角为 ( ) (A )030 (B )060 (C )090 (D )01206、设q p ,是简单命题,则“p 且q 为假”是“p 或q 为假”的 ( ) A 、必要不充分条件 B 、充分不必要条件C 、充要条件D 、既不充分也不必要条件7、8、9、 10、 11、12、二:填空题(本大题共4小题,每小题4分,共16分。

把答案填在题中横线上)13)不等式01>-xx的解集为 14)已知⎩⎨⎧≥+-≤-05302y x y x ,则22-+y x 的最大值是15)设常数0a >,42ax⎛ ⎝展开式中3x 的系数为32,则2l i m ()nn a a a →∞++⋅⋅⋅+=_____。

16、三:解答题(本大题共6小题,共74分。

解答应写出文字说明、证明过程或演算步骤) 17、18、19、(12分)20、(12分)21、[参考答案]一:选择题(本大题共12小题,每小题5分,共60分。

2012深圳杯数学建模夏令营D题

2012深圳杯数学建模夏令营D题

程序问题一程序1A=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',1,'b2:c611'); B=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',2,'b2:c788'); C=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',3,'b2:c271'); D=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',4,'b2:c213'); E=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',5,'b2:c96'); F=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',6,'b2:c35'); G=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',7,'b2:c21'); H=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',8,'b2:c7');I=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',9,'b2:c11'); J=xlsread('F:\数学建模第二期培训\第三题\新建文件夹 (2)\data',10,'b2:c30'); Ax=A(:,1);Ay=A(:,2);Bx=B(:,1);By=B(:,2);Cx=C(:,1);Cy=C(:,2);Dx=D(:,1);Dy=D(:,2);Ex=E(:,1);Ey=E(:,2);Fx=F(:,1);Fy=F(:,2);Gx=G(:,1);Gy=G(:,2);Hx=H(:,1);Hy=H(:,2);Ix=I(:,1);Iy=I(:,2);Jx=J(:,1);Jy=J(:,2);plot(Ax,Ay,'b+',Bx,By,'rh',Cx,Cy,'g*',Dx,Dy,'cd',Ex,Ey,'mo',Fx,Fy,'yp ',Gx,Gy,'kx',Hx,Hy,'b+',Ix,Iy,'b+',Jx,Jy,'rx')legend('A型孔','B型孔','C型孔','D型孔','E型孔','F型孔','G型孔','H型孔','I 型孔','J型孔')title('各种类型孔的分布图')xlabel('x')ylabel('y')grid on %画出分格线plot(Ax,Ay,'b+',Bx,By,'rh',Cx,Cy,'g*',Dx,Dy,'cd',Ex,Ey,'mo',Fx,Fy,'yp ',Gx,Gy,'kx',Hx,Hy,'b+',Ix,Iy,'b+',Jx,Jy,'rx')legend('A型孔','B型孔','C型孔','D型孔','E型孔','F型孔','G型孔','H型孔','I 型孔','J型孔')title('各种类型孔的分布图')xlabel('x')ylabel('y')grid on %画出分格线程序2m=10;Alpha=1;Beta=5;Rho=0.1;NC_max=200;Qx100;%为使程序运行速度更快,取蚂蚁数为10function[R_best,L_best,L_ave,Shortest_Route,Shortest_Length]=yiqunsuanfa(C1,N C_max,m,Alpha,Beta,Rho,QX)%% 主要符号说明%% C1 n个城市的坐标,n×2的矩阵%% NC_max 最大迭代次数%% m 蚂蚁个数%% Alpha 表征信息素重要程度的参数%% Beta 表征启发式因子重要程度的参数%% Rho 信息素蒸发系数%% QX 信息素增加强度系数%% R_best 各代最佳路线%% L_best 各代最佳路线的长度%%第一步:变量初始化n=size(C1,1);%*表示问题的规模(城市个数)D=zeros(n,n);%D表示完全图的赋权邻接矩阵for i=1:nfor j=1:nif i~=jD(i,j)=((C1(i,1)-C1(j,1))^2+(C1(i,2)-C1(j,2))^2)^0.5;elseD(i,j)=eps;endD(j,i)=D(i,j);endendEta=1./D; %Eta为启发因子,这里设为距离的倒数Tau=ones(n,n);%Tau为信息素矩阵Tabu=zeros(m,n);%存储并记录路径的生成NC=1;%迭代计数器R_best=zeros(NC_max,n);%各代最佳路线L_best=inf.*ones(NC_max,1);%各代最佳路线的长度L_ave=zeros(NC_max,1);%各代路线的平均长度while NC<=NC_max%停止条件之一:达到最大迭代次数%%第二步:将m只蚂蚁放到n个城市上Randpos=[];for i=1:(ceil(m/n))Randpos=[Randpos,randperm(n)];endTabu(:,1)=(Randpos(1,1:m))';%%第三步:m只蚂蚁按概率函数选择下一座城市,完成各自的周游for j=2:nfor i=1:mvisited=Tabu(i,1:(j-1));%已访问的城市J=zeros(1,(n-j+1));%待访问的城市P=J;%待访问城市的选择概率分布Jc=1;for k=1:nif length(find(visited==k))==0J(Jc)=k;Jc=Jc+1;endend%下面计算待选城市的概率分布for k=1:length(J)P(k)=(Tau(visited(end),J(k))^Alpha)*(Eta(visited(end),J(k))^Beta);endP=P/(sum(P));%按概率原则选取下一个城市Pcum=cumsum(P);Select=find(Pcum>=rand);to_visit=J(Select(1));Tabu(i,j)=to_visit;endendif NC>=2Tabu(1,:)=R_best(NC-1,:);end%第四步:记录本次迭代最佳路线L=zeros(m,1);for i=1:mR=Tabu(i,:);for j=1:(n-1)L(i)=L(i)+D(R(j),R(j+1));endL(i)=L(i)+D(R(1),R(n));endL_best(NC)=min(L);pos=find(L==L_best(NC));R_best(NC,:)=Tabu(pos(1),:);L_ave(NC)=mean(L);NC=NC+1%第五步:更新信息素Delta_Tau=zeros(n,n);for i=1:mfor j=1:(n-1)Delta_Tau(Tabu(i,j),Tabu(i,j+1))=Delta_Tau(Tabu(i,j),Tabu(i,j+1))+QX/ L(i);endDelta_Tau(Tabu(i,n),Tabu(i,1))=Delta_Tau(Tabu(i,n),Tabu(i,1))+QX/L(i) ;endTau=(1-Rho).*Tau+Delta_Tau;%%第六步:禁忌表清零Tabu=zeros(m,n);end%%第七步:输出结果Pos=find(L_best==min(L_best));Shortest_Route=R_best(Pos(1),:);Shortest_Length=L_best(Pos(1));subplot(1,2,1)DrawRoute(C1,Shortest_Route)subplot(1,2,2)plot(L_best)hold onplot(L_ave)程序3function DrawRoute(C1,R)N=length(R);scatter(C1(:,1),C1(:,2));hold onplot([C1(R(1),1),C1(R(N),1)],[C1(R(1),2),C1(R(N),2)])hold onfor ii=2:Nplot([C1(R(ii-1),1),C1(R(ii),1)],[C1(R(ii-1),2),C1(R(ii),2)])hold onendtitle('旅行商问题优化结果')程序4%求总的路程路线和时间function [Tabu montime]=sj5(M,X,Y,p,Ax,Ay,Bx,By,Cx,Cy,Dx,Dy,Ex,Ey,Fx,Fy,Gx,Gy,Hx,Hy,Ix, Iy,Jx,Jy)n=size(M,1);Tabu=1;tk=3;jk=1;for i=1:n-1tempt=[];for k=1:nif isempty(find(Tabu==k))&&length(find(Tabu==p(2, k)))>0, tempt=[tempt k];endendif length(tempt)==0temptendsum1=inf;for kk=temptif sum1>X(tk,kk)ik=kk;sum1=X(tk,ik);endendTabu=[Tabu ik];tk=ik;endTabu(1)=1;time=0;mon=0;tx=Tabu(1:end-1);ty=Tabu(2:end);for i=1:n-1mon=mon+X(tx(i),ty(i));if p(1,tx(i))~=p(1,ty(i))zty=abs(p(1,tx(i))-p(1,ty(i)));if zty<4st=zty;elsest=8-zty;endif Y(tx(i),ty(i))<18*st;time=time+18*zty;elsetime=time+Y(tx(i),ty(i));endelsetime=time+Y(tx(i),ty(i));endendDrawRoute(M,Tabu,Ax,Ay,Bx,By,Cx,Cy,Dx,Dy,Ex,Ey,Fx,Fy,Gx,Gy,Hx,Hy,Ix,I y,Jx,Jy)。

2013年全国大学生数学建模竞赛专科组D题

2013年全国大学生数学建模竞赛专科组D题

公共自行车服务系统运行规律研究-----以浙江省温州市鹿城区为例摘要:根据浙江省温州市鹿城区公共自行车使用数据,基于公共自行车服务模式和使用规则,对公共自行车服务系统的运行规律从各个站点运行规律、借车人借车规律和高峰日系统的具体运行规律三个方面进行研究。

并基于所发现的规律对目前公共自行车服务系统的设置进行了评价,同时给出了提高系统效率的建议。

首先分别统计各站点20天中每天及累计的借车频次和还车频次。

结果显示,几乎所有站点在20天中的累计借车频次排序和还车频次排序大致相同。

另外也对每次用车时长的分布情况也进行了分析,发现用车时长主要集中在30分钟之内,约占用车量的90%;而用车时长在60分钟之内的用车量占总用车量的99%。

这与鹿城区公共自行车租赁收费标准(1小时内免费,超过1小时收费)相吻合。

以上现象也间接说明了政府制定的收费机制的有效性。

然后统计了20天中各天使用公共自行车的不同借车人数量,并统计了数据中出现过的每张借车卡累计借车次数的分布情况。

结果表明,20天中借车人数呈现明显的周期性波动,并且波动周期大约为7天。

而且波动量非常大,从最高峰第20天约20000个借车人到最低峰不到5000人。

还发现每张借车卡累计借车次数的分布情况如下:20天中55%的借车人借车次数不超过10次,35%的在10-30次之间,9%的在30次以上。

最后基于用车次数最多的第20天的数据对系统进行更为详细的分析。

发现借还车站点之间最长距离为32号站到45号站的距离,最短为73号站到115号站,99号站到150号站的距离。

借还车频次最高的站点分别为42号站点(770次)和56号站点(743次)。

两个站点借还车主要时间段相同,98.6%的借车在7:00到21:00之间,99.9%的还车在7:00到21:30之间。

发现两个站点运行规律的不同点有:56号站点比42号站点在早晨上班时的用车高峰更明显。

56号站点比45号站点用车时长在30分钟之内的多3.3%,说明还车高峰期人们会较快使用完自行车。

13年研究生数学建模D题

13年研究生数学建模D题
IAQI P IAQPHi IAQPL 0 (CP BPLo ) IAQPL 0 BPHi BPLo
(5.1)
式中: IAQI P ——污染物项目 P 的空气质量分指数;
CP ——污染物项目 P 的质量浓度值; BPHi ——文献[2]中表 1 中与 CP 相近的污染物浓度限值的最高值; BPLo ——文献[2]中表 1 中与 CP 相近的污染物浓度限值的最低值;
相关系数ρXY 取值在-1 到 1 之间,ρXY = 0 时,称 X,Y 不相关;|ρXY| = 1 时,称 X,Y 完全相关,此时,X,Y 之间具有线性函数关系; |ρXY| < 1 时,X 的变动引起 Y 的部分变动, ρXY 的绝对值越大, X 的变动就会引起 Y 的变动越大, 当|ρXY|> 0.8 时称为高度相关,当|ρXY| < 0.3 时,称为低度相关,其他为中 度相关。 由表 5.1 可以看出,在 10 个监测点的相关系数的均值中,PM2.5 的浓度含 量与 CO 的浓度含量的相关系数为 0.801, 因此两者高度相关; 与 SO2 、 NO2、 PM10 的相关系数分别为 0.664、0.697、0.762,因此 PM2.5 与三者分别中度相关; 与 O3 相关系数为-0.327,因此定义两者低度相关。ຫໍສະໝຸດ 4.4 对问题 4 的分析
灰箱模型只能针对已发生的类似情况进行预测,因此针对突发状况,即某 监测点的 PM2.5 浓度突然增至数倍,我们可以将此突发情况简化成已知边界条 件的点源扩散模型,将 PM2.5 监测点浓度最高点作为点源,突发状况发生前各 监测点的浓度作为初始条件,并考虑当天的天气影响因素。最后按环保部新修 订的《环境空气质量标准》给出重度污染和可能安全区域。
二、问题的假设

2013年全国研究生数学建模竞赛D题

2013年全国研究生数学建模竞赛D题

2013年全国研究生数学建模竞赛D 题空气中PM2.5问题的研究大气为地球上生命的繁衍与人类的发展提供了理想的环境。

它的状态和变化,直接影响着人类的生产、生活和生存。

空气质量问题始终是政府、环境保护部门和全国人民关注的热点问题。

2013年7月12日《中国新闻网》记者周锐报道:“2013年初以来,中国发生大范围持续雾霾天气。

据统计,受影响雾霾区域包括华北平原、黄淮、江淮、江汉、江南、华南北部等地区,受影响面积约占国土面积的1/4,受影响人口约6亿人”(中国国家发展和改革委员会(发改委)2013年7月11日公布在官方网站上的一份报告披露了上述信息,中新社北京7月11日电)。

对空气质量监测,预报和控制等问题,国家和地方政府均制定了相应政策、法规和管理办法。

2012年2月29日,环境保护部公布了新修订的《环境空气质量标准》 (GB3095—2012)[1],本次修订的主要内容:调整了环境空气功能区分类,将三类区并入二类区;增设了颗粒物(粒径小于等于2.5μm)浓度限值和臭氧8小时平均浓度限值;调整了颗粒物(粒径小于等于10μm)、二氧化氮、铅和苯并(a)芘等的浓度限值;调整了数据统计的有效性规定。

与新标准同步还实施了《环境空气质量指数(AQI)技术规定(试行) 》 (HJ633—2012)[2]。

新标准将分期实施,京津冀、长三角、珠三角等重点区域以及直辖市和省会城市已率先开始实施并发布AQI(Air Quality Index);今年113个环境保护重点城市和国家环保模范城市也已经实施;到2015年所有地级以上城市将开始实施;2016年1月1日,将在全国实施新标准。

上述规定中,启用空气质量指数AQI 作为空气质量监测指标,以代替原来的空气质量监测指标――空气污染指数API (Air Pollution Index)。

原监测指标API 为无量纲指数,它的分项监测指标为3个基本指标(二氧化硫2SO 、二氧化氮2NO 和可吸入颗粒物PM10)。

2013数模D题-附件2统计数据

2013数模D题-附件2统计数据

区名年年平均降水量年平均最高气温年平均最低气温年平均风速A2002 1.0520.38.44 6.21 A2003 2.2618.628.47 5.98 A2004 1.8420.118.56 6.14 A2005 1.5819.848.53 5.58 A2006 1.3420.319.33 5.43 A2007 1.3920.239.12 5.02 A2008 1.4320.038.85 4.9 A2009 2.0519.858.67 4.79 A2010 1.8219.558.83 4.76 A2011 1.5719.198.83 4.54区名年年平均降水量年平均最高气温年平均最低气温年平均风速B2002 1.4420.8210.5 6.1 B2003 2.8318.8510.24 5.28 B2004 2.2920.7710.75 5.33 B2005 2.3719.9210.61 5.18 B2006 1.7220.8911.71 5.42 B2007 1.9620.7911.99 5.08 B2008 1.4320.4411.69 4.95 B2009 1.820.3511.54 4.96 B2010 2.1920.4411.72 5.03 B2011 1.7119.5211.26 4.6区名年年平均降水量年平均最高气温年平均最低气温年平均风速C2002 1.5120.3210.07 5.96 C2003 3.1318.3710.01 6.28 C2004 2.119.810.57 6.22 C2005 1.5619.2810.33 5.78 C2006 1.4420.2411.36 5.56 C2007 1.4620.0511.18 5.46 C2008 1.5319.7910.4 5.39 C2009 1.6519.6910.62 5.17 C2010 1.5519.5710.49 5.3 C2011 2.419.0410.4 5.01区名年年平均降水量年平均最高气温年平均最低气温年平均风速D2002 2.421.1411.48 4.5 D2003 3.1319.5210.81 4.38 D2004 2.321.0811.27 4.34 D2005 2.9720.2410.92 4.07 D2006 1.8821.4111.97 4.14 D2007 2.0721.311.99 3.8 D2008 2.2920.811.53 3.64 D2009 2.4120.4211.91 3.65 D2010 3.3120.8111.89 3.81 D2011 2.120.5811.89 3.67区名年年平均降水量年平均最高气温年平均最低气温年平均风速E2002 1.6520.6910.47 4.28 E2003 3.7318.9110.11 4.28E2004 2.8920.4510.39 4.6 E2005 2.3119.710.04 4.3 E2006 1.7920.5610.9 4.14 E2007 2.2820.4110.84 4.08 E2008 2.2419.710.16 4.11 E2009 2.042010.17 3.98 E2010 2.2620.2610.33 4.14 E2011 2.3319.599.86 3.75区名年年平均降水量年平均最高气温年平均最低气温年平均风速F2002120.810.3 5.3 F2003 2.519.210 5.2 F2004220.510.4 5.3 F2005 1.719.910.45 F2006 1.520.911.35 F2007 1.220.911.2 4.7 F2008 1.620.510.8 4.9 F2009 1.520.310.7 4.9 F2010 1.820.410.9 4.8 F2011 1.719.610.4 4.6区名年年平均降水量年平均最高气温年平均最低气温年平均风速G2002 3.621.112.1 5.4 G2003419.611.2 5.3 G2004 2.921.412.1 5.5 G2005 3.920.111.8 4.9 G2006 2.921.212.5 4.9 G2007 4.221.412.9 4.7 G2008 3.720.712.2 4.7 G2009 2.620.412.2 4.7 G2010320.612.4 4.6 G2011 2.120.211.7 4.3区名年年平均降水量年平均最高气温年平均最低气温年平均风速H2002 1.6421.1411 5.1 H2003 2.6119.4910.29 5.03 H2004 2.120.9710.76 5.2 H2005 2.220.110.6 4.8 H2006221.211.4 4.7 H2007 1.821.211.6 4.6 H2008220.911.1 4.9 H2009 2.320.811.1 4.8 H2010 1.820.911.1 4.9 H2011 2.120.110.8 4.6区名年年平均降水量年平均最高气温年平均最低气温年平均风速I2002 1.921.411.6 3.9 I2003 3.819.511.2 3.4 I2004 2.921.412 3.2 I2005 2.820.411.7 2.7 I2006 2.321.412.5 2.8I2007 2.821.312.63 I2008 2.320.811.93 I2009 2.221.211.9 3.1 I2010 2.321.412 3.3 I2011 2.120.811.8 2.5区名年年平均降水量年平均最高气温年平均最低气温年平均风速J2002 2.521.111.8 4.6 J2003 4.119.511.1 4.5 J2004 2.221.311.6 4.6 J2005420.111.24 J2006 3.121.312 3.9 J2007 3.321.112.2 3.5 J2008 2.920.511.6 3.6 J2009 2.520.711.9 3.6 J2010 2.420.912 3.6 J2011220.411.7 3.1冰雹次数221110冰雹次数1111120冰雹次数1111111冰雹次数1111121冰雹次数111211冰雹次数11111112冰雹次数111111冰雹次数1121111冰雹次数1111210冰雹次数11111111。

2013全国数学建模竞赛题目A-B

2013全国数学建模竞赛题目A-B

2013高教社杯全国大学生数学建模竞赛题目(请先阅读“全国大学生数学建模竞赛论文格式规范”)A题车道被占用对城市道路通行能力的影响车道被占用是指因交通事故、路边停车、占道施工等因素,导致车道或道路横断面通行能力在单位时间内降低的现象。

由于城市道路具有交通流密度大、连续性强等特点,一条车道被占用,也可能降低路段所有车道的通行能力,即使时间短,也可能引起车辆排队,出现交通阻塞。

如处理不当,甚至出现区域性拥堵。

车道被占用的情况种类繁多、复杂,正确估算车道被占用对城市道路通行能力的影响程度,将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车站等提供理论依据。

视频1(附件1)和视频2(附件2)中的两个交通事故处于同一路段的同一横断面,且完全占用两条车道。

请研究以下问题:1.根据视频1(附件1),描述视频中交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。

2.根据问题1所得结论,结合视频2(附件2),分析说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异。

3.构建数学模型,分析视频1(附件1)中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。

4.假如视频1(附件1)中的交通事故所处横断面距离上游路口变为140米,路段下游方向需求不变,路段上游车流量为1500pcu/h,事故发生时车辆初始排队长度为零,且事故持续不撤离。

请估算,从事故发生开始,经过多长时间,车辆排队长度将到达上游路口。

附件1:视频1附件2:视频2附件3:视频1中交通事故位置示意图附件4:上游路口交通组织方案图附件5:上游路口信号配时方案图注:只考虑四轮及以上机动车、电瓶车的交通流量,且换算成标准车当量数。

附件3视频1中交通事故位置示意图附件4附件5上游路口信号配时方案本题附件1、2的数据量较大,请竞赛开始后从竞赛合作网站“中国大学生在线”网站下载:试题专题页面:/service/jianmo/index.shtml试题下载地址:/service/jianmo/sxjmtmhb/2013/0525/969401.shtml2013高教社杯全国大学生数学建模竞赛题目(请先阅读“全国大学生数学建模竞赛论文格式规范”)B题碎纸片的拼接复原破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。

2013年夏令营深圳杯数学建模D题

2013年夏令营深圳杯数学建模D题

2013深圳夏令营数学建模承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完明白, 在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题.我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出.我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性.如有违反竞赛规则的行为,我们将受到严肃处理.我们参赛选择的题号是(从题中选择一项填写): D 题所属学校:贵州民族大学参赛队员:1.姓名:2.姓名:3.姓名:指导教师或指导教师组负责人(打印并签名):2013深圳夏令营数学建模编号专用页赛区评阅编号(由赛区组委会评阅前进行编号):全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):自然灾害保险问题的研究摘要本文讨论了如何建立合理的自然灾害保险的问题,并建立了关于P省农业保险公司盈亏的简洁数学模型。

引入。

原理对于问题一,对于问题二,对于问题三,关键词:费率、线性回归、农业保险一、问题重述根据2013年3月5日《环球时报》转摘美国《商业周报》的相关报道,“在2012年全世界发生的10大自然灾害中,有4场是发生在中国。

包括3场严重的夏季洪涝灾和席卷苏鲁冀等沿海地区的台风‘达维’造成的灾害。

另外,还有很多地区遭受了严重干旱、冰雹等自然灾害,共造成290亿美元的损失,但通过投保由保险公司赔付的比例仅占总损失的4%左右,这个比例相对美国的自然灾害保险赔付率相差甚远。

”另据报道:“2013年3月20日发生在广东、广西等省部分地区的一场大风和冰雹灾害,造成直接经济损失达13亿多元。

”这个事实警示我们,中国需要重视和加强自然灾害保险的研究和实践,特别是针对严重自然灾害的保险体系建设和对策方案的研究,推动由政府主导的自然灾害政策性保险方案的实施。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

自然灾害保险问题的研究摘要我国是农业大国,又是世界上遭受自然灾害损失最为严重的国家之一。

近10年来,自然灾害给我国造成的经济损失每年都在1000亿元以上。

自然灾害对农业经济发展的影响非常严重。

但与国际上大灾风险主要通过保险机制来分担化解的做法不同,我国自然灾害损失的救助工作主要依靠国家财政援助和生产自救进行,有关自然灾害风险防范的保险体系尚未真正建立。

因此,必需改革目前的保险体制,探索建立巨灾保险救助和通过资产证券化等非传统风险转移方式分散农业巨灾风险的新途径,有效地提升保险在国家灾害救助体系中的积极作用,因此我们分析了近几年天气,各地区的农作物种植面积,受灾,成灾,绝收面积的有关数据,得出了自然灾害的变化趋势,通过Excel,matlab等软件建立了几个模型以及分析出了受灾面积的函数y=-879.8x+2E+6,R*R=0.089,成灾面积y=-132.6X+21663,R*R+0.003绝收面积的函数y=-328.1X+66308,R*R=0.307并且还分析了出了降水量,风速,冰雹在近几年的变化趋势,为今后的预防工作和提出更加合理的保险险种方案做出了充分的准备。

关键词:自然灾害、保险险种、灾害变化趋势、土地种植面积、模型的建立一、问题重述根据2013年3月5日《环球时报》转摘美国《商业周报》的相关报道,“在2012年全世界发生的10大自然灾害中,有4场是发生在中国。

包括3场严重的夏季洪涝灾和席卷苏鲁冀等沿海地区的台风‘达维’造成的灾害。

另外,还有很多地区遭受了严重干旱、冰雹等自然灾害,共造成290亿美元的损失,但通过投保由保险公司赔付的比例仅占总损失的4%左右,这个比例相对美国的自然灾害保险赔付率相差甚远。

”另据报道:“2013年3月20日发生在广东、广西等省部分地区的一场大风和冰雹灾害,造成直接经济损失达13亿多元。

”这个事实警示我们,中国需要重视和加强自然灾害保险的研究和实践,特别是针对严重自然灾害的保险体系建设和对策方案的研究,推动由政府主导的自然灾害政策性保险方案的实施。

农业灾害保险是国家政策性保险之一,即政府为保障国家农业生产的发展,基于商业保险的原理并给予政策扶持的一类保险产品。

农业灾害保险也是针对自然灾害,保障农业生产的重要措施之一,是现代农业金融服务的重要组成部分,它与现代农业技术、现代农业信息化及市场建设共同构成整个农业现代化体系。

农业灾害保险险种是一种准公共产品,基于投保人、保险公司和政府三方面的利益,按照公平合理的定价原则设计,由保险公司经营的保险产品,三方各承担不同的责任、义务和风险。

农业灾害保险分种植业保险和养殖业保险两大类,现有几十个险种,因不同地区的气象条件和作物种类不同,其险种和设置方案都不尽相同。

农业灾害保险除遵循保险的共同原理外,有其自身的特点。

比如,其损失规律有别于人寿保险和通常的财产保险(如汽车险)等。

政府作为投保人和承保人之外的第三方介入以体现对国家安全和救灾的责任。

附件1给出了P省种植业现行的部分险种方案,请你们从实际出发,查阅和参考附件中的数据资料,通过分析建模,研究解决下面的问题:(1)对附件2中的数据做必要的统计分析,研究P省现有农业灾害保险险种方案可能存在的风险,并分析其方案是否存在不合理性。

(2)针对P省的具体情况,选取其中部分农业灾害保险险种,设计更实际可行的农业灾害保险的险种方案,包括标的、保险金、保费、费率、赔付率、政府补贴率等;并对方案的有效性(即保险公司和投保人的风险大小)及可行性做出定量分析。

(3)将你们的模型推广应用。

根据某省(市、区)的实际情况(或参见附件3),查阅相关资料,提出相应的农业灾害保险的险种方案,并对可能存在的风险做出分析;针对其它方面的自然灾害保险问题进行研究。

(4)结合你们的模型结果,从地方政府、保险公司和投保人三个方面,提出有利于自然灾害保险长远发展的对策方案,希望能用定量依据或方法说明其对策方案的可行性和有效性,并给政府相关部门写一篇建议书。

附件1:2012年P省政策性农业保险统颁条款(部分);附件2:P省10地区的2002~2011年的主要气象数据;附件3:全国各省(市、区)的1998~2011年的农作物受灾统计数据。

二、问题分析我国是世界上自然灾害发生频繁,据联合国统计资料显示,自20世纪以来中国是继美国、日本之后世界上自然灾害最严重的国家之一,人员伤亡和财产损失严重,给国民经济和人民生活造成巨大的损失。

特使2008年的南方冰雪灾害和汶川地震给我们留下了惨痛的记忆。

同时也暴露出我国巨灾保险制度的不足与缺失。

自然灾害导致的巨灾风险是保险公司需要面对的重大风险事件,巨灾事件带来的损失巨大一般保险公司无法承担,这样给保险公司的经营带来了严峻的挑战,有效进行风险监管,提高偿付能力已经成为保险业生存和科学发展的内在要求。

保险公司应该如何加强风险监管以维持自身经营的稳定性,同时保障受灾人民的生活?这是一个复杂的问题,需要保险公司和政府双方共同合作推进巨灾保险的制度完善。

针对问题一我们需要进行数据的统计与分析并且预测将来自然灾害的变化趋势,于是我们需要对附表中所给的数据进行了分类统计,得出未来几年各自然灾害的变化趋势,以及主要集中在哪些地区,各种自然灾害所占比例,它们主要集中在几月份,各地区的受灾状况,以及全国自然灾害的发展趋势,然后与P省现有农业灾害保险险种方案进行对比,提出合理化建议。

对于问题二我们选取P省其中苹果种植农业灾害保险险种进行改进,使它更具有合理性有效性。

针对问题三我们根据黑龙江省的实际情况提出了相应的农业灾害保险的险种方案,并对其中的风险进行了分析。

针对问题四从地方政府、保险公司和投保人三个方面,提出有利于自然灾害保险长远发展的对策方案三、模型的建立与求解一、问题一的分析(一)通过对附表二的数据分析我们得出暴雨的年降水量的变化趋势以及其函数,风速每年的变化趋势及函数并运matlab软件画出了他们的曲线图1、年降水量变化2、年风速变化由以上几个图我们可以知道暴雨的降水量在呈下降的趋势,风速也逐渐下降但他们在时间分布的并不均匀,没有规律,所以灾害的突发性挺强的,而且还可以看出降水量最高曾达到过250~300mm,最大风速15~20m/s所以说只能是总体上呈下降的趋势但某一年实际情况并不确定,这给农民带来了很大的风险,也就加重了政府和保险公司的责任,尤其对于现在P省现有农业灾害保险险种方案来说如果重大自然灾害给农民带来的损失是不可估量的,那么政府将承担大部分补偿工作,这样的设计具有非常大的不合理性,给国家带来很大的负担,例如P省的苹果西瓜等保险金额比较高的作物受到非常大的损害,那么政府补贴也将占到80%。

(二)选取西瓜种植灾害保险险种根据P省十个地区最近十年的天气变化数据我们画出了暴雨、冰雹、六级以上的风主要集中的月份,我们需要根据它们做好防范的准备,要求保险公司在灾害集中的月份把工作深入到农民家里去与农民一起做好预防灾害的准备工作,从而把损失减到最小,一方面减轻了农民的负担,另一方面也减轻了政府的负担。

根据上图我们可以得知暴雨在近十年来的分布的月份比较分散从三月份到十一月份均有分布,所以我们选取了保费比较高的西瓜种植保险方案进行改进。

如下:西瓜种植保险条款保险对象第一条凡本省行政区域内的农户、农民专业合作组织或集体经济组织种植西瓜的面积达5亩(含)以上,且生长管理正常的,均可向保险公司投保本保险。

在保险合同订立时,上述种植者或组织简称为被保险人。

投保本保险时,投保人必须将被保险人所有种植西瓜的地块全部投保,不得选择性投保。

保险金额、保险费率和保险费第四条每亩保险金额1000元,保险费率7%,每亩保险费70元。

品种保险金额(元/亩)费率保险费总保险费(元)政府补贴(元)农户交纳(元)西瓜1000 7% 70 42 28赔付率已开花正在形成果实的地块,在保险期限内遭受雹灾,保险公司在每亩有效保险金额的50%以内酌情赔付。

全部损失:果实成熟度在60%以下的保险西瓜,因雹灾全部绝产,按每亩有效保险金额的70%计算赔偿;成熟度在60%以上的保险西瓜,因雹灾全部绝产,按每亩有效保险金额的100%计算赔偿;部分损失:果实成熟度在60%以下的保险西瓜,因雹灾部分绝产,结合损失率与每亩有效保险金额的70%计算赔偿;成熟度在60%以上的保险西瓜,因雹灾部分绝产,结合损失率与每亩有效保险金额计算赔偿。

轻微损失:果实成熟度在60%以下的保险西瓜,因雹灾造成轻微损失,按西瓜每亩损失程度及数量计算赔付,但最高不超过每亩有效保险金额的70%;成熟度在60%以上的保险西瓜,因雹灾造成轻微损失,按西瓜每亩损失程度及数量计算赔付,但最高不超过每亩有效保险金额的80%。

已采摘部分西瓜的瓜园,在计算赔款时按照采摘部分占保险金额的比例相应扣除。

已采摘90%到完毕的瓜园,保险公司不再承担保险责任。

(三)浙江农业灾害保险从以上两个图中我们可以看出黑龙江省土地种植面积是相对比较大的,但同时通过我们对附表三的分析也得出它也是这十几年来受灾面积最大的省,如果不制定出合理的保险计划,这都将是政府人民的一大负担。

为此我们设计了更实际可行的计划如下:一、保险方案实行5倍赔付封顶方案,种养两业分别计算:当年农业保险累计赔付在在2至3倍部分的责任,由政府与共保体按1:1承担;在3至5倍部分的责任,由政府与共保体按2:1承担;赔付率达到500%时封顶,特殊情况需要赔付时由政府承担。

另外,浙江省政府建立了种植险巨灾风险准备金,每年按照当年种植业保险保费25%提取。

当种植业保险参保品种全年赔款总额在种植业保险保费1.3倍以下,由共保体全额赔付;当种植业保险参保品种全年赔款总额超过种植业保险保费1.3倍时,巨灾风险准备金负责承担1.3-2倍部分,超过2倍以上部分,按政府超赔方案进行分摊。

近十年来受灾的变化趋势图,二、存在的风险问题:1、保险公司或地方政府等单一主体无力单独承担农业巨灾风险。

农业保险危险单位大、承保周期长,存在灾害损失范围广、程度大的巨灾风险。

随着农业保险参保面积剧增,经办机构承担的风险不断积累,一旦发生大面积农业灾害,将造成经办机构收不抵支,极大影响经办机构的经营稳定性。

与此同时,经办机构还面临着国内外商业再保险承保能力不足,分保成本相对较高的难题。

另外,由于大多数省份缺乏农业保险风险分散机制,地方政府往往不得不承诺分担农业巨灾风险,会给地方财政带来较大的支出困境。

若真正发生农业巨灾,赔付封顶的政策可能会很难执行或引起诉讼。

2、部分已有地方性农业保险巨灾基金难以化解较大的巨灾风险现有的地方性农业保险巨灾基金大多规模有限且过于分散,往往仅由经办机构的经营盈余累积而建立,缺乏原始启动资金,且不利于实现灾害风险在全国范围内的分散,从而限制了巨灾风险的分散渠道和范围,难以完全化解突发的农业巨灾风险。

相关文档
最新文档