倒谱推导分析

根据对语音产生机理的研究可知,语音信号X(n)可看作声门激励信号和声道冲激响应信号的卷积
在这里插入图片描述

为了便于处理各卷积信号,系统通常采用同态处理的方法进行解卷积,即将卷积关系变为求和关系,分离参与卷积的各个信号。

一般的同态系统可分解为三个部分。

第一个子系统完成卷积信号转化为加性信号的运算
倒谱推导分析_第1张图片
第二个子系统进行线性处理得到

第三个子系统恢复卷积新号
倒谱推导分析_第2张图片
复倒谱:
在这里插入图片描述
倒谱:
在这里插入图片描述

clc;clear;
clc;clear;
[x,fs]=audioread('C3_4_y_1.wav');%读取信号
N=1024;%fft点数
x=x(1:1000);%取其中一段信号
len=length(x);
t=(0:len-1)/fs;%计算时间
y=fft(x,N);%根据公式计算倒谱
yy=real(ifft(log(abs(y))));
yy=yy(1:1000)';

subplot(311)
plot(t,x);
xlabel('时间/s');ylabel('幅值');
axis([0 max(t) -1 1]);
title('原信号');
subplot(312)
plot(t,yy)
xlabel('倒谱率/s');ylabel('幅值');
axis([0 t(512) -0.2 0.2]);
title('倒谱信号');

yc=cceps(x);
yn=icceps(yc);
subplot(313)
plot(t,yn);
xlabel('时间/s');ylabel('幅值');
axis([0 max(t) -1 1]);
title('还原信号');

倒谱推导分析_第3张图片

你可能感兴趣的:(matlab,语音)