上机2：描述统计和指数.doc

描述统计和指数 分组频数与命令的编制 x=randn(50,1);　%生成50行1列的标准正态分布的随机数，存入x中； edges=[-inf -1 -0.5 0 0.5 1 inf]　%对x中数据分组时的各组组限值； freq=histc(x,edges) % 以edges为组限对x中的数据进行分组，返回各组频数存入freq中。 打开一个空白m文件，将1.2节中的3条语句复制到文件中，将文件保存为mycomd.m。在命令窗口中运行mycomd，即相当于运行这3条语句。 分组频率与函数的编制 打开一个空白m文件，将下列语句复制到文件中，将文件保存为freq_rate.m: function f=freq_rate(x,edges) freq=histc(x,edges); total_freq=sum(freq); f=freq/total_freq; 然后在命令窗口中运行: x=randn(50,1); edges=[-inf -1 -0.5 0 0.5 1 inf]; freq_rate(x,edges)。 直方图 x=randn(50,1);hist(x); %默认分为10组 hist(x,5) %分为5组 [n,xout] = hist(x,5) %返回各组的频数和组中值. hist(x,[-1 -0.5 0 0.5 1]) %指明各组组中值 折线图 freq_rate=[8;15;23;34;15;5] cumsum(freq_rate) %向下累计 flipud(cumsum(flipud(freq_rate)))% 向上累计 plot(x,y) 以x为横坐标，y为纵坐标，绘制曲线。 plot(75:50:325, freq_rate,.-) .-中的’.’表示每个数据点用一个点来表示，-表示各点用实线连接 plot(-3:0.1:3, normpdf(-3:0.1:3,0,1),g*--) 画标准正态分布的概率密度曲线。g表示绿色, ’*’表示每个数据点用一个星号来表示，--表示各点用虚线连接 常用的点类型、线类型和颜色如下表所示： 练习：运行如下绘制图形的例子。 clear; x=-3:0.1:3; % 生成向量x。或用: x=linspace(-3,3,100) y1=normpdf(x,0,1); %计算标准正态分布的概率密度函数值，放入y1。 y2=normcdf(x,0,1); %计算标准正态分布的累积分布函数值，放入y2。 figure; % 创建一个新的图形窗口。 subplot(1,2,1); % 定义第一个子图区域。 plot(x,y1); % 用实线画曲线。 subplot(1,2,2); % 定义第二个子图区域。 plot(x,y2, .); % 用‘.’号画曲线。 上述命令画出了标准正态分布的概率密度函数和累积分布函数曲线。 饼图 pie(freq_rate) pie3(freq_rate) explode=[0 1 0 1 0 0] pie3(freq_rate,explode) 柱状图 bar(freq_rate) bar3(freq_rate) 平均数、众数、中位数、标准差 mode(x) %众数 mean(x) 　%平均数 x=[75 125 175 225 275] f=[5 6 20 4 5] x_bar=dot(x,f)/sum(f)　　 %加权平均： 练习：编制一个求加权平均数的函数。 median(x)　%中位数 prod(x)^(1/length(x)) %x中数据的几何平均 geomean(x)　% x中数据的几何平均 prod(x.^f)^(1/sum(f)) %加权几何平均 range(x)　　%全距 mad(x) %平均差 std(x) %样本标准差 std(x,1) %总体标准差 var(x) %样本方差 var(x,1) %总体方差 物量指数和价格指数计算 例4-2-1: Pn=[10 100 3]; P0=[10 110 4]; P1=[12 120 5]; Q0=[100 50 1000]; Q1=[130 60 800]; KQ_bar_1=dot(Q1,P0)/dot(Q0,P0) KQ_bar_2=dot(Q1,P1)/dot(Q0,P1) KQ_bar_3=dot(Q1,Pn)/dot(Q0,Pn) KP_bar_1=dot(P1,Q0)/dot(P0,Q0) KP_bar_2=dot(P1,Q1)/dot(P0,Q1) 练习：计算课本第77页例4-3-1。 作业（将相应Matlab语句抄在纸上上交） (1)编制一个求加权标准差的函数