概率统计计算matlab.doc

《统计计算与软件》期末考查报告 年级：2013 专业：数学与应用数学 成绩： 1、生成4行5列的N(1,4)随机数矩阵，（1）计算出各行各列的均值、方差、标准差、4阶矩，以及任意两行两列的协方差和相关系数；（2）采用for循环语句编写程序计算出该矩阵的第二行样本的偏斜度以及第三列样本的峰度，并用matlab命令进行验证程序计算结果的正确性。 解： A=normrnd(1,2,4,5) A = 0.6702 -0.7273 0.9863 0.5488 2.1051 2.2554 1.1547 4.0653 3.2347 3.2012 3.1865 -1.4282 -0.5393 -1.1781 4.0884 3.2185 -1.2270 1.7428 1.0651 1.1719 B=mean(A) B = 2.3327 -0.5570 1.5637 0.9176 2.6416 B=mean(A,2) B = 0.7166 2.7823 0.8259 1.1943 B=mean(A,3) 这个没必要算可以看到这个B就是原始矩阵A B = 0.6702 -0.7273 0.9863 0.5488 2.1051 2.2554 1.1547 4.0653 3.2347 3.2012 3.1865 -1.4282 -0.5393 -1.1781 4.0884 3.2185 -1.2270 1.7428 1.0651 1.1719 C=var(A) 这个命令算出来的方差是无偏的，除以的是（n-1,命令var（算出来的是有偏的列向量方差，除以的是n。 C = 1.4278 1.3890 3.6820 3.3063 1.6182 C=var(A,1,2) 算的行向量的有偏方差 C = 0.8231 0.9906 5.4356 2.0557 D=var(A,1,2) D = 0.8231 0.9906 5.4356 2.0557 F=std(A) 算的列向量的无偏标准差 F = 1.1949 1.1785 1.9189 1.8183 1.2721 G=std(A,1,2)算的行向量有偏标准差，1表示有偏，2表示行向量 G = 0.9073 0.9953 2.3314 1.4338 H=moment(A,4) H = 2.1964 2.3406 14.7080 12.0338 2.3073 I=moment(A,4,2) I = 1.6138 1.9752 37.9551 10.2502 J=cov(A(1,:),A(3,:)) 为什么只算第一行和第三行的协方差 J = 1.0289 1.9108 1.9108 6.7944 1.0289为向量A(1,:)和它自身的协方差也就是第一行向量的方差；同样6.7944表示A(3,:)的方差。 A(1,:)和的协方差，矩阵A(1,:)的协方差，由于cov(a,b)=cov(b,a),所以这个数也为1.9108. 利用cov命令可算的各列向量之间的协方差，cov（A 要算各行向量的协方差，只需用A’命令，将行向量变为了列向量，再利用命令cov（’）算得各行之间的协方差 cov(A) ans = 1.4278 -0.3955 -0.2901 -0.4082 0.2607 -0.3955 1.3889 2.0308 1.9186 0.2578 -0.2901 2.0308 3.6820 3.4810 -0.5321 -0.4082 1.9186 3.4810 3.3062 -0.5848 0.2607 0.2578 -0.5321 -0.5848 1.6181 cov(A) ans = 1.0290 0.8066 1