clc;
clear?all;
[filenamepathnamefilterindex]?=?uigetfile（‘*.txt‘?‘請讀入男生訓(xùn)練集‘）;
fileAddrs?=?[pathnamefilename];
[MH?MW]=textread（fileAddrs‘%f%f‘）;
[filenamepathnamefilterindex]?=?uigetfile（‘*.txt‘?‘請讀入女生訓(xùn)練集‘）;
fileAddrs?=?[pathnamefilename];
[FH?FW]=textread（fileAddrs‘%f%f‘）;
FA=[FH?FW];FA=FA‘;
MA=[MH?MW];MA=MA‘;
for?k=1:50
???NT（:k）=FA（:k）;
end
??for?k=51:100
???NT（:k）=MA（:k-50）;
??end?
??%C=2
??z1=NT（:1）;z2=NT（:51）;
??%
?for?k=1:100
?????a=z1b=z2
?????G1=zeros（2100）;G2=zeros（2100）;
?????for?i=1:100
??????d1=sqrt（（NT（1i）-z1（11））^2+（NT（2i）-z1（21））^2）;
??????d2=sqrt（（NT（1i）-z2（11））^2+（NT（2i）-z2（21））^2）;
??????if?d1??????????G1（:i）=NT（:i）;
??????else
??????????G2（:i）=NT（:i）;
??????end
??end
??G1（:find（sum（abs（G1）1）==0））=[];
??G2（:find（sum（abs（G2）1）==0））=[];
??z1=mean（G12）;
??z2=mean（G22）;
????if?isequal（az1）==1&&isequal（z2b）==1
?????break
??end
?end
?z1z2G1G2size（G12）size（G22）
?for?i=1:size（G12）
????x=G1（1i）;
????y=G1（2i）;
????plot（xy‘G.‘）;
????hold?on
?end
?for?i=1:size（G22）
????x=G2（1i）;
????y=G2（2i）;
????plot（xy‘R*‘）;
????hold?on
?end
plot（z1（11）z1（21）‘B+‘）
hold?on
plot（z2（11）z2（21）‘B+‘）
hold?on

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1255??2012-05-09?22:09??c.m

?????文件??????66204??2012-05-04?11:17??2?-?副本?（2）.jpg

?????文件??????31956??2012-05-04?11:16??2?-?副本?（1）.jpg

?????文件??????45398??2012-05-04?18:57??3.tif

?????文件???????1578??2012-05-09?22:09??c3.m

-----------?---------??----------?-----??----

???????????????146391????????????????????5