多媒体语音信号处理LTP长时预测器

前人栽树后人乘凉：记得关注！
感谢室友！感谢同学分享！
题目：

代码：

clc;
clear all;

%读取音频文件
[s,fs] = audioread(‘F:\JI course\第二学期\ECTE801Multimedia Signal Processing\ECTE801多媒体信号处理\ECTE801多媒体信号处理\Speech+Files-20200213\Male16k-2017.wav’);

fs = 16000; % 设定采样频率
tf = 0.02; % 设定合适的帧长度 0.1 - 0.3
N = fs * tf; % 每一帧信号的长度

%求LPC Frame 的个数
LPC_num_frames = ceil(length(s) / N);

%防止无法操作补0
pad = LPC_num_frames * N - length(s);
s = [s;zeros(pad,1)];

buffer = reshape(s, N, LPC_num_frames);
p = 16;
LPC_Fres = []; % LPC分析滤波器线性预测的余量信号

%模块一：LP分析滤波器。
for i = 1 : LPC_num_frames
[a, g] = lpc(buffer(:, i), p);
z = zeros(p, 1);
[LPC_res, Z] = filter(a, 1, buffer(:, i), z);