(完整word版)java求矩陣的特征值和特征向量(AHP層次分析法計算權(quán)重)(附源代碼)

（完整word版）java求矩陣的特征值和特征向量（AHP層次分析法計算權(quán)重）（附源代碼）java求矩陣的特征值和特征向量（AHP層次分析法計算權(quán)重）（附源代碼）這幾天做一個項目，需要用到求矩陣的特征值特征向量.我C++學(xué)的不好，所以就去網(wǎng)站找了很多java的源代碼，來實(shí)現(xiàn)這個功能。很多都不完善,甚至是不準(zhǔn)確。所以自己參考寫了一個。這個用于我一個朋友的畢業(yè)設(shè)計.結(jié)果肯定正確。話不多說，貼源代碼!importjava。math。BigDecimal;importjava.util。Arrays；/**AHP層次分析法計算權(quán)重＊@sineejdk1。6木@author劉興＊@version1.0@date2012。05.25**/publicclassAHPComputeWeight{/***@paramargs*/publicstaticvoidmain(String[]args){/木木a為N*N矩陣*///double】][]a={{1,1,1},{1,1,1},{1,1,1}};double]][]a二{{1,3,5},{2,3,1,},{4,7,3}}；//double】][]a={{1,1/5,1/3},{5,1,1},{3,1,1}}；//double[]】]a={{1,1/2,2,1},{2,1,3,4},{1/2,1/3,1,1},{1,1/4,1,1}}//double[][]a={{1,0.5,0.5},{2,1,1},{2,1,1}}；//double】][]a={{1,1/4,1/3,1},{4,1,3,5},{3,1/3,1,4},{1,1/5,1/41}}；//double[][]a={{1,2,3,5},{0。5,1,2,3},{0。33,0.5,1,2},{0。2,0.33,0.5,1}};intN=a[0].length；double】]weight=newdouble】N]；AHPComputeWeightinstance=AHPComputeWeight。getlnstance();instance。weight(a,weight,N);System。out.println(Arrays.toString(weight));}//單例privatestaticfinalAHPComputeWeightacw=newAHPComputeWeight()；//平均隨機(jī)一致性指針privatedouble[]Rl={0。00,0。00,0.58,0.90,1.12,1.21,1。32,1。411.45,1.49};//隨機(jī)一致性比率(完整word版)java求矩陣的特征值和特征向量(AHP層次分析法計算權(quán)重)(附源代碼)privatedoubleCR=0.0;//最大特征值privatedoublelamta=0。0;/＊*＊私有構(gòu)造*/privateAHPComputeWeight(){}/*＊*返回單例*木@return＊/publicstaticAHPComputeWeightgetInstance(){returnacw；}/*＊*計算權(quán)重＊＊@parama＊@paramweight＊@paramN＊/publicvoidweight(double[][]a，double[]weight,intN){//初始向量Wkdouble[]w0=newdouble[N];for(inti=0;i〈N;i++){w0[i]=1。0/N；}//一般向量W(k+1)double[]w1=newdouble[N];//W(k+1)的歸一化向量double[]w2=newdouble[N]；doublesum=1.0；doubled=1.0;//誤差doubledelt=0.00001;（完整word版）java求矩陣的特征值和特征向量（AHP層次分析法計算權(quán)重）（附源代碼）while（d〉delt）{d=0。0;sum=0;//獲取向量intindex=0;for（intj=0；j〈N；j++）{doublet=0.0；for（intl=0；l<N;l++）t+=a[j][l]*w0[l]；//w1[j]=a[j][0]＊w0[0]+a[j][1]＊w0[1]+a[j][2]＊w0[2]；w1[j]=t；sum+=w1[j];}//向量歸一化for（intk=0；k〈N；k++）{w2[k]=w1[k]/sum;//最大差值d=Math.max（Math.abs（w2[k]-w0[k]），d）；//用于下次迭代使用w0[k]=w2[k]；}}//計算矩陣最大特征值lamta.Cl,RIlamta=0.0；for（intk=0；k<N;k++）{lamta+=w1[k]/（N*w0[k]）;}doubleCI=（lamta—N）/（N-1）;if（RI[N-1]!=0）{CR=CI/RI[N—1]}//四舍五入處理lamta=round（lamta,3）；CI=Math.abs（round（CI,3））；CR=Math.abs（round（CR,3））;for（inti=0；i<N；i++）{w0[i]=round（w0[i],4）；w1[i]=round（w1[i],4）;（完整word版）java求矩陣的特征值和特征向量（AHP層次分析法計算權(quán)重）（附源代碼）w2[i]=round（w2[i],4）；}//控制臺打印輸出System。out.println（"lamta二”+lamta）;System。out.println（”CI二"+Cl）;System。out。println（"CR二”+CR）;//控制臺打印權(quán)重System.out。println(”w0[]二"”；TOC\o"1-5"\h\zfor(inti二0；i〈N;i++){System。out。print(w0[i]+””)；}System。out.println(”")；System.out。println(”w1[]=");for(inti=0；i〈N；i++;{System.out。print(w1[i]+”");}System。out。println("");System。out。println("w2[]=");for(inti=0;i<N；i++){weight[i]=w2[i];System.out。print(w2[i]+”")；}System.out.println(”";；}/*＊*四舍五入＊＊@paramv*@paramscale木@return＊/publicdoubleround(doublev,intscale){if(scale<0){thrownewIllegalArgumentException(”Thescalemustbeapositiveintegerorzero”;;}BigDecimalb=newBigDecimal(Double.toString(v))；BigDecimalone=newBigDecimal("1")；returnb.div