使用MATLAB对单程票卡分布的研究.docVIP

PAGE \* MERGEFORMAT3 使用MATLAB对单程票卡生态环境的研究 前提 单程票系统稳定；（根据5年运营数据可得到） 系统内票卡总量基本不变；（常识） 每天流失量等于更新量； 每天更新量基本为常数。 假设 设备系统总量为1，即100%； 票卡每日使用率，为p； 使用时票卡的流失率为u； 则，每日补充系统的u*p的新卡，系统稳定。 设计 *假设卡片能用N次，即卡片的最高使用次数不超过N； 则可用N+1维向量V对票卡分布进行描述： R=（r1，r2，r3，…，rN）； 其中ri表示使用次数为i的卡在系统中所占的比例，显然Σri=100%； 依上述假设(二)可知： 设第a天为Ra，第a+1天为R a+1 则： R a+1（1）=（1-p）*Ra（1）+ u*p； 当b1: R a+1（b）=（1-p）*Ra（b）+(1-u)*p*Ra（b-1）; 令：V=(1,R); 构造变换矩阵T ,实现上述递推关系。 T= 1 u*p 0 0 …. 10 1-p 1-u*p 0 …. -1 0 0 1-p 1-u*p … -10 0 0 1-p … -1……… ………… ….0 0 0 0 … 0 其中，u为流失率；p为使用率。 则： 递推公式为 Vi+1=Vi*T  调试 程序设N=300；u=0.02；p=0.25 程序模拟天数为1500天，数据已保持稳定。  源码 N=300; Day=2000; min=0.000001; u=0.02;% 使用丢失率 p=0.25;%票卡使用率，重要!!!!!!!!!!!!!!不是每天每张票卡均被使用一次 % u*p为票卡率，亦为系统补充率!!!!!!!!!!!!! V=zeros(1,N); V(1)=1; F=V; T=zeros(N,N); for i=1:N-1 T(i,i+1)=(1-u)*p; T(i+1,i+1)=1-p; end %初始化 in=1; for i=2:Day V=V*T; V(1)=u*p+(1-p)*in; in=V(1); for j=2:N %·防止计算零误差 if V(j)min V(j)=0; end end if iN V(i)=1-sum(V(1:i-1)); if V(i-1)=V(i) %防止计算零误差导致数据震荡；此时数据应为下降趋势，由此判断 V(i)=V(i-1)*(1-u*p); end else V(N)=0; end F=[F;V]; end save ticket F;%保存计算文件。