估计出租车的总数。假如你在马路上见到连续经过的10辆出租车的车牌号如下:0357,0476,0156,0410,0243,0674,0823,0520,0580,0789.试用不同的方法来估算你所在的城市的出租车总数。要求:(1)尽可能多的用不同的方法估算;(2)通过适当的假设和数值模拟计算对不同方法的准确性和稳定性进行比较。
A=[0357,0476,0156,0410,0243,0674,0823,0520,0580,0789];
B=sort(A); %把数按从小到大的顺序排列
a=sum(A)/10; %求A的平均数
x1=2*a-1 %平均值模型
x2=2*median(B)-1 %中位数模型
x3=B(1,10)+B(1,1)-1 %两短间隔对称模型
x4=(1+1/10)*B(1,10)-1 %平均间隔模型
x5=(1+1/(2*10-1))*(B(1,10)-1/(2*10)) %区间均分模型
A=randi(950,100,6); %950个数里取100个样本,每个样本6个数据。
m(1:100,1)=mean(A,2); %计算m矩阵的第一列存放这100个样本的每一个样本的平均值
m(1:100,2)=median(A,2); %计算m矩阵的第二列存放这100个样本的每一个样本的中位数
m(1:100,3)=max(A,[],2); %计算m矩阵的第三列存放这100个样本的每一个样本的最大值
m(1:100,4)=min(A,[],2); %计算m矩阵的第四列存放这100个样本的每一个样本的最小值
%计算预估值x0。x(:,1)分别存放模型i的预估值,利用五个模型的公式计算可得。
x(:,1)=2.*m(:,1); %平均值模型
x(:,2)=2.*m(:,2); %中位数模型
x(:,3)=m(:,3)+m(:,4)-1; %两端间隔对称模型
x(:,4)=(1+1/6).*m(:,3)-1; %平均间隔模型
x(:,5)=(1+1/12).*m(:,3); %区间均分模型
%计算每个模型的平均值,平均值误差,标准差
p(1,:)=mean(x,1);%存放每个模型的平均值
p(2,:)=p(1,:)-950; %每个模型的平均值误差
p(3,:)=std(x,0,1) %每个模型的标准差