%filename:manwoman.m
%different man from woman.
%===========================================================
clear;
% [x,fs,bits]=wavread(‘d.wav‘, [1 5000]); % 读入声音文件(*.wav)
[x,fs,bits]=wavread(‘d.wav‘, [5000 10944]); % 读入声音文件(*.wav)
sound(x,fs,bits); % 数据通过声卡转换为声音
%===========================================================
% pause;
data1=x(:,1);
n=0:length(x)-1;
time=n/fs;
subplot(3,1,1) % 绘制2行1列的第1张子图
plot(time,data1) % 以时间为横轴，数据为纵轴作图
xlabel(‘Time (sec.)‘) % 标注横坐标
ylabel(‘Signal Level (Volts)‘) % 标注纵坐标
grid on % 添加网格
% pause;
% 对采集数据作滤波处理
blocksize =length(data1); % 计算窗函数长度
window = hanning(blocksize); % 计算汉宁窗函数(此函数为MATLAB自带)
data2=window.*data1; % 对数据先作加窗处理
% wp=[70,400]; %100Hz--400Hz
% wp=wp*2/fs;
wp=3000*2/fs;
[b,a]=ellip(4,0.1,20,wp); % 构造椭圆滤波器
data=filter(b,a,data2); % 求加窗处理后的数据经过滤波器的响应
subplot(3,1,2) % 绘制2行1列的第1张子图
plot(time,data) % 以时间为横轴，数据为纵轴作图
xlabel(‘Time (sec.)‘) % 标注横坐标
%ylabel(‘Signal Level (Volts)‘) % 标注纵坐标
grid on % 添加网格
% pause;
[xmax,index]=max(data1);
timewin=floor(0.015*fs);
xwin=data1(index-timewin:index+timewin);
[y,lags]=xcov(xwin);
subplot(3,1,3)
plot(lags,y)
grid on
ylen=length(y);
halflen=(ylen+1)/2 +30;
yy=y(halflen: ylen);
[ymax,maxindex] = max(yy);
fmax=fs/(maxindex+30);
disp([‘Maximum occurred at ‘, num2str(fmax), ‘ Hz‘])
f=fmax;
justify=f; %Hz
if justify<180 %limitHz以下判断是男的
display ‘男‘
else
display ‘女‘
end
原文：/zhang-pengcheng/p/4476533.html
