%

% pr10_3_1

clear all; clc; close all;

filedir=[];                             % 设置路径

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

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

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

x=x-mean(x);                            % 消除直流分量

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

xl=length(x);                           % 数据长度

time=(0:xl-1)/fs;                       % 计算出时间刻度

p=12;                                   % LPC的阶数为12

wlen=200; inc=80;                       % 帧长和帧移

msoverlap = wlen - inc;                 % 每帧重叠部分的长度

y=enframe(x,wlen,inc)';                 % 分帧

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

% 语音分析:求每一帧的LPC系数和预测误差

for i=1 : fn

u=y(:,i);                           % 取来一帧

A=lpc(u,p);                         % LPC求得系数

aCoeff(:,i)=A;                      % 存放在aCoeff数组中

errSig = filter(A,1,u);             % 计算预测误差序列

resid(:,i) = errSig;                % 存放在resid数组中

end

%%%%%%%%%%%%%%%%%%%%%

%这个地方添加了程序

%对预测系数和预测误差进行上采样

aCoeff=upsample(aCoeff,2);

resid=upsample(resid,2);

%当采样率是8k的时，将参数进行修改，以扩展到16k

fn=size(aCoeff,2);

inc=inc*2;

wlen=wlen*2;

xl=xl*2;

fs=fs*2;

%%%%%%%%%%%%%%%%%%%%%

% 语音合成:求每一帧的合成语音叠接成连续语音信号

for i=1:fn

A = aCoeff(:,i);                    % 取得该帧的预测系数

residFrame = resid(:,i);            % 取得该帧的预测误差

synFrame = filter(1, A', residFrame); % 预测误差激励,合成语音

outspeech((i-1)*inc+1:i*inc)=synFrame(1:inc);  % 重叠存储法存放数据

% 如果是最后一帧,把inc后的数据补上

if i==fn

outspeech(fn*inc+1:(fn-1)*inc+wlen)=synFrame(inc+1:wlen);

end

end;

ol=length(outspeech);

if ol

outspeech=[outspeech zeros(1,xl-ol)];

end

% 发声

% wavplay(x,fs);

% pause(1)

wavplay(outspeech,fs);

wavwrite(outspeech,fs,'bandwidth_audio.wav');

% 作图

% subplot 211; plot(time,x,'k');

% xlabel(['时间/s' 10 '(a)']); ylabel('幅值'); ylim([-1 1.1]);

% title('原始语音信号')

% subplot 212; plot(time,outspeech,'k');

% xlabel(['时间/s' 10 '(b)']); ylabel('幅值'); ylim([-1 1.1]);

% title('合成的语音信号')