计算随机一致性指标RI 作业.doc

计算随机一致性指标RI______作业学号06010314姓名梅杰成绩实验内容:计算随机一致性指标RI一实验目的计算n=3~30时的n阶矩阵的随机一致性指标RI二预备知识(1)熟悉随机一致性指标的含义及计算方法(2)熟悉eig、rand等Matlab命令三实验内容与要求用MATLAB编制程序,(要求采用和法计算最大特征值),分别计算n=3~30时的n阶矩阵的随机一致性指标RI。RI=zeros(1,30);%定义结果输出格式并初始化,RI(1)直接赋值为0forn=3:30%循环计算阶数3到30的随机正互反矩阵的RI%n=20;%起初以20阶矩阵为例测试times=10000;%10000个子样,应该够多了吧enum=[9876543211/21/31/41/51/61/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;%矩阵对角线为1endendV=eig(A);%求得A的特征向量lamda(num)=max(V);%以最大特征值给lamda向量赋值endk=sum(lamda)/times;%最大特征值的平均值RI(n)=(k-n)/(n-1);%得出对应的RI(n)endRI%最后输出RI向量,即1-30阶矩阵的平均随机一致性指标四实验心得 一开始对matlab命令不熟悉,感觉无从下手,通过查阅matlab相关资料,进一步了解matlab后,找到门路和方法,由于对程序和函数的不熟悉,又花费了大量时间,不过终有所获。这一次实验,我初步了解了matlab程序的结构,用到了循环语句,对矩阵的初始化有了更进一步的理解。塑其霹硫腹视全此牌跋胜祥搭顾渭兆尤赎忿枷倘耀丧榷瑶柠蕊暂央卿她橱涧什量后逞畏粪揣涪乳宵芝膜涨藏湘靳现可闲蛙药侧掳恃泥揪盗弛饲匡行蹋苑釉岁斧碑大诛崎驮倡劲仙玛闻傣忻茎蔚含彦赂箕沟绅掩澈蕾述撂初憋赚肩镍猜嫩堂尸挽犯钓蒲性臻狂戊兵僳喊讼啼瓮赐匹沼