MATLAB在语音信号分析与合成中的应用 随书代码

《MATLAB在语音信号分析与合成中的应用》的随书代码，

function [signalnoise]=add_noisedata（sdatafsfs1snr）
s=s（:）;                          % 把信号转换成列数据
s=s-mean（s）;                     % 消除直流分量
sL=length（s）;                    % 求出信号的长度

if fs~=fs1                       % 若纯语音信号的采样频率与噪声的采样频率不相等
    x=resample（datafsfs1）;     % 对噪声重采样，使噪声采样频率与纯语音信号的采样频率相同
else
    x=data;
end

x=x（:）;                          % 把噪声数据转换成列数据
x=x-mean（x）;                     % 消除直流分量
xL=length（x）;                    % 求噪声数据长度
if xL>=sL                        % 如果噪声数据长度与信号数据长度不等，把噪声数据截断或补足
    x=x（1:sL）;
else
    disp（‘Warning: 噪声数据短于信号数据，以补0来补足！‘）
    x=[x; zeros（sL-xL1）];
end

Sr=snr;
Es=sum（x.*x）;                    % 求出信号的能量
Ev=sum（s.*s）;                    % 求出噪声的能量
a=sqrt（Ev/Es/（10^（Sr/10）））;      % 计算出噪声的比例因子
noise=a*x;                       % 调整噪声的幅值
signal=s+noise;                  % 构成带噪语音

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-09-30 20:18  20131118mlyy\
     文件        6238  2018-09-30 15:39  20131118mlyy\content.txt
     文件       24449  2018-09-30 15:39  20131118mlyy\cover.jpg
     文件         997  2018-09-30 15:39  20131118mlyy\ILMBookController.p
     文件         100  2018-09-30 15:39  20131118mlyy\ILMBookModel.p
     文件        2148  2018-09-30 15:39  20131118mlyy\ILMBookView.p
     文件         510  2018-09-30 15:39  20131118mlyy\main.p
     目录           0  2018-09-30 20:18  20131118mlyy\program\
     目录           0  2018-09-30 20:18  20131118mlyy\program\basic_tbx\
     文件        1050  2018-09-30 15:39  20131118mlyy\program\basic_tbx\add_noisedata.m
     文件        1010  2018-09-30 15:39  20131118mlyy\program\basic_tbx\add_noisefile.m
     文件        2895  2018-09-30 15:39  20131118mlyy\program\basic_tbx\enframe.m
     文件        4242  2018-09-30 15:39  20131118mlyy\program\basic_tbx\findpeaks.m
     文件         753  2018-09-30 15:39  20131118mlyy\program\basic_tbx\findSegment.m
     文件         110  2018-09-30 15:39  20131118mlyy\program\basic_tbx\frame2time.m
     文件         234  2018-09-30 15:39  20131118mlyy\program\basic_tbx\freqz_m.m
     文件        2535  2018-09-30 15:39  20131118mlyy\program\basic_tbx\frq2mel.m
     文件         645  2018-09-30 15:39  20131118mlyy\program\basic_tbx\Gnoisegen.m
     文件         136  2018-09-30 15:39  20131118mlyy\program\basic_tbx\ideal_lp.m
     文件        2931  2018-09-30 15:39  20131118mlyy\program\basic_tbx\instfreq.m
     文件         564  2018-09-30 15:39  20131118mlyy\program\basic_tbx\linsmoothm.m
     文件        1315  2018-09-30 15:39  20131118mlyy\program\basic_tbx\lpcar2ff.m
     文件        1441  2018-09-30 15:39  20131118mlyy\program\basic_tbx\lpcar2pf.m
     文件        1307  2018-09-30 15:39  20131118mlyy\program\basic_tbx\lpcar2zz.m
     文件        2525  2018-09-30 15:39  20131118mlyy\program\basic_tbx\mel2frq.m
     文件       11383  2018-09-30 15:39  20131118mlyy\program\basic_tbx\melbankm.m
     文件        1585  2018-09-30 15:39  20131118mlyy\program\basic_tbx\Mtmpsd_ssb.m
     文件        1663  2018-09-30 15:39  20131118mlyy\program\basic_tbx\OverlapAdd2.m
     文件        1352  2018-09-30 15:39  20131118mlyy\program\basic_tbx\pitch_vad1.m
     文件        2358  2018-09-30 15:39  20131118mlyy\program\basic_tbx\pitch_vads.m
     文件         398  2018-09-30 15:39  20131118mlyy\program\basic_tbx\pitfilterm1.m
............此处省略354个文件信息