#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
#define sr 44100.0

void main () 
{	
	int i,samp; //variabili numero intero a 32 (4 byte)
	double x; //numero con virgola a DOppia precisione (8 Byte)
	short int *seno; //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)


	x = pi2/(double)tablen; 
	printf ("frequenza: "); 
	scanf ("%f", &frequenza);
	printf ("durata: ");
	scanf ("%f", &secondi);
	samp = (int)(sr*secondi); //numeri campioni dell'onda

	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);

				}
	
	FILE *fp; //serve per aprire il file wave 
	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]); //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);
}