#include?
#include?
#include?
#include

int?number=18;
int?numbers=520000;
int??kk=0;
long?counter=0;
long?N_sample;
short?*lpdata;

#define?mfcc_dims??13
#define?pi?3.1415926
#define?SEQU_length??1280

void?waveread（char?*filename）;
double?mfcc（double?dataR[]double?dataI[]double?mel[13]int?N）;

struct?MFCC_STRUCT
{
	int?frames;
	double?**data;
};

void?waveread（char?*filename）
{
???long?count;
???char?cbuf[5];
???FILE?*fp;
???if（（fp=fopen（filename“rb+“））==NULL）?
???{
	???exit（0）;
???}
???counter=0;

???fread（cbufsizeof（char）4fp）;
???cbuf[4]=‘\0‘;
???if（strcmp（cbuf“RIFF“）!=0）
	???exit（1）;
???fseek（fpsizeof（char）*4SEEK_CUR）;
???fread（cbufsizeof（char）4fp）;
???cbuf[4]=‘\0‘;
???if（strcmp（cbuf“WAVE“）!=0）?
	???exit（1）;
???fread（cbuf41fp）;
???cbuf[4]=‘\0‘;
???if（strcmp（cbuf“fmt?“）!=0）
	???exit（1）;
???fread（&count41fp）;
???fseek（fpcount+4SEEK_CUR）;
???
???fread（&count41fp）;
???N_sample=count/2;

??if（N_sample>numbers）
??{
????counter=numbers;
??}
???else?
???{
????counter=N_sample;
???}
????lpdata=（short?*）malloc（counter*sizeof（short））;
?
????fread（lpdata2counterfp）;
????fclose（fp）;
	
}

void?swap?（double?a?double?b）
{
???double?t;
???t?=?a;
???a?=?b;
???b?=?t;
}

void?bitrp?（double?xreal?[]?double?ximag?[]?int?n）
{

???int?i?j?a?b?p;

???for?（i?=?1?p?=?0;?i????{
?????p?++;
???}
?????for?（i?=?0;?i?	?{
???????a?=?i;
???????b?=?0;
?????????for?（j?=?0;?j?		?{
????????????b?=?（b?<????????????a?>>=?1;
		?}
???????????if?（?b?>?i）
		???{
??????????????swap?（xreal?[i]?xreal?[b]）;
??????????????swap?（ximag?[i]?ximag?[b]）;
		???}?
	?}
}


double?mfcc（double?dataR[]double?dataI[]double?mel[13]int?N）
{
??int?imkjtindex1index2;
??int?n=256;
??double?wreal[512]wimag[512]trealtimagurealuimagarg;
?
???int?Fi_upFi_down;
???double?delt_f?=?8000.0?/?N;
???double?MFCC_Y[20]?MFCC_k?=?0;
???int?delt_F_up?delt_F_down;
???int?F[22]?=?{0100200300400500600700800900100011491320151617412000229726393031348240004000};
???double?cos_par[20]?=?{0.0654498469?0.1963495408?0.3272492347?0.4581489286?0.5890486225
0.7199483164?0.8508480103?0.9817477042?1.1126473981?1.2435470920
1.3744467859?1.5053464798?1.6362461737?1.7671458676?1.8980455615
2.0289452554?2.1598449493?2.2907446432?2.4216443371?2.5525440310};
???
???mel[0]=0;
??for?（i=0;i??{
???mel[0]+=dataR[i]*dataR[i];
??}
???mel[0]=log（mel[0]）;

???for（i=0;i?????dataR[i]=dataR[i]*（0.54-0.46*cos（2*pi*i/（N-1）））;



??bitrp?（dataR?dataI?n）;

??arg?=?-?2?*?3.14159?/?n;
??treal?=?cos?（arg）;
??timag?=?sin?（arg）;
??wreal?[0]?=?1.0;
??wimag?[0]?=?0.0;
??for?（j?=?1;?j???{
???wreal[j]=（double）cos（（j-1）*arg）;?
???wimag[j]=（double）sin（（j-1）*arg）;?
??}
??for?（m?=?2;?m?<=?n;?m?*=?2）
??{
??for?（k

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

????I.A....?????31418??2006-10-22?09:57??dtw_train\0_10.wav

????I.A....?????25318??2006-10-22?09:56??dtw_train\0_11.wav

????I.A....?????26134??2006-10-22?09:56??dtw_train\0_12.wav

????I.A....?????17988??2006-10-22?10:02??dtw_train\0_20.wav

????I.A....?????16480??2006-10-22?10:02??dtw_train\0_21.wav

????I.A....?????15722??2006-10-22?10:02??dtw_train\0_22.wav

????I.A....?????18738??2006-10-22?10:07??dtw_train\0_30.wav

????I.A....?????18188??2006-10-22?10:06??dtw_train\0_31.wav

????I.A....?????25418??2006-10-22?10:06??dtw_train\0_32.wav

????I.A....?????18014??2006-10-22?16:22??dtw_train\0_40.wav

????I.A....?????16712??2006-10-22?16:22??dtw_train\0_41.wav

????I.A....?????17982??2006-10-22?16:23??dtw_train\0_42.wav

????I.A....?????17556??2006-10-22?20:57??dtw_train\0_50.wav

????I.A....?????16376??2006-10-22?20:57??dtw_train\0_51.wav

????I.A....?????17312??2006-10-22?20:58??dtw_train\0_52.wav

????I.A....?????17322??2006-10-22?22:00??dtw_train\0_60.wav

????I.A....?????17198??2006-10-22?22:01??dtw_train\0_61.wav

????I.A....?????16012??2006-10-22?22:01??dtw_train\0_62.wav

????I.A....?????29424??2006-10-22?10:12??dtw_train\1_10.wav

????I.A....?????27978??2006-10-22?10:11??dtw_train\1_11.wav

????I.A....?????29228??2006-10-22?10:11??dtw_train\1_12.wav

????I.A....?????13292??2006-10-22?10:14??dtw_train\1_20.wav

????I.A....?????12474??2006-10-22?10:14??dtw_train\1_21.wav

????I.A....?????11960??2006-10-22?10:14??dtw_train\1_22.wav

????I.A....?????20822??2006-10-22?10:17??dtw_train\1_30.wav

????I.A....?????16834??2006-10-22?10:17??dtw_train\1_31.wav

????I.A....?????17894??2006-10-22?10:17??dtw_train\1_32.wav

????I.A....?????10930??2006-10-22?16:17??dtw_train\1_40.wav

????I.A....?????14132??2006-10-22?16:18??dtw_train\1_41.wav

????I.A....?????13286??2006-10-22?16:18??dtw_train\1_42.wav

............此處省略182個文件信息