我在MATLAB中遇到以下错误:
???下标索引必须是实数正整数或 逻辑值。
==>中的错误Lloyd_Max at 74 D(w_count)= mean((x - 中心(XQ))^ 2);
这是我的代码:
function [ xq,centers,D ] = Lloyd_Max( x,N,min_value,max_value )
%LLOYD_MAX Summary of this function goes here
% Detailed explanation goes here
x = x';
temp = (max_value - min_value)/2^N;
count=1;
for j=0:temp:((max_value - min_value)-temp),
centers(count) = (j + j + temp )/2;
count = count + 1;
end
for i=1:length(centers),
k(i) = centers(i);
end
w_count = 0;
while((w_count < 2) || (D(w_count) - D(w_count - 1) > 1e-6))
w_count = w_count + 1;
count1 = 2;
for i=2:(count-1),
T(i) = (k(i-1) + k(i))/2;
count1 = count1 +1 ;
end
T(1) = min_value;
T(count1) = max_value;
index = 1;
for j=2:count1,
tempc = 0;
tempk = 0;
for k=1:10000,
if(x(k) >= T(j-1) && x(k) < T(j))
tempk = tempk + x(k);
tempc = tempc + 1;
end
end
k(index) = tempk;
k_count(index) = tempc;
index = index + 1;
end
for i=1:length(k),
k(i) = k(i)/k_count(i);
end
for i=1:10000,
if (x(i) > max_value)
xq(i) = max_value;
elseif (x(i) < min_value)
xq(i) = min_value;
else
xq(i) = x(i);
end
end
for i=1:10000,
cnt = 1;
for l=2:count1,
if(xq(i) > T(l-1) && xq(i) <= T(l))
xq(i) = cnt;
end
cnt = cnt +1 ;
end
end
D(w_count) = mean((x - centers(xq)).^2);
end
end
我打电话给它并有这些输入:
M = 10000
t=(randn(M,1)+sqrt(-1)*randn(M,1))./sqrt(2);
A= abs(t).^2;
[xq,centers,D] = Lloyd_Max( A,2,0,4 );
我试着评论while和D,结果: 我得到了xq和中心都正常,xq在1-4范围内,中心1-4个索引和0.5-3.5范围。
我不知道这里出了什么问题......请帮助我。
提前感谢!
神秘莫测!
谢谢大家的帮助! 我只是推出了for循环:
for i=1:10000,
if (x(i) > max_value)
xq(i) = max_value;
elseif (x(i) < min_value)
xq(i) = min_value;
else
xq(i) = x(i);
end
end
它就像魅力......这个循环再次启动了数组。对不起。再次感谢你!
答案 0 :(得分:1)
在您的函数中间某处有一个作业xq(i) = x(i)
,但是您从A
从x
计算A
的外部t
传递randn
由xq
采样,因此您不能保证{{1}}是一个整数。
答案 1 :(得分:0)
我不确定你的目标是什么,但你的向量xq
不包含整数,它包含双精度数。如果要像使用centers(xq)
一样使用索引向量,则向量的所有元素都必须是整数。
经过一点点检查,看起来xq是x值,你应该找到一些方法将它们映射到它们所属的最接近的单元格的整数(我猜'中心'代表单元格的中心?)< / p>