GCWG=[23.75 	23.81 	23.69 	23.73 	23.75 	23.79 	23.81 	23.88 	23.79 
24.53 	24.80 	24.65 	25.69 	25.81 	25.77 	25.87 	25.30 	25.42 
26.20 	26.58 	26.37 	28.39 	28.78 	28.79 	29.02 	28.44 	28.21 
28.17 	28.97 	28.54 	31.32 	31.82 	31.66 	31.78 	31.75 	31.58 
30.04 	31.25 	30.71 	32.67 	32.56 	32.71 	32.66 	32.78 	32.62 
31.92 	32.48 	32.42 	32.01 	32.12 	32.02 	32.12 	32.15 	32.30 ];
for i=1:5
    GPW(i,:)=(GCWG(i,:)+GCWG(i+1,:))/2;
end
GCOG=[12.48 	13.40 	13.65 	11.80 	12.28 	10.78 	10.95 	10.05 	10.75 
13.05 	14.35 	15.08 	12.28 	12.90 	12.25 	11.45 	10.95 	12.50 
13.35 	15.18 	16.08 	12.48 	14.02 	13.35 	11.75 	12.01 	14.35 
13.81 	15.88 	17.55 	13.05 	15.00 	14.35 	12.40 	12.65 	15.60 
14.00 	16.50 	18.72 	13.15 	15.90 	15.00 	12.75 	13.01 	16.95 
14.15 	17.51 	19.34 	13.25 	16.35 	15.30 	12.78 	13.38 	17.80 ];
for i=1:5
    GPO(i,:)=(GCOG(i,:)+GCOG(i+1,:))/2;
end
DD=(GPW-GPO)./GPO*100;

load matlab.mat %装载双轴拉伸数据
l=[4.62	 4.52	4.48	6.38	6.35	6.06	8.00	7.80	7.66]/100;%转化成米单位
t=[0.00115,0.00111,0.00108,0.00128,0.00123,0.00116,0.00143,0.00129,0.0013];
%%%%%%%%%第一只袜子AB段 对应论文里是sample A的 part1 后面以此类推
TempAB1=find (AB1(:,1)>=20 & AB1(:,1)<=30);%寻找形变在20%到30%之间的力
ABT1=AB1(TempAB1,2)/(t(1)*l(1))*9.8;%计算应力
YAB1=AB1(TempAB1,1);%导出拉伸
P=polyfit(YAB1,ABT1,1);
ME(1,1)=P(1);

%%%%%%%%%%%%%%%%%AB2
TempAB2=find (AB2(:,1)>=20 & AB2(:,1)<=30);%寻找形变在20%到30%之间的力
ABT2=AB2(TempAB2,2)/(t(2)*l(2))*9.8;%计算应力
YAB2=AB2(TempAB2,1);%导出拉伸
P=polyfit(YAB2,ABT2,1);
ME(1,2)=P(1);%最终的E

%%%%%%%%%%%%%%%%%AB3
TempAB3=find (AB3(:,1)>=20 & AB3(:,1)<=30);%寻找形变在20%到30%之间的力
ABT3=AB3(TempAB3,2)/(t(3)*l(3))*9.8;%计算应力
YAB3=AB3(TempAB3,1);%导出拉伸
P=polyfit(YAB3,ABT3,1);
ME(1,3)=P(1);

%%%%%%%%%%%%%%%%%AB4
TempAB4=find (AB4(:,1)>=20 & AB4(:,1)<=30);%寻找形变在20%到30%之间的力
ABT4=AB4(TempAB4,2)/(t(4)*l(4))*9.8;%计算应力
YAB4=AB4(TempAB4,1);%导出拉伸
P=polyfit(YAB4,ABT4,1);
ME(1,4)=P(1);

%%%%%%%%%%%%%%%%%AB5
TempAB5=find (AB5(:,1)>=20 & AB5(:,1)<=30);%寻找形变在20%到30%之间的力
ABT5=AB5(TempAB5,2)/(t(5)*l(5))*9.8;%计算应力
YAB5=AB5(TempAB5,1);%导出拉伸
P=polyfit(YAB5,ABT5,1);
ME(1,5)=P(1);

