基于制定合理医疗保障分配制度的数学模型.docx

1、基于制定合理医疗保障分配制度的数学模型基于制定合理医疗保障分配制度的数学模型摘要：随着社会的发展，各企业越来越注重员工的福利措施，医疗保障基金就是其中的一项。本文针对某一拥有四个子公司的集团关于医疗保障基金额度的分配问题，运用了科学的方法，建立了两个符合实际情况的预测模型，求解出了最佳的分配方案。在问题的分析中首先通过过去24年的医疗费用走势图预测2004年四个子公司的医疗保障费用将超过80万元。由此通过建立模型验证并求出四子公司的分配方案。模型一：通过Matlab软件基于多项式拟和原理，建立了1980-2003年关系的预测模型，求得各个子公司在1980-2003年所需保障费用的大体发展趋势与

2、走向，到2004年各个子公司所需保障费用分别为19.9764万元、17.1314 万元、24.5819万元、25.5532万元，具有较好的中短期预测效果。模型二：由于方差能够反映一组数据波动的大小，通过Matlab软件计算出各子公司在1980年消费指数条件下的医疗费用方差。运用式子预测出各子公司2004年在1980年消费指数条件下的医疗费用。再通过2004年的通货膨胀指数计算出到2004年各个子公司所需保障费用分别为20.7189万元、17.7369万元、24.1846万元、25.8031万元。这和模型一的结果是非常接近的，由此得出模型一是合理的。模型一和模型二的结果都超过80万元，于是80万

3、元显然不能满足四个子公司的需求。为了使四子公司得到公平的分配，这里运用权重比例分配法，合理分配了80万元，分配方案为18.31万元、15.71万元、22.54万元、23.43万元。本模型在计算与统计过程中，引入了多项式拟和和方差估计两种方法，在实际工程中，可以对商品的销统计、人口的统计及其他的预测控制带来帮助。本文建立的分配方案模型不仅可以用于基金分配，还可以用于材料分配、人力资源分配等问题。关键词：多项式拟和 方差估计 合理分配一.问题的重述： 传统的合作医疗保障制度随着社会的发展，日益不能满足各项需求，因此根据已有数据，运用数学建模的方法，制定出合理的、现代医疗保障分配制度是一个重要问题。

4、某集团下分别有A、B、C、D四个子公司，各子公司财务分别独立核算，每个子公司都实施了对雇员的医疗保障计划，由各子公司自行承担雇员的全部医疗费用。为进一步规范各个子公司的医疗保障计划，集团董事会规定，在2003年底,各个子公司均需以银行活期存款的方式，设立医疗保障基金，基金专门用于支付2004年度雇员的医疗费用。 并规定每个子公司的医疗保障基金只能用于支付本子公司雇员。董事会综合考虑了各种因素，确定本集团设立的2004年度医疗保障基金的总额度为80万元，这一额度在四个子公司之间分配。对于各子公司，如果2004年度总的医疗费用支出低于该子公司的医疗保障基金的额度， 则雇员可以及时得到医疗方面的保障

5、。 而如果总的医疗费用超过了医疗保障基金的额度， 则子公司需要通过其他渠道来筹措超出部分的额度。这会导致某些雇员无法及时报销医疗费用。并且每个子公司的雇员人数以及每一年龄段的雇员比例，在各年度都保持相对稳定，四个子公司的个年度的医疗费用支出见表1；相关年度的通货膨胀指数见表2；已知2004年银行活期存款利率为1%由此确定80万元医疗保障基金在四个子公司之间的分配方案表1：年度 公司A 公司B 公司C 公司D1980 8.28 8.81 8.02 10.501981 8.76 9.31 8.36 10.761982 9.29 10.41 9.20 11.341983 10.73 11.61 10

