- 在线时间
- 79 小时
- 最后登录
- 2017-11-5
- 注册时间
- 2014-12-3
- 听众数
- 11
- 收听数
- 0
- 能力
- 0 分
- 体力
- 2 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 68
- 相册
- 0
- 日志
- 0
- 记录
- 2
- 帖子
- 127
- 主题
- 9
- 精华
- 0
- 分享
- 0
- 好友
- 38
升级   66.32% TA的每日心情 | 开心 2016-10-1 10:51 |
---|
签到天数: 53 天 [LV.5]常住居民I
 群组: 2015美赛优秀论文解析 群组: 2015美赛备战交流群组 |
1体力
clc;
close;
clear all;
% 控制输出结果精度
format short;
% 原始数据
x=[308.58 310 295 346 367
195.4 189.9 187.2 205 222.7
24.6 21 12.2 15.1 14.57
20 25.6 23.3 29.2 30
19.98 19 22.3 23.5 27.655
170 174 197 216.4 235.8
57.55 70.74 76.8 80.7 89.85
88.56 70 85.38 99.83 103.4
11.19 13.28 16.82 18.9 22.8
4.03 4.26 4.34 5.06 5.78
13.7 15.6 13.77 11.98 13.95];
n1 = size(x,1);
% 数据标准化处理
for i = 1:n1
x(i,:) = x(i,:)/x(i,1);
end
%保存中间变量,亦可省略此步,将原始数据赋予变量data
data = x;
% 分离参考数列(母因素)
consult=data(6:n1,:);
m1=size(consult,1);
% 分离比较数列(子因素)
compare=data(1:6,:);
m2=size(compare,1);
for i=1:m1
for j=1:m1
t(j,:)=compare(j,:)-consult(i,:);
end
min_min=min(min(abs(t')));
max_max=max(max(abs(t')));
% 通常分辨率都是取0.5
resolution=0.5;
%计算关联系数
coefficient=(min_min+resolution*max_max)./(abs(t)+resolution*max_max);
% 计算关联度
corr_degree=sum(coefficient')/size(coefficient,2);
r(i,:)=corr_degree;
end
% 输出关联度值并绘制柱形图
r
bar(r,0.90);
axis tight;
% 从左至右依次是 固定资产投资,工业投资,农业投资,科技投资,交通投资
legend('固定投资','工业投资','农业投资','科技投资','交通投资');
% 以下程序是为了给x轴添加汉字标签
% 其基本原理是先去掉x轴上的固有标签,然后用文本标注x轴
% 去掉x轴上的默认标签
set(gca,'XTickLabel','');
% 设定X轴刻度的位置,这里有6个母因素
n=6;
% 这里注意 :x_range范围如果是[1 n]会导致部门柱形条不能显示出来
% 所以范围要缩一点
x_value = 1:1:n;
x_range = [0.6 n+.4];
% 获取当前图形的句柄
set(gca,'XTick',x_value,'XLim',x_range);
% 在x轴上标记6个母因数
profits={'国民收入','工业收入','农业收入','商业收入','交通收入','建筑业收入'};
y_range = ylim;
% 用文本标注母因数名称
handle_date = test(x_value,y_range(1)*ones(1,n)+.018,profits(1:1:n));
% 把文本的字体设置合适的格式和大小并旋转一定的角度
set(handle_date,'HorizontalAlignment','right','VerticalAlignment','top','Rotation',35,...
'fontname','Arial','fontsize',10.5);
% y轴标记
ylabel('y');
title('投资对收入的作用');
|
zan
|