数模中的表和图.doc

例1证卷交易所的证卷交易价格的综合指标称为该交易所的“交易指数”,记为。若当日交易指数为,上一日交易指数为,则称为该交易指数的涨跌,这是一个连续型随机变量。 上海证卷交易所2005年前50个交易日的交易指数的涨跌观测值如下： 13.93 -6.92 -6.13 -14.79 -15.70 -2.83 -11.01 -4.28 -9.03 -0.87 5.70 -21.92 -0.48 -17.80 -5.87 8.20 -2.67 -28.87 -1.23 1.26 19.61 -11.98 7.46 -0.73 -5.27 -4.47 -4.61 1.20 6.18 53.50 -5.51 -30.70 2.84 -12.70 7.70 3.89 16.37 39.08 16.66 -12.15 -15.22 -19.30 -0.06 2.01 -15.64 7.28 13.64 -8.07 6.50 21.75 根据交易指数的涨跌的观测数据，编制频数分布表。基本步骤如下： 首先，将上述数据载入Matlab系统的工作内存 clear data=[13.93 -6.92 -6.13 -14.79 -15.70 -2.83 -11.01 -4.28 -9.03 -0.87,5.70 -21.92 -0.48 -17.80 -5.87 8.20 -2.67 -28.87 -1.23 1.26 19.61 -11.98 7.46 -0.73 -5.27 -4.47 - 4.61 1.20 6.18 53.50 -5.51 -30.70 2.84 -12.70 7.70 3.89 16.37 39.08 16.66 -12.15 -15.22 -19.30 -0.06 2.01 -15.64 7.28 13.64 -8.07 6.50 21.75];; 【说明】指令clear的功能是清空Matlab系统的工作内存，每一项新的工作之前运行这一指令，有助于防止运算错误的发生。 编制频数分布表。基本步骤如下： ⑴ 数据分组 ①确定数据组个数　根据样本容量，确定分组数，推荐公式。k=ceil(1+3.322*log10(numel(data))) k = 7 【说明】这里用了三个Matlab系统函数：log10()求常用对数，numel()求矩阵的元素个数，ceil()向正无穷大圆整数据。另外，这里还涉及到了Matlab的数值计算问题，稍后予以讲解。 分组数也可参考下表确定： 样本容量 分组数 ＜50 5～6 50～100 6～10 100～250 7～12 ＞250 12～20 ② 计算极差 计算公式为。 ③ 确定组距 计算公式为，一般总取为数据的最小测量单位的整数倍。 ④ 确定各组端点 计算公式为,。的确定方法，一般地取比数据的最小值小半个测量单位。 ⑵ 统计各组频数 各组频数就是数据落入各个小组中的个数，记为。 上述计算极差、确定组距、确定各组端点和统计各组频数的工作均可由Matlab系统函数hist完成。hist的输入参数通常有两个，第一个是数据向量，第二个是小组个数，hist的输出参数有两个，第一个返回各组的数据频数，第二个返回各个数据组的区间位置值（组中值）。 [ni,ak]=hist(data,k) ni = 4 13 17 11 2 1 1 ak = Columns 1 through 4 -24.6857 -12.6571 -0.6286 11.4000 Columns 5 through 7 23.4286 35.4571 47.4857 ⑶ 计算频率 ① 计算各组频率　计算公式为。 fi=ni/length(data) fi = Columns 1 through 4 0.0816 0.2653 0.3469 0.2245 Columns 5 through 7 0.0408 0.0204 0.0204 ② 计算各组累积频率　计算公式为。 mfi=cumsum(fi) mfi = Columns 1 through 4 0.0816 0.3469 0.6939 0.9184 Columns 5 through 7 0.9592 0.9796 1.0000 【说明】cumsum()是Matlab系统函数，功能是对向量的元素累积求和 ⑷ 编制频数分布表 将上述计算结果列表如下，即得 表1-4　1995年前50个交易日上证指数涨跌分组频数分布