Matlab错误:下标索引必须是实数正整数或逻辑

时间:2013-01-15 06:20:34

标签: matlab

我在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

它就像魅力......这个循环再次启动了数组。对不起。再次感谢你!

2 个答案:

答案 0 :(得分:1)

在您的函数中间某处有一个作业xq(i) = x(i),但是您从Ax计算A的外部t传递randnxq采样,因此您不能保证{{1}}是一个整数。

答案 1 :(得分:0)

我不确定你的目标是什么,但你的向量xq不包含整数,它包含双精度数。如果要像使用centers(xq)一样使用索引向量,则向量的所有元素都必须是整数。

经过一点点检查,看起来xq是x值,你应该找到一些方法将它们映射到它们所属的最接近的单元格的整数(我猜'中心'代表单元格的中心?)< / p>