?

%求回歸方程 clc, clear %x1,...,x6的數(shù)據(jù)和權(quán)重數(shù)據(jù) % aw = xlsread('E:\a-建模\第六輪\結(jié)果\附件2.汽泵前置泵部分監(jiān)測數(shù)據(jù).xls', 2,'J3:P3603'); aw = xlsread('E:\a-建模\第六輪\結(jié)果\附件3.汽輪機本體部分監(jiān)測數(shù)據(jù).xlsx', 2,'A17:X5938'); % aw = xlsread('E:\a-建模\第六輪\結(jié)果\附件4.一次風(fēng)機部分監(jiān)測數(shù)據(jù).xls', 2,'L21:P3412'); %把x1,...,x6的數(shù)據(jù)和權(quán)重數(shù)據(jù)保存在純文本文件zhb.txt中 w=aw(end,:); %提取權(quán)重向量 a=aw([1:end-1],:); %提取指標(biāo)數(shù)據(jù) a(:,[1,3])=-a(:,[1,3]); %把成本型指標(biāo)轉(zhuǎn)換成效益型指標(biāo),針對附件三 a(:,:) = -a(:,:); ra=tiedrank(a) %對每個指標(biāo)值分別編秩，即對a的每一列分別編秩 [n,m]=size(ra); %計算矩陣sa的維數(shù) RSR=mean(ra,2)/n %計算秩和比 W=repmat(w,[n,1]); WRSR=m*mean(ra.*W,2)/n %計算加權(quá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 %計算標(biāo)準(zhǔn)正態(tài)分布的p分位數(shù)+5 X=[ones(n,1),Probit']; %構(gòu)造一元線性回歸分析的數(shù)據(jù)矩陣 [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

　　

?

總結(jié)

以上是生活随笔為你收集整理的3-matlab求秩和比代码的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。