MATLAB自定义量化函数为何仅返回一个输出？

核心问题排查与解决 #

1. 函数调用方式错误（最直接原因） #

如果调用函数时只写了一个接收变量，比如：

Xq = cuantificacion(x, xmax, xmin, b);

MATLAB只会返回第一个输出参数Xq。必须用两个变量接收所有返回值：

[Xq, SNq] = cuantificacion(x, xmax, xmin, b);

2. 代码内部逻辑错误（导致SNq计算异常） #

你的代码里还有几个会影响SNq正常输出的问题：

• 循环逻辑错误：MATLAB是1-based索引，循环变量i应从1开始；且直接对整个向量x操作会重复累加，需遍历每个元素。另外，量化噪声功率应该计算误差的平方和，而非误差和：

  for i = 1:N
      p = p + x(i)^2;
      q = q + (Xq(i) - x(i))^2;
  end
  

• 对数函数误用：MATLAB中log是自然对数，计算信噪比需用常用对数log10，否则结果数值会偏差：

  SNq = 10*log10(Px/Pq);
  

• 交互式输入不合理：函数内部用input获取样本数N，会导致调用时需手动输入，且易与输入向量x的长度不匹配。建议直接从x的长度获取：

  N = length(x);
  

• 量化边界处理疏漏：abs(x)<xmax的判断应包含等于的情况，否则等于xmax的元素会被错误归类：

  if abs(x(i)) <= xmax
  

修正后的完整代码 #

function [Xq, SNq] = cuantificacion(x, xmax, xmin, b)
    N = length(x); % 直接从输入向量获取样本数
    L = 2^b;
    delta = (xmax - xmin)/L;
    
    Xq = zeros(size(x));
    % 逐元素处理量化逻辑
    for i = 1:N
        if abs(x(i)) <= xmax
            Xq(i) = (fix(abs(x(i))/delta) + 1/2)*delta*sign(x(i));
        else
            Xq(i) = ((L-1)/2)*delta*sign(x(i));
        end
    end
    
    % 计算信号功率与量化噪声功率
    p = 0;
    q = 0;
    for i = 1:N
        p = p + x(i)^2;
        q = q + (Xq(i) - x(i))^2;
    end
    Px = p/N;
    Pq = q/N;
    SNq = 10*log10(Px/Pq);
end

内容的提问来源于stack exchange，提问作者user18559097