%%%%%%%%%%%%%%%%%AB6
TempAB6=find (AB6(:,1)>=20 & AB6(:,1)<=30);%寻找形变在20%到30%之间的力
ABT6=AB6(TempAB6,2)/(t(6)*l(6))*9.8;%计算应力
YAB6=AB6(TempAB6,1);%导出拉伸
P=polyfit(YAB6,ABT6,1);
ME(1,6)=P(1);

%%%%%%%%%%%%%%%%%AB7
TempAB7=find (AB7(:,1)>=20 & AB7(:,1)<=30);%寻找形变在20%到30%之间的力
ABT7=AB7(TempAB7,2)/(t(7)*l(7))*9.8;%计算应力
YAB7=AB7(TempAB7,1);%导出拉伸
P=polyfit(YAB7,ABT7,1);
ME(1,7)=P(1);

%%%%%%%%%%%%%%%%%AB8
TempAB8=find (AB8(:,1)>=20 & AB8(:,1)<=30);%寻找形变在20%到30%之间的力
ABT8=AB8(TempAB8,2)/(t(8)*l(8))*9.8;%计算应力
YAB8=AB8(TempAB8,1);%导出拉伸
P=polyfit(YAB8,ABT8,1);
ME(1,8)=P(1);

%%%%%%%%%%%%%%%%%AB9
TempAB9=find (AB9(:,1)>=20 & AB9(:,1)<=30);%寻找形变在20%到30%之间的力
ABT9=AB9(TempAB9,2)/(t(9)*l(9))*9.8;%计算应力
YAB9=AB9(TempAB9,1);%导出拉伸
P=polyfit(YAB9,ABT9,1);
ME(1,9)=P(1);

%%%%%%BC1
TempBC1=find (BC1(:,1)>=20 & BC1(:,1)<=30);%寻找形变在20%到30%之间的力
BCT1=BC1(TempBC1,2)/(t(1)*l(1))*9.8;%计算应力
YBC1=BC1(TempBC1,1);%导出拉伸
P=polyfit(YBC1,BCT1,1);
ME(2,1)=P(1);

%%%%%%%%BC2
TempBC2=find (BC2(:,1)>=20 & BC2(:,1)<=30);%寻找形变在20%到30%之间的力
BCT2=BC2(TempBC2,2)/(t(2)*l(2))*9.8;%计算应力
YBC2=BC2(TempBC2,1);%导出拉伸
P=polyfit(YBC2,BCT2,1);
ME(2,2)=P(1);

%%%%%%%%BC3
TempBC3=find (BC3(:,1)>=20 & BC3(:,1)<=30);%寻找形变在20%到30%之间的力
BCT3=BC3(TempBC3,2)/(t(3)*l(3))*9.8;%计算应力
YBC3=BC3(TempBC3,1);%导出拉伸
P=polyfit(YBC3,BCT3,1);
ME(2,3)=P(1);

%%%%%%%%BC4
TempBC4=find (BC4(:,1)>=20 & BC4(:,1)<=30);%寻找形变在20%到30%之间的力
BCT4=BC4(TempBC4,2)/(t(4)*l(4))*9.8;%计算应力
YBC4=BC4(TempBC4,1);%导出拉伸
P=polyfit(YBC4,BCT4,1);
ME(2,4)=P(1);

%%%%%%%%BC5
TempBC5=find (BC5(:,1)>=20 & BC5(:,1)<=30);%寻找形变在20%到30%之间的力
BCT5=BC5(TempBC5,2)/(t(5)*l(5))*9.8;%计算应力
YBC5=BC5(TempBC5,1);%导出拉伸
P=polyfit(YBC5,BCT5,1);
ME(2,5)=P(1);

%%%%%%%%BC6
TempBC6=find (BC6(:,1)>=20 & BC6(:,1)<=30);%寻找形变在20%到30%之间的力
BCT6=BC6(TempBC6,2)/(t(6)*l(6))*9.8;%计算应力
YBC6=BC6(TempBC6,1);%导出拉伸
P=polyfit(YBC6,BCT6,1);
ME(2,6)=P(1);

