% (c31.m)%M is the number of coefficients to be computedT0=pi;N0=256;T=T0/N0;M=10;t=0:T:T*(N0-1);t=t';f=exp(-t/2);f(1)=0.604;Dn=fft(f/N0)[Dnangle,Dnmag]=cart2pol(real(Dn),imag(Dn));k=0:length(Dn)-1;k=k';subplot(211),stem(k,Dnmag)subplot(212), stem(k,Dnangle)%Computation of trigonometric Fourier coefficients:%c31;clgC0=Dnmag(1); Cn=2*Dnmag(2:M);Amplitudes=[C0;Cn]Angles=Dnangle(1:M);Angles=Angles*(180/pi);disp('Amplitudes Angles')[Amplitudes Angles]k=0:length(Amplitudes)-1; k=k';subplot(211),stem(k,Amplitudes)subplot(212), stem(k,Angles)