Question

During lab 4, we have seen numerical implementation of Fourier Series for periodic signals. As first part of this assignment, you need to write a Matlab function that would take an array representing a single period of a signal (x), corresponding time array (t), and return the Fourier Series coefficients (Ck) in exponential form. The function should also be able to take two (2) optional input arguments: number of Fourier coefficients (Nk) and plot option (p). Use the template ‘fourier_series_exp.m’ for this problem.

(Template "fourier_series_exp.m")

function [Ck]=fourier_series_exp(x,t,Nk,p)

% Ck = exponential fourier series cofficient
% x = single period of a signal
% t = time corrosponding to 'x'
% Nk = (optional input) number of exponential terms
% p = plotting option ; p=0, no plots, p = 1 plot Ck vs k and reconstructed signal
% dT = t(2)-t(1) = temporal resolution of signal (x)
% T = peiod of signal 'x'
% w0= angular frequency of signal 'x'

dT=t(2)-t(1);
T= dT*length(t);
w0=2*pi/T;   

% Check the number of inputs, 'nargin' returns number of input arguments
if nargin <2
error('Not enough input argument!')
elseif nargin == 2
Nk=101; % you can set any default value you like
p=0; % not plots
elseif nargin ==3
p=0; % not plots
end
k=-floor(Nk/2):floor(Nk/2); % if Nk=11, k=-5:5; if Nk=12, k=-6:6

%% evaluate Ck
%
% % % write this code segment
%
%% plot spectrum and reconstructed signal   
if p==1
% plot abs(Ck) vs k and angle(Ck) vs k
%
% % % write this code segment
%
  
% plot 3 cycles of the signal 'x' and the reconstructed signal
%
% % % write this code segment
%
end

end

Answer 1

function [Ck]=fourier_series_exp(x,t,Nk,p)

% Ck = exponential fourier series cofficient
% x = single period of a signal
% t = time corrosponding to 'x'
% Nk = (optional input) number of exponential terms
% p = plotting option ; p=0, no plots, p = 1 plot Ck vs k and reconstructed signal
% dT = t(2)-t(1) = temporal resolution of signal (x)
% T = peiod of signal 'x'
% w0= angular frequency of signal 'x'
Fs=100;

dT=t(2)-t(1);
T= dT*length(t);
w0=2*pi/T;   

% Check the number of inputs, 'nargin' returns number of input arguments
if nargin <2
error('Not enough input argument!')
elseif nargin == 2
Nk=101; % you can set any default value you like
p=0; % not plots
elseif nargin ==3
p=0; % not plots
end
k=-floor(Nk/2):floor(Nk/2); % if Nk=11, k=-5:5; if Nk=12, k=-6:6

%% evaluate Ck
%
% % % write this code segment
%Ck = fft(x);

L=length(x);
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(x,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2);

%
%% plot spectrum and reconstructed signal   
if p==1
% plot abs(Ck) vs k and angle(Ck) vs k
%
% % % write this code segment
% Plot single-sided amplitude spectrum.
subplot(3,1,1); plot(x); title('Original Signal - Single Period');
subplot(3,1,2); plot(f,2*abs(Y(1:NFFT/2)));
str = strcat('FFT at N = ',num2str(NFFT),' Points'); title(str);
xlabel('Frequency (Hz)'); ylabel('|Y(f)|');

%Reconstructed Signal
y =ifft(Y);
subplot(3,1,3); plot(y); title('Reconstructed Signal');
%
end

end

Original Signal - Single Period FFT at N 16 Points 05 10 15 20 25 30 35 45 50 Frequency (Hz) Reconstructed Signal 0.05 0.05