%%%%%%%%BC7
TempBC7=find (BC7(:,1)>=20 & BC7(:,1)<=30);%寻找形变在20%到30%之间的力
BCT7=BC7(TempBC7,2)/(t(7)*l(7))*9.8;%计算应力
YBC7=BC7(TempBC7,1);%导出拉伸
P=polyfit(YBC7,BCT7,1);
ME(2,7)=P(1);

%%%%%%%%BC8
TempBC8=find (BC8(:,1)>=20 & BC8(:,1)<=30);%寻找形变在20%到30%之间的力
BCT8=BC8(TempBC8,2)/(t(8)*l(8))*9.8;%计算应力
YBC8=BC8(TempBC8,1);%导出拉伸
P=polyfit(YBC8,BCT8,1);
ME(2,8)=P(1);

%%%%%%%%BC9
TempBC9=find (BC9(:,1)>=20 & BC9(:,1)<=30);%寻找形变在20%到30%之间的力
BCT9=BC9(TempBC9,2)/(t(9)*l(9))*9.8;%计算应力
YBC9=BC9(TempBC9,1);%导出拉伸
P=polyfit(YBC9,BCT9,1);
ME(2,9)=P(1);

%%%%%%CD1
TempCD1=find (CD1(:,1)>=20 & CD1(:,1)<=30);%寻找形变在20%到30%之间的力
CDT1=CD1(TempCD1,2)/(t(1)*l(1))*9.8;%计算应力
YCD1=CD1(TempCD1,1);%导出拉伸
P=polyfit(YCD1,CDT1,1);
ME(3,1)=P(1);

%%%%%%CD2
TempCD2=find (CD2(:,1)>=20 & CD2(:,1)<=30);%寻找形变在20%到30%之间的力
CDT2=CD2(TempCD2,2)/(t(2)*l(2))*9.8;%计算应力
YCD2=CD2(TempCD2,1);%导出拉伸
P=polyfit(YCD2,CDT2,1);
ME(3,2)=P(1);

%%%%%%CD3
TempCD3=find (CD3(:,1)>=20 & CD3(:,1)<=30);%寻找形变在20%到30%之间的力
CDT3=CD3(TempCD3,2)/(t(3)*l(3))*9.8;%计算应力
YCD3=CD3(TempCD3,1);%导出拉伸
P=polyfit(YCD3,CDT3,1);
ME(3,3)=P(1);

%%%%%%CD4
TempCD4=find (CD4(:,1)>=20 & CD4(:,1)<=30);%寻找形变在20%到30%之间的力
CDT4=CD4(TempCD4,2)/(t(4)*l(4))*9.8;%计算应力
YCD4=CD4(TempCD4,1);%导出拉伸
P=polyfit(YCD4,CDT4,1);
ME(3,4)=P(1);

%%%%%%CD5
TempCD5=find (CD5(:,1)>=20 & CD5(:,1)<=30);%寻找形变在20%到30%之间的力
CDT5=CD5(TempCD5,2)/(t(5)*l(5))*9.8;%计算应力
YCD5=CD5(TempCD5,1);%导出拉伸
P=polyfit(YCD5,CDT5,1);
ME(3,5)=P(1);

%%%%%%CD6
TempCD6=find (CD6(:,1)>=20 & CD6(:,1)<=30);%寻找形变在20%到30%之间的力
CDT6=CD6(TempCD6,2)/(t(6)*l(6))*9.8;%计算应力
YCD6=CD6(TempCD6,1);%导出拉伸
P=polyfit(YCD6,CDT6,1);
ME(3,6)=P(1);

%%%%%%CD7
TempCD7=find (CD7(:,1)>=20 & CD7(:,1)<=30);%寻找形变在20%到30%之间的力
CDT7=CD7(TempCD7,2)/(t(7)*l(7))*9.8;%计算应力
YCD7=CD7(TempCD7,1);%导出拉伸
P=polyfit(YCD7,CDT7,1);
ME(3,7)=P(1);

