?

%求回歸方程 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求秩和比代码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。