zhengshu=x(1:4,1:1000)
futou1=x(1:4,1001)
futou2=x(1:4,1002)
futou3=x(1:4,1003)
futou1=futou1/sqrt(sum(futou1.*futou1));
fu12=futou1'*futou2
futou1mo=sqrt(sum(futou1.*futou1));
fu12=fu12/futou1mo
futou2=futou2-fu12*futou1
futou2=futou2/sqrt(sum(futou2.*futou2));
w=futou2'*futou1
fu13=futou1'*futou3
fu13=fu13/futou1mo
futou3=futou3-fu13*futou1
fu23=futou2'*futou3
fu23=fu23/futou1mo
futou3=futou3-fu23*futou2
futou3=futou3/sqrt(sum(futou3.*futou3));