function?ds（filer）

[yfs]=audioread（file）;
y=y（:1）;
yy=downsample（yr）;?????%降采樣
t=interp（yyr）;?????????%插值

n=0:length（t）-1;????????%建立一個信號等長的序列
time=n/fs;??????????????%建立時間序列，作為橫坐標
subplot（211）;?????????%圖1：時域波形圖
plot（timet）;???????????%作圖
title（‘音頻信號波形圖‘）??%標題
xlabel（‘時間/s‘）;???????%標注橫坐標
ylabel（‘幅值‘）;?????????%標注縱坐標
grid?on;????????????????%打開網格線
%=======頻域圖像======
N=length（t）;????????????%取信號矩陣的長度
Y1=fft（tN）;????????????%N點傅里葉變換
mag=abs（Y1）;????????????%取模
f=n*fs/N;???????????????%頻率序列
subplot（212）;?????????%圖2：頻譜圖
???????????%取1/2作圖
plot（f（1:fix（N/2））mag（1:fix（N/2）））;
title（‘音頻信號fft頻譜圖‘）;%標題
xlabel（‘頻率/Hz‘）;???????%標注橫坐標
ylabel（‘幅度‘）;??????????%標注縱坐標
grid?on;?????????????????%打開網格線
%======基因頻率提取======
[~index]=max（t）;
timewin=floor（0.015*fs）;
xwin=t（index-timewin:index+timewin）;
[y~]=xcov（xwin）;
ylen=length（y）;
halflen=（ylen+1）/2?+30;
yy=y（halflen:?ylen）;
[~maxindex]?=?max（yy）;
fmax=fs/（maxindex+30）;
disp（[file‘基音頻率為?‘?num2str（fmax）?‘?Hz‘]）
%======通過基因頻率判斷男女聲======
if?fmax<250;
????disp（[file‘?是男聲文件‘]）;
else
????disp（[file‘?是女聲文件‘]）;
end;
sound（tfs）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1296??2015-05-24?13:59??ds.m
?????文件????????1294??2015-05-19?20:59??man_or_woman.m