Modélise votre fichier wav stéréo en abadie.jo

Pour soutenir l'équipe des développeurs de la carte audio linux
et faire avancer le proget du couplage à un casque neuronal,
don possible à partir de 1 euro



/* Imaginer et creer par Abadie joris*/
/* en utilisant un concept de Andre pierre jocelyn*/
/* abadie.jo en binaire 4 octets, 2 octets pour l'amplitude, un octets pour nbdepoints et un octet pour jo le nombre sans dimesion caractéristique de la forme*/
/*Cette oeuvre est protégée et n'est divulguée que pour utilisation sous linux dans les conditions Creative common ci-dessous*/
/* Cette oeuvre est mise à disposition selon les termes de la Licence Creative Commons Paternité - Pas d'Utilisation Commerciale - Partage des Conditions Initiales à l'Identique 2.0 France. */
/*Aucune garantie n'est donner sur les effets de cet executable*/
/* pour compiler : gcc modeliseur-stereo.c -lm -o modeliseur-stereo*/

#include 
#include 
#include 
#include 
#define total 258
#define total1 258
#define filtrepts 3

#define efface "joris.wav"
#define cecci1 "21abadie.jo"
#define cecci2 "22abadie.jo"

int jo;
int amplitudemax;
int amplitudemin;
int positionmonte;
int positiondescent;
int position;
int Mise_en_memoire1 [69000];
int Mise_en_memoire2 [69000];
int pourtableau1;
int pourtableau2;
int amplitudemaxc2;
int amplitudeminc2;
int positionmontec2;
int positiondescentc2;
int positionc2;
int Mise_en_memoire1c2 [69000];
int Mise_en_memoire2c2 [69000];
int pourtableau1c2;
int pourtableau2c2;

int main(int argc, char *argv[])

