The Program

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#define MAX 300

main()

{

int r, i, s, t, max_coeff;

float scale;

float h1[MAX], h2[MAX], g[MAX];

scanf("%d", &r);

if (r<= 0 || r > MAX) {

printf("Invalid resolution!\n");

exit(1);

}

scale = 2 * M_PI/r;

for(i = 0; i < r; ++i)

scanf("%f", &g[i]);

for(s = 0; s <= max_coeff; ++s) {

h1[s] = 0;

h2[s] = 0;

for(t = 0; t < r; ++t) {

h1[s] += g[t] * cos(scale*s*t);

h2[s] += g[t] * sin(scale*s*t);

}

}

printf("\n\nSpectrum:\n");

for(i = 1; i <= max_coeff; ++i)

printf("%f\n",

sqrt(h1[i]*h1[i] + h2[i]*h2[i]) );

}

This is the program pg241.c

Previous slide

Next slide

Back to first slide

View graphic version