0% found this document useful (0 votes)
9 views3 pages

IDTFT-Based FIR LPF Design Guide

The document details the design of an Ideal Low Pass Filter (LPF) using the Inverse Discrete Time Fourier Transform (IDTFT) approach. It includes symbolic integration and differentiation, evaluation of the IDTFT expression, and the calculation of filter coefficients from the impulse response. Additionally, it presents the magnitude response of the FIR-LPF and compares it with elementary cascade LPF filters.

Uploaded by

mn22eeb0a09
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
9 views3 pages

IDTFT-Based FIR LPF Design Guide

The document details the design of an Ideal Low Pass Filter (LPF) using the Inverse Discrete Time Fourier Transform (IDTFT) approach. It includes symbolic integration and differentiation, evaluation of the IDTFT expression, and the calculation of filter coefficients from the impulse response. Additionally, it presents the magnitude response of the FIR-LPF and compares it with elementary cascade LPF filters.

Uploaded by

mn22eeb0a09
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

FIR Ideal Filter design using Inverse Discrete Time Fourier

Transform (IDTFT) Approach


(A) Symbolic integration and differentiation

%%%% Symbolic Integration and Differentiation


clear all
clc
%%%% Symbolic Integration
syms x a
f=exp(2*x);
y=int(f,[0 1])

f1=exp(2*x*a);
y1=int(f1,[0 1])

%%%% Symbolic differentiation


f2=x^2+2*x+1;
y2=diff(f2)

%%%% Substitute in the above symbolic expression


y3=subs(y2,1)

(B) Design the Ideal LPF with following specifications

From the Inverse Discrete time Fourier Transform (IDTFT) approach


%%% LPF design using Inverse Discrete Time Fourier Transform (IDTFT) approach %%%%
clear all
clc

%%%% Step1: Evaluate the IDTFT Expression using symbolic integration


syms w n
f= exp(i*w*n); %%% Here i = complex no, w = frequency, n = samples

y=(0.5/pi)*int(f,w, [-pi/2,pi/2]); %%%%%%%% IDTFT evaluation with the lower and


upper limits from -pi/2 to pi/2, respectively

%%%% Here intergration is done by using the command 'int'

%%%% The below expression 'out' is obtainted by typing 'y' in the command
%%%% window and then copy-paste it here (DONT WRITE/TYPE IT)
out=-(exp(-(pi*n*1i)/2)*(exp(pi*n*1i)-1)*5734161139222659i)/(36028797018963968*n);

%%%% Step 2: Obtain the Impulse response plot to get the filter coefficient %%%%

%%%% For n = 0, the above expression becomes indeterminate form, so apply the L-
Hopitals rule
%%% Based on the rule: Evaluate the ratio of differentiation ([Link] 'n') between
numerator and denominator
outnum=-(exp(-(pi*n*i)/2)*(exp(pi*n*i)-1)*5734161139222659i); %%% Numerator
expression
outden=(36028797018963968*n); %%% Denominator expression
num1=diff(outnum); %%%% Differentiate the numerator term using 'diff' command
den1=diff(outden); %%%% Differentiate the denominator term using 'diff' command

%%%% Substitute the value of n= 0 , in the differentiation ratio using 'subs'


command
%n=0;
h0=subs(num1,0)/subs(den1,0); %%% Symbolic expression to evaluate the output of h0
from the command window
h0= 0.5; %%%% h0 value can be obtained, from the command window, by typing h0
after executing the above line
ii=1;

%%% Plot the Impulse response


for n= -15:1:15
if n==0
h=h0;
stem(n,h,'r','LineWidth',2);hold on; grid on
else
h=-(exp(-(pi*n*1i)/2)*(exp(pi*n*1i) -
1)*5734161139222659i)/(36028797018963968*n); %%% Integration output
stem(n,h,'r','LineWidth',2);hold on; grid on
end
end
xlabel ('No of Samples');ylabel ('Impulse response');

%%% Step 3: To get the filter coefficients from the truncated impulse response
plot
%%%%%% Visually find the length to be 11 (-5 to 5) becz most of the information in
the impulse response plot is concentrated with the length of 11
%%%% Take only positive values of n ( due to the symmetry of an impulse response)
clear n %%%% Command to delete the symbolic variable n

ii=1;
for n=-5:1:5
if n== 0
y=h0;
else
y=-(exp(-(pi*n*1i)/2)*(exp(pi*n*1i) -
1)*5734161139222659i)/(36028797018963968*n);
end
f_coeff(ii)=real(y);
ii=ii+1;
end
Filter_Coefficients= f_coeff;

%%%%% Step 4: Apply the Z-transform with above calculate filter coefficient and
then multiply the transfer function
%%%%% TF with higher negative degree term (here z^-5) to get the causal TF
%%%% Obtained Digital LPF transfer function using IDTFT approach

b=[0.06366 0 -0.106 0 0.3183 0.5 0.3183 0 -0.106 0 0.06363];a=1;


[H,w] = freqz(b,a); %%%% Magnitude response

figure(2)
plot(w/pi,abs(H),'b','LineWidth',2); grid on; hold on
xlabel ('Normalized Frequency (w/ \pi)');ylabel ('Magnitude in absolute value');
title('Magnitude response of FIR-LPF: IDTFT Approach ')

