数字变声器的lpc函数问题！！！

Question

[s,fs] =audioread('88.mp3');

s = s/max(s);

L = length(s);

FN = floor(L/FL)-2;

for n = 3:FN

     s_w = s(n*FL-WL+1:n*FL).*hw;
    [A, E] = lpc(s_w, P);
    代码到这里运行时，A这个1X11行矩阵会出现后面十个数没有定义的情况！！！！！
    原因好像是音频文件前部分没有声音，所以s前部分取值为零，s_w部分值也为零，这个怎么解决啊！
    如果自己定义A矩阵的话音频输出就会有问题！！！

Answer 1

以下回答参考 皆我百晓生、券券喵儿 等免费微信小程序相关内容作答，并由本人整理回复。

这个问题是由于音频文件中前几帧的采样值为空导致的。在进行低通滤波之前，我们需要确保音频数据中的所有帧都有至少一个非零值。

你可以通过以下方式解决这个问题：

1. 将音频文件分割成小块，然后对每个小块应用低通滤波。
2. 对每个小块应用高斯窗函数进行加权平均，以消除其前几帧的噪声影响。
3. 将处理过的每一帧合并到一个大数组中。

以下是一个示例代码，使用MATLAB和numpy库来实现上述步骤：

% 加载音频文件
[s, fs] = audioread('88.mp3');
s = s / max(abs(s)); % 预加重

% 获取音频长度和采样率
L = length(s);
Fs = fs; % 采样率

% 定义低通滤波参数
P = 4; % 高通截止频率（Hz）
WL = 500; % 滤波器宽度（Hz）

% 创建一个长度等于音频长度的空矩阵
A = zeros(1, L); 

% 进行低通滤波
for n = 1:L-1
    s_w = s(n*WL+1 : (n+1)*WL); % 取出当前帧的所有样本点
    A(n) = sum(exp(-abs(s_w)/2 * (P^2 - WL^2))); % 计算加权平均
end

% 输出结果
disp(A)

请注意，这段代码只是一个基本示例，实际应用可能需要根据具体情况进行调整。例如，你可能需要修改高通截止频率或滤波器宽度等参数，或者将音频文件分割成更小的块。