小波变换法求取一帧语音信号的基音周期检测

%小波法
%读音频文件----这部分运行有问题，运行显示：Index exceeds matrix dimensions（超出矩阵维度？）
clear all;
close all;
clc;
[filenamepathname]=uigetfile（‘*.wav‘‘请选择语音文件:‘）;
[soundfsnbits]=wavread（[pathnamefilename]）;%sound 为采样数据，fs采样率。nbits量化比特数
%假设声音是双声道，这里只取单声道作分析，sound（:1）为左声道，右声道为sound（:2） 
a=size（sound（:1））;
Tn=（a（1）-1）/fs;     
t=0:1/fs:Tn;
figure（1）
subplot（311）;
plot（tsound（:1））xlabel（‘t（s）‘）ylabel（‘sound‘）
title（‘录制的语音段‘）; 
I=sound（:1）;
s=I（11401:11700）;%可以根据需要选择不同的语音段
subplot（312）;
  plot（s）;
 title（‘选取的语音段‘）; 
 %小波去噪
[cl]=wavedec（s6‘db6