%%%% Elementary Cascade LPF_ FIR Filter: Frequency response %%%%%


b1=[1 1];a1=2; %%%%% Elementary Digital LPF
[h1,w] = freqz(b1,a1);
b2=conv([1 1],[1 1]);a2=4; %%%%% Cascade Two LPF
[h2,w] = freqz(b2,a2);
b3=conv([1 1],conv([1 1],[1 1]));a3=8; %%%%% Cascade Three LPF
[h3,w] = freqz(b3,a3);
b4=conv([1 1],conv([1 1],conv([1 1],[1 1])));a4=16; %%%%% Cascade Four LPF
[h4,w] = freqz(b4,a4);
%figure(3)
plot(w/pi,abs(h1),'g-.','LineWidth',2); grid on; hold on
plot(w/pi,abs(h2),'r--','LineWidth',2); grid on; hold on
plot(w/pi,abs(h3),'m--','LineWidth',2); grid on; hold on
plot(w/pi,abs(h4),'b-.','LineWidth',2); grid on; hold on
xlabel ('Normalized Frequency (w/ \pi)');ylabel ('Magnitude in absolute value');
legend('LPF_IDTFT','Single LPF','Two LPF','Three LPF','Four LPF')

Common questions

Powered by AI

The design of an ideal low-pass filter (LPF) using the Inverse Discrete Time Fourier Transform (IDTFT) approach involves several key steps: First, symbolically integrate the IDTFT expression over a given frequency range to obtain an impulse response expression. Next, compute the impulse response at the specific point where it is indeterminate by applying L'Hopital's rule. After that, plot the impulse response to determine the filter coefficients, focusing on a symmetrical range [-5,5] where most of the information is concentrated. The numerical values of the impulse response (focusing only on positive values due to symmetry) provide the filter coefficients. Finally, apply the Z-transform to these coefficients and multiply the resulting transfer function by a term representing delay to produce a causal transfer function .

Cascading LPF sections involves combining multiple low-pass filter stages to achieve greater attenuation of unwanted frequencies. Each additional section effectively increases the filter's roll-off rate and improves stop-band attenuation. The document describes constructing increasingly complex filters by convolving the coefficients of simpler filters to form composite transfer functions (e.g., progressing from single LPF, to two, three, and four LPF cascades). This technique enhances the filter's overall performance by creating steeper transition bands with better stop-band characteristics .

The symmetry property of the impulse response is crucial because it implies that the filter can efficiently process signals using only the positive values of n, given that the negative side is simply the mirror image of the positive side. This simplifies the computation and facilitates the construction of the filter coefficients, allowing for more efficient implementation of the low-pass filter .

The magnitude response plot is significant in evaluating the performance of an FIR LPF designed using the IDTFT approach. It shows how the filter attenuates or passes signals across different frequencies, which is essential to ensure that low frequencies are passed while higher frequencies are attenuated. By plotting the magnitude response, one can assess the passband, stopband, and transition band characteristics, thereby verifying the effectiveness of the filter design and making adjustments if necessary .

L'Hopital's rule is applied in the computation of the impulse response when the expression becomes indeterminate for n=0. The impulse response derived from the IDTFT is initially indeterminate at n = 0, so the rule is employed to evaluate the limit by calculating the derivatives of the numerator and the denominator separately with respect to n. This operation allows the determination of a finite value for the impulse response at this point, which is found to be h0 = 0.5 .

In plotting the magnitude response of FIR filters, frequency normalization serves to scale the frequency axis relative to the Nyquist frequency, represented as a fraction of π. This normalization enables the magnitude response to be more easily interpreted across different sampling rates and filter configurations, facilitating comparison and analysis of the filter's performance irrespective of the specific sample rate .

The application of the Z-transform is essential for obtaining a causal transfer function as it converts the sequence of filter coefficients from the impulse response into a polynomial representation in the z-domain. This transformation accommodates delays in signal processing by introducing appropriate powers of z (such as z^-5), which ensures causality by accounting for the system's requirement to only output responses based on past and present input samples, not future ones .

The filter coefficients derived from the impulse response are used directly to form the numerator of the transfer function for the digital low-pass filter. In the IDTFT approach, these coefficients represent the weights applied to different delays in the signal and constitute the polynomial in the Z-domain. The resulting transfer function is then delayed by multiplying with an additional factor such as z^-5 to ensure causality in practical systems .

Symbolic computation enhances the FIR filter design process by allowing exact analytical integration and differentiation, leading to precise expressions for the impulse response. This approach facilitates handling complex mathematical operations inherent in IDTFT, such as integration over infinite limits and differentiation of indeterminate forms, thereby yielding accurate filter coefficients crucial for high-performance filter design .

Cascading multiple LPF sections enhances filter efficiency by significantly improving the steepness of the filter’s roll-off and increasing the attenuation rate in the stop band. The document indicates that composite filters with cascaded sections exhibit better performance metrics due to aggregated filtering effects, which can more effectively distinguish and suppress undesired higher frequencies while maintaining low-frequency passband fidelity, thus facilitating more stringent filtering requirements without a disproportionate increase in complexity .

You might also like