% 概率密度函数的初始聚类
% 平均分段
for k = 1:K
T = size(samples(k).data,1);
samples(k).segment=floor([1:T/N:T T+1]);
end
%对属于每个状态的向量进行K均值聚类,得到连续混合正态分布
for i = 1:N
%把相同聚类和相同状态的向量组合到一个向量中
vector = [];
for k = 1:K
seg1 = samples(k).segment(i);
seg2 = samples(k).segment(i+1)-1;
vector = [vector ; samples(k).data(seg1:seg2,:)];
end
mix(i) = getmix(vector, M(i));
end
我写的和你一样啊,但是floor([1:T/N:T T+1]);是由N+1项的,下面循环for i = 1:N-1应该是for i = 1:N。。。