matlab语音信号基频参数提取

%aa=wavread('tezhengtiqu1.wav')
[x,fs,bits]=wavread('抢地主0.wav'); % 读入声音文件（*.wav）
data1=x(:,1);
n=0:length(x)-1;
time=n/fs;
% 对采集数据作滤波处理
blocksize =length(data1); % 计算窗函数长度
window = hanning(blocksize); % 计算汉宁窗函数(此函数为MATLAB自带)
data2=window.*data1; % 对数据先作加窗处理
wp=3000*2/fs;
[b,a]=ellip(4,0.1,20,wp); % 构造椭圆滤波器
data=filter(b,a,data2);
%plot(data);% 求加窗处理后的数据经过滤波器的响应
[xmax,index]=max(data1);
timewin=floor(0.015*fs);
xwin=data1(index-timewin:index+timewin);
[y,lags]=xcov(xwin);
ylen=length(y);
halflen=(ylen+1)/2 +30;
yy=y(halflen: ylen);
%画基频图
plot(yy);
[ymax,maxindex] = max(yy);
fmax=fs/(maxindex+30);
%fmaxce2;
fmaxce2=fmax;
%fmax=roundn(fmax,-3);
fmax21=roundn((fmax-5),3);
fmax22=roundn((fmax+5),3);
%fmax2= [fmax21,'~~',fmax22];

projfile=fullfile(matlabroot,'toolbox','yuyin101','PID.pjt')
projpath=fileparts(projfile)
open(cc,projfile)
visible(cc,1)
cd(cc,projpath)
build(cc,'all',60)
load(cc,'Timer.out',30)
open(cc,'main.c','text')
activate(cc,'main.c','text')
insert(cc,'main.c',61)
halt(cc)
restart(cc)
run(cc,'runtohalt',50)
ddatV=read(cc,address(cc,'JP1'),'int16',1)
JP1=int16(fmaxce2);
ddatV=JP1;
%idatV=read(cc,address(cc,'fOut'),'single',256)
%write(cc,address(cc,'JP1'),fmaxce2));
run(cc,'runtohalt')
%ddatV=read(cc,address(cc,'fIn'),'single',256)
%idattV=read(cc,address(cc,'fOut'),'single',1000)
%plot(abs(idattV))