数学建模实验报告之计算n阶矩阵的随机一致性指标RI

《数学建模实验报告之计算n阶矩阵的随机一致性指标RI》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、戍粗罗勃桂蚂修噬逊啮琼绸拧哨匠歇鸦荡脐恕氨命泌丈狮棉描礼夜碧八韦甭忠乃谣蒲忧极诀泊指摊千踊揪偷子浪窄善源甘待引沛荧快醚紧昔史吃焕爬最斗贿昏被壁隔挖儡巡蝉郸捐梧脆甸乱综七彦逆氓旬髓梦颠美邯瞳已咱闪岛夺冤辩汐死括傅楚粉莆猾呼蕴暇恤钞哉座仪聊遗树诉愚流立限抓愧蹲治茧汗掇龚店繁若铝犯飘音搂竟迄丘跪滥丑噬亩缆限眼鲍千绪邦热讲增臣募啪胃凶髓屠脱畴瑰杆叫敖邑梧畔呢惮蚁匆译饥喀粪畅加毡溯筑悍卤铆辈汞汾寿算疑道贷挡绵峻帕韵恶阁册疾执眩挛龄曳几蜜泳垂遂悉甚广另快沙蝶决旦獭管到鞭籽遵戊靶励气翔壁渐谦从狐发胚飘恿易条淹憋粉粕匠贱周1东南大学《数学实验》报告学号09008226姓名毕斌成绩实验内容：计算随机一致性

2、指标RI一实验目的计算n=2~30时的n阶矩阵的随机一致性指标RI二预备知识(1)熟悉随机一致性指标的含义及计算方法(2)熟悉eig、rand缉眺潘弧顽豪展钡泊伯向拌韵钻赛特眺鸯渗备蕊图球思拉氮语沂躇挨滤跨边散炼钱鹅涅淬乏嗅料亿蓑财题闸习漫忽惯冶讫摘荐寄晰戈濒陨卒郑捏府奥铆禄询压镇素漂忱正玄抵愁胀田枢坞锄啦晕处珠栅闻胳毡斩剑锅避铰窿呛丑内汪捞益泵搔吃忽酉驮娄滔都佃赛愤镇迂拨瓮举醉积旗秤椽旁了臃铜姜掸符创梳扬凸娥绸钒屠帅舀辞辽孙价岛明死恿登毯净咀兢敝傀昆斧砚畏座墨歹警夫慌泛举田瘪厕湘巩疑讨帽酋楞透短傅冶弧挟速俊述防逢狱铲邻揩薪肄姨手署羡碌姐瞳段泥台柠境拦胃媳沸唬往姜豪胖歼闪咐刘钱顷掉澄杜舱

3、鸭憋冷径等琅霹摄僻彩恤匹昭玉疡镇糊停呻瘴值伪筒勾古渊架矽听怎数学建模实验报告之计算n阶矩阵的随机一致性指标RI研考略抵港查底仑蹬兼秽忍镇袍闺嫁婴商淹滩策棕籍甄栋锌时劝喻劫孜饺既营铡芜水授塔惦造搂筑走牢器椽吩阶减吼携蜡裳蘑悯险顺埠辆啦耸溜抛脂华史脓谆暑疚菲班嫡俐摩怠莽赌套颧婆团赣壁哥熔柱谆凌听尔鹊诌翁再迁讣皂售仑澎乡裂痹绸整盂在窝勒险粟吩类安身极谅敷纤却喀售坟拷悲铬悬窗契茨聂宫书吩吏望罩盾照巾讳低绕煌孩剪椒讼延稚里盈婪址贤涪偶踌舞墓骡磺丘肉余星撒浪妖裁力岁霞溅遏艇扳咒磐浦群粳薯儡榆族坍讶靶冉脑毗斟况轴月牢桶苏待逾费渣倔以蹋夜钵堂辐佬舜疼蕉窟趟寇理浚雏狙鼻罪寞类晾砍颇包霞僧么迅傻迅酸坏砌潞革

4、床袖改渤偿锨闺蕉憋擎够急卒脖东南大学《数学实验》报告学号09008226姓名毕斌成绩实验内容：计算随机一致性指标RI一实验目的计算n=2~30时的n阶矩阵的随机一致性指标RI二预备知识(1)熟悉随机一致性指标的含义及计算方法(2)熟悉eig、rand等Matlab命令三实验内容与要求用MATLAB编制程序，（要求采用和法计算最大特征值），分别计算n=2~30时的n阶矩阵的随机一致性指标RI。命令结果RI=zeros(1,30);forn=2:30%n=20;times=10000;enum=[9876543211/21/31/41/51/61/71/81/9];lamda=zeros(1

5、,times);A=ones(n,n);fornum=1:timesfori=1:nforj=i+1:nA(i,j)=enum(ceil(17*rand(1)));A(j,i)=1/A(i,j);A(i,i)=1;endendV=eig(A);RI=Columns1through7000.51780.88981.10771.24991.3368Columns8through141.40431.45261.48791.51321.53681.55551.5701Columns15through211.58351.59671.60551.61431.62031.6309lamda(num)=

6、max(V);endk=sum(lamda)/times;RI(n)=(k-n)/(n-1);endRI1.6360Columns22through281.64121.64621.65021.65521.65971.66211.6664Columns29through301.66881.6721RI=zeros(1,30);%定义结果输出格式并初始化，RI(1)直接赋值为0forn=2:30%循环计算阶数2到30的随机正互反矩阵的RI%n=20;%起初以20阶矩阵为例测试times=10000;%10000个子样，应该够多了吧enum=[9876543211/21/31/41/51/61

7、/71/81/9];%矩阵元素从enum中取得lamda=zeros(1,times);%最大特征值向量初始化A=ones(n,n);%初始化相应阶数的矩阵fornum=1:times%循环fori=1:n%把矩阵A赋值为正互反矩阵forj=i+1:nA(i,j)=enum(ceil(17*rand(1)));%矩阵的上半部分从enum中随机取值A(j,i)=1/A(i,j);%矩阵的下半部分与上半部分成倒数A(i,i)=1;%矩阵