#include <stdlib.h> //standard libreria. include significa che i contenuti dei files .h vengono inclusi nel programma
#include <stdio.h> //standard input output
#include <time.h> //funzione del tempo di sistema
#include "wavfile.h" //header contenente intestazione file wave. quando ci sono le virgolette il file di intestazione sta nella stessa cartella del sorgente
#include <string.h>
#include <math.h>
#include <malloc.h>

#define tablen 4096 // simboli di costanti
#define pi2 6.283185307179586476925286766559

void main () 
{	
	float sr=44100.0;
	int i,samp; //variabili numero intero a 32 (4 byte)
	double x; //numero con virgola a DOppia precisione (8 Byte)
	short int *seno,*buffer; //puntatore a intero a 16 bit
	double sintable [tablen]; //alloco tablen (4096 double per 8 byte (4096* byte))
	struct WavHdr header; // riserva uno spazio di memoria come è stato definito
	float secondi,frequenza;//con virgola semplice (4 byte)
	FILE *fp; //serve per aprire il file wave

	x = pi2/(double)tablen; 
	printf ("frequenza: "); 
	scanf ("%f", &frequenza);
	//printf ("durata: ");
	//scanf ("%f", &secondi);




	fp=fopen ("Flauto_eb.wav","rb"); //apre il file .wav
	if (fp==NULL) {printf ("errore in lettura"); exit(1);}
	fread (&header,sizeof(struct WavHdr),1,fp);// legge l'intestazione del file .wav
	secondi=((float)header.samp_num/(float)header.samp_rate)/2;//calcolo della durata del file Flauto_eb.wav
	samp = (int)((float)header.samp_rate*secondi); //numeri campioni dell'onda
	buffer= (short int*)calloc (samp,sizeof(short int));
	fread(buffer,sizeof(short int), samp, fp);// lettura dei campioni sonori
	fclose(fp);
	seno= (short int*)calloc (samp,sizeof(short int)); //calloc permette una allocazione dinamica,
	//(array è di misura variabile)



				for	(i=0; i<tablen; i++) 
				{
					sintable[i] = sin (x*i);

				}
	
	 
	int z =0; //fase istantanea
	int step = (int)(frequenza*(float)tablen/sr);


// il modo per ottenere la frequenza desiderata è quello di far leggere la tabella con un passo variabile: più è grande il passo, più velocemente viene letta la tabella, quindi maggiore sarà la frequenza. 
	
 
				for (i=0; i<samp; i++)  
				{ 
					z = z % tablen; // % è il resto della divisione fra z (la fase) e tablen: la fase non deve superare la lunghezza della tabella (tablen -1)
					seno[i]=(short int)(32000.0*sintable[z]*((float)buffer[i]/32000.0));//si normalizza a 1
					//seno [i] = (short int) (32000.0*sintable [z]); //riscalamento onda a 16 bit
					z= z+ step; //aggiornamento della fase istantanea
				} 

	fp = fopen ("sinusoide.wav","wb"); // apertura in scrittura del file wave

	strncpy(header.riff, "RIFF",4); //record relativo all'intestazione del file wave contenente l'onda creata
	header.file_dim = samp*2+44;
	strncpy(header.wave_format,"WAVEfmt ",8);          
	header.reserved1 = 16;   
	header.reserved2 = 1;   
	header.reserved3 = 1;
	header. samp_rate = 44100;
	header.fc =  44100*2;
	header.stereo_mono = 2;
	header.bit_num = 16;
	strncpy(header.dati,"data",4);                 
	header.samp_num= samp*2;
	
	fwrite ((void *)&header,sizeof(header),1,fp);// scrittura su file wave
	fwrite ((void *)seno,sizeof(short int),samp,fp); //scrittura dei campioni
	fclose (fp);
}