Question

我在MatLab中遇到一些问题。

我试图用给定的汉明权重r来计算长度为m的二进制字的数量，其中只包含多达k个连续的零。汉明权重是二进制字中非零项的数量。

我根据论文"Constant-Weight and Constant-Charge Binary Run-Length Limited Codes" (Kurmaev)实施了以下代码。

%% Different cases for the weight r
if (r==0)
    if (m<=k)
        numberOfBinaryStrings = 1;
    else
        numberOfBinaryStrings = 0;
    end
else

%% Computation of the number of binary strings
% Determination of the sum bound
bound = min(m,k);

d = k+1;
tmp = 0;

for j=0:bound
    for s=0:r-1
        if (m-j-1-s*d < r-1)
            bin1 = 0;
        else
            bin1 = nchoosek(m-j-1-s*d,r-1);
        end

        if (m-j-1-(k+1)-s*d < r-1)
            bin2 = 0;
        else
            bin2 = nchoosek(m-j-1-(k+1)-s*d,r-1);
        end

        tmp = tmp + (-1)^s * nchoosek(r-1,s) * ( bin1 - bin2 );
    end
end
numberOfBinaryStrings = tmp;

端

代码适用于给定的k和低字长以及汉明权重r。在某些参数，特别是大参数，我得到负面结果，这不应该是。我已经尝试用gammaln-function替换nchoosek函数以避免溢出。但在那里我也得到了负面结果。

你有什么想法，我能做什么？谢谢！

Answer 1

对于m足够小，你可以：

生成长度为m的所有二进制单词作为矩阵，每行包含一个单词。这可以通过conversion to binary轻松完成。
仅保留体重r的人。每个单词的权重可以通过summing沿着第二维的上述矩阵获得。
计算没有k+1个连续零的单词数。这些单词可以通过用k+1个向量观察被否定词的2D-convolution输出来检测。

代码：

m = 7;
r = 4;
k = 2;
h = dec2bin(0:2^m-1) - '0'; % step 1
h = h(sum(h,2)==r, :); % step 2
result = sum(all(conv2(1-h, ones(1,k+1), 'valid') < k+1, 2)); % step 3

Answer 2

使用http://www.mathworks.com/matlabcentral/fileexchange/22725-variable-precision-integer-arithmetic在Matlab中获取任意大小的整数。这应该消除任何可能的溢出问题。

如果这不能解决您的问题，那么您在某处会出现逻辑错误。在那种情况下，我建议制作另一种解决方案，并试图产生一个不同的最小例子。然后试着找出原因。

这是Python中的另一个解决方案，您可以将其用于比较。

#! /usr/bin/env python

stored_word_counts = {}

def word_counts(word_length, ones, max_zero_run):
    if 0 == ones:
        if max_zero_run < word_length:
            # Impossible.
            return 0
        else:
            # String of all zeros.
            return 1
    key = (word_length, ones, max_zero_run)
    if key not in stored_word_counts:
        # We will try all places we can put a 1.
        mid = (ones+1)/2 # Cut in half, round up.
        ones_before = mid - 1
        ones_after = ones - mid
        stored_word_counts[key] = sum(
                word_counts(pos, ones_before, max_zero_run) * word_counts(word_length - pos - 1, ones_after, max_zero_run)
                    for pos in xrange(ones_before, word_length - ones_after)
                )
    return stored_word_counts[key]

print(word_counts(50, 20, 5)) # Change this line as needed.

确定长度为m和汉明权重r的二进制字符串的数量，它可以包含连续的零到k

2 个答案: