语音信号线性预测分析(9)

clear;close all clc

fid=fopen('sx86.txt','r'); p1=fscanf(fid,'%f') fclose(fid);

p=filter([1-0.68],1,p1);%预加重滤波 x1=fra(320,160,p); for i=60:79 x=x1(i,:) a3=lpc(x,15);

a4=a3(:);%将线性预测系数付给矩阵a4

lsf=a_lsf_conversion(a4)%调用函数a_lsf_conversion实现从LPC系数到lsf参数的转换

isf=lsf;

isf(16)=0.5*acos(a4(16,1));%isp的最后一个参数取为a的最后一个参数，isf的最后一个参数取为0.5*acos(a4(16,1)）

isp=cos(isf);

hold on%让连续20帧isp及isf绘制在一个图形中 figure(1); for j=1:16

isf2(i-59,j)=isf(j); end figure(2); plot(isp) end

%EOF a_isf_lpc_conversion_20frame.m 运行结果：

连续20帧语音信号的ISP轨迹图