6、.51 12.891984 10.88 11.39 10.70 13.201985 11.34 12.53 11.20 14.241986 11.97 13.58 13.01 14.901987 12.02 13.70 13.24 15.281988 12.16 13.32 13.82 15.201989 12.83 14.32 14.74 16.961990 13.90 15.84 17.33 19.231991 14.71 14.67 18.50 20.991992 16.11 14.99 17.72 23.221993 16.40 14.56 18.45 23.191994 17.07

7、14.55 19.77 24.041995 16.96 14.80 19.56 23.721996 16.88 15.41 19.70 24.881997 17.20 15.76 20.52 27.341998 19.87 16.76 22.51 28.121999 20.19 17.68 23.10 28.382000 20.00 17.33 23.24 28.812001 19.81 17.03 23.22 26.712002 19.40 16.95 23.76 19.822003 20.48 16.66 24.50 20.48在下图中，轴从1到24分别按顺序表示从1980年到2003年的

8、各年份， 轴表示相应的支出金额（单位：万元）。图形如下：表2：日期 通货膨胀指数 日期 通货膨胀指数1980年1月1日 100 1992年1月1日 1811981年1月1日 105 1993年1月1日 1821982年1月1日 115 1994年1月1日 1891983年1月1日 128 1995年1月1日 1911984年1月1日 130 1996年1月1日 1971985年1月1日 140 1997年1月1日 20001986年1月1日 150 1998年1月1日 2101987年1月1日 151 1999年1月1日 2111988年1月1日 152 2000年1月1日 2151989年1

9、月1日 160 2001年1月1日 2171990年1月1日 175 2002年1月1日 2191991年1月1日 180 2003年1月1日 226二问题的分析：从医疗保障基金的性质来看，其目的是为员工的医疗费用进行补偿，由于集团是由四个子公司组成，为了在四个公司之间的金额分配更合理，更公平，从而达到真正为公司员工服务的目的，集团就应当统筹安排基金总额及分配到各子公司的具体金额。如果基金总额不足以让每个子公司都能及时给雇员医疗保障，那么集团不能让某些子公司满足而另外的相差甚远，只能通过合理的分配让每个子公司自行筹措一部分，这样才能体现出公平性。如果基金总额能够使每个子公司雇员得到及时报销，则

10、只需满足每子公司有足够的报销费用，这样才能体现出合理性。.根据医疗保障发放的实际情况，医疗费用报销可随到随报、按月报、按季度报以及年终报销等制度。为了使利息最高，采取年终报销最好。由于2003年度各子公司医疗费用支出总额超出了计划设立的2004年保障基金本息和（20.48+16.66+|24.50+20.48=82.1280）显然各子公司必须自行筹措部分支出. 每个子公司的雇员人数以及每一年龄段的雇员比例，在各年度保持相对稳定，这说明各子公司每年在该方面的支出也是相对稳定的。我们可以把数据转化为同一年度进行验证.三模型的假设：1 各子公司各年度医疗费用的变化主要是因为通货膨胀引起的。2 雇员是

11、在年终报销医疗费用。3 2004年社会安定，全国GDP增长保持稳定，全国无重大灾害，疫情发生。4 银行利率变化的影响不考虑。四符号说明：Z 基金总额度(80万)Z 基金存入银行后一年的本息和（80.8万） A公司第i年所分配到的基金额度 B公司第i年所分配到的基金额度 C公司第i年所分配到的基金额度 D公司第i年所分配到的基金额度五模型的建立与求解：由于每个公司的雇员人数以及每一年龄段的雇员比例在各年度保持相对稳定，我们把各公司各年的医疗费用支出转换为1980年的医疗费用支，由此预测出2004年时医疗费用支，最后再乘以2004年的通货膨胀指数，得到2004年实际各公司的医疗费用支费用。方法一：

