global M;
E=eye(8,8);%生成8*8單位陣
for i=2:9
    N(i-1,:)=M(i,1:8)-M(1,1:8); %2-9行各指標(biāo)對(duì)第一行的差值
    c1(i-1)=M(i,9)/M(i,10)-M(1,9)/M(1,10);%2-9行負(fù)荷率對(duì)第一行的差值
    c2(i-1)=M(i,10)-M(1,10);%2-9行煤耗對(duì)第一行的差值
end 
for j=1:9
    t(j)=M(j,9)/M(j,10);  %負(fù)荷率
end
b1=c1';b2=c2';
X=N\b1; %各指標(biāo)對(duì)負(fù)荷率的影響系數(shù)
Y=N\b2; %各指標(biāo)對(duì)煤耗的影響系數(shù)
AA1=diag(N(1,:));%生成對(duì)角陣，對(duì)角元素為第2行與第1行各指標(biāo)差值
AA2=diag(N(2,:));%生成對(duì)角陣，對(duì)角元素為第3行與第1行各指標(biāo)差值
AA3=diag(N(3,:));%生成對(duì)角陣，對(duì)角元素為第4行與第1行各指標(biāo)差值
AA4=diag(N(4,:));%生成對(duì)角陣，對(duì)角元素為第5行與第1行各指標(biāo)差值
AA5=diag(N(5,:));%生成對(duì)角陣，對(duì)角元素為第6行與第1行各指標(biāo)差值
AA6=diag(N(6,:));%生成對(duì)角陣，對(duì)角元素為第7行與第1行各指標(biāo)差值
AA7=diag(N(7,:));%生成對(duì)角陣，對(duì)角元素為第8行與第1行各指標(biāo)差值
AA8=diag(N(8,:));%生成對(duì)角陣，對(duì)角元素為第9行與第1行各指標(biāo)差值
AAA=[AA1^3,AA1^2,AA1^1,E;
    AA2^3,AA2^2,AA2^1,E;
    AA3^3,AA3^2,AA3^1,E;
    AA4^3,AA4^2,AA4^1,E;
    AA5^3,AA5^2,AA5^1,E;
    AA6^3,AA6^2,AA6^1,E;
    AA7^3,AA7^2,AA7^1,E;
    AA8^3,AA8^2,AA8^1,E];%系數(shù)矩陣
BBB=[E,E,E,E,E,E,E,E];
qwe=pinv(AAA');
XXX=BBB*qwe;
XA=XXX(1:8,:);
XB=XXX(9:16,:);
XC=XXX(17:24,:);
XD=XXX(25:32,:);
XF=XXX(33:40,:);
XG=XXX(41:48,:);
XH=XXX(49:56,:);
oy=[14.84-15.5,-2,-0.75,1,1.13,-3,0.8,-5];
AA=diag(oy);
oz=oy*(XA*AA^6+XB*AA^5+XC*AA^4+XD*AA^3+XF*AA^2+XG*AA+XH)
T=oz*X