/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 13.04.2 ] */

/* [wxMaxima: input   start ] */
kill(all)$
N:8$ k:makelist(k-N/2,k,0,N);T:1$ f0:1/T$ 
x:1$ a:1/2;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
/* Reference Signal */
R(t):=if abs(t)<= a/2 then 1 else 0$
wxplot2d(R(t),[t,-T/2,T/2],[y,-0.2,1.2],[xlabel,"t"],[ylabel,"Amplitude"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
/* Calculation x(t) -> X(f) */
A:exp(-%i*2*%pi*f0*k*t);
C:expand(integrate(A*x,t,-a/2,a/2))/T;
wxplot2d([discrete,k/T,realpart(C)],[style,[points,2,1,1]],[xlabel,"f[Hz]"],[ylabel,"Ck"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
/* Calculation X(f) -> x(t) */
B:exp(%i*2*%pi*f0*k*t);
xx:B.C;
wxplot2d([realpart(xx),R(t)],[t,-T/2,T/2],[ylabel,"Amplitude"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wxplot2d(realpart(xx),[t,-4*T/2,4*T/2],[ylabel,"Amplitude"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
/* Increase number of Ck */
kill(all)$
N:32$ k:makelist(k-N/2,k,0,N);T:1$ f0:1/T$ 
x:1$ a:1/2;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
A:exp(-%i*2*%pi*f0*k*t);
C:expand(integrate(A.x,t,-a/2,a/2))/T;
wxplot2d([discrete,k/T,realpart(C)],[style,[points,2,1,1]],[xlabel,"f[Hz]"],[ylabel,"Ck"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
B:exp(%i*2*%pi*f0*k*t);
xx:B.C;
wxplot2d(realpart(xx),[t,-T/2,T/2],[ylabel,"Amplitude"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
/* Increase Interval T with number of Ck */
kill(all)$
N:32$ k:makelist(k-N/2,k,0,N);T:4$ f0:1/T$ 
x:1$ a:1/2;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
A:exp(-%i*2*%pi*f0*k*t);
C:expand(integrate(A.x,t,-a/2,a/2))/T;
wxplot2d([discrete,k/T,realpart(C)],[style,[points,2,1,1]],[xlabel,"f[Hz]"],[ylabel,"Ck"]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
B:exp(%i*2*%pi*f0*k*t);
xx:B.C;
wxplot2d(realpart(xx),[t,-T/2,T/2],[ylabel,"Amplitude"]);
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$