小波变换提取基音频率 Matlab

%
% pr3_4_1 
clear all; clc; close all;

[xfs]=wavread（‘awav.wav‘）;             % 读入语音数据
N=length（x）;                            % 信号长度
x=x-mean（x）;                            % 消除直流分量
J=2;                                    % 设小波变换级数为J
[CL] = wavedec（xJ‘db1‘）;             % 对时间序列进行一维多分辨分解
CaLen=N/2.^J;                           % 估计近似部分的系数长度
Ca=C（1:CaLen）;                          % 取近似部分的系数
Ca=（Ca-min（Ca））./（max（Ca）-min（Ca））;     % 对近似部分系数做规正处理
for i=1:CaLe