1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70
|
void FrequencyTransposition(double *signal_reel, float tk, int N)
{
int i;
float fp;
double signal_re[N];
double cosinus[N];
double sinus[N];
double signal_re_transpose[N];
double signal_im_transpose[N];
FILE * FicSortie;
printf("\nN=%d",N);
printf("\ntk=%f",tk);
for(i=0;i<N;i++)
{
signal_re[i]=signal_reel[i]; /*Récupération du signal partie réelle*/
}
printf("\n\nFréquence de la porteuse en MHz :");
scanf("%f",&fp);
printf("\n\nSignal partie réelle transposé par une porteuse de %lf MHz :",fp);
/*Partie reelle*/
for(i=0;i<N;i++)
{
cosinus[i] = cos(2*3.1416*fp*i*tk);
signal_re_transpose[i]=signal_re[i]*cosinus[i]; /*Multiplication du signal reel par cos(2*pi*fp*tk)*/
}
/*Partie imaginaire*/
for(i=0;i<N;i++)
{
sinus[i] = sin(2*3.1416*fp*i*tk);
signal_im_transpose[i]=signal_re[i]*sinus[i]; /*Multiplication du signal par sin(2*pi*fp*tk)*/
}
FicSortie= fopen("F:/Users/garnier/Matlab/imag.txt","w");
for (i=0;i<N;++i)
{fprintf(FicSortie,"%lf\n ",signal_im_transpose[i]);
}
fclose(FicSortie);
FicSortie= fopen("F:/Users/garnier/Matlab/reel.txt","w");
for (i=0;i<N;++i)
{fprintf(FicSortie,"%lf\n ",signal_re_transpose[i]);
printf(" %lf ",signal_re_transpose[i]);
}
fclose(FicSortie);
return;
} |
Partager