2015-03-25 40 views
0

我有3组数据(AA,B1和C1),具有不同的长度和大小。我的代码的目的是能够计算数据之间的相似性得分,例如AA和B1,AA和C1,B1和C1之间的相似性得分。所以下面是我的代码,应该能够计算相似性分数,但循环中存在一些问题。对于每对不同长度的数据,只选择最高值。输出应该是AA-B1:0.2226,AA-C1:0.2037和B1-C1:0.1111,它们代表每对的相似性分数。Matlab循环计数相似度分数问题

实际上,我的代码的输出,特别是max_val{i}是根据配对的大小设定的。例如,大小为1 x 2的配对在max_val输出中应该有一个值,而不是三个。谢谢。

A1={[4,3,4,3,3]}; 
A2={[3,1,2,4]}; 
A3={[1,2,4]}; 
AA=[A1,A2,A3]; 
B1={[2,2,4,4]}; 
C1={[4,4,4,3,2,2]}; 

set={[AA],[B1],[C1]}; 
comb_set=nchoosek(set,2); %combinations of two sets 

for h=1:size(comb_set,1) 
comb_pair=comb_set(h,:)'; 
sets=comb_pair; 

cat=horzcat(sets{:}); 
c=reshape(repmat(sets{1},numel(sets{2}),1),numel(sets{1})*numel(sets{2}),1); 
d=repmat(sets{2}(:),length(sets{1}),1); 
pairs=[c d]; 
ind=cellfun(@numel,pairs(:,1)) > cellfun(@numel,pairs(:,2)); 
pairs(ind,[1 2]) = pairs(ind,[2 1]) %possible pairs of the row of subset 
p=cell(size(pairs,1),1); 

for i=1:size(pairs,1) 
%the two vectors 
[a,b]=deal(pairs{i,:}); 
%sliding window indices, and compute the sum 
idx=hankel(1:numel(a),numel(a):numel(b)); 
count_minus{i}=bsxfun(@minus,b(idx.'),a); %count minus between pairs 
count_total{i}=numel(a)+numel(b); %count total 
count_intersect{i}=sum(count_minus{i}'==0)'; %count no. of intersection 
union{i}=count_total{i}-count_intersect{i}; %union 
subset{i}=count_intersect{i}./union{i}; %subset each pair similarity score 
max_val{i}=max(subset{i}) %maximum similarity score 
bsum=cellfun(@(x) sum(x),max_val); 
total{i}=sum(bsum~=0); 
average=sum(bsum)/total{i} 
end 
end 

回答

1

你有2个为ih循环和内循环使用max_val{i}。这意味着相同的max_val单元结构将用于h的每个值 - 在此例中为size(comb_set,1)=3。对于h的每次迭代,上一次迭代中定义的max_val{i}将被覆盖。由于在第一次和第二次运行size(pairs,1)=3该单元将有一个长度为3.在最后一个size(pairs,1)=1。单元格的长度仍然为3,但您只会覆盖第一个元素 - 您可以查看您提供的代码的输出,并看到最后两个元素与前一个迭代中的最后两个元素相等。您需要以不同的方式定义max_val。例如,您可以定义一个数组以将最大值保留在i循环中,并将其写入max_values{h}

当你有这样的错误时,更容易生成一个最小的工作示例并检查出现了什么问题。这段代码很难阅读,也因为缺少缩进。 CTRL + A和CTRL +我让你的生活更轻松。