%%%%%%CD8
TempCD8=find (CD8(:,1)>=20 & CD8(:,1)<=30);%寻找形变在20%到30%之间的力
CDT8=CD8(TempCD8,2)/(t(8)*l(8))*9.8;%计算应力
YCD8=CD8(TempCD8,1);%导出拉伸
P=polyfit(YCD8,CDT8,1);
ME(3,8)=P(1);

%%%%%%CD9
TempCD9=find (CD9(:,1)>=20 & CD9(:,1)<=30);%寻找形变在20%到30%之间的力
CDT9=CD9(TempCD9,2)/(t(9)*l(9))*9.8;%计算应力
YCD9=CD9(TempCD9,1);%导出拉伸
P=polyfit(YCD9,CDT9,1);
ME(3,9)=P(1);

%%%%%%DE1
TempDE1=find (DE1(:,1)>=20 & DE1(:,1)<=30);%寻找形变在20%到30%之间的力
DET1=DE1(TempDE1,2)/(t(1)*l(1))*9.8;%计算应力
YDE1=DE1(TempDE1,1);%导出拉伸
P=polyfit(YDE1,DET1,1);
ME(4,1)=P(1);

%%%%%%DE2
TempDE2=find (DE2(:,1)>=20 & DE2(:,1)<=30);%寻找形变在20%到30%之间的力
DET2=DE2(TempDE2,2)/(t(2)*l(2))*9.8;%计算应力
YDE2=DE2(TempDE2,1);%导出拉伸
P=polyfit(YDE2,DET2,1);
ME(4,2)=P(1);

%%%%%%DE3
TempDE3=find (DE3(:,1)>=20 & DE3(:,1)<=30);%寻找形变在20%到30%之间的力
DET3=DE3(TempDE3,2)/(t(3)*l(3))*9.8;%计算应力
YDE3=DE3(TempDE3,1);%导出拉伸
P=polyfit(YDE3,DET3,1);
ME(4,3)=P(1);

%%%%%%DE4
TempDE4=find (DE4(:,1)>=20 & DE4(:,1)<=30);%寻找形变在20%到30%之间的力
DET4=DE4(TempDE4,2)/(t(4)*l(4))*9.8;%计算应力
YDE4=DE4(TempDE4,1);%导出拉伸
P=polyfit(YDE4,DET4,1);
ME(4,4)=P(1);

%%%%%%DE5
TempDE5=find (DE5(:,1)>=20 & DE5(:,1)<=30);%寻找形变在20%到30%之间的力
DET5=DE5(TempDE5,2)/(t(5)*l(5))*9.8;%计算应力
YDE5=DE5(TempDE5,1);%导出拉伸
P=polyfit(YDE5,DET5,1);
ME(4,5)=P(1);

%%%%%%DE6
TempDE6=find (DE6(:,1)>=20 & DE6(:,1)<=30);%寻找形变在20%到30%之间的力
DET6=DE6(TempDE6,2)/(t(6)*l(6))*9.8;%计算应力
YDE6=DE6(TempDE6,1);%导出拉伸
P=polyfit(YDE6,DET6,1);
ME(4,6)=P(1);

%%%%%%DE7
TempDE7=find (DE7(:,1)>=20 & DE7(:,1)<=30);%寻找形变在20%到30%之间的力
DET7=DE7(TempDE7,2)/(t(7)*l(7))*9.8;%计算应力
YDE7=DE7(TempDE7,1);%导出拉伸
P=polyfit(YDE7,DET7,1);
ME(4,7)=P(1);

%%%%%%DE8
TempDE8=find (DE8(:,1)>=20 & DE8(:,1)<=30);%寻找形变在20%到30%之间的力
DET8=DE8(TempDE8,2)/(t(8)*l(8))*9.8;%计算应力
YDE8=DE8(TempDE8,1);%导出拉伸
P=polyfit(YDE8,DET8,1);
ME(4,8)=P(1);

