新建文件夹（3）（1）.zip

新建文件夹(3)(1).zip

% 
clear all; clc; close all;

fle=‘snn27.wav‘;                            % 指定文件名
[xxfs]=wavread（fle）;                       % 读入一帧语音信号
u=filter（[1 -.99]1xx）;                    % 预加重
wlen=length（u）;                             % 帧长
p=12;                                       % LPC阶数
a=lpc（up）;                                 % 求出LPC系数
U=lpcar2pf（a255）;                          % 由LPC系数求出功率谱曲线
freq=（0:256）*fs/512;                        % 频率刻度
df=fs/512;                                  % 频率分辨率
U_log=10*log10（U）;                          % 功率谱分贝值
subplot 211; plot（u‘k‘）;                   % 作图
axis（[0 wlen -0.5 0.5]）;
title（‘预加重波形‘）;
xlabel（‘样点数‘）; ylabel（‘幅值‘）
subplot 212; plot（freqU_log‘k‘）;
title（‘声道传递函数功率谱曲线‘）;
xlabel（‘频率/Hz‘）; ylabel（‘幅值/dB‘）;

n_frmnt=4;                                  % 取四个共振峰
const=fs/（2*pi）;                            % 常数  
rts=roots（a）;                               % 求根
k=1;                                        % 初始化
yf = [];
bandw=[];
for i=1:length（a）-1                     
    re=real（rts（i））;                        % 取根之实部
    im=imag（rts（i））;                        % 取根之虚部
    formn=const*atan2（imre）;               % 计算共振峰频率
    bw=-2*const*log（abs（rts（i）））;           % 计算带宽
    
    if formn>150 & bw <700 & formn        yf（k）=formn;
        bandw（k）=bw;
        k=k+1;
    end
end

[y ind]=sort（yf）;                          % 排序
bw=bandw（ind）;
F = [NaN NaN NaN NaN];                      % 初始化
Bw = [NaN NaN NaN NaN];
F（1:min（n_frmntlength（y））） = y（1:min（n_frmntlength（y）））;   % 输出最多四个
Bw（1:min（n_frmntlength（y））） = bw（1:min（n_frmntlength（y）））; % 输出最多四个
F0 = F（:）;                                  % 按列输出
Bw = Bw（:）;
p1=length（F0）;                              % 在共振峰处画线
for k=1 : p1
    m=floor（F0（k）/df）;
    P（k）=U_log（m+1）;
    line（[F0（k） F0（k）][-10 P（k）]‘color‘‘k‘‘linestyle‘‘-.‘）;
end
fprintf（‘F0=%5.2f   %5.2f   %5.2f   %5.2f\n‘F0）;
fprintf（‘Bw=%5.2f   %5.2f   %5.2f   %5.2f\n‘Bw）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-10-10 15:39  源码\
     目录           0  2016-04-15 20:40  源码\9550\
     目录           0  2016-04-15 19:53  源码\9550\pro\
     目录           0  2016-04-15 20:17  源码\9550\pro\basic_tbx\
     文件        1050  2013-10-23 04:28  源码\9550\pro\basic_tbx\add_noisedata.m
     文件        1010  2013-10-23 04:31  源码\9550\pro\basic_tbx\add_noisefile.m
     文件        2895  2009-11-02 10:08  源码\9550\pro\basic_tbx\enframe.m
     文件        4242  2010-10-28 23:41  源码\9550\pro\basic_tbx\findpeaks.m
     文件         753  2013-06-12 09:17  源码\9550\pro\basic_tbx\findSegment.m
     文件         110  2012-06-12 21:40  源码\9550\pro\basic_tbx\frame2time.m
     文件         234  2004-05-27 03:14  源码\9550\pro\basic_tbx\freqz_m.m
     文件        2535  2010-08-01 21:41  源码\9550\pro\basic_tbx\frq2mel.m
     文件         645  2013-05-02 05:20  源码\9550\pro\basic_tbx\Gnoisegen.m
     文件         136  2004-05-27 03:05  源码\9550\pro\basic_tbx\ideal_lp.m
     文件        2931  1999-10-30 05:01  源码\9550\pro\basic_tbx\instfreq.m
     文件         564  2013-10-24 03:54  源码\9550\pro\basic_tbx\linsmoothm.m
     文件        1315  2007-05-04 21:01  源码\9550\pro\basic_tbx\lpcar2ff.m
     文件        1441  2007-05-04 21:01  源码\9550\pro\basic_tbx\lpcar2pf.m
     文件        1307  2007-05-04 21:01  源码\9550\pro\basic_tbx\lpcar2zz.m
     文件        2525  2010-08-01 21:41  源码\9550\pro\basic_tbx\mel2frq.m
     文件       11383  2010-01-03 09:58  源码\9550\pro\basic_tbx\melbankm.m
     文件        1585  2013-10-23 23:52  源码\9550\pro\basic_tbx\Mtmpsd_ssb.m
     文件        1663  2008-01-14 10:30  源码\9550\pro\basic_tbx\OverlapAdd2.m
     文件        1352  2013-10-24 03:58  源码\9550\pro\basic_tbx\pitch_vad1.m
     文件        2358  2013-08-02 10:28  源码\9550\pro\basic_tbx\pitch_vads.m
     文件         398  2013-01-03 06:33  源码\9550\pro\basic_tbx\pitfilterm1.m
     文件         363  2013-06-08 23:35  源码\9550\pro\basic_tbx\polydetrend.m
     文件        1888  2009-06-04 00:57  源码\9550\pro\basic_tbx\rfft.m
     文件        2206  2016-04-15 19:43  源码\9550\pro\LPC.m
     文件         684  2012-11-23 00:30  源码\9550\pro\snn27.wav
     文件      268288  2018-10-10 18:36  源码\9550\基于LPC法的共振峰检测20160415.doc
............此处省略6个文件信息