最小距离法分类

function unsupCla（inphotooutphotoweishu）
pn1=imread（inphoto）;
pn1=double（pn1）;
[mnz]=size（pn1）;
num=weishu;
yanse（1:121:3） = [02550; 00255; 2552550; 2550255; 0255255; 125125125; 12500; 01250; 00125; 1251250; 1250125; 0125125];;
%%
pn2（1:m1:n1:3）=0;              %定义存放分类后的文件矩阵
c_center（1:num1:z）=0;
for p=1:z
    mx=max（max（pn1（1:m1:np）））;
    mn=min（min（pn1（1:m1:np）））;
    for q=1:5
         c_center（pq）=（mx*q-mn*q）/（num+1）+mn;      %即分为五类
    end
end
new_c_center（1:z1:num）=0;     %定义存放新的类中心
A=ones（znum）;   %定义存放各象元到各类中心的光谱距离
B=ones（1num）;
%%
threod1=24;   %设置秩代最大次数
threod2=1;    %循环时相邻两次类中心最小差值
flag=0;       %赋予秩代次数初始值
%%
while （flag    lenth=B;  %定义各分类象元的个数
    he=A;     %定义存放各分类象素之和的矩阵
    for i=1:m
        for j=1:n
            for k=1:num
            %计算各象元分别到各类的光谱欧式距离
            distance=0;
            for dis=1:z
            distance=distance+（pn1（ijdis）-c_center（disk））^