秩和比综合评价模型

秩和比综合评价模型，数学建模中的常用算法。适合广大初学者

clcclear
aw=load（‘shuju.txt‘）;%该数据横向是指标，纵向是指标每年的数据
a=aw（[1:end-1]:）;%取矩阵1到n-1行
w=aw（end:）;%取矩阵最后一行
a（:[26]）=-a（:[26]）;
%这里的2，3，4，5，6指标为负向指标，也就是指标值越小越好，为与正向指标同步，因此这里乘以个负号
ra=tiedrank（a）;
%[R tie] = tie（X）计算向量X中值的秩，如果有任何X值被束缚，tie计算它们的平均秩。
%返回值是对非参数测试信号秩和ranksum所要求的关系的调整，以及对Spearman秩相关的计算。
%例子：从最小到最大，两个20的值是2和3，所以它们都是2。5（平均2和3）:
%tiedrank（[10 20 30 40 20]）
%ans =1.0000    2.5000    4.0000    5.0000    2.5000 
[nm]=size（ra）;%求矩阵维度
RSR=mean（ra2）/n;
%mean求数组的平均数或者均值
%mean（A2）返回值为该矩阵的各行向量的均值（每行的平均值）
W=repmat（w[n1]）;%用于复制平铺矩阵，相当于讲w矩阵看成一个元素，形成n×1维的矩阵并用W矩阵记录
WRSR=sum（ra.*W2）/n;
[sWRSRind]=sort（WRSR）;%sort为排序函数
p=[1:n]/n;
p（end）=1-1/（4*n）;
Probit=norminv（p01）+5;
%norminv函数，我只看了第一种使用方式， X = NORMINV（PAS） ，其中，P取概率。A取均值。S取方差。
%然后返回值X就是指满足均值为A，方差为S的高斯分布的累计概率密度值。即F（a）=P，返回值就是a。
%至于F的含义：对连续函数，所有小于等于a的值，其出现概率的和为F（a）=P（xX=[ones（n1）Probit‘];
[ababintrrintstats]=regress（sWRSRX）;
%[b，bint，r，rint，statsl=regess（y，x，alpha）其中因变量数据向量y和自变量数据矩阵x按以下排列方式输入      
%对一元线性回归，取k=1即可。alpha为显著性水平（缺省时设定为0.05），
%输出向量b，bint为回归系数估计值和它们的置信区间，r，rint为残差及其置信区间，
%stats是用于检验回归模型的统计量，有三个数值%第一个是R2，其中R是相关系数，第二个是F统计量值，第三个是与统计量F对应的概率P，当P<α
WRSRfit=ab（1）+ab（2）*Probit;
y=[1983:1992]‘;
xlswrite（‘ex147.xlsx‘[y（ind）ra（ind:）sWRSR]1）;
xlswrite（‘ex147.xlsx‘[y（ind）ones（n1）[1:n]‘p‘Probit‘WRSRfit‘[n:-1:1]‘]2）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       9689  2018-04-01 14:42  秩和比综合评价\ex147.xlsx

     文件       1833  2018-04-01 14:41  秩和比综合评价\main.m

     文件         90  2018-04-03 22:25  秩和比综合评价\read me.txt

     文件        313  2018-03-14 09:11  秩和比综合评价\shuju.txt

     文件    1310779  2017-09-01 15:24  秩和比综合评价\基于秩和比法的经济发展综合评价方_省略_2014年湖南省宏观经济数据为例_罗冬良.caj

     文件     227143  2017-09-01 15:24  秩和比综合评价\基于秩和比综合评价法的继电保护可靠性分析_张智锐.caj

     文件     569966  2017-09-01 15:24  秩和比综合评价\秩和比法在碘缺乏病防治工作综合评价中的应用_李珊珊.caj

     目录          0  2018-04-03 22:25  秩和比综合评价

----------- ---------  ---------- -----  ----

              2119813                    8