系统仿真（一）：到达教学楼时间拟合

到达教学楼时间拟合

一、数据收集

到达教学楼时间统计范围：7：50-8：40。
统计数据方法：每分钟一组。
数据记录方式：表格，共51组原始数据。原始数据如图1所示，直方图如图2所示。

图1 到达教学楼时间记录表图

图2 到达教学楼时间统计直方图

数据优化：通过观察数据，由于样本量较小，时间间隔过小，出现个别偶然性数据噪声。现优化51组数据原1分钟每组合并为3分钟每组，为方便数据统计，设8：00为原点0，原7：50-8：40范围由（-10，40）调整为（-4，11），优化后数据如表1所示，直方图如图3所示，共16组。

表1 优化后数量表

图3 优化后数据直方图

二、数据拟合

初步观察图3所示得统计数据绘制分布直方图，判断到达教学楼不同时间的人数符合高斯分布，故采用高斯分布模型拟合数据。
采用Matlab中拟合函数fit( xx, yy,‘gauss1’ )。
并求得参数样本均值为3.5，样本方差为4.7610。
画图，拟合曲线和原始数据图像如图3所示。

图3 原始数据和拟合曲线图

三、拟合检验

使用cdf方法检验,调用matlab中normcdf()函数。
kstest(函数)检验是否符合置信区间为95%的正态分布。
首先，由于正态分布的均值、标准差未知，因此先计算样本的均值和标准差，再做正态拟合。通过样本资料的计算建立假设：H0：样本数据服从均值为5.483，标准差为3.207的正态分布，H1：样本数据不服从均值为5.483，标准差为3.207的正态分布。取alpha=0.05，完成cdf的检验，最后检验是否符合置信区间为95%的正态分布，符合则输出符合，反之不符合。
检验结果如图4所示，即服从。

图4 KS检验数据是否符合高斯分布

四、总结

根据统计和结果分析，由原始数据分析得同学们到达教学楼时间均值为8:16，众值为8:21，由优化后数据分析到达时间服从。

附录

clear;clc;close all
format compact
%正态分布的拟合
num = 16;
filename='time.xlsx'; %读取文件地址
% x=xlsread(filename,"Sheet1",'c2:c17'); %参数依次是：文件，sheet名称，读取的范围
% y=xlsread(filename,"Sheet1","d2:d17");
x=xlsread(filename,"Sheet1",'a2:a52'); %参数依次是：文件，sheet名称，读取的范围
y=xlsread(filename,"Sheet1","b2:b52");
xx = x(:);
yy = y(:);
bar(xx,yy)
example_mean = mean(xx);%样本均值
example_std = std(xx);%标准差
% Set up fittype and options.
[fitresult, gof] = fit( xx, yy,'gauss1' );%直接调用 高斯函数模型
fitresult;
% Plot fit with data.
figure;
plot(fitresult)
hold on
plot(xx, yy,'b*');
legend('拟合曲线','原始数据', 'Location', 'NorthEast');
title(['正态分布拟合,num=',num2str(num)])
xlabel('x');
ylabel('y');
grid on
saveas(gcf,'pic.png')
% 输出拟合参数
a = fitresult.a1;
wa = fitresult.b1;
xc = fitresult.c1;alpha = 0.05;
[mu,sigma]=normfit(xx);%求正态分布的参数值
p1=normcdf(xx,mu,sigma);%cdf检验
[H1,S1]=kstest(xx,[xx,p1],alpha);%检验是否符合置信区间为95%的正态分布
if H1 == 0disp("服从正态分布");
elsedisp("不服从");
end

系统仿真（一）：到达教学楼时间拟合相关推荐

wpf datagrid 计算两个列的差值_天体到达指定方位角的时间计算
接上一篇天体的升起.中天.落下时间 - 且听风吟的文章 - 知乎,天体的升起落下和中天是满足一些特殊条件,对于一般情况,如指定一个方位角,要求其到达的时间,我们也可以计算出来. 可以先判断一个天体的方 ...
如何在cad中导入谷歌地图_如何在Google地图中设置出发和到达时间
如何在cad中导入谷歌地图 Unless you work from home, Google Maps is fantastic at knowing where you need to be, a ...
三中新教学楼信息系统集成方案
第 1 章用户需求分析信息化是我国现代化建设的组成部分,可以大大提高我国综合国力,它关系到各级政府部门,各经济.科技.文化机构,对整个社会的方方面面都会产生重要影响.目前,由国家各部委单位主持的国 ...
二、四维拟合股市理论的时空哲学
笔者将其归为股市时空哲学,意思很明显,这个拟合理论有时候可用,有时候不可用.书中举例都是可用的例子,但是对于中枢行情,要谨慎这个V字时空转化理论无法拟合. 这幅图是空间转化为时间发生共振的两个例子. ...
日期类的时间从为什么是从1970年1月1日（格林威治时间）
I suspect that Java was born and raised on a UNIX system. UNIX considers the epoch (when did time be ...
任务的插入时间复杂度优化到 O(1)，Timing Wheel时间轮是怎么做到的?
点击上方蓝色"方志朋",选择"设为星标" 回复"666"获取独家整理的学习资料! 在kafka中,有许多请求并不是立即返回,而且处理完一些异 ...
winform 让他间隔一段时间执行事件且只执行一次_Redis 事件机制详解
点击上方"程序员历小冰",选择"置顶或者星标" 你的关注意义重大! Redis 采用事件驱动机制来处理大量的网络IO.它并没有使用 libevent 或者 li ...
Setup 和Hold （建立时间和保持时间）解析
本文是第一篇,有空会继续更新.(转载请注明出处!!!) STA分析是基于同步电路设计模型的,在数据输入端,假设外部也是同时钟的寄存器的输出并且经过若干组合逻辑进入本级,而输出也被认为是驱动后一级的同时 ...
事务、事件（文件、时间、调度和执行）、复制、分片（范围、哈希）、简单的论坛系统分析
1. 事务一个事务包含了多个命令,服务器在执行事务期间,不会改去执行其它客户端的命令请求. 事务中的多个命令被一次性发送给服务器,而不是一条一条发送,这种方式被称为流水线. 它可以减少客户端与服务器 ...

系统仿真（一）：到达教学楼时间拟合

一、数据收集

二、数据拟合

三、拟合检验

四、总结

附录

系统仿真（一）：到达教学楼时间拟合相关推荐

最新文章

热门文章

系统仿真（一）：到达教学楼时间拟合

一、数据收集

二、数据拟合

三、 拟合检验

四、总结

附录

系统仿真（一）：到达教学楼时间拟合相关推荐

最新文章

热门文章

三、拟合检验