Pasted

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "complex.h"
 
#define FS 32.78 *pow(10,6)
 
 
void gen_sig(float *sig, int size);
 
void g_sig(struct complex *data, struct complex32 *data32, struct complex16 *data16,int size);
 
int main(){
    const int size = 307200;
    printf("%d\n", size);
    float sig[size];
    struct complex *data;
    struct complex32 *data32;
    struct complex16 *data16;
 
 
 
 
 
    //data   = malloc(sizeof(struct complex));
   // data32 = malloc(sizeof(struct complex32));
    //data16 = malloc(sizeof(struct complex16));
 
    printf("FS: %f\n" ,FS);
    if ((data=malloc(sizeof(struct complex)*(size_t)size))==NULL)
    {
      fprintf(stderr, "Out of memory!\n");
      exit(1);
    }
 
  if ((data16=malloc(sizeof(struct complex16)*(size_t)size))==NULL)
    {
      fprintf(stderr, "Out of memory!\n");
      exit(1);
    }
 
  if ((data32=malloc(sizeof(struct complex32)*(size_t)size))==NULL)
    {
      fprintf(stderr, "Out of memory!\n");
      exit(1);
    }
 
 
    g_sig(data,data32, data16, size);
    
return 0;
}
 
void g_sig(struct complex *data, struct complex32 *data32, struct complex16 *data16,int size){
 
    int write = 0;
    int j;
    // Generating a table got 100 MHz channel
    int cp0 = 160; // First cyclic prefix
    int cp  = 144; // Normal cyclic prefix
    int cpe = 512; // Extended cyclic prefix
    int deltaF = 200; //Hz
    int ofdm_size = 2048;
    for (j=0;j<size;j++){
        data[j].r = cos(20*pow(10,6)*j/FS);//Is is Mhz on this line;
        data[j].i = sin(20*pow(10,6)*j/FS);
    }
    //printf("%d",(1<<15)); 
    //Truncate data to fixed-point representation
    // I've used 16 and 32 bits representation, 
    for (j=0; j<size;j++){
        data16[j].r = (short)data[j].r*(1<<15)-1;
        data16[j].i = (short)data[j].i*(1<<15)-1;
        data32[j].r = (int)data[j].r*(1<<15)-1;
        data32[j].i = (int)data[j].i*(1<<15)-1;
        printf("%f\n %f\n", data16[j].r, data16[j].i);
                        
    }
}
 
 
 
 
 
void gen_sig(float *sig, int size){
    int write = 0;
    // Generating a table got 100 MHz channel
    int cp0 = 160; // First cyclic prefix
    int cp  = 144; // Normal cyclic prefix
    int cpe = 512; // Extended cyclic prefix
    int deltaF = 200; //Hz
    int ofdm_size = 2048;
    
    if (write){ //write to file.txt
        FILE *f;
        f = fopen("file.txt", "w"); 
        if(f==NULL){
            exit(1); 
        }
        fprintf(f, "x[%d]={ ", size);
            
        for (int i=1;i<size; i+=2){
            sig[i-1] = sin(20*i/FS);
            sig[i]   = cos(20*i/FS);
            fprintf(f, "%f, ", sig[i-1]);
        }
        fprintf(f,"}; \n");
    }   
    for (int i=1;i<size; i+=2){
        sig[i-1] = sin(20*i/FS);
        sig[i]   = cos(20*i/FS);
    }
}