{

/*j ouvre mon fichier de donner et deux fichiers de sortie*/
FILE* oups = NULL;
FILE* joris = NULL;
FILE* jorisc = NULL;

int abadie = 0;
int ji;

oups = fopen(efface, "rb");
joris = fopen(cecci1, "wb");
jorisc = fopen(cecci2, "wb");


	if ((oups != NULL)&&(joris != NULL)&&(jorisc != NULL))
	{
int amplitude = 0;
double ludion11 = 65500;/**amplitude haute***/
double ludion21 = 65500;
double ludion12 = 0;/**amplitude basse***/
double ludion22 = 0;

int sens = 1;
int sensc = 1;

/*debut effet memoire*/
double tamponcomparateur = 0; /* tamponcomparateur est une mise en memoire de la derniere valeur de l'amplitude canal1*/
double tamponcomparateur1 = 0;
double tamponcomparateur2 = 0;
double tamponcomparateur3 = 0;
double tamponcomparateur4 = 0;
double tamponcomparateur5 = 0;

double tamponcomparateurc = 0; /* tamponcomparateurc est une mise en memoire de la derniere valeur de l'amplitude canal2*/
double tamponcomparateurc1 = 0;
double tamponcomparateurc2 = 0;
double tamponcomparateurc3 = 0;
double tamponcomparateurc4 = 0;
double tamponcomparateurc5 = 0;

double nbdepoints1 = 0;/*canal 1*/
double nbdepoints2 = 0;
double nbdepointsc1 = 0;/*canal 2*/
double nbdepointsc2 = 0;

int canal = 1;
int purge = 0;
int assemble = 1;
positionmonte = 0;
positiondescent = 0;
position = 1;
int trouve = 1;
int trouvet = 1;

abadie = fgetc(oups); /* abadie vient chercher tout ce qui bouge*/
		while (abadie != EOF)
		{
if ( trouve < 7 )/*je détruis ici les premiers, */ 
{

switch (trouve) {

case 1 :

if ( abadie == 0x64 )
{
trouve = 2;
}
break;

case 2 :

if ( abadie == 0x61 )
{
trouve = 3;
}
break;

case 3 :

if ( abadie == 0x74 )
{
trouve = 4;
trouvet = 2;
}
break;

case 4 :
if ( (abadie == 0x61) && (trouvet == 2))
{
trouve = 7;
}
break;

default:
trouve = 1;
trouvet = 1;
}

}


			else
			{

unsigned char tab[total] =
{0x80,0x81,0x82,0x83,0x84,0x85,0x86,0x87,0x88,0x89,0x8A,0x8B,0x8C,0x8D,0x8E,0x8F,0x90,0x91,0x92,0x93,0x94,0x95,0x96,0x97,0x98,0x99,0x9A,0x9B,0x9C,0x9D,0x9E,0x9F,0xA0,0xA1,0xA2,0xA3,0xA4,0xA5,0xA6,0xA7,0xA8,0xA9,0xAA,0xAB,0xAC,0xAD,0xAE,0xAF,0xB0,0xB1,0xB2,0xB3,0xB4,0xB5,0xB6,0xB7,0xB8,0xB9,0xBA,0xBB,0xBC,0xBD,0xBE,0xBF,0xC0,0xC1,0xC2,0xC3,0xC4,0xC5,0xC6,0xC7,0xC8,0xC9,0xCA,0xCB,0xCC,0xCD,0xCE,0xCF,0xD0,0xD1,0xD2,0xD3,0xD4,0xD5,0xD6,0xD7,0xD8,0xD9,0xDA,0xDB,0xDC,0xDD,0xDE,0xDF,0xE0,0xE1,0xE2,0xE3,0xE4,0xE5,0xE6,0xE7,0xE8,0xE9,0xEA,0xEB,0xEC,0xED,0xEE,0xEF,0xF0,0xF1,0xF2,0xF3,0xF4,0xF5,0xF6,0xF7,0xF8,0xF9,0xFA,0xFB,0xFC,0xFD,0xFE,0xFF,0x00,0x01,0x02,0x03,0x04,0x05,0x06,0x07,0x08,0x09,0x0A,0x0B,0x0C,0x0D,0x0E,0x0F,0x10,0x11,0x12,0x13,0x14,0x15,0x16,0x17,0x18,0x19,0x1A,0x1B,0x1C,0x1D,0x1E,0x1F,0x20,0x21,0x22,0x23,0x24,0x25,0x26,0x27,0x28,0x29,0x2A,0x2B,0x2C,0x2D,0x2E,0x2F,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x3A,0x3B,0x3C,0x3D,0x3E,0x3F,0x40,0x41,0x42,0x43,0x44,0x45,0x46,0x47,0x48,0x49,0x4A,0x4B,
0x4C,0x4D,0x4E,0x4F,0x50,0x51,0x52,0x53,0x54,0x55,0x56,0x57,0x58,0x59,0x5A,0x5B,0x5C,0x5D,0x5E,0x5F,0x60,0x61,0x62,0x63,0x64,0x65,0x66,0x67,0x68,0x69,0x6A,0x6B,0x6C,0x6D,0x6E,0x6F,0x70,0x71,0x72,0x73,0x74,0x75,0x76,0x77,0x78,0x79,0x7A,0x7B,0x7C,0x7D,0x7E,0x7F};


unsigned char tab1[total1] = 
{0x00,0x01,0x02,0x03,0x04,0x05,0x06,0x07,0x08,0x09,0x0A,0x0B,0x0C,0x0D,0x0E,0x0F,0x10,0x11,0x12,0x13,0x14,0x15,0x16,0x17,0x18,0x19,0x1A,0x1B,0x1C,0x1D,0x1E,0x1F,0x20,0x21,0x22,0x23,0x24,0x25,0x26,0x27,0x28,0x29,0x2A,0x2B,0x2C,0x2D,0x2E,0x2F,0x30,0x31,0x32,0x33,0x34,0x35,0x36,0x37,0x38,0x39,0x3A,0x3B,0x3C,0x3D,0x3E,0x3F,0x40,0x41,0x42,0x43,0x44,0x45,0x46,0x47,0x48,0x49,0x4A,0x4B,0x4C,0x4D,0x4E,0x4F,0x50,0x51,0x52,0x53,0x54,0x55,0x56,0x57,0x58,0x59,0x5A,0x5B,0x5C,0x5D,0x5E,0x5F,0x60,0x61,0x62,0x63,0x64,0x65,0x66,0x67,0x68,0x69,0x6A,0x6B,0x6C,0x6D,0x6E,0x6F,0x70,0x71,0x72,0x73,0x74,0x75,0x76,0x77,0x78,0x79,0x7A,0x7B,0x7C,0x7D,0x7E,0x7F,0x80,0x81,0x82,0x83,0x84,0x85,0x86,0x87,0x88,0x89,0x8A,0x8B,0x8C,0x8D,0x8E,0x8F,0x90,0x91,0x92,0x93,0x94,0x95,0x96,0x97,0x98,0x99,0x9A,0x9B,0x9C,0x9D,0x9E,0x9F,0xA0,0xA1,0xA2,0xA3,0xA4,0xA5,0xA6,0xA7,0xA8,0xA9,0xAA,0xAB,0xAC,0xAD,0xAE,0xAF,0xB0,0xB1,0xB2,0xB3,0xB4,0xB5,0xB6,0xB7,0xB8,0xB9,0xBA,0xBB,0xBC,0xBD,0xBE,0xBF,0xC0,0xC1,0xC2,0xC3,0xC4,0xC5,0xC6,0xC7,0xC8,0xC9,0xCA,0xCB,
0xCC,0xCD,0xCE,0xCF,0xD0,0xD1,0xD2,0xD3,0xD4,0xD5,0xD6,0xD7,0xD8,0xD9,0xDA,0xDB,0xDC,0xDD,0xDE,0xDF,0xE0,0xE1,0xE2,0xE3,0xE4,0xE5,0xE6,0xE7,0xE8,0xE9,0xEA,0xEB,0xEC,0xED,0xEE,0xEF,0xF0,0xF1,0xF2,0xF3,0xF4,0xF5,0xF6,0xF7,0xF8,0xF9,0xFA,0xFB,0xFC,0xFD,0xFE,0xFF};



			double force2;
			int force;
			int force1;
			int groupe1;
			int groupe2;
amplitudemax = 100000;
amplitudemin = 1;
amplitudemaxc2 = 100000;
amplitudeminc2 = 1;

				for (ji=0; ji<255; ji++)
				{
					if ( abadie == tab[ji])
					{
					force1 = ji;
					}
				}

/***debut assemble***/
				if (assemble == 1)
				{
				groupe1 = force1;
				assemble++;
				}

				else
				{
				groupe2 = force1;
				force2 = (double)(groupe1 + 256 * groupe2);
				force= (int)(force2);

					switch ( canal)
					{
					case 1:
					{
/*sens 1 canal 1 je monte*/

						switch (sens)
						{
						case 1:

						/*indique la sortie sens 2 donc j entre en sens 1*/
						{
						double variationAmpl;
						variationAmpl = abs(tamponcomparateur2 - ludion12);
/*

							if ((force < tamponcomparateur)&&( force < tamponcomparateur1 )&&( force < tamponcomparateur2 )&&(variationAmpl > 80 ))*/
if ((force < tamponcomparateur)&&(variationAmpl > 80 ))/*verifie que j ai atteinds le pic*/
							{


								amplitude = (int)tamponcomparateur2;


/*début analyse jo*/

amplitudemax = amplitude;
float difference;
/*la valeur maximal est amplitude puisque j'arrive a la fin d'un front montant*/
float repere;

difference = (float) (amplitudemax - amplitudemin);

repere = ((2 - 1.63212) * difference) + amplitudemin;

for (positionmonte = 0; repere > Mise_en_memoire1 [positionmonte]; positionmonte++) ;

jo = (int)((nbdepoints1 - (2 * positionmonte))*50/nbdepoints1);


if (jo < 7){
jo = 7;
}
/*fin de l'analyse jo**/



char binnb1;
char binnb2;
int pts1;
pts1 = nbdepoints1 * 0.00390625;
int pts2;
pts2 = nbdepoints1 - (256*pts1);
binnb1 = tab1[pts2];
binnb2 = tab1[pts1];
char binamp;
int binamp1;
int binamp2;
char binamp3;
binamp1 = amplitude * 0.00390625;
binamp = tab1[binamp1];
binamp2 = amplitude - (256*binamp1);
binamp3 = tab1[binamp2];
char binjo;
binjo = tab1[jo];

fwrite(&binamp3, sizeof(char), 1, joris);
fwrite(&binamp, sizeof(char), 1, joris);
fwrite(&binnb1, sizeof(char), 1, joris);
fwrite(&binnb2, sizeof(char), 1, joris);
fwrite(&binjo, sizeof(char), 1, joris);

								nbdepoints1 = 0;
								nbdepoints2 = 1;
								sens++; /* j envoie la main a la descente*/
								ludion11 = tamponcomparateur2; /*filtre amplitude*/
								positionmonte = 0;
								pourtableau1 = 0;



								}


	

								else
								{
								nbdepoints1++;
								pourtableau1++;
								Mise_en_memoire1 [pourtableau1] = force;
								}

						}
						break;


/* sens 2 ************************descente******canal1***********************************************/
						default: 
						{
						double variationAmpl;
						variationAmpl = abs(ludion11 - tamponcomparateur2);
/*
							if (( force > tamponcomparateur )&&( force > tamponcomparateur1 )&&( force > tamponcomparateur2 )&&(variationAmpl > 80 )) */
							if (( force > tamponcomparateur )&&(variationAmpl > 80 )) /* je change de sens j'envoie sens 2 + filtre*/

							{

								amplitude = (int)tamponcomparateur2;



/*début analyse jo*/
amplitudemin = amplitude;
float difference;
/*la valeur min est amplitude puisque j'arrive a la fin d'un front descendant*/
float repere;


difference = (float) (amplitudemax - amplitudemin);

repere = ((2 - 1.63212) * difference) + amplitudemin;

for (positiondescent = (nbdepoints2 -1); repere > Mise_en_memoire2 [positiondescent]; positiondescent--) ;

position = nbdepoints2 -1 -positiondescent;

jo = (int)((nbdepoints2 - (2 * position))*50/nbdepoints2);

if (jo < 7){
jo = 7;
}/*corection faible nb de points*/
/*fin de l'analyse jo**/


char binnb1;
char binnb2;
int pts1;
pts1 = nbdepoints2 * 0.00390625;
int pts2;
pts2 = nbdepoints2 - (256*pts1);
binnb1 = tab1[pts2];
binnb2 = tab1[pts1];
char binamp;
int binamp1;
int binamp2;
char binamp3;
binamp1 = amplitude * 0.00390625;
binamp = tab1[binamp1];
binamp2 = amplitude - (256*binamp1);
binamp3 = tab1[binamp2];
char binjo;
binjo = tab1[jo];

fwrite(&binamp3, sizeof(char), 1, joris);/**amplitude poids faible*/
fwrite(&binamp, sizeof(char), 1, joris);/**amplitude poids fort*/
fwrite(&binnb1, sizeof(char), 1, joris);/**durée poids faible**/
fwrite(&binnb2, sizeof(char), 1, joris);/**durée poids fort**/
fwrite(&binjo, sizeof(char), 1, joris);/**jo*/

								sens = 1; /* je prepare changement de sens*/
								nbdepoints2 = 0;
								nbdepoints1 = 1;
								ludion12 = tamponcomparateur2; /*filtre amplitude*/
								pourtableau2 = 0;
								positiondescent = 0;
								}


								else
								{
								nbdepoints2++;
								pourtableau2++;
								Mise_en_memoire2 [pourtableau2] = force;

								}

							}
						}
/*****pas touche *****/

/* ici je mets en memoire pour traitement apres lecture suivante*/

					tamponcomparateur = (double)force;
					tamponcomparateur1 = tamponcomparateur;
					tamponcomparateur2 = tamponcomparateur1;
					tamponcomparateur3 = tamponcomparateur2;
					tamponcomparateur4 = tamponcomparateur3;
					tamponcomparateur5 = tamponcomparateur4;

					canal++;
					}/***fin canal 1***/
					
					break;
/************************************canal2***************************************************************************/

					default: /*****debut canal 2********/
					{
/*sens 1*/
						switch (sensc)
						{
						case 1:

						{
						double variationAmpl;
						variationAmpl = abs(tamponcomparateurc2 - ludion22);
					/*		if ((force < tamponcomparateurc)&&( force < tamponcomparateurc1 )&&( force < tamponcomparateurc2 )&&(variationAmpl > 80 ))*/


							if ((force < tamponcomparateurc)&&(variationAmpl > 80 ))/*verifie que j ai atteinds le pic*/
							{


								amplitude = (int)tamponcomparateurc2;

/*début analyse jo*/

amplitudemaxc2 = amplitude;
float difference;
/*la valeur maximal est amplitude puisque j'arrive a la fin d'un front montant*/
float repere;

difference = (float) (amplitudemaxc2 - amplitudeminc2);

repere = ((2 - 1.63212) * difference) + amplitudeminc2;

for (positionmontec2 = 0; repere > Mise_en_memoire1c2 [positionmontec2]; positionmontec2++) ;

jo = (int)((nbdepointsc1 - (2 * positionmontec2))*50/nbdepointsc1);


if (jo < 7){
jo = 7;
}
/*fin de l'analyse jo**/


char binnb1;
char binnb2;
int pts1;
pts1 = nbdepointsc1 * 0.00390625;
int pts2;
pts2 = nbdepointsc1 - (256*pts1);
binnb1 = tab1[pts2];
binnb2 = tab1[pts1];
char binamp;
int binamp1;
int binamp2;
char binamp3;
binamp1 = amplitude * 0.00390625;
binamp = tab1[binamp1];
binamp2 = amplitude - (256*binamp1);
binamp3 = tab1[binamp2];
char binjo;
binjo = tab1[jo];

fwrite(&binamp3, sizeof(char), 1, jorisc);/**amplitude poids faible*/
fwrite(&binamp, sizeof(char), 1, jorisc);/**amplitude poids fort*/
fwrite(&binnb1, sizeof(char), 1, jorisc);/**durée poids faible**/
fwrite(&binnb2, sizeof(char), 1, jorisc);/**durée poids fort**/
fwrite(&binjo, sizeof(char), 1, jorisc);/**jo*/




								nbdepointsc1 = 0;
								nbdepointsc2 = 1;
								sensc++; /* j envoie la main a la descente*/
								ludion21 = tamponcomparateurc2; /*filtre amplitude*/
								positionmontec2 = 0;
								pourtableau1c2 = 0;
								}



	
								else
								{
								nbdepointsc1++;
								pourtableau1c2++;
								Mise_en_memoire1c2 [pourtableau1c2] = force;
								}

						}
						break;
/* sens 2 ****************descente********************************************************************************************/
						default: 

						{
						double variationAmpl;
						variationAmpl = abs(ludion21 - tamponcomparateurc2);
					/*		if (( force > tamponcomparateurc )&&( force > tamponcomparateurc1 )&&( force > tamponcomparateurc2 )&&(variationAmpl > 80 ))*/

							if (( force > tamponcomparateurc )&&(variationAmpl > 80 )) /* je change de sens j'envoie sens 2 + filtre*/
							{

								amplitude = (int)tamponcomparateurc2;
/*début analyse jo*/
amplitudeminc2 = amplitude;
float difference;
/*la valeur min est amplitude puisque j'arrive a la fin d'un front descendant*/
float repere;

difference = (float) (amplitudemaxc2 - amplitudeminc2);

repere = ((2 - 1.63212) * difference) + amplitudeminc2;

for (positiondescentc2 = (nbdepointsc2 -1); repere > Mise_en_memoire2c2 [positiondescentc2]; positiondescentc2--) ;

position = nbdepointsc2 -1 -positiondescentc2;

jo = (int)((nbdepointsc2 - (2 * position))*50/nbdepointsc2);

if (jo < 7){
jo = 7;
}/*corection faible nb de points*/
/*fin de l'analyse jo**/

char binnb1;
char binnb2;
int pts1;
pts1 = nbdepointsc2 * 0.00390625;
int pts2;
pts2 = nbdepointsc2 - (256*pts1);
binnb1 = tab1[pts2];
binnb2 = tab1[pts1];
char binamp;
int binamp1;
int binamp2;
char binamp3;
binamp1 = amplitude * 0.00390625;
binamp = tab1[binamp1];
binamp2 = amplitude - (256*binamp1);
binamp3 = tab1[binamp2];
char binjo;
binjo = tab1[jo];

fwrite(&binamp3, sizeof(char), 1, jorisc);/**amplitude poids faible*/
fwrite(&binamp, sizeof(char), 1, jorisc);/**amplitude poids fort*/
fwrite(&binnb1, sizeof(char), 1, jorisc);/**durée poids faible**/
fwrite(&binnb2, sizeof(char), 1, jorisc);/**durée poids fort**/
fwrite(&binjo, sizeof(char), 1, jorisc);/**jo*/


								sensc = 1; /* je prepare le passage a la monter*/
								nbdepointsc2 = 0;
								nbdepointsc1 = 1;
								ludion22 = tamponcomparateurc2; /*filtre amplitude*/
								pourtableau2c2 = 0;
								positiondescentc2 = 0;
								}



								else
								{
								nbdepointsc2++;
								pourtableau2c2++;
								Mise_en_memoire2c2 [pourtableau2c2] = force;
								}

						}/**fin sensc***/
						}/***fin switch sensc****/
					tamponcomparateurc = (double)force;/* ici je mets en memoire pour traitement apres lecture suivante*/
					tamponcomparateurc1 = tamponcomparateurc;
					tamponcomparateurc2 = tamponcomparateurc1;
					tamponcomparateurc3 = tamponcomparateurc2;
					tamponcomparateurc4 = tamponcomparateurc3;
					tamponcomparateurc5 = tamponcomparateurc4;

					canal = 1;
					}/***fin canal 2***/
					}/***fin switch canal****/
/**************************fin canal2****************************/
				assemble = 1;
				}/****fin assemble***/
			}/**fin choix purge***/
		abadie = fgetc(oups); /* abadie va chercher la suivante*/
		}/***fin de la boucle***/
	fclose(jorisc), jorisc = NULL;;/*****je ferme jorisc *****/
	fclose(joris), joris = NULL;;/*****je ferme joris *****/
	fclose(oups), oups = NULL;;
	}

return 0;
}