mx=[4.10 2.36 2.31 2.14 3.35 1.66 1.83 2.65 1.20 2.12 1.29 1.52 1.89 1.13 1.55 2.32 3.52 2.48 2.20 1.27 3.00 1.46 1.13 1.20 1.42 3.80 1.33 2.15 1.84 1.17 1.74 1.63 3.25 1.41 1.53 5.00 4.00 1.56 2.55 2.74 1.46 3.20 1.68 1.81 2.85 4.22 1.92 1.14 3.80 3.08 4.00 3.45 3.10 3.28 1.47 3.65 3.50 1.45 2.27 3.00 3.35 2.40 1.26 3.62 2.80 3.15 3.41 3.25 2.82 1.12 1.56 2.10 3.40 1.12 1.92 3.05 1.45 4.40 3.35 3.35 3.82 1.61 1.70 2.37 3.30 3.58 2.27 4.15 1.55 2.80 3.15 4.15 2.00 1.68 3.55 3.90 2.67 9.25 1.94 2.72 1.40 3.15 2.80 2.50 2.95 2.75 2.85 2.51 1.75 3.88 1.63 2.20 1.12 1.51 1.70 3.65 1.21 2.90 2.28 3.45 1.39 2.74 2.35 1.20 2.50 3.02 1.30 3.90 1.24 4.55 1.61 4.25 1.61 3.05 1.77 1.97 1.94 1.97 1.66 3.20 4.06 2.44 2.11 1.29 5.50 4.05 3.58 2.12 3.15 2.20 3.62 3.40 1.76 4.00 1.21 3.60 2.60 1.13 3.25 2.50 3.30 1.87 1.50 3.27 5.25 2.95 7.40 2.67 2.08 3.20 1.14 1.93 3.00 2.60 2.12 5.20 2.90 3.02 5.85 3.50 2.13 2.53 3.45 1.13 2.10 1.30 1.74 3.04 1.81 5.10 3.65 1.97 8.35 3.00 2.10 3.15 1.59 1.11 4.85 2.70 4.20 1.62 1.34 1.93 3.03 3.25 1.85 3.25 4.30 3.45 2.65 2.10 2.10 3.53 1.58 3.75 1.60 2.15 5.10 3.70 3.30 3.70 4.80 1.47 1.19 1.39 1.44 3.70 1.56 3.20 1.97 1.08 1.82 1.43 1.52 2.56 2.78 3.93 3.50 3.35 1.38 3.60 2.58 3.30 1.49 1.69 3.00 2.15 3.15 3.85 1.45 3.50 1.69 2.53 3.50 1.57 1.34 1.32 5.90 1.65 2.25 2.73 3.05 1.56 3.40 1.42 2.92 1.16 2.67 2.20 1.88 2.86 3.05 3.75 2.85 3.16 2.17 2.13 3.30 1.66 1.26 1.66 2.65 2.10 5.70 2.15 3.10 2.90 3.40 2.78 3.30 1.73 1.97 3.05 5.30 2.90 1.43 2.07 1.39 5.70 1.21 3.10 5.05 2.98]; %给定已知的序列
for i=1:304
X(i)=mx(i);
end
%计算直径D
D=zeros(304,304);
for i=1:303
for j=2:304
if(isx=0;
d=0;
for t=i:j
sx=sx+X(t);
end
meanx=sx/(j-i+1);
for k=i:j
d=d+(X(k)-meanx).^2;
end
D(j:i)=d;
end
end
end
D;
%计算分类损失函数
Lb=zeros(304,303);
J=zeros(304,303);
for k=2:303
for l=3:304
if(kmin=1000;
if(k==2) %k=2时的最小损失函数
for j=2:1
temp=D(j-1,1)+D(l,j);
if(tempmin=temp;
minj=j;
end
end
else %k>2时的最小损失函数
for j=k:l
temp=Lb(j-1,k-1)+D(l,j);
if(tempmin=temp;
minj=j;
end
end
end
Lb(l,k)=min;
J(l,k)=minj;
end
end
end
Lb;
J;
%绘制分类损失函数随k变化趋势图
y=zeros(1,38);
for i=1:38
y(i)=Lb(304,i+1);
end
y;
k=2:39;
plot(k,y);hold on
plot(k,y,'r*')
xlabel('分类数')
ylabel('L[b(n,k)]')
函数或变量 'minj' 无法识别。