suiji6.m

这是数模中有关层次分析的源码

n=15000;
a=[];
for i=1:n
    m=floor(rand(1,15)*9+1);
    for l=1:15
        if floor(rand*2+1)==1
            m(l)=1/m(l);
        end
    end
    m=log(m);
    p=[(m(1)+m(6)-m(2))^2 (m(1)+m(7)-m(3))^2 (m(1)+m(8)-m(4))^2 (m(1)+m(9)-m(5))^2 (m(2)+m(10)-m(3))^2 (m(2)+m(11)-m(4))^2 (m(2)+m(12)-m(5))^2 (m(3)+m(13)-m(4))^2 (m(3)+m(14)-m(5))^2 (m(4)+m(15)-m(5))^2];
    s=sum(p);
    a=[a s];
end
b=sum(a)/n