12、拟合法1，采用拟合方法对各年医疗支出费用（转化成1980年支出费用）进行预测时：拟合方法对各年医疗支出费用源程序见（A）经过实验,用7次方为最高次数拟合效果很好: 此时,通货膨胀指数函数为227.0469.取为2272，采用拟合方法对通货膨胀指数进行预测：拟合方法对通货膨胀指数进行预测源程序见（B）最后求得各公司的医疗费用如下：A公司预测所得到的2004年的医疗费用为19.9764万元B公司预测所得到的2004年的医疗费用为17.1314 万元C公司预测所得到的2004年的医疗费用为24.5819万元D公司预测所得到的2004年的医疗费用为25.5532万元总费用为：87.2429万元方法二：

13、方差法将各子公司每年的医疗费用开支转化为1980年的费用. 如下表：(单位:万元)：年份 A公司 B公司 C公司 D公司1980 8.2800 8.8100 8.0200 10.50001981 8.3429 8.8667 7.9619 10.24761982 8.0783 9.0522 8.0000 9.86091983 8.3828 9.0703 8.2109 10.07031984 8.3692 8.7615 8.2308 10.15381985 8.1000 8.9500 8.0000 10.17141986 7.9800 9.0533 8.6733 9.93331987 7.9603

14、 9.0728 8.7682 10.11921988 8.0000 8.7632 9.0921 10.00001989 8.0187 8.9500 9.2125 10.60001990 7.9429 9.0514 9.9029 10.98861991 8.1722 8.1500 10.2778 11.66111992 8.9006 8.2818 9.7901 12.82871993 9.0110 8.0000 10.1374 12.74181994 9.0317 7.6984 10.4603 12.71961995 8.8796 7.7487 10.2408 12.41881996 8.568

15、5 7.8223 10.0000 12.62941997 8.6000 7.8800 10.2600 13.67001998 9.4619 7.9810 10.7190 13.39051999 9.5687 8.3791 10.9479 13.45022000 9.3023 8.0605 10.8093 13.40002001 9.1290 7.8479 10.7005 12.30882002 8.8584 7.7397 10.8493 9.05022003 9.0619 7.3717 10.8407 9.0619由于方差能够反映一组数据波动的大小，其表达式为，对于一组离散且无明显分布概率的数

16、据，还可表示为。通过如下Matlab的程序可以得到A,B,C,D四个子公司在1980年消费指数条件下的医疗费用方差份分别为：，。在下图是1980年到2003年各年份各子公司按1980年消费指数的医疗支出费用的走势图。（单位：万元）从此图中我们不难看出A,B,C三个子公司的医疗费用的波动是比较小的，同上面所作的方差一致。且从图中可以明显看出2003年A,C两子公司的医疗费用在2002年的基础上是增加的，而B公司则刚好相反。于是我们可以预测2004年A,B,C三个子公司在1980年物价情况下的医疗费用符合关系式：。由此我们求出：（万元），（万元），（万元）。接下来我们再对D公司进行预测，由于无论是

17、从方差还是从走势图都显示出D公司的医疗费用波动是非常大的，所以D公司2004年的医疗费用不能用上面的方法预测。又从走势图中可以看出D公司的医疗费用在1990到2000年间很高，而在2002到2003年又降下来并趋近于平稳。由此我们可以这样认为：该公司在1990到2000年间，雇员的健康程度逐年下滑（暂且认为年龄增长），而在2000到2003年由于老雇员的退出和年轻雇员的引进，年医疗费用就降下来并趋于稳定。因此可以用1990到2003年的平均值作为2004年的医疗费用：。由方案一预测出的2004年通货膨胀指数2.26乘以上面算出的结果就得到了2004年各公司应分配的医疗费用分别是:20.7189

