提高班第三次作业.doc

数学建模提高班第三次作业 【作业1】 全国重点水泥企业经济效益综合评价例。 利用主成分综合评价全国重点水泥企业的经济效益。原始数据（数据来自1984年中国统计年鉴）见表5-10。 1.利用一般的Matlab命令将数据标准化； 2. 利用Matlab命令计算出相关矩阵的方差及特征值； 3.写出主成分的线性表达式； 4.利用统计工具箱中命令得到如上结果； 5.利用SPSS进行主成分分析； 对得到的结果作出解释。 解： 利用一般的Matlab命令将数据标准化 根据题意利用MATLAB编程如下： load dsc.txt; [n,p]=size(dsc); MEAN=mean( dsc) ; %求各变量的均值 STD=std( dsc); %求各变量的标准差 MEAN=ones(n,p)*diag(MEAN); STD=ones(n,p)*diag(STD); x=(dsc-MEAN)./STD %原始数据标准化 求解得标准化后的数据为： x = -0.3767 -0.3579 -0.1136 -0.3667 -0.0388 -1.1777 -1.3633 -0.2771 0.0882 -0.2939 0.0508 -0.2873 0.3562 -1.1777 -0.4073 2.2531 -0.6046 -0.6229 0.0568 -0.5810 -0.4277 0.2944 1.6425 -0.7741 -1.8223 -1.7643 -1.6884 -1.7200 -1.1725 -0.1472 1.0792 -0.5482 1.5891 1.9944 1.2163 2.1414 1.2879 0.8833 -1.8658 0.6039 2.1248 0.9100 2.2981 1.1681 0.7553 -1.9138 -0.4190 1.6432 0.9641 0.5177 0.6216 0.5421 0.8587 0 0.2635 1.2591 -0.6277 -0.9322 -0.8054 -0.7824 -2.8714 1.9138 0.2565 -0.2771 0.1066 -0.1974 0.1210 -0.1911 0.3219 1.1777 1.0442 -0.0964 0.3083 0.3104 0.9802 0.4389 -0.0857 -0.1472 0.0975 -0.3223 -0.3660 -0.2021 -0.4453 -0.2843 0.1121 -0.7361 -1.0711 -0.7967 -0.5076 -0.2567 -0.9264 -0.4629 0.4385 -0.2944 -0.4471 -0.6160 -0.3998 -0.1222 -0.0164 -0.1454 -0.2095 0.8833 0.6795 -1.2711 0.4792 1.8441 -0.2316 1.4787 0.8449 0 -0.5452 -0.6838 -0.9556 -0.8270 -1.1177 -0.9481 -0.1698 0.4416 1.0559 -0.0964 利用Matlab命令计算出相关矩阵的方差及特征值 根据题意利用MATLAB编程如下： R=cov(x) %注释： 由于数据已经过标准化处理，故x的协方差矩阵等于其相关系数矩阵，即R=corrcoef(x). [V,D]=eig(R); %注释： 函数eig的功能是对矩阵R进行正交对角化变换，矩阵D是以R的特征值为对角元的对角矩阵（对角元按从小到大的顺序排列），矩阵V是正交变换矩阵。 DD=[ ]; %将特征值对角矩阵D改写为列向量DD for i=p:-1:1 %此处要注意eig函数的输出D中特征值的排列顺序 DD=[DD;D(i,i)]; end O