原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学建模提高班第三次作业
【作业 1】
全国重点水泥企业经济效益综合评价例。
利用主成分综合评价全国重点水泥企业的经济效益。 原始数据 （数据来自 1984 年中国统计年鉴）
见表 5-10 。
1. 利用一般的 Matlab 命令将数据标准化；
2. 利用 Matlab 命令计算出相关矩阵的方差及特征值；
3. 写出主成分的线性表达式；
4. 利用统计工具箱中命令得到如上结果；
5. 利用 SPSS进行主成分分析；
6. 对得到的结果作出解释。
解：
1. 利用一般的 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
2. 利用 Matlab 命令计算出相关矩阵的方差及特征值
根据题意利用 MATLAB编程如下：
R=cov(x) % 注释： 由于数据已经过标准化处理，故 x 的协方差矩阵等于其相关系数矩阵，
即 R=corrcoef(x).
[V,D]=eig(R); % 注释： 函数 eig 的功能是对矩阵 R 进行正交对角化变换，矩阵 D 是以 R 的
特征值为对角元的对角矩阵（对角元按从小到大的顺序排列） ，矩阵 V 是正交变换
文档评论（0）