%求回归方程
clc, clear
%x1,...,x6的数据和权重数据
% aw = xlsread('E:\a-建模\第六轮\结果\附件2.汽泵前置泵部分监测数据.xls', 2,'J3:P3603');
aw = xlsread('E:\a-建模\第六轮\结果\附件3.汽轮机本体部分监测数据.xlsx', 2,'A17:X5938');
% aw = xlsread('E:\a-建模\第六轮\结果\附件4.一次风机部分监测数据.xls', 2,'L21:P3412');
%把x1,...,x6的数据和权重数据保存在纯文本文件zhb.txt中
w=aw(end,:); %提取权重向量
a=aw([1:end-1],:); %提取指标数据
a(:,[1,3])=-a(:,[1,3]); %把成本型指标转换成效益型指标,针对附件三
a(:,:) = -a(:,:);
ra=tiedrank(a) %对每个指标值分别编秩,即对a的每一列分别编秩
[n,m]=size(ra); %计算矩阵sa的维数
RSR=mean(ra,2)/n  %计算秩和比
W=repmat(w,[n,1]);
WRSR=m*mean(ra.*W,2)/n  %计算加权秩和比
[sWRSR, ind] = sort(WRSR);
p=[1:n]/n;    %计算累积频率
p(end)=1-1/(4*n) %修正最后一个累积频率,最后一个累积频率按1-1/(4*n)估计
Probit=norminv(p,0,1)+5  %计算标准正态分布的p分位数+5
X=[ones(n,1),Probit'];  %构造一元线性回归分析的数据矩阵
[ab,abint,r,rint,stats]=regress(sWRSR,X)  %一元线性回归分析WRSRfit=ab(1)+ab(2)*Probit  %计算WRNR的估计值
% [sWRSRfit,ind]=sort(WRSRfit,'descend')  %对WRNR的估计值按照从大到小排序
myp = Probit';
myw = sort(WRSR);%方程检验
yuanshi = [0.754084649 0.829721651 0.734280548 0.733794539 0.744495448 0.744385963 0.72802354  0.661016138 0.89500779  0.895390987];
de = [0.701423111  0.783380847 0.678977312 0.67850058  0.689637931 0.689145903 0.674474351 0.630014879 0.875795584 0.879279344];
can = de - yuanshi;
xiangdui = (can) ./ yuanshi

  

3-matlab求秩和比代码相关推荐

  1. matlab求系统根轨迹代码_广州数控GSK RH06六轴焊接机器人的轨迹规划

    技术参数 运动范围 任务点位置 问题:如实现上图中的 运动轨迹(以基坐标系为参照),请给出相关的轨迹规划算法及其实现方式(结合运动学及动力学方程进行相关的轨迹点选取及工程实现),其中各个点的坐标为A= ...

  2. matlab求系统根轨迹代码_根轨迹法、PID参数整定和matlab指令计算

    收获 (1)理解根轨迹的概念及其在控制系统设计中的作用: (2)手绘根轨迹草图,以及如何使用极端及绘制根轨迹: (3)熟悉在反馈控制系统中应用广泛的关键部件:PID控制器: (4)理解根轨迹在参数设计 ...

  3. matlab求系统根轨迹代码_要想正确画出根轨迹,先搞清楚这8大法则再说!

    话说宝刀君这几年给学弟学妹们辅导自控时,学弟学妹们面对根轨迹,经常是一脸懵逼. 常常是拿到题后,零极点都标出来了,法则也写出来了,但就是根轨迹画不出来. 有时好不容易画出来了,返回去和法则一对照,忍不 ...

  4. 3、利用matlab求f(x)的一阶导函数(完整代码)

    已知,利用matlab求f(x)的一阶导函数 在 MATLAB 中,可以使用符号计算工具箱 Symbolic Math Toolbox 来求 $f(x)$ 的一阶导函数.具体步骤如下: 打开 MATL ...

  5. matlab 求矩阵秩,求矩阵秩的两种方法及MATLAB的应用

    摘    要: 高等代数是一门逻辑思维比较强和理论知识比较深的学科, 它具有丰富的数学知识, 涉及许多重要的数学思想, 其在数学领域的应用很广泛, 如行列式.矩阵的相关计算和求解线性方程组的解方面的应 ...

  6. matlab路面功率谱密度的代码,[转载]转:matlab求功率谱密度代码实例

    (如果不先对数据都减去一个均值,那么接近0Hz部分的功率会非常非常大~这是有害的,一般做数据处理的时候如果只关注交流量的话,就要先使得整个序列的均值为0,方法就是给每个数减去该序列的一个均值2012- ...

  7. matlab求阶乘函数的代码,用matlab编写阶乘函数,调用该函数生成1!,……10!的阶乘表...

    由热心网友 zzfg511bhk 提供的答案1: for i=1:10 A(i,:)=[i fac(i)] end 其中 fac(n)为求阶乘的函数保存在fac.m文件中 答:Matlab 2015R ...

  8. 在matlab中可以计算式子的最大值吗,matlab求最大值

    matlab遗传算法求最大值问题例题完整代码_IT认证_资格考试/认证_教育专区.matlab遗传算法求最大值问题例题完整代码 本文件的目的是减少您打字的烦恼 %主程序:用...... 11MATLA ...

  9. 精馏塔matlab,MATLAB图解精馏塔理论塔板数程序代码

    <MATLAB图解精馏塔理论塔板数程序代码>由会员分享,可在线阅读,更多相关<MATLAB图解精馏塔理论塔板数程序代码(6页珍藏版)>请在人人文库网上搜索. 1.MATLAB图 ...

  10. matlab 求n 的和,MATLAB求1的阶乘加到n的阶乘和 不要现有的函数,要自己编写出来的...

    点击查看MATLAB求1的阶乘加到n的阶乘和 不要现有的函数,要自己编写出来的具体信息 答:myfactorial = @(n) factorial(n) myfactorial = @(n) pro ...

最新文章

  1. Java项目:学生选课系统(java+javaweb+jdbc)
  2. 新书《活用UML-需求分析高手》 序
  3. 嵌入式 linux下利用backtrace追踪函数调用堆栈以及定位段错误
  4. linux libodbc.so.1,CentOS6.0虚拟机上安装nginx启动的错误---缺少libpcre.so.1共享库
  5. 避免使用FOR –反假战役
  6. 必须要知道Java如何取得当前路径
  7. 跟随冒号畅游编程世界
  8. 《上市公司信息披露电子化规范》简介
  9. SubSonic 零代码的数据访问
  10. 如何给main传参数
  11. 什么是3D打印?游戏建模具体怎么做?哪个更有发展
  12. JavaScript时间格式化工具函数
  13. 使用注册表清理远程桌面连接历史记录
  14. 泛微E8、E9二次开发、泛微开发获取流程文档主、明细表单值,提供泛微ecology8二次开发完整项目下载,泛微把流程文档内容推送HR、ERP、SAP操作,泛微与ERP、SAP、HR集成
  15. nginx在Linux系统安装
  16. 转载:【菜鸟专用】使用LaTeX轻松撰写精美个人简历
  17. 订单信息表和订单明细表
  18. 逻辑思维训练——假设法
  19. 3d Max人物动画学习笔记(二)蒙皮
  20. 【小学信息技术教资面试】教案模板

热门文章

  1. 有限域f9的特征是多少_密码学数学基本第十一讲有限域.ppt
  2. Python分析身份证所在地(省、市、区),生日,生肖,星座,代码,性别,校验码。
  3. uniapp开发App如何引入阿里巴巴矢量库图标
  4. (2)安装宝塔与docker及docker镜像下载加速
  5. 海思hitool工具使用
  6. 如何解决Maven导入Oracle驱动出现Missing artifact com.oracle:ojdbc14:jar:10.2.0.2.0,Mvaen:com.oracle:ojdbc14:10.
  7. linux的源码安装步骤(以安装nginx为例)
  8. 12款黑科技神器,务必低调使用
  9. css实现页面标签的跳转
  10. tar.bz2 解压命令