کد:
clear all;
clc;
X = xlsread('input.xlsx',1,'B1');
Y = xlsread('input.xlsx',1,'B2');
Z = xlsread('input.xlsx',1,'B3');
TXY = xlsread('input.xlsx',1,'B4');
TXZ = xlsread('input.xlsx',1,'B5');
TYZ = xlsread('input.xlsx',1,'B6');
I1=X+Y+Z;
I2=X*Y+X*Z+Y*Z-TXY^2-TXZ^2-TYZ^2;
I4=[X TXY TXZ;TXY Y TYZ;TXZ TYZ Z];
I3=det(I4);
a=I1;
b=I2;
c=I3;
p=[1 -a b -c];
r=roots(p);
k=r(1);
l=r(2);
m=r(3);
w=[k l m];
Q=sort(w,'descend')
for i=1:3
a(i)=det([(Y-Q(i)) TYZ ; TYZ (Z-Q(i)) ]);
b(i)=-det([TXY TYZ ;TXZ (Z-Q(i)) ]);
c(i)=det([TXY (Y-Q(i));TXZ TYZ ]);
k(i)=1/((a(i))^2+(b(i))^2+(c(i))^2)^.5;
ll(i)=a(i)*k(i);
mm(i)=b(i)*k(i);
nn(i)=c(i)*k(i);
aaaa(i)=(ll(i))^2+(mm(i))^2+(nn(i))^2;
end
F=[ll(1) mm(1) nn(1);ll(2) mm(2) nn(2);ll(3) mm(3) nn(3)]
G=['l(1)' 'm(1)' 'n(1)';'l(2)' 'm(2)' 'n(2)';'l(3)' 'm(3)' 'n(3)']
for nnn=1:9
string=strcat(G(nnn),'=',num2str((F(nnn))));
disp(string);
end
clear all;
clc;
X = xlsread('input.xlsx',1,'B1');
Y = xlsread('input.xlsx',1,'B2');
Z = xlsread('input.xlsx',1,'B3');
TXY = xlsread('input.xlsx',1,'B4');
TXZ = xlsread('input.xlsx',1,'B5');
TYZ = xlsread('input.xlsx',1,'B6');
I1=X+Y+Z;
I2=X*Y+X*Z+Y*Z-TXY^2-TXZ^2-TYZ^2;
I4=[X TXY TXZ;TXY Y TYZ;TXZ TYZ Z];
I3=det(I4);
a=I1;
b=I2;
c=I3;
p=[1 -a b -c];
r=roots(p);
k=r(1);
l=r(2);
m=r(3);
w=[k l m];
Q=sort(w,'descend')
for i=1:3
a(i)=det([(Y-Q(i)) TYZ ; TYZ (Z-Q(i)) ]);
b(i)=-det([TXY TYZ ;TXZ (Z-Q(i)) ]);
c(i)=det([TXY (Y-Q(i));TXZ TYZ ]);
k(i)=1/((a(i))^2+(b(i))^2+(c(i))^2)^.5;
ll(i)=a(i)*k(i);
mm(i)=b(i)*k(i);
nn(i)=c(i)*k(i);
aaaa(i)=(ll(i))^2+(mm(i))^2+(nn(i))^2;
end
F=[ll(1) mm(1) nn(1);ll(2) mm(2) nn(2);ll(3) mm(3) nn(3)]
G=['l(1)' 'm(1)' 'n(1)';'l(2)' 'm(2)' 'n(2)';'l(3)' 'm(3)' 'n(3)']
for nnn=1:9
string=strcat(G(nnn),'=',num2str((F(nnn))));
disp(string);
end