基音周期预测.pdfVIP

数字信号处理实验 专业班级 组别 成员 实验内容： 编程求解出各自声音信号的基音周期。 1 程序代码（或者软件流程图等） (1) function zhouqi=jiyinzhouqi(filename) %帧长和帧位移是重要的参数，位移是帧长的 0~1/2 %短时自相关分析 %filename 语音文件 *.wav %zhouqi 基音周期 , 以毫秒为单位表示 [signal,fs]=wavread(filename); % 用于得到声音文件的数据和采样率 shift=0.02; % 每次移动 20 毫秒 shift=round(fs*shift); % 帧移 n1=fix(fs*0.01)+1; % 分析起点 0.01ms, 帧长 20ms n2=fix(fs*0.03)+1; shift_count=fix((length(signal)-n1)/shift); value =zeros(1,shift_count); % 存放每次移位后的帧的基音周期 for ii=1:shift_count % 分析次数 if n2length(signal) % 防止溢出 data=signal(n1:n2); % 加窗，提取一帧数据 N=n2-n1+1; % 每一帧的长度 R=zeros(1,N); % 创建一个一行 N 列的矩阵 for k=1:N-1 % 求自相关序列 for jj=1:N-k R(k)=R(k)+data(jj)*data(jj+k); % 矩阵用于储存每次自相关的结果； end end value(ii)=find_maxn(R); % 调用基音周期分析函数，求最大值所 对应的位置，即基音周期 n1=n1+shift; % 移动帧，计算下一帧的基音周期 第 - 1 - 页 数字信号处理实验 n2=n2+shift; end end figure(1) stem(value); % 画出基因周期走势图 axis([0 length(value) 0 1000]) aver=mean(value); % 基音周期的平均值 , 未去除野点 value=value(logical(abs(value-aver)=aver/5));% 找出偏移均值超出均值的 1/5 的基音周期，将其去除 len= length(value); % 去除大野点后剩余的基音点数 for jj=1:3:len/3 % 中值平滑，滑动窗口宽度 3，精度为 中值 1/4 （剔除野点） average=(value(jj)+value(jj+1)+value(jj+2))/3; for kk=1:3 if abs((value(jj-1+kk))-average)average/4 value(jj-1+kk)=0; % 将野点置零，同时数组长度减一 end end end value=value(( value~=0)); % 出去所有野点后的基音周期数组 len= length(value); % 去除野点以后的基音点数 fig