数学建模之出租车数目估计

数学建模之出租车数目估计

代码

m=100;%样本数目
n=20;%样本容量
x=zeros(5,m);
y=zeros(5,20);
for w=1:20
    n=w*20;
for k=1:m
        pe=randperm(1000);
        temp=pe(1:n);
        x(1,k)=2*mean(temp)-1;
        x(2,k)=2*median(temp)-1;
        x(3,k)=max(temp)+min(temp)-1;
        x(4,k)=(1+1/n)*max(temp)-1;
        x(5,k)=(1+1/(2*n-1))*(max(temp)-1/2*n);
end
  y(:,w)=mean(x,2);
end
for i=1:5
    plot(y(i,:));
    hold on;
end
ylabel('模型估计值');
xlabel('样本容量');
title('估计量收敛曲线');
legend('模型一','模型二','模型三','模型四','模型五');
 
% error=abs(y-1000)
% st=std(x,1,2)
% kx=500:20:1500;
% for i=1:5
% figure(i)
% hist(x(i,:),20,'red');
% name=['第',num2str(i),'个模型预估车辆数目']
% xlabel(name);
% ylabel('频数统计');
% title('统计直方图');
% end