18、万元，17.7369万元，24.1846万元和25.8031万元。这个结果和方案一的结果是非常接近的，于是我们可以认为方案一是合理的。六.模型的评价与推广：1.通过多种方法对此模型求解结果的比较可以看出采用拟合方法对各年医疗支出费用的预测是合理的，符合公司的实际情况，可以在公司推广和应用。2.在预测2004年各个子公司A、B、C、D的医疗费用分配问题时，对近24年来这些大量波动且随机的数据，我们采用了拟合和方差两种方案，并对结果的比较，选取最优方案，保证了模型的正确性和准确性。3.在对题中所给的大量数据进行处理时，为了避免大量的计算，我们利用了Matlab语言编程。这种做法值得推广到更大的数据

19、处理问题中。七.参考文献：1.XXXXX2.XXXXX八.文中涉及的各源程序如下：（A）拟合方法对各年医疗支出费用源程序：A=load(e:7.txt);f=A;b=f(4,:);t=1980:2003;ni=polyfit(t,b,5);nit1=1980:2004;y1=polyval(ni,t1);plot(t,b,k*,t1,y1,k);title(A公司医疗支出费用预测图表);xlabel(年份);ylabel(A公司医疗支出费用);y1（B）拟合方法对通货膨胀指数进行预测源程序：A=load(e:6.txt);b=A;t=1980:2003;ni=polyfit(t,b,7);ni

20、t1=1980:2004;y1=polyval(ni,t1);plot(t,b,k*,t1,y1,k);title(通货膨胀指数预测图表);xlabel(年份);ylabel(通货膨胀指数);y1（C）方差程序：n1=;n2=;n3=;n4=;a=;k1=0;k2=0;k3=0;k4=0;b1=0;b2=0;b3=0;b4=0;x1=0;x2=0;x3=0;x4=0;c=8.2800 8.8100 8.0200 10.5000 8.3429 8.8667 7.9619 10.2476 8.0783 9.0522 8.0000 9.8609 9.0932 9.8390 8.9068 10.923

21、7 8.3692 8.7615 8.2308 10.1538 8.1000 8.9500 8.0000 10.1714 7.9800 9.0533 8.6733 9.9333 7.9603 9.0728 8.7682 10.1192 8.0000 8.7632 8.7368 10.0000 8.0188 8.9500 9.2125 10.6000 7.9429 9.0514 9.9029 10.9886 8.1722 8.1500 10.2778 11.6611 8.9006 8.2818 9.7901 12.8287 9.0110 8.0000 10.1374 12.7418 9.0317

22、7.6984 10.4603 12.7196 8.8796 7.7487 10.2408 12.4188 8.5685 7.8223 10.0000 12.6294 8.6000 7.8800 10.2600 13.6700 9.4619 7.9810 10.7190 13.3905 9.5687 8.3791 10.9479 13.4502 9.3023 8.0605 10.8093 13.4000 9.1290 7.8479 10.7005 12.3088 8.8584 7.7397 10.8493 9.0502 9.0619 7.3717 10.8407 9.0619;for i=1:2

23、4 n1(i,1)=c(i,1)-8.612975; n2(i,2)=c(i,2)-8.422133; n3(i,3)=c(i,3)-9.601929; n4(i,4)=c(i,4)-11.36790;endfor k1=n1(:,1).2; k2=n2(:,2).2; k3=n3(:,3).2; k4=n4(:,4).2;endfor i=1:24 b1=b1+k1(i,1); b2=b2+k2(i,1); b3=b3+k3(i,1); b4=b4+k4(i,1);endfor x1=b1/24;x2=b2/24;x3=b3/24;x4=b4/24;endx1x2x3x4（D）网络神经源程序

24、：input=1:24;h=load(e:1.txt);target=h;net=newff(1 24,5,1,tansig,purelin,trainlm);net.trainparam.epochs=100;net.trainparam.goal=0.001;net=train(net,input,target);figure(1)plot(input,target,*);title()xlabel(input);ylabel(target);pause;hold on;xtest=(1:0.01:25);ytest=sim(net,xtest);plot(xtest,ytest,r);title();w1=net.iw(1)w2=net.iw(2)b1=net.b(1)b2=net.b(2)ytest