%%%%%%DE9
TempDE9=find (DE9(:,1)>=20 & DE9(:,1)<=30);%寻找形变在20%到30%之间的力
DET9=DE9(TempDE9,2)/(t(9)*l(9))*9.8;%计算应力
YDE9=DE9(TempDE9,1);%导出拉伸
P=polyfit(YDE9,DET9,1);
ME(4,9)=P(1);

%%%%%%EF1
TempEF1=find (EF1(:,1)>=20 & EF1(:,1)<=30);%寻找形变在20%到30%之间的力
EFT1=EF1(TempEF1,2)/(t(1)*l(1))*9.8;%计算应力
YEF1=EF1(TempEF1,1);%导出拉伸
P=polyfit(YEF1,EFT1,1);
ME(5,1)=P(1);

%%%%%%EF2
TempEF2=find (EF2(:,1)>=20 & EF2(:,1)<=30);%寻找形变在20%到30%之间的力
EFT2=EF2(TempEF2,2)/(t(2)*l(2))*9.8;%计算应力
YEF2=EF2(TempEF2,1);%导出拉伸
P=polyfit(YEF2,EFT2,1);
ME(5,2)=P(1);

%%%%%%EF3
TempEF3=find (EF3(:,1)>=20 & EF3(:,1)<=30);%寻找形变在20%到30%之间的力
EFT3=EF3(TempEF3,2)/(t(3)*l(3))*9.8;%计算应力
YEF3=EF3(TempEF3,1);%导出拉伸
P=polyfit(YEF3,EFT3,1);
ME(5,3)=P(1);

%%%%%%EF4
TempEF4=find (EF4(:,1)>=20 & EF4(:,1)<=30);%寻找形变在20%到30%之间的力
EFT4=EF4(TempEF4,2)/(t(4)*l(4))*9.8;%计算应力
YEF4=EF4(TempEF4,1);%导出拉伸
P=polyfit(YEF4,EFT4,1);
ME(5,4)=P(1);

%%%%%%EF5
TempEF5=find (EF5(:,1)>=20 & EF5(:,1)<=30);%寻找形变在20%到30%之间的力
EFT5=EF5(TempEF5,2)/(t(5)*l(5))*9.8;%计算应力
YEF5=EF5(TempEF5,1);%导出拉伸
P=polyfit(YEF5,EFT5,1);
ME(5,5)=P(1);

%%%%%%EF6
TempEF6=find (EF6(:,1)>=20 & EF6(:,1)<=30);%寻找形变在20%到30%之间的力
EFT6=EF6(TempEF6,2)/(t(6)*l(6))*9.8;%计算应力
YEF6=EF6(TempEF6,1);%导出拉伸
P=polyfit(YEF6,EFT6,1);
ME(5,6)=P(1);

%%%%%%EF7
TempEF7=find (EF7(:,1)>=20 & EF7(:,1)<=30);%寻找形变在20%到30%之间的力
EFT7=EF7(TempEF7,2)/(t(7)*l(7))*9.8;%计算应力
YEF7=EF7(TempEF7,1);%导出拉伸
P=polyfit(YEF7,EFT7,1);
ME(5,7)=P(1);

%%%%%%EF8
TempEF8=find (EF8(:,1)>=20 & EF8(:,1)<=30);%寻找形变在20%到30%之间的力
EFT8=EF8(TempEF8,2)/(t(8)*l(8))*9.8;%计算应力
YEF8=EF8(TempEF8,1);%导出拉伸
P=polyfit(YEF8,EFT8,1);
ME(5,8)=P(1);

%%%%%%EF9
TempEF9=find (EF9(:,1)>=20 & EF9(:,1)<=30);%寻找形变在20%到30%之间的力
EFT9=EF9(TempEF9,2)/(t(9)*l(9))*9.8;%计算应力
YEF9=EF9(TempEF9,1);%导出拉伸
P=polyfit(YEF9,EFT9,1);
ME(5,9)=P(1);

ME1=sort(7*ME,'descend');
for i=1:9
    Cylinder(:,i)=ME1(:,i).*DD(:,i)*t(i)./(GPW(:,i)*0.01)/1000;%转化成kpa
end