我在实现一篇论文里面的算法

有几个地方不太明白，想请教请教

1、使用覆盖窗口方法计算短时能量函数。    是不是就是书上的计算短时能量的那种方法呀？

2、一个序列加窗过后的因果序列部分具有最小相位特点。      这里的因果序列部分怎么取呀，我是直接加窗过后就接着下一步运算了

3、将语音信号进行量化,取倒数过程变换成映射过程 ,即若原来为最大值的将其变成最小值, 原来为次最大值则变成次最小值。     这个是怎么实现呀？

4、我再做语音信号处理时，想对一个序列(x)加窗

要求是加窗的长度是序列长度(N)的1/4到1/8

我如果直接  y=x.*hamming(N/4)'      *这里的N/4是整数了

会显示矩阵维度不一样

于是，我就想着先让窗的长度为N/4，然后在窗的后面补0

即win1=hamming(N/4)';

a=zeros(1,N-N/4);

win=[win1,a];

y=x.*win; 不知道这样对不对

我现在写完代码结果是不正确的，觉得可能是这些地方出错了，求指导呀

这是我目前的代码

filedir=[];                     %设置路径

filename='SP_OYCY_2_02.wav';        %设置文件名

fle=[filedir filename];         %构成完整的路径和文件名

[x,fs]=audioread(fle);          %读入数据文件

x = x / max(abs(x));            %幅值归一化

NN=length(x);                    %求数据长度

time=(0:NN-1)/fs;                % 设置时间

wlen=1100;inc=440;                %确定帧长和帧移

IS=0.1;NIS=fix((IS*fs-wlen)/inc+1); %计算NIS

overlap=wlen-inc;

win1=boxcar(wlen);               %给出窗

y=enframe(x,win1,inc)';         %信号分帧

fn=size(y,2);                   %确定帧数

frameTime=frame2time(fn, wlen, inc, fs);% 计算各帧对应的时间

amp1=sum(y.^2);                 %求短时平均能量

ampm1=multimidfilter(amp1,10);   %中值平滑

for i=1:fn;

amp2(i)=ampm1(fn+1-i);

end

amp3=[amp2,ampm1];                 %合并

E=1./(amp3).^0.001;                 %求解函数

N1=length(E);                      %确定离散时间傅里叶变换长度

cn1=ifft(E,N1);

win2=hamming(N1/6)';

a=zeros(1,N1-N1/6);

win=[win2,a];

cn=cn1.*win;              %加窗的因果序列

N=N1;n=0:N-1;y=cn;          %绘制相位谱

Y=fft(y,N);

A=abs(Y);f=n*fs/N;

ph=2*angle(Y(1:N/2));

ph=ph*180/pi;

angX=ph(1:N/2);

subplot 311; plot(time,x,'k');                  %画出语音波形图

title('纯语音波形');ylabel('幅值');

axis([0 2.5 -1 1]);

subplot 312;plot(frameTime,ampm1,'k');          %画出短时能量图

xlabel('时间/s');ylabel('幅值');title('短时能量');

subplot 313;plot(frameTime,ph(1:N/2));           %画出相位谱图

xlabel('频率/hz'),ylabel('相角'),title('相位谱');

grid on;

for i=2:N1/2;                           %求最小相位延迟函数

ak(i)=angX(i)-angX(i-1);

end

figure;plot(frameTime,ak);title('最小群延迟函数');

QQ截图.png

(28.5 KB, 下载次数: 0)

2019-3-10 19:18 上传

我得出来的最小相位群延迟函数图，本来峰峰值处应该对应音节端点的