《Matlab语音信号分析与合成》第2版程序

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  2016-10-09 15:28  程序1\
     目录           0  2016-10-09 15:28  程序1\basic_tbx\
     文件        1050  2013-10-22 15:28  程序1\basic_tbx\add_noisedata.m
     文件        1010  2013-10-22 15:31  程序1\basic_tbx\add_noisefile.m
     文件        2895  2009-11-01 21:08  程序1\basic_tbx\enframe.m
     文件        4242  2010-10-28 10:41  程序1\basic_tbx\findpeaks.m
     文件         753  2013-06-11 20:17  程序1\basic_tbx\findSegment.m
     文件         110  2012-06-12 08:40  程序1\basic_tbx\frame2time.m
     文件         234  2004-05-26 14:14  程序1\basic_tbx\freqz_m.m
     文件        2535  2010-08-01 08:41  程序1\basic_tbx\frq2mel.m
     文件         645  2013-05-01 16:20  程序1\basic_tbx\Gnoisegen.m
     文件         136  2004-05-26 14:05  程序1\basic_tbx\ideal_lp.m
     文件        2931  1999-10-29 16:01  程序1\basic_tbx\instfreq.m
     文件         564  2013-10-23 14:54  程序1\basic_tbx\linsmoothm.m
     文件        1315  2007-05-04 08:01  程序1\basic_tbx\lpcar2ff.m
     文件        1441  2007-05-04 08:01  程序1\basic_tbx\lpcar2pf.m
     文件        1307  2007-05-04 08:01  程序1\basic_tbx\lpcar2zz.m
     文件        2525  2010-08-01 08:41  程序1\basic_tbx\mel2frq.m
     文件       11383  2010-01-02 20:58  程序1\basic_tbx\melbankm.m
     文件        1585  2013-10-23 10:52  程序1\basic_tbx\Mtmpsd_ssb.m
     文件        1663  2008-01-13 21:30  程序1\basic_tbx\OverlapAdd2.m
     文件        1352  2013-10-23 14:58  程序1\basic_tbx\pitch_vad1.m
     文件        2358  2013-08-01 21:28  程序1\basic_tbx\pitch_vads.m
     文件         398  2013-01-02 17:33  程序1\basic_tbx\pitfilterm1.m
     文件         363  2013-06-08 10:35  程序1\basic_tbx\polydetrend.m
     文件        1888  2009-06-03 11:57  程序1\basic_tbx\rfft.m
     文件        1182  2013-10-23 10:16  程序1\basic_tbx\simplesubspec.m
     文件         400  2013-05-01 06:07  程序1\basic_tbx\SNR_singlech.m
     文件        1650  2008-11-05 15:59  程序1\basic_tbx\SpecColorMap.m
     文件         655  2013-06-08 05:38  程序1\basic_tbx\stftms.m
     目录           0  2016-10-09 15:28  程序1\EMD\
............此处省略348个文件信息