2016-8-15 第一次训练 代码

求上路的实际出租车数量

clcclear allx0=[1:24];y0=[26340   21119   16365   12522   12520   20249   34261   47568   54608   51322   45526   44621   45330   43347   40866   41749   42999   41767   42325   43186   38832   36115   34679   33191];y0 = y0 ./ 65000; %上路专车实际比例a = polyfit(x0,y0,6) %6阶多项式拟合上路专车比例y1=polyval(a,x0); %拟合值得到上路出租车比例plot(x0, y0, x0, y1);y1 = y1.*66648 %求得24h上路出租车数

结果

a =
0.0000 -0.0000 0.0012 -0.0219 0.1950 -0.6948 0.9962

y1 =
1.0e+04 *
Columns 1 through 11
3.1700 1.5288 1.0955 1.3930 2.0671 2.8661 3.6235 4.2410 4.6735 4.9153 4.9886
Columns 12 through 22
4.9320 4.7921 4.6154 4.4424 4.3027 4.2120 4.1702 4.1615 4.1552 4.1090 3.9727
Columns 23 through 24
3.6946 3.2279

上路专车比例以及拟合图像

排队论计算平均等待时间

simu1.m

function [i,t]=simu1i=1;e=0;w=0;x(i)=random('exp',20);c(i)=x(i);b(i)=x(i);while(b(i)<=480)   y(i)=10*rand(1)+10;   e(i)=b(i)+y(i);   w=w+b(i)-c(i);   i=i+1;   x(i)=random('exp',20);   c(i)=c(i-1)+x(i);   b(i)=max(c(i),e(i-1));endi=i-1;t=w/i;%fprintf('一天完成服务人数：%f位 \n', i);%fprintf('每位顾客平均等待时间：%f分钟 \n', t);

simu2.m

function simu2M(1)=0;T(1)=0;for i=1:1000    [M(i),T(i)]=simu1;endmean_M=mean(M);mean_T=mean(T);fprintf('模拟1000次，平均每天完成服务人数：%f位 \n', mean_M);fprintf('模拟1000次，每位顾客平均等待时间：%f分钟 \n', mean_T);

结果

模拟1000次，平均每天完成服务人数：22.911000位
模拟1000次，每位顾客平均等待时间：14.996549分钟

AHP代码

clcclear allA=[1 1/4 1/5 3 ;4 1 1/3 7 ;5 3 1 6 ;1/3 1/7 1/6 1][n,n]=size(A);x=ones(n,100);y=ones(n,100);m=zeros(1,100);m(1)=max(x(:,1));y(:,1)=x(:,1);x(:,2)=A*y(:,1);m(2)=max(x(:,2));y(:,2)=x(:,2)/m(2);p=0.0001;i=2;k=abs(m(2)-m(1));while  k>p  i=i+1;  x(:,i)=A*y(:,i-1);  m(i)=max(x(:,i));  y(:,i)=x(:,i)/m(i);  k=abs(m(i)-m(i-1));enda=sum(y(:,i));w=y(:,i)/a;t=m(i);disp(w);disp(t);         %以下是一致性检验CI=(t-n)/(n-1);RI=[0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59];CR=CI/RI(n);if CR<0.10    disp('此矩阵的一致性可以接受!');    disp('CI=');disp(CI);    disp('CR=');disp(CR);end