%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function?[yfsposition]?=?audioSeg（filenmclipLen）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%This?is?main?audio?stream?segment?functionit?can?classfication?the?audio?file?to?different?segments.
%
%input:
%?????filenm?=?the?flie?which?need?to?be?classification;
%?????clipLen?=?segment?clip?lengthunit?is?framedefault?value?is?100;
%output:
%?????y?=?the?points?vector?of?the?file;
%?????fs?=?sample?rate;
%?????position?=?type?change?points;
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?Written?by?Qian?Yong?Gao???????????????????%%
%%?Data:June?232006??????????????????????????%%
%%?Update:June?232006?by?Qian?Yong?Gao???????%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if?nargin?????clipLen?=?100;
end
??????????
[fidmsg]?=?fopen（filenm‘rb‘）;%open?file
error（msg）;
fseek（fid8-1）;%file?point?rewind
header?=?fread（fid4‘uchar‘）;

%?makesure?‘*wav‘?file?format
if?header‘?~=?‘WAVE‘
????error（‘This?file?is?not?a?wavfile‘）;
else
????fseek（fid80）;
????iPCM?=?fread（fid1‘ushort‘）;%PCM?code?format
????iMono?=?fread（fid1‘ushort‘）;%Mono?channel?
????if?（iPCM?~=?1）?&?（iMono?~=?1）
????????error（‘This?wavefile?is?not?Mono?PCM?type‘）;
????end
end
fs?=?fread（fid1‘ushort‘）;%sample?rate
fseek（fid80）;
nBit?=?fread（fid1‘ushort‘）;%each?sample?bits
nByte?=?nBit/8;%each?sample?bytes;
header?=?fread（fid4‘uchar‘）;
if?header‘?~=?‘data‘
????error（‘The?file?is?corrupt‘）;
end????
if?nByte?==?1
????nData?=?fread（fid1‘ulong‘）;?%file?samples?number
????rid?=?‘uchar‘;
else
????nData?=?fread（fid1‘ulong‘）?/?2;
????rid?=?‘short‘;
end????

position?=?0;%beginning?position?is?‘0‘
nRead?=?0;?%n?samples?read
nframe?=?0;?%n?frame?read
segLen?=?0;?%n?clip?read
y?=?[];%signal?read?
frameLen?=?floor（0.02?*?fs）;?%frame?length
frameMov?=?ceil（0.01?*?fs）;?%frame?step
overlap?=?frameLen?-?frameMov;%frame?overlap

while?nRead?<=?nData?-?frameLen
????y1?=?fread（fidframeLenrid）;%read?a?frame
????nRead?=?nRead?+?frameMov;%read?samples?number?increase?frameMov
????if?fid?==?‘uchar‘
????????y1?=?（y1?-?128）?/?128?;
????????fseek（fid-overlap0）;%8?bitrewind?file?pointer?overlap?bites
????else?
????????y1?=?y1?/?32768;
????????fseek（fid-2?*?overlap0）;%16?bitrewind?file?pointer?2*overlap?bites
????end????
????nframe?=?nframe?+?1;%increase?read?frame?number
????y?=?[y;y1‘];%read?signal?????
????%when?read?frame?is?clipLen‘s?integer?timesnClip?increase?1
????segLen?=?segLen?+?1;
????if?segLen?==?3*clipLen
????????clipMat?=?y（[nframe-3*clipLen+1:nframe]:）;%clip?matrix?which?wait?for?classifcation
????????f?=?fft（clipMat‘.*（hamming（frameLen）*ones（13*clipLen））1024）;
????????f?=?abs（f）;
????????mfc?=?melfc（f3030）;?
????????segN?=?BeyasInfoCrit（mfcframeLenframeMov0.6）;
????????if?（segN>=clipLen）?&?（segN<=2*clipLen）
????????????position?=?[position;nframe-3*clipLen+segN];
????????????segLen?=?3?*?clipLen?-?segN;
????????else
???????

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????3420??2009-12-08?16:47??音頻分割\audioSeg.m

?????文件???????1516??2009-12-08?16:27??音頻分割\BeyasInfoCrit.m

?????文件???????6706??2009-12-08?16:43??音頻分割\mfcc.m

?????目錄??????????0??2009-12-22?16:31??音頻分割

-----------?---------??----------?-----??----

????????????????11642????????????????????4