DSP Lab Manual 2018-19 PrintCopy MANU

DSP LAB MANUAL Dept of ECE

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

LABORATORY MANUAL
DIGITAL SIGNAL PROCESSING
III B.TECH -II Semester (ECE)

A.Y 2019-2020
1
VISION AND MISSION OF THE INSTITUTE
INSTITUTE VISION:
To establish as an ideal academic institutions in the service of the nation, the world
and the humanity by graduating talented engineers to be ethically strong, globally
competent by conducting high quality research, developing breakthrough technologies,
and disseminating and preserving technical knowledge.
INSTITUTE MISSION:
To fulfill the promised vision through the following strategic characteristics and
aspirations:
1. An atmosphere that facilitates personal commitment to the educational success of

students in an environment that values diversity and community.
2. Prudent and accountable resource management;
3. Undergraduate programs that integrate global awareness, communication skills and
team building across the curriculum;
4. Leadership and service to meet society’s needs;
5. Education and research partnerships with colleges, universities, and industries to
graduate education and training that prepares students for interdisciplinary
engineering research and advanced problem solving;
6. Highly successful alumni who contribute to the profession in the global society.
2
VISION AND MISSION OF THE DEPARTMENT
DEPARTMENT VISION:
Imparting quality technical education through research, innovation and team work for a
lasting technology development in the area of Electronics and Communication
Engineering.
MISSION:
To develop a strong centre of excellence for education and research with excellent
infrastructure and well qualified faculties to instill in them a passion for knowledge.
To achieve the Mission the department will
1. Establish a unique learning environment to enable the students to face the challenges
of the Electronics and Communication Engineering field.
2. Promote the establishment of centre of excellence in niche technology areas to nurture
the spirit of innovation and creativity among faculty and students.
3. Provide ethical and value based education by promoting activities addressing the
societal needs.
4. Enable students to develop skills to solve complex technological problems of current
times and also provide a framework for promoting collaborative and multidisciplinary
activities.
PROGRAMME EDUCATIONAL OBJECTIVES

1. PEO 1: have successful careers in Industry.
2. PEO 2: show excellence in higher studies/ Research.
3. PEO 3: Show good competency towards Entrepreneurship.
3
PROGRAM OUTCOMES
a An ability to apply knowledge of Science, Mathematics, Engineering & Computing
fundamentals for the solutions of Complex Engineering problems
b An ability to identify, formulates, research literature and analyze complex engineering
problems using first principles of mathematics and engineering sciences.
c An ability to design solutions to complex process or program to meet desired needs.
d Ability to use research-based knowledge and research methods including design of
experiments to provide valid conclusions.
e An ability to use appropriate techniques, skills and tools necessary for computing practice.
f Ability to apply reasoning informed by the contextual knowledge to assess social issues,
consequences & responsibilities relevant to the professional engineering practice.
g Ability to understand the impact of engineering solutions in a global, economic,
environmental, and societal context with sustainability.
h An understanding of professional, ethical, Social issues and responsibilities.
i An ability to function as an individual, and as a member or leader in diverse teams and in

multidisciplinary settings.
j An ability to communicate effectively on complex engineering activities within the
engineering community.
k Ability to demonstrate and understanding of the engineering and management principles as
a member and leader in a team.
l Ability to engage in independent and lifelong learning in the context of technological
change.
PROGRAM SPECIFIC OUTCOMES
PSO1 Analyze and design analog & digital circuits or systems for a given specification and function.
Implement functional blocks of hardware-software co-designs for signal processing and
PSO2
communication applications.
4
COURSE OUTCOMES FOR DIGITAL SIGNAL PROCESSING LAB
1. Generate different types of Signals using MATLAB and DSP

Processor.
2. Plot DFT, FFT, and phase and Magnitude spectrums of various
Digital Signal Processing
Lab signals.
3. Implement low pass and high pass FIR and IIR Filters.
4. Generate DTMF Signals
5
DIGITAL SIGNAL PROCESSING LAB
List of Experiments:
1) To find DFT/IDFT of given DT signal

2) To find frequency response of a given system (Transfer function /Differential equation
form).
3) Generation of sinusoidal waveform /signal based on recursive difference equations
4) . Implementation of FFT of given sequence.
5) Determination of power spectrum of a given signal.
6) . Implementation of LP FIR filters for a given sequence.
7) Implementation of HP IIR filters for a given sequence.
8) Generation of DTMF signals.
9) Implementation of Decimation process.
10) Implementation of Interpolation process.
11) Implementation of I/D sampling rate converters.
12) Impulse response of first order and second order systems.
6
INSTRUCTIONS TO THE STUDENT FJHINSTRUCTIOTO

THE STUDENTTHE STUDENT
1. Students are required to attend all labs.
2. Students will work individually in hardware laboratories and in computer laboratories.
3. While coming to the lab bring the lab manual cum observation book, record etc.
4. Take only the lab manual, calculator (if needed) and a pen or pencil to the work area.
5. Before coming to the lab, prepare the pre-lab questions. Read through the lab
experiment to familiarize yourself with the components and assembly sequence.
6. Utilize 3 hours time properly to perform the experiment (both in software and
hardware) and note down the readings properly. Do the calculations, draw the graph
and take signature from the instructor.
7. If the experiment is not completed in the prescribed time, the pending work has to be
done in the leisure hour or extended hours.
8. You have to submit the completed record book according to the deadlines set up by
your instructor.
9. For practical subjects there shall be a continuous evaluation during the semester for 25
sessional marks and 50 end examination marks.
10. Of the 25 marks for internal, 15 marks shall be awarded for day-to-day work and 10
marks to be awarded by conducting an internal laboratory test.
INDEX
SL.NO. EXPERIMENT NAME PAGENO.
1 Generation of sinusoidal waveform /signal based on recursive difference 9

equations.
7
2 Histogram of White Gaussian Noise and Uniformly Distributed Noise. 13
3 To find DFT/IDFT of given DT signal. 18
4 To find frequency response of a given system (Transfer function 25

/Differential equation form).
5 Obtain the Fourier series coeffients by formula and using FFt and compare for half 30
sinewave
6 Implementation of FFT of given sequence. 34

.
7 Determination of power spectrum of a given signal. 39
8 Implementation of LP FIR filters for a given sequence. 43
9 Implementation of HP IIR filters for a given sequence. 50
10 Generation of Narrow band signal through filtering. 54
11 Generation of DTMF signals. 59

.
12 Implementation of Decimation process. 66
13 Implementation of Interpolation process. 70
14 Implementation of I/D sampling rate converters. 74

15 Impulse response of first order and second order systems. 78
EXPERMENT NO-1
Generating of a sinusoidal signal based on recursive difference equation.
AIM: -
8
To write a MATLAB program for generating a sinusoidal signal based on recursive difference
equation.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
A recursive system is one whose present output depends on present and past values. A non
recursive system output depends only on the present input. A recursive discrete system can be
used to generate a sinusoidal signal. A signal in general phenomena is a physical quantity that
changes with an independent variable and carries information. In standard, it is defined as a
mathematical function of one or more independent variable. There are certain parameters that
are common to all the signals and those parameters are called basic characteristics of the
signal. They are amplitude, frequency, phase, energy, power, and bandwidth. The process of
finding the above six parameters of the given signal is called Signal Analysis. This can be
accomplished by representing the signal in either Time-domain or Frequency-domain. Those
parameters which are difficult to be found in time domain can be determined in frequency
domain.
PROCEDURE:-
 Open MATLAB
 Open new M-file
 Type the program
 Save in current directory
 Compile and Run the program
 For the output see command window\ Figure window
PROGRAM:-
clc;
clear all;
close all;
f_hz=400;
9
fs=8000;
f0=2*pi*f_hz/fs;
a1=2*cos(f0); a2=-1;
b1=sin(f0);
xnm1=1; ynm1=0; ynm2=0;
y1=zeros(1,200);
for i=1:200
y1(i)=b1*xnm1+a1*ynm1+a2*ynm2;
ynm2=ynm1;
ynm1=y1(i);
end
plot(1:length(y1),y1);
xlabel('Time in ms-----> ');
ylabel('Amplitude in volts-----> ');
title(['Sinusoidal Signal of frequency ', num2str(f_hz), ' Hz']);
RESULT:-
Sinusoidal Signal of frequency 400 Hz
7
5
Amplitude in volts----->
-1
0 20 40 60 80 100 120 140 160 180 200
Time in ms----->
Real Life Applications
Trigonometry is especially important in architecture because it allows the architect to

calculate distances and forces related to diagonal elements., for example on bridges and tall
structures, the diagonal has to be strong and accurate to keep the structure
standing. Architecture Space flight relies on calvulations and conversions to polat coordiates.
10
because they help model orbital motions. Polar coordinates express a position on a two-
dimensional plane using an angle from a fixed direction and a distance from a fixed point.
Polar coordinates can be converted to Cartesian coordinates- the coordinate plane that we are
used to seeing, and have been seeing since elementary. Polar coordinates can be converted to
the Cartesian coordinates (x,y) by using sine and cosine functions. By multiplying the polar
coordinates by cosine, the x coordinate can be obtained. By multiplying the polar coordinates
by sine the y coordinate can be found. If a trumpet sounds at 440 Hz, at various amplitudes,
the summation of sine waves or in other words Fourier series will be 440 Hz, 880 Hz, 1, 320
Hz, 1,760Hz.
As we know sound travels in waves and frequencies.A French scientist and mathematician by
the name of Jean Baptiste Fourier proved that any waveform that repeats itself after a period
of time (such as a musical sound) can be expressed as the sum of an infinite set of sine curves.
As we know sound travels in waves and frequencies. Ballisic trajectories rely on sin/cos.
EXERCISE PROGRAMS:-
1. Write program to get Discrete time Sinusoidal Signal?
2. Write program to get Fourier Transform of Sinusoidal Signal?
3. Write program to get Inverse Fourier Transform of Sinusoidal Signal?
4. Write Program for the following Function
Y=exp (-2*∏*f*t)+exp(-8*∏*f*)
Y= ((exp (-1.56∏f)*Sin(2∏f)+cos(2∏f)?
5. Write a mat lab program for generating u(n)-u(n-1)?
6. Write program to get Discrete time co-Sinusoidal Signal?
7. Write program to get Discrete time saw tooth Signal?
8. Write program to get Discrete time triangular Signal?
9. Write program to get addition of two sinusoidal sequences?
10. Write program to get exponential sequence?
11. Write program to get Discrete time Co-Sinusoidal Signal?
12. Write program to get Fourier Transform of Co-Sinusoidal Signal?
13. Write program to get Inverse Fourier Transform of Co-Sinusoidal Signal?
14. Write program to get exponential decaying sequence?
15. Write program to get exponential growing sequence?
17. Write program to get complex exponential sequence?
18. Write a mat lab program for generating sinc function?
19. Write a mat lab program for generating signum function?
20. Write a program to generate negative ramp signal
11
VIVA QUESTIONS:-
1. Define Signal?
2. Define determistic and Random Signal?
3. Define Delta Function?
4. What is Signal Modeling?
5. Define Periodic and a periodic Signal?
6. Represent impulse signal in terms of unit step signal
7. Define Unit step signal
8. Define ramp signal
9. Define continuous time signal
10. Define discrete time signal.
11. Define impulse signal.
12. Define ramp signal
13. Define unit step signal
14. Define exponent ional signal
15. Define sinusoidal signal
16. Define C.T.S
17. Define D.T.S.
18. Compare C.T.S & D.T.S
19. Define Stem, Plot, Plot3,fplot, ezplot, linspace, flyplr, grid,mesh and legend.
20. Draw the C.T.S & D.T.S diagrams.
21. Define signal and signal processing
22. Differentiate digital and analog signals?
23. How the DSP processor will differ from conventional processors?
24. Expand the abbreviation TMS320C 5X/6X
25. What kind of processor is DSP processor?
26. What are the main building
27. blocks of DSP processor?
28. What is the main function of MAC unit?
29. Explain VLIW architecture?
30. What is meant by circular buffer?
EXPERMENT NO:2
HISTOGRAM OF GAUSSIAN NOISE
AIM: -
12
To write a MATLAB program for computation of statistical parameters of Gaussian

distribution and uniform distribution.
SOFTWARE REQURIED:-
1.MATLAB R2010a.
2.Windows XP SP2.
THEORY:-
In random variable or random process if all the samples are clustered at center value,
then it is Gaussian distribution. Instead, if all the samples are equally repeated then it is called
uniform distribution. The histogram is analogous to probability density function. It plots the
sample values with their iteration number. The Gaussian and uniform distributions differ in
statistical averages such as mean, variance, standard deviation, skew, kurtosis, etc.,.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
clc;
clear all;
close all;
x=randn(1,2000);
13
y=rand(1,2000);
hist(x);
xlabel('Sample value');
ylabel('Probability density');
title('Gaussian pdf');
figure;
hist(y);
xlabel('Sample value');
ylabel('Probability density');
title('Uniform pdf');
mx1=mean(x); display(mx1);
my1=mean(y); display(my1);
ux2=var(x); display(ux2);
uy2=var(y); display(uy2);
ux3=skewness(x); display(ux3);
uy3=skewness(y); display(uy3);
ux4=kurtosis(x); display(ux4);
uy4=kurtosis(y); display(uy4);
OUTPUT:-
mx1 =
0.0220
my1 =
0.5042
ux2 =
1.0159
uy2 =
0.0832
ux3 =
-0.0329
uy3 =
-0.0071
ux4 =
3.1472
uy4 =
1.8301
14
Gaussian pdf
600
500
400
Probability density
300
200
100
0
-4 -3 -2 -1 0 1 2 3 4
sample value
Uniform pdf
250
200
Probability density
150
100
50
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Sample value
VIVA QUESTIONS:-
1.What is a noise ?
2.What is gaussion noise?
3.What is correlation ?
4.State Paeseval’s energy theorem for a periodic signal?
5.What is Signum function?
6. How many types of correlation are there?
7. how many types of noises are there?
15
8.what is auto correlation?

9.what are the types of noise?
10.what is ramp signal?
11.what is signal
12.what is white noise?
13.what is periodic signal?
14.what is non periodic signal?
15.what is system bandwidth?
16.what is signal bandwidth?
17.what is causal system?
18.what is non causal system?
19.what is periodic signal?
21.what is signum function?
22.what is sinc function?
23.what are the types of correlation?
24.what is ramp function?
25.what is correlation ?
26.application of periodic function?
27.what is gaussian noise?
28.what is continuous time signal?
29.what is descrete time signal?
30.what is plot?
31.what is stem?
32.what is time delay?
33.what is xalbel?
34.what is y label?
35.what is psd?
36.what is random signal?
37.what is signal?
38.what are the types of signals?
39.what is bandwidth?
40.what is correlation?
41.what are the correlation types?
42.what is y label?
43.what is mean square?
44.what is periodi?
46..what is gaussian noise?
47.what is continuous time signal?
48.what is descrete time signal?
49.what is plot?
50.what is time delay?
16
EXERCISE PROGRAMS
1.Write a MATLAB program to plot the co-sinusoidal signal.

2. Write a MATLAB program to subtract two sinusoidal signals.
3. Write a MATLAB program to subtract and multiply twosinusoidal signals.
4. Write a MATLAB program to right shift the signal to 5 times of the original signal.
5. Write a MATLAB program to left shift the signal to 8 times of the original signal.
6. Write a MATLAB program to add two different signals with 2 <t<5
7. Write a MATLAB program to shift a positive time line signal to negative timeline signal.
8. Write a MATLAB program to subtract co-sinusoidal signals.
9. Write a MATLAB program to subtract two sinusoidal signals
10. Write a MATLAB program to division and multiply two co-sinusoidal signals.
11. Write a MATLAB program to generate time scaling of a sequence.
12. Write a MATLAB program to generate time shifting of a sequence.
13. Write a MATLAB program to generate time folding of a sequence.
14. Write a MATLAB program to generate amplitude scaling of a sequence with amplitude 5.
15. Write a MATLAB program to generate time scaling of a sequence with time 2sec.
16. Write a MATLAB program to add two different signals with 4 <t<8
17. Write a MATLAB program to shift a negative time line signal to positive timeline signal.
18. Write a MATLAB program to subtract sinusoidal signals.
19. Write a MATLAB program to subtract and divide two sinusoidal signals
20. Write a MATLAB program to add and multiply two co-sinusoidal signals.
21. Write a MATLAB program to generate time scaling of a signal.
22. Write a MATLAB program to generate time shifting of a signal.
23. Write a MATLAB program to generate time folding of a signal.
24. Write a MATLAB program to generate amplitude scaling of a sequence with amplitude 10
25. Write a MATLAB program to generate time scaling of a sequence with time 5sec.
Realtime Applications:
• Stream processing
• Block processing
• Vector processing
• Digital processing applications
17
EXPERMENT NO-3
DFT/IDFT OF A SEQUENCE WITHOUT USING THE INBUILT FUNCTIONS
AIM:-
To find the DFT& IDFT of a sequence without using the inbuilt functions.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
Given a sequence of N samples f(n), indexed by n = 0..N-1, the Discrete Fourier

Transform (DFT) is defined as F(k), where k=0..N-1:
F(k) are often called the 'Fourier Coefficients' or 'Harmonics'.
The sequence f(n) can be calculated from F(k) using the Inverse Discrete Fourier Transform
(IDFT):
In general, both f(n) and F(k) are complex.
The DFT is the most important discrete transform, used to perform Fourier analysis in
many practical applications.[1] In digital signal processing, the function is any quantity or
signal that varies over time, such as the pressure of a sound wave, a radio signal, or daily
temperature readings, sampled over a finite time interval (often defined by a window
function). In image processing, the samples can be the values of pixels along a row or column
of a raster image. The DFT is also used to efficiently solve partial differential equations, and
to perform other operations such as convolutions or multiplying large integers.
18
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
%DFT%
clc;
clear all;
close all;
a=input ('enter the input sequence');
N=length(a);
disp('length of input sequence is ');
N
for k=1:N;
x(k)=0;
for i=1:N;
x(k)=x(k)+a(i)*exp((-j*pi*2/N)*((i-1)*(k-1)));
end;
end;
k=1:N;
disp('the output is');
x(k)
subplot(2,1,1);
stem(k,abs(x(k)));
grid;
xlabel ('discrete frequency');
ylabel('magnitude');
title('magnitude response of dft');
subplot(2,1,2);
stem(angle(x(k))*180/(pi));
grid;
xlabel('discrete frequency');
ylabel('phase angle');
title('phase response of dft');
19
%IDFT%
clc;
clear all;
close all;
a=input('enter the input sequence');
disp('the length of input sequence is');
N=length(a);
N
for n=1:N;
x(n)=0;
for k=1:N;
x(n)=x(n)+a(k)*exp((j*pi*2*(n-1)*(k-1)/N));
end;
end;
n=1:N;
x=1/N*x(n);
disp('the output is');
x(n)
stem(n,abs(x));
grid;
xlabel('discrete time');
title('magnitude response of the idft');
grid;
20
OUTPUT:-
enter the input sequence[1 2 3 4]
length of input sequence is
N=
the output is
ans =
10.0000 -2.0000 + 2.0000i -2.0000 - 0.0000i -2.0000 - 2.0000i
enter the input sequence[10 -2+2j -2 -2-2j]

the length of input sequence is
N= 4
the output is
1.0000 2.0000 + 0.0000i 3.0000 - 0.0000i 4.0000 - 0.0000i

magnitude response of the idft
4
3.5
2.5
magnitude
1.5
0.5
0
1 1.5 2 2.5 3 3.5 4
discrete time
21
RESULT:-
DFT&IDFT of a given discrete time signal are executed using mat lab software.
EXERCISE PROGRAM:-
1. Write a matlab program to find the circular convolution of two sequences?

2. Write a matlab program to find the circular convolution of
x1(n)={2,3,-1,2};x2(n)={-1,2,-1,2}?
x1(n)={1,-1,2,3}; x2(n)={2,0,1,1}?
x1(n)={1,1,-1,2}; x2(n)={0,1,2,3}?
5. Write a matlab program to find the DFT of x(n) ={1 1 1 1 0 0 0 0}?
6. Write a matlab program to find the DFT of x(n) ={1 2 1 2}?
7. Write a matlab program to find the DFT of x(n) ={1 0 -1 0}?
8. Write a matlab program to find the IDFT of X(k) ={1,1,-2j,-1,1+2j }?
9. Write a matlab program to find the IDFT of X(k) ={1 0 1 0}?
10. Write a matlab program to find the to compare circular and linear convolution of two
sequences?
13. Write a matlab program to find the IDFT of X(k) ={0 -1 0 -1}?
x1(n)={1,0,-1,0}; x2(n)={0,1,0,1}?
x1(n)={1,1,-1,-1}; x2(n)={1 2 3 4}?
18. Write a matlab program to find the DFT of x(n) ={1 0 -1 0}?
19. Write a matlab program to find the IDFT of X(k) ={1,1+2j,-1,1-2j }?
20. Write a matlab program to find the IDFT of X(k) ={0 0 1 1}?
22
VIVA QUESTIONS:-
1. Define Symmetric and Anti-Symmetric Signals?

2. Define Continuous and Discrete Time Signals?
3. What are the Different types of representation of discrete time signals?
4. What are the Different types of Operation performed on signals?
5. Define DFT
6. How DFT can be calculated in matrix form?
7. Write the standard formula of DFT
8. What is the relation between DFT and DFS
9. What is the relation between DFT and Z-Transform
10. What is the relation between DFT and Laplasetransform.
11. Whether DFT is a linear transform?
12. What is the difference between circular convolution & linear convolution?
13. Can you implement linear convolution using circular convolution?
14. How FFT algorithms are classified?
15. How to calculate output of DFT using MATLAB?
16. Where DFT is used?
17. What is the difference between DFT and IDFT?
18. How to compute maximum length N for a circular convolution using DFT and IDFT.
19. Explain the function of twiddle factor?
20. Give the practical application DFT & IDFT?
21. Explain the significance of convolution.
22. Define linear convolution.
23. Why linear convolution is called as a periodic convolution?
24. Why zero padding is used in linear convolution?
25. What are the four steps to find linear convolution?
26. What is the length of the resultant sequence in linear convolution?
27. How linear convolution will be used in calculation of LTI system response?
28. List few applications of linear convolution in LTI system design.
29. Give the properties of linear convolution.
30. How the linear convolution will be used to calculate the DFT of a signal
The Discrete Fourier Transform (DFT) is one of the most important tools in
Digital Signal Processing. First, the DFT can calculate a signal's frequency spectrum. This is a
direct examination of information encoded in the frequency, phase, and amplitude of the
component sinusoids. For example, human speech and hearing use signals with this type of
encoding. Second, the DFT can find a system's frequency response from the system's impulse
response, and vice versa. This allows systems to be analyzed in the frequency domain, just as
convolution allows systems to be analyzed in the time domain. Third, the DFT can be used as
an intermediate step in more elaborate signal processing techniques. The classic example of
this is FFT convolution, an algorithm for convolving signals that is hundreds of times faster
than conventional methods.
23
EXPERMENT NO-4
TO FIND FREQUENCY RESPONSE OF GIVEN SYSTEM GIVEN IN TRANSFER

FUNCTION/ DIFFERENTIAL EQUATION FORM.
AIM: -
To study frequency response of second order system using MATLAB.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
Second order systems are the systems or networks which contain two or
more storage elements and have describing equations that are second order
differential equations.
The frequency response of second order filters is characterised by three

filter parameters: the gain k, the corner frequency and the quality factor Q.
24
A second order filter is a circuit that has a transfer function of the form:
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
%frequency response of differential equation %
clc ;
clear all;
b=[1,4];
a=[1,-5];
w=-2*pi:pi/8:2*pi;
[h]=freqz(b,a,w);
subplot(2,1,1);
stem(w,abs(h));
xlabel('freq/w');
grid;
title('magnitude response of diffenrtial equataion');
subplot(2,1,2);
stem(w,angle(h));
25
xlabel('freq/w');
ylabel('phase in rad');
grid;
title('phase response of diffenrtial equataion');
OUTPUT:-
26
magnitude response of diffenrtial equataion

1.5
magnitude 1
0.5
0
-8 -6 -4 -2 0 2 4 6 8
freq/w
phase response of diffenrtial equataion
4
2
phase in rad
-2
-4
-8 -6 -4 -2 0 2 4 6 8
freq/w
RESULT:-
Hence the frequency response is executed by using MATLAB.
EXERCISE PROGRAM:-
27
1. Write a matlab program to find the frequency response of the following difference
equation y(n)-7y(n-1)+9y(n-2)=3x(n)-2x(n-1)?
equation 3y(n)+5y(n-1)=9x(n)?
equation 9 y(n)-2y(n-1)+7y(n-2)-3y(n-3)=6x(n)+x(n-1)?
equation 8y(n)+6y(n-1)=4x(n)+2x(n-1)?
equation 3y(n)-8y(n-1)+9y(n-2)=9x(n)+5x(n-1) ?
equation 6y(n)-5y(n-1)=9x(n)+5x(n-1) -7x(n-2)?
equation 9y(n)-8y(n-1)+2y(n-2)=9x(n)-3x(n-1) ?
equation 2y(n)-8y(n-1)=9x(n)+5x(n-1) ?
equation 9y(n)-8y(n-1)+9y(n-2)=9x(n)+5x(n-1)-x(n-2) ?
equation 3y(n)-8y(n-1)=7x(n)-3x(n-1) ?
equation 11y(n-2)-5y(n-4)=6x(n)+5x(n-3) -9x(n-7)?
equation 9y(n-1)-8y(n-5)+2y(n-7)=5x(n-1)-1x(n-3) ?
equation y(n-3)-8y(n-7)=5x(n-1)+3x(n-5) ?
equation 9y(n-2)-8y(n-4)+9y(n-6)=9x(n-1)+9x(n-3)-x(n-5) ?
equation 3y(n+2)-8y(n+4)=7x(n-5)-3x(n-9) ?
equation 11y(n-2)-5y(n-4)=6x(n-2)+5x(n-3) -9x(n-7)?
equation 9y(n-1)-8y(n-4)+2y(n-7)=5x(n-1)-1x(n-3) ?
equation y(n-2)-8y(n-7)=5x(n-1)+3x(n-5) ?
equation 9y(n-2)-8y(n-4)+9y(n-6)=9x(n+1)+9x(n-3)-x(n-5) ?
equation 3y(n+2)-8y(n+4)=7x(n+5)-3x(n-9) ?
VIVA QUESTIONS:-
28
1. What is the commend to find phase angle?

2. What is the commend to find frequency response?
3. What is transition band?
4. What is the formula for Z-transform?
5. What is the relationship b/w impulse response& frequency response?
7. What is the relation between DFT and Z-Transform.
8. Give any example of a FIR interpolator?
9. Define impulse response.
10. Give me one example for impulse response.
11. Write the Formula for impulse response.
12. What are major role in order & length?
13. Define power spectral Density?
14. What is the need for spectral estimation?
15. Determine the power spectrum density?
16. What is the relation between auto correlation & spectral density?
17. Give the estimation of auto correlation function & power density for random Signals?
18. Explain power spectrum estimation using the Bartlett window?
19. Give the formula for PSD.
20. What is filter?
21. Give mathematical definition of circular convolution.
22. Why circular convolution is called as periodic convolution?
23. Difference between linear convolution and circular convolution.
24. Explain the circular shift.
25. How circular convolution is used to calculate the Z
26. transform of a signal?
27. List few Applications of circular convolution
28. What are the different methods used to calculate circular convolution?
29. Explain properties of circular convolution?
30. Explain modulo N operation.
Realtime Applications
In electronics this stimulus would be an input signal. In the audible range it is usually
referred to in connection with electronic amplifiers, microphones and loudspeakers. Radio
spectrum frequency response can refer to measurements of coaxial cable, twisted-pair cable,
video switching equipment, wireless communications devices, and antenna systems.
Infrasonic frequency response measurements include earthquakes and electroencephalography
(brain waves).Frequency response requirements differ depending on the application. In high
fidelity audio, an amplifier requires a frequency response of at least 20–20,000 Hz, with a
tolerance as tight as ±0.1 dB in the mid-range frequencies around 1000 Hz, however, in
telephony, a frequency response of 400–4,000 Hz, with a tolerance of ±1 dB is sufficient for
intelligibility of speech.
Experiment 5
29
Fourier Series Coefficients
Aim: To write a MATLAB program for obtaining Fourier series coefficients by formula and
using Fast Fourier Transform algorithm, to verify the same in Code composer studio, and to
plot all the corresponding graphs.
Software Required: MATLAB 7.0 with signal processing, communications, control systems,
and statistics tool boxes.
Theory:
The fast Fourier transform represents a given signal in frequency domain by reducing the
number of complex multiplication to (N/2)log2N and complex additions to Nlog2N in DFT
algorithm. The time taken by fft to find Fourier transform is lesser than dft. It is implemented
in decimation in time (DIT FFT) and decimation in frequency (DIF FFT) by using radix based
butterfly structures.
MATLAB Program:
lc;
close all;
clear acll;
n =0:0.075:1;
x=cos(2*pi*0.5*n);
N=length(x);
for k=1:N
y=0;
for n=1:N
y=y+x(n).*exp(-j*2*pi*(k-1)*(n-1)/N);
end
C1(k)=y/N ;
end
display(C1);
C2=ifft(x);
m1=abs(C1);
m2=abs(C2);
display(C1);
subplot(3,1,1);
stem(x);
xlabel('Discrete Time')
ylabel('Amplitude')
title('Time Domain Signal')
subplot(3,1,2)
stem(m1);
xlabel('k')
ylabel('ck')
title('Fourier Series Coefficients using Formula');
30
subplot(3,1,3);
stem(m2);
xlabel('k')
ylabel('Ck')
title('Fourier Series Coefficients using FFT');
Output:
C1 =
Columns 1 through 4
0.0238 0.0893 - 0.4314i 0.0747 - 0.1591i 0.0727 - 0.0923i
Columns 5 through 8
0.0721 - 0.0579i 0.0718 - 0.0348i 0.0717 - 0.0164i 0.0717 + 0.0000i
Columns 9 through 12
0.0717 + 0.0164i 0.0718 + 0.0348i 0.0721 + 0.0579i 0.0727 + 0.0923i
0.0747 + 0.1591i 0.0893 + 0.4314i
C1 =
Columns 1 through 4
0.0238 0.0893 - 0.4314i 0.0747 - 0.1591i 0.0727 - 0.0923i
Columns 5 through 8
0.0721 - 0.0579i 0.0718 - 0.0348i 0.0717 - 0.0164i 0.0717 + 0.0000i
0.0717 + 0.0164i 0.0718 + 0.0348i 0.0721 + 0.0579i 0.0727 + 0.0923i
0.0747 + 0.1591i 0.0893 + 0.4314i
Time Domain Signal

1
Amplitude
-1
0 2 4 6 8 10 12 14
Discrete Time
Fourier Series Coefficients using Formula
0.5
ck
0
0 2 4 6 8 10 12 14
k
Fourier Series Coefficients using FFT
0.5
Ck
0
0 2 4 6 8 10 12 14
k
Result: The MATLAB program for obtaining Fourier series coefficients by formula and
using Fast Fourier Transform algorithm has been written, the same is verified.
31
EXERCISE PROGRAM:-
1. Determine the Fourier series coefficients ak for x1(t) shown below. x1(t)= x1(t + 10) ?
2. Determine the Fourier series coefficients bk for x2(t) shown below. x2(t)= x2(t + 10) ?
3. Determine the Fourier series coefficients ck for x3(t) shown below. x3(t)= x3(t + 10)?
4. Determine the Fourier series coefficients for a cos wave by using formula method?
5. Determine the Fourier series coefficients for a cosine wave using FFT ?
6. Determine the Fourier series coefficients of the following signal
x(t) = 1+cos(2πt)
x(t) = [1+cos(2πt)] [sin(10πt+π/6)] ?
7. Find the Fourier series coefficients for (periodic extension of)

f(t) = 1 − t 2 , t ∈ [−1, 1] ?
8. . Find the cosine Fourier series coefficients for (periodic extension of)
f(t) = ½ 1, t ∈ [0, 1); 0, t ∈ [1, 4). ?
9. Find the Fourier series for (periodic extension of) f(t) = ½ 1, t ∈ [0, 2); −1, t ∈ [2, 4). ?
10. Determine the Fourier series coefficients for a given even signal f(t)=2cosπt ?
11. Determine the Fourier series coefficients for a odd function l f(t) = sin t ?
12. Determine the Fourier series coefficients for a Saw-tooth Wave ?
13. Consider the function below {

1-x 0 < x ≤1; 0 1<x≤2 }
And find out fourier series coefficients?
14. Determine the Fourier series coefficients ck for xn(t) shown below. xn(t)= xn(t + 50)?
x(t) = 1+cos(2πt)
x(t) = [1+sin(2πt)] [cos(10πt+π/6)] ?
16. Find the Fourier series coefficients for (periodic extension of)
f(t) = 1+t 2 , t ∈ [1, 2] ?
17. Find the cosine Fourier series coefficients for (periodic extension of)
f(t) = 2/31, t ∈ [0, 1); 0, t ∈ [1, 4). ?
18. Determine the Fourier series coefficients ck for xn(t) shown below.
xn(t)= xn(t+ 350)?

x(t) = 1+cos(5πt)
x(t) = [1+sin(5πt)] [cos(14πt+π/9)] ?
x(t) = 1+cos(2πt)
x(t) = [1+tan(2πt)] [cot(10πt+π/6)] ?
VIVA QUESTIONS:-
32
1. What is ment by fourier series?

2. Difine Fourier transforms?
3. Give the differences between fourier series and fourier transforms?
4. What is the formula for Fourier series?
5. What is the formula for Fourier transforms?
6. What is the formula for Z-transforms?
7. What is the formula for Laplace transforms?
8. Give the differences between Z-transform and laplace transforms?
9. Applications of Fourier Series and Fourier Transforms?
10. What are the Fourier Coefficents?
11. Do exponential fourier series also have fourier coefficients to be evaluated.?
12. What is the polar form of the fourier series?
14. What is a line spectrum?
15. Fourier series is not true in case of discrete time signals.
16. What is the disadvantage of exponential Fourier series?
17. Fourier series uses which domain representation of signals?.
18. How does Fourier series make it easier to represent periodic signals?
19. If transfer function of a system is H(z) = 6 + z-1 + z-2 then system is ?
23. Explain modulo N operation.
Fourier series are of great importance in both theoretical and applied mathematics. For
orthonormal families of complexvalued functions {φn}, Fourier Series are sums of the φn that
can approximate periodic, complexvalued functions with arbitrary precision. This paper will
focus on the Fourier Series of the complex exponentials. Of the many possible methods of
estimating complexvalued functions, Fourier series are especially attractive because uniform
convergence of the Fourier series (as more terms are added) is guaranteed for continuous,
bounded functions. Furthermore, the Fourier coefficients are designed to minimize the square
of the error from the actual function. Finally, complex exponentials are relatively simple to
deal with and ubiquitous in physical phenomena. This paper first defines generalized Fourier
series, with an emphasis on the series with complex exponentials. Then, important properties.
EXPERMENT NO-6
33
IMPLEMENTATION OF FFT OF GIVEN SEQUENCE
AIM: -
Implementation of FFT of given sequence.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
A fast Fourier transform (FFT) is an algorithm to compute the discrete Fourier

transform (DFT) and its inverse. Fourier analysis converts time (or space) to frequency and
vice versa; an FFT rapidly computes such transformations by factorizing the DFT matrix into
a product of sparse (mostly zero) factors.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
%fft%
clc;
clear all;
close all;
xn=input('enter the input sequence');
N=input('enter the number of samples');
n=0:1:N-1;
xk=fft(xn,N);
k=0:1:N-1;
subplot(2,1,1);
stem(k,abs(xk));
xlabel('frq/w');
title('magnitude response of fft');
subplot(2,1,2);
stem(k,angle(xk));
xlabel('frq/w');
ylabel('phase');
34
title('phase response of fft');
OUTPUT:-
enter the input sequence[1 2 3 4]

enter the number of samples8
magnitude response of fft
10
magnitude
0
0 1 2 3 4 5 6 7
frq/w
phase response of fft
4
2
phase
-2
-4
0 1 2 3 4 5 6 7
frq/w
RESULT:-
Hence the FFT of a given input sequence is performed & executed by using
MATLAB.
EXERCISE PROGRAM:-
35
1. Write a matlab program to find the cross correlation using FFT?

5. Write a matlab program to find the IDFT of X(k) ={1,1+j,-2j,1+2j ,-j, +j}?
6. Write a matlab program to find the IDFT of X(k) ={1,0,-2j,-1,+2j,-7j }?
7. Write a matlab program to find the IDFT of X(k) ={1,1,-2j,-1,1+2j }?
8. Write a matlab program to find the IDFT of X(k) ={1,1+j,-2j,-1+-j,1+2j }?
11. Write a matlab program to find the IDFT of X(k) ={1+j,0,1-2j,-1,1+2j,1-7j }?
12. Write a matlab program to find the IDFT of X(k) ={1,2+5j,2-2j,-1,5+2j }?
13. Write a matlab program to find the IDFT of X(k) ={5+6j,1+6j,5-2j,-1+9j,1-2j }?
16. Write a matlab program to find the IDFT of X(k) ={2+j,0,7+2j,-1-5j,9+2j,6-7j }?
17. Write a matlab program to find the IDFT of X(k) ={1,2+5j,2-2j,-1,5+2j }?
18. Write a matlab program to find the IDFT of X(k) ={6j,1+6j,5-2j,-1-9j,1-2j }?
19. Write a matlab program to find the DFT of x(n) ={1 0 1 0 1 0 0 0 0 }?
36
VIVA QUESTIONS:-
1. Whether linear convolution equation is a difference equation?

2. Whether DFT is a linear transform?
3. What is the difference between circular convolution & linear convolution?
4. Can you implement linear convolution using circular convolution?
5. How FFT algorithms are classified?
9. How to compute maximum length N for a circular convolution using DFT and IDFT.
10. Explain the function of twiddle factor?
11. Give the practical application DFT & IDFT?
12. Explain the role of DFT & IDFT when the signal converted from the time domain to
frequency domain?
13. Differentiate between time variant and time invariant system. If x 1(n)={1,2,3,4}
and x 2(n)={1,2,3} Find the convolution using tabular representation.
14. Draw all elementary standard discrete time signals.
15. Differentiate between causal and Non causal system.
16. If x 1(n)={1,2,3,4} and x 2(n)={5,6,7,8} Find the circular representation for the above
sequences.
17. How can you compute Fourier transform form Z-transform ?
18. FFT is in complex domain how to use it in real life signals optimally?
19. What is the difference between FFT and IFFT?
20. Explain using convolution the effects of taking an FFT of a sample with no windowing
.
21. What is a filter?
22. Differentiate analog filter and digital filter
23. Define FIR filter
24. What are the differences between recursive and non recursive systems?
25. List a few Applications of FIR filters.
26. Explain advantages of FIR filters over IIR filters.
27. Explain limitations of FIR filters.
28. What is the different method to design FIR filters?
29. Explain different window functions.
30. Differentiate rectangular, triangular and Kaiser windows.
37
Real Time Applications
The Discrete Fourier Transform (DFT) is one of the most important tools in Digital Signal
Processing. First, the DFT can calculate a signal's frequency spectrum. This is a direct
examination of information encoded in the frequency, phase, and amplitude of the component
sinusoids. For example, human speech and hearing use signals with this type of encoding.
Second, the DFT can find a system's frequency response from the system's impulse response,
and vice versa. This allows systems to be analyzed in the frequency domain, just as
convolution allows systems to be analyzed in the time domain. Third, the DFT can be used as
an intermediate step in more elaborate signal processing techniques. The classic example of
this is FFT convolution, an algorithm for convolving signals that is hundreds of times faster
than conventional methods.
38
EXPERMENT NO-7
DETERMINATION OF POWER SPECTRUM
AIM: - To obtain power spectrum of given signal using MATLAB.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
THEORY:-
The power spectrum of a time-series x(t) describes how the variance of the data x(t) is
distributed over the frequency components into which x(t) may be decomposed. This
distribution of the variance may be described either by a measure µ or by a statistical
cumulative distribution function S(f) = the power contributed by frequencies from 0 up to f.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
% power spectrum %
clc;
clear all;
close all;
f1=input('enter the first frequencey f1=');
f2=input('enter the second frequencey f2=');
fs=input('enter the sampling frequencey fs=');
t=0:1/fs:1;
x=2*sin(2*pi*f1*t)+3*sin(2*pi*f2*t)+rand(size(t));
psd1=abs(fft(x).^2);
subplot(2,1,1);
plot(t*fs,10*log(psd1));
xlabel('frequency');
title('psd using square magnitude method');
psd2=abs(fft(xcorr(x),length(t)));
subplot(2,1,2);
plot(t*fs,10*log(psd2));
39
xlabel('frequency');
title('psd using auto corelation method');
OUTPUT:-
enter the first frequencey f1=200

enter the second frequencey f2=400
enter the sampling frequencey fs=1000
psd using square magnitude method

150
magnitude
100
50
0
0 100 200 300 400 500 600 700 800 900 1000
frequency
psd using auto corelation method
140
magnitude
120
100
80
0 100 200 300 400 500 600 700 800 900 1000
frequency
RESULT:-
Hence the power spectral density is performed & executed by using MATLAB.
40
EXERCISE PROGRAM:-
1. Write a matlab program for power spectrum estimate using Welch method?
2. Write a matlab program to plot the frequency response of a first order system?
3. Write a matlab program to plot the frequency response of the system?
4. Write a matlab program to generate the periodic sequence?
5. Write a matlab program to generate the aperiodic sequence?
6. Write a matlab program to demonstrate the property of digital frequency?
7. Write a matlab program to illustrate the concept of aliasing?
8. Write a matlab program to plot magnitude and phase response of first order lowpass
filter?
9. Write a matlab program to plot magnitude and phase response of first order highpass
filter?
10. Write a matlab program to plot magnitude and phase response of second order
bandpass filter?
bandstop filter?
12. Write a matlab program to plot magnitude and phase response of second order lowpass
filter?
highpass filter?
14. Write a matlab program to plot magnitude and phase response of first order bandpass
filter?
15. Write a matlab program to plot magnitude and phase response of first order bandstop
filter?
16. Write a matlab program to plot the frequency response of a first order system?
17. Write a matlab program to plot the frequency response of the system?
18. Write a matlab program to generate the periodic sinusoidal sequence?
19. Write a matlab program to generate the a periodic sinusoidal sequence?
20. Write a matlab program to demonstrate the property of digital signal?
VIVA QUESTIONS:-
1. Give the expressions for finding the Average power of a signal/sequence?

2. Give the expressions for finding the energy of a signal/sequence?
3. What is power spectrum?
4. Why there are two peaks in the magnitude spectrum of sine wave?
5. What is spectrogram? Which built in function is used to solve a given difference
equation?
6. What is frequency response? Give equation for first order system andsecond order
system?
7. .What is an LTI system?
8. What is steady state response?
9. Suppose we have a system with transfer function H(z) = 1 / ((z –1)*(z –0.9)). Is the
system stable or unstable?
41
10. What is Auto Regressive Model? How is the order of auto regressive model is
decided?
11. Differentiate between linear and circular convolution.
12. Determine the unit step response of the linear time invariant system with impulse
response h(n)=a nu(n) a<1&-a<1
13. Determine the range of values of the parameter a for which linear time invariant
system with impulse response h(n)=a u(n) is stable.
14. Consider the special case of a finite duration sequence given as X(n)={2 4 0 3},
resolve the sequence x(n) into a sum of weighted sequences.
15. Describe impulse response of a function?
16. Where to use command filter or impz, and what is the difference between these two?
20. Explain the function of twiddle factor? List some a dvantages of digital filters over
analog filters.
21. Write some differences between FIR and IIR filters
22. What are the different methods to design IIR filters.
23. Why IIR filters are not reliable What are different applications of IIR filters
24. What are advantages of IIR filters.
25. What are disadvantages of IIR filters.
26. Differentiate Butterworth and Chebyshev approximations.
27. What is meant by impulse response.
28. What is the importance of impulse response to calculate the o/p response of the filter?
29. Describe impulse response of a function?
30. Where to use command filter or impz, and what is the difference between these two?
Applications:
It and the cumulative spectral density are very useful for identifying periodic components in
time series. For instance, say you want to study weather cycles or even something like daily
sales figures data for cycles. Cycles will show up pretty nicely as spikes in the PSD or CSD.
You can also use the PSD or CSD to help diagnose whether your model fits well. If the spikes
at identified frequencies are gone that’s a sign that you have identified the periodic
components.
42
EXPERMENT NO-8
IMPLEMENTATION OF LP FIR FILTERS FOR GIVEN SEQUENCE/SIGNAL.
AIM: -
Implementation of Low Pass FIR filters for given sequence.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2.Windows XP SP2.
THEORY:-
A Finite Impulse Response (FIR) filter is a discrete linear time-invariant system

whose output is based on the weighted summation of a finite number of past inputs. An FIR
transversal filter structure can be obtained directly from the equation for discrete-time
convolution.
In this equation, x(k) and y(n) represent the input to and output from the filter at time n. h(n-
k) is the transversal filter coefficients at time n. These coefficients are generated by using
FDS (Filter Design Software or Digital filter design package).
FIR – filter is a finite impulse response filter. Order of the filter should be specified.
Infinite response is truncated to get finite impulse response. Placing a window of finite length
does this. Types of windows available are Rectangular, Barlett, Hamming, Hanning,
Blackmann window etc. This FIR filter is an all zero filter.
PROCEDURE:-
 Open MATLAB
 Open new M-file
43
PROGRAM:-
% LOW PASS FILTER %
clc;
clear all;
close all;
rp=input('enter the pass band ripple:rp=');
rs=input('enter the stop band ripple:rs=');
fp=input('enter the pass band freq :fp=');
fs=input('enter the stop band freq:fp=');
f=input('enter the sampling freq:f=');
wp=2*fp/f;
ws=2*fs/f;
num=-20*log10(sqrt(rp*rs))-13;
den=14.6*(fs-fp)/f;
n=ceil(num/den);
n1=n+1;
if(rem(n,2)~=0);
n1=n;
n=n-1;
end;
c=input('enter the type of window function 1.rectangular 2.trangular 3.kaiser:n=');
if(c==1);
y=rectwin(n1);
disp('rectangular window filter response');
end;
if(c==2);
y=triang(n1);
disp('triangular window filter response');
end;
if(c==3);
y=kaiser(n1);
disp('kaiser window filter response');
end;
% low pass filter %
b=fir1(n,wp,y);
[n,o]=freqz(b,1,256);
m=20*log10(abs(n));
plot(o/pi,m);
xlabel('normalised freq output i');
ylabel('gain in db');
title('low pass filter');
44
OUTPUT:-
enter the pass band ripple:rp=0.02

enter the stop band ripple:rs=0.01
enter the pass band freq :fp=1000
enter the stop band freq:fp=1500
enter the sampling freq:f=10000
enter the type of window function 1.rectangular 2.trangular 3.kaiser:n=1
rectangular window filter response
low pass filter

10
-10
-20
-30
gain in db
-40
-50
-60
-70
-80
-90
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
normalised freq output i

triangular window filter response
45
low pass filter

0
-5
-10
-15
gain in db
-20
-25
-30
-35
-40
-45
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

kaiser window filter response
46
low pass filter

10
-10
-20
-30
gain in db
-40
-50
-60
-70
-80
-90
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
RESULT:-
The implementation of Butterworth Low passes FIR Filters Completed.
47
EXERCISE PROGRAM:-
1. Write a matlab program to design FIR low pass filter using hamming/hanning
window?
2. Write a matlab program to design FIR low pass filter using hamming& blackman
window?
3. Write a matlab program to design FIR low pass filter with cutoff frequency .5pi using
frequency sampling method?
4. Write a matlab program to plot the frequency response of low pass filter using Kaiser
window for different values of beeta?
5. Write a matlab program to design FIR low pass filter for the given specifications using
Kaiser window?
6. Write a matlab program to design a 25-tap Hilbert transformer using Bartlett and
hamming windows and plot their frequency response?
7. Write a matlab program to design a 25-tap lowpass filter with cutoff frequency .5pi
radians using rectangular and hamming windows and plot their frequency response?
8. Write a matlab program to design a 25-tap highpass filter with cutoff frequency .5pi
radians using rectangular and Blackman windows and plot their frequency response?
9. Write a matlab program to design a 25-tap bandpass filter with cutoff frequency .25pi
& .75pi radians using rectangular and hamming windows and plot their frequency
response?
10. Write a matlab program to design a 25-tap bandstop filter with cutoff frequency .25pi
& .75pi radians using rectangular and hamming windows and plot their frequency
response?
11. Write a matlab program to design a 25-tap Hilbert transformer using Bartlett window
and plot the frequency response?
12. Write a matlab program to design a 25-tap Hilbert transformer using hamming window
and plot the frequency response?
13. Write a matlab program to design FIR low pass filter using hamming window?
14. Write a matlab program to design FIR low pass filter using blackman window?
15. Write a matlab program to design FIR high pass filter using blackman window?
16. Write a matlab program to design FIR high pass filter with cutoff frequency .5pi using
frequency sampling method?
17. Write a matlab program to plot the frequency response of high pass filter using Kaiser
window for different values of beeta?
18. Write a matlab program to design FIR high pass filter for the given specifications
using Kaiser window?
19. Write a matlab program to design a 25-tap Hilbert transformer using Bartlett and
hamming windows and plot their frequency response?
20. Write a matlab program to design a 25-tap high pass filter with cutoff frequency .5pi
radians using rectangular and hamming windows and plot their frequency response?
48
VIVA QUESTIONS:-
1. Give the expression for finding the magnitude & phase response of FIR filter?
2. Compare different windows & their characteristics?
3. How FIR filters are designed using rectangular window?
4. What is Gibbs phenomenon? How it can be reduced?
5. Compare FIR&IIR filters?
6. What is "decimation-in-time" versus "decimation-in-frequency"?
7. What is "bit reversal"?
8. What do you mean by phase spectrum and magnitude spectrum/ give comparison?
9. How do you reduce spectral leakage?
10. What do you mean by spectral resolution?
11. What is FIR and IIR filter define, and distinguish between these two?
12. What is window method? How you will design an FIR filter using window method?
13. What are low-pass and band-pass filter and what is the difference between these two?
14. What is the matlab command for Hamming window? Explain.
15. What do you mean by built in function ‘abs’ and where it is used?
16. Explain how the FIR filter are stable?
17. Why is the impulse response "finite"?
18. What does "FIR" mean?
19. What are the advantages of FIR Filters (compared to IIR filters)?
20. What are the disadvantages of FIR Filters (compared to IIR filters)?
21. Define Gibb’s phenomena
22. What is meant by ringing effect?
23. Why do we need to represent a signal in frequency domain?
24. Why Fourier Series converges only for periodic signals?
25. How the ringing effect can be rectified?
26. What are the different forms of Fourier series?
27. What are the limitations of Fourier series?
28. Write few applications of Fourier series.
29. Explain the MATLAB functions ‘tic’ and ‘toc’
30. Explain the MATLAB function ‘legend’
Real time Applications
The calculation of the digital filter function (DSP) is done by the CPU. For that a FIR
algorithm of Infineon's TriLib is used. The coefficients for the filter can be calculated with an
appropriate software program. For this application note the program "ScopeFIR" was used.
Finally the GPTA of the TC1775 generates a PWM signal with the digital filtered.
49
EXPERMENT NO 9
IMPLEMENTATION OF HP IIR FILTERS FOR A GIVEN SEQUENCE/SIGNAL.
AIM: To implement the analog & digital High Pass filter.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
Infinite impulse response (IIR) is a property applying to many linear time-invariant

systems. Common examples of linear time-invariant systems are most electronic and digital
filters. Systems with this property are known as IIR systems or IIR filters, and are
distinguished by having an impulse response which does not become exactly zero past a
certain point, but continues indefinitely. This is in contrast to a finite impulse response in
which the impulse response h(t) does become exactly zero at times t > T for some finite T,
thus being of finite duration.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
clc;
clear all;
close all;
disp('enter the sepecifications of iir filter');
rp=input('enter the pass band ripple:rp=');
rs=input('enter the stop band ripple:rs=');
wp=input('enter the pass band freq:wp=');
ws=input('enter the stop band freq:ws=');
fs=input('enter the sampling freq fs=');
w1=2*wp/fs;
w2=2*ws/fs;
[N,wc]=buttord(w1,w2,rp,rs,'s');
disp('freq resp of iir high pass filter is:');
[b,a]=butter(N,wc,'high','s');
w=0:0.001:pi;
50
[n,omega]=freqs(b,a,w);
m=20*log10(abs(n));
subplot(2,1,1);
plot(omega/pi,m);
xlabel('normalised freq');
ylabel('gain');
title('magnitude response');
subplot(2,1,2);
plot(angle(n));
xlabel('normalised freq');
ylabel('phase');
title('phase response');
OUTPUT:-
enter the sepecifications of iir filter

enter the stop band ripple:rs=60
enter the pass band freq:wp=1500
enter the stop band freq:ws=3000
enter the sampling freq fs=7000
freq resp of iir high pass filter is:
51
magnitude response
500
0
gain
-500
-1000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
normalised freq
phase response
4
2
phase
-2
-4
0 500 1000 1500 2000 2500 3000 3500
normalised freq
RESULT:-
Thus IIR lowpass filter is designed using MATLAB.
EXERCISE PROGRAM:-
1. Write a matlab program to generate IIR chebyshev analog highpass filter?

2. Write a matlab program to design a Butterworth highpass filter for the specifications?
3. Write a matlab program to design a Butterworth bandpass filter for the specifications?
5. Write a matlab program to design a Butterworth bandreject filter for the
specifications?
6. Write a matlab program to design a chebyshev -II high pass filter for the
specifications?
7. Write a matlab program to design a chebyshev -I bandreject filter for the
specifications?
8. Write a matlab program to design a chebyshev -II bandreject filter for the
specifications?
9. Write a matlab program to design a chebyshev -I highpass filter for the specifications?
10. Write a matlab program to design a chebyshev -II highpass filter for the specifications?
11. Write a matlab program to design a chebyshev -I high pass filter for the specifications?
12. Write a matlab program to generate IIR chebyshev analog lowpass filter?
13. Write a matlab program to design a Butterworth lowpass filter for the specifications?
14. Write a matlab program to design a Butterworth bandpass filter for the specifications?
52
16. Write a matlab program to design a Butterworth bandreject filter for the
specifications?
17. Write a matlab program to design a chebyshev -I lowpass filter for the specifications?
18. Write a matlab program to design a chebyshev -II lowpass filter for the specifications?
19. Write a matlab program to design a chebyshev -I bandpass filter for the specifications?
20. Write a matlab program to design a chebyshev -II bandpass filter for the
specifications?
VIVA QUESTIONS:-
1. What are the steps in designing the IIR filters?

2. State the disadvantages of impulse invariant transformation?
3. Why impulse invariant transformation is not suitable for design of high pass filters?
4. What is frequency relationship for bilinear transformation?
5. What is the frequency relationship for bilinear transformation?
6. Why interpolate,i needed for any signal/sequence?
7. What is the "interpolation factor"?
8. Which signals can be interpolated?
9. Can interpolate will happens in multiple stages? If yes give reason?
21. What are the characteristics of Butterworth filters?
22. How the s-plane is mapped to z-plane in impulse invariant transformation?
23. How is the non-periodic nature of the input signal handled?
24. If a have two vectors how will i check the orthogonality of those vectors.
25. What is the importance of decimation for a given signal/sequence?
26. What do you mean Aliasing? What is the condition to avoid aliasing for sampling?
27. How does poly phasefiltering save computations in a decimation filter?
28. Give any practical application of decimation?
29. Which signals can be downsampled?
30. What happens if I violate the Nyquist criteria in down sampling or decimating?
Real time Applications

IIR filter structure which is composed of cascaded sections of second order Direct Form I
filters that use magnitude truncation. IIR filters are used in Small monitor loudspeaker. IIR
filters are used electronic crossover for a 3-way loudspeaker.
53
EXPERMENT NO-10
NARROWBAND SIGNAL THROUGH FILTERING
AIM: -
Generation of Sine Wave & Illustration of the Sampling Process in the Time Domain.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
Sinusoidal Signal Generation
The sine wave or sinusoid is a mathematical function that describes a smooth repetitive
oscillation. It occurs often in pure mathematics, as well as physics, signal processing,
electrical engineering and many other fields. Its most basic form as a function of time (t)
where:
• A, the amplitude, is the peak deviation of the function from its center position.
• ω, the angular frequency, specifies how many oscillations occur in a unit time
interval, in radians per second
• φ, the phase, specifies where in its cycle the oscillation begins at t = 0.
A sampled sinusoid may be written as:
PROCEDURE:-
 Open MATLAB
 Open new M-file
54
PROGRAM:-
clc;
t=0:0.0005:1;
a = 10
f = 13;
xa=a*sin(2*pi*f*t)
subplot(2,1,1)
plot(t,xa);
grid
xlabel('Time,msec')
;ylabel('Amplitude'
);
title('Continuous-time signal axis([01-
10.2 10.2])
subplot(2,1,2
); T = 0.01;
n = 0:T:1;
xs=a*sin(2*pi*f*n)
; k = 0:length(n)-1;
stem(k,xs);
grid
xlabel('Time index n');
ylabel('Amplitude');
title('Discrete-time signal
x[n]'); axis([0 (length(n)-1)
10.2 10.2])
55
OUTPUT:-
RESULT:-
Sinusoidal signal is generated by using MATLAB.
56
EXERCISE PROGRAM:-
1. Write program to get Discrete time Sinusoidal Signal?

4. Write a matlab program for generating u(n)-u(n-1)?
5. Write program to get Discrete time co-Sinusoidal Signal?
6. Write program to get Discrete time saw tooth Signal?
7. Write program to get Discrete time triangular Signal?
8. Write program to get addition of two sinusoidal sequences?
9. Write program to get exponential sequence?
10. Write program to get Fourier Transform of Co-Sinusoidal Signal?
11. Write program to get Inverse Fourier Transform of Co-Sinusoidal Signal?
13. Write program to get exponential growing sequence?
14. Write program to get addition of two Co-sinusoidal sequences?
15. Write program to get continues time Square Signal?
16. Write program to get continues time Sinusoidal Signal?
19. Write a matlab program for generating u(n)+u(n-2)?
20. Write program to get continues time co-Sinusoidal Signal?
VIVA QUESTIONS:-
1. Define sinusoidal signal?

2. Define C.T.S?
3. Define D.T.S?
4. Compare C.T.S & D.T.S?
5. Draw the C.T.S & D.T.S diagrams?
7. What is filter?
10. Which built in function is used to solve a given difference equation?
11. What is frequency response? Give equation for first order system andsecond order
system?
57
12. What is an LTI system?

13. What is steady state response?
14. What is Auto Regressive Model? How is the order of auto regressive model is
decided?
15.
16. Differentiate between linear and circular convolution.
17. Determine the unit step response of the linear time invariant system with impulse
response h(n)=a nu(n) a<1&-a<1
18. Determine the range of values of the parameter a for which linear time invariant
system with impulse response h(n)=a u(n) is stable.
21. Can IIR filters be Linear phase? how to make it linear Phase?
22. What is the special about minimum phase filter?
24. In signal processing, why we are much more interested in orthogonal transform?
25. What are the filter specifications required to design the analog filters?
26. What is meant by frequency response of filter?
27. What is meant by magnitude response?
28. What is meant by phase response?
29. Differentiate ideal filter and practical filter responses.
30. What are the different types of analog filter approximations?
Detection of QRS complexes in ECG signals is required to determine heart rate, and it
is an important step in the study of cardiac disorders. ECG signals are usually affected by
noise of low and high frequency. To improve the accuracy of QRS detectors several methods
have been proposed to filter out the noise and detect the characteristic pattern of QRS
complex. Most of the existing methods are at a disadvantage from relatively high
computational complexity or high resource needs making them less optimized for its
implementation on portable embedded systems, wearable devices or ultra-low power chips.
We present a new method to detect the QRS signal in a simple way with minimal
computational cost and resource needs using a novel non-linear filter.
58
EXPERMENT NO-11
DTMF SIGNAL GENERATION
AIM: -
The objective of this program is To Generate Dual Tone Multiple Frequency (DTMF)
Signals.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
Dual Tone Multiple Frequency (DTMF) Signals.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
% Dual Tone Multiple Frequency (DTMF) Signals.
clc;
clearall;
closeall;
number=input('enter a phone number with no

spaces','s'); %number=1;
fs=8192; % fs is the sampling Frequency
T=0.5; % T stores how for how long a tone will be played
x= 2*pi*[697 770 852 941];
59
y= 2*pi*[1209 1336 1477

1633]; t=[0:1/fs:T]'
tx=[sin(x(1)*t),sin(x(2)*t),sin(x(3)*t),sin(x(4)*t)]/2;
ty=[sin(y(1)*t),sin(y(2)*t),sin(y(3)*t),sin(y(4)*t)
]/2; for k=1:length(number)
switch
number(k)
case '1'
tone =
tx(:,1)+ty(:,1);
sound(tone);
stem(tone);
case '2'
tone =
tx(:,1)+ty(:,2);
sound(tone);
stem(tone);
case '3'
tone =
tx(:,1)+ty(:,3);
sound(tone);
stem(tone);
case '4'
tone = tx(:,2)+ty(:,1); sound(tone); stem(tone); otherwise

disp('invalid number');
end
pause(2.70
)
60
end;
OUTPUT:-
Input: 01234
61
62
RESULT:-
Dual Tone Multiple Frequency (DTMF) Signals are generated by using MAT LAB.
EXERCISE PROGRAM:-
1. Write a matlab program to generate a sine wave with amplitude = 3, frequency 20Hz?
2. Write a matlab program to generate a cos wave with amplitude = 3, frequency 20Hz?
3. Write a matlab program to generate a triangular wave with amplitude = 8, frequency
10Hz?
4. Write a matlab program to generate a square wave with amplitude = 2, frequency
10kHz?
5. Write a matlab program to generate a sinc wave with amplitude = -8, frequency5Khz?
6. Write a matlab program to generate a sine wave with amplitude = 7, frequency 29Hz.
7. Write a matlab program to generate a cos wave with amplitude = 9, frequency 50Hz.
100Hz.
63
10kHz.
10. Write a matlab program to generate a sinc wave with amplitude = 5, frequency5Khz.
11. Write a matlab program to generate a sine wave with amplitude = 17, frequency
29kHz.
12. Write a matlab program to generate a cos wave with amplitude = 19,
frequency600kHz.
100Hz.
14. Write a matlab program to generate a sawtooth wave with amplitude = 20, frequency
15kHz.
15. Write a matlab program to generate a sinc wave with amplitude = 8, frequency85Khz.
20Hz?
1kHz?
18. Write a matlab program to generate a sinc wave with amplitude = frequency5Khz?
19. Write a matlab program to generate a sine wave with amplitude = 7, frequency 29Hz?
20. Write a matlab program to generate a cos wave with amplitude = 9, frequency 50Hz?
VIVA QUESTIONS:-
1. Define Signal?
2. Define determistic and Random Signal?
3. Define Delta Function?
4. What is Signal Modeling?
5. Define Periodic and a periodic Signal?
11.What is the need for spectral estimation?
12.Determine the power spectrum density?
17. What is filter?
21. Define sinusoidal signal
22. Define C.T.S .
23. Define D.T.S.
24. Compare C.T.S & D.T.S.
64
27. Define signal and signal processing.
28. Differentiate digital and analog signals?
29. How the DSP processor will differ from conventional processors?
30. Expand the abbreviation TMS320C 5X/6X
In earlier days, our telephone systems were operated by manually in a telephone

exchange room. The callers will pick up the phone and giving instruction to the operator to
connect their destination line. The DTMF technology provides ultimate solutions for the
telephone industries which is used to switch two lines automatically. The DTMF stands for
‘Dual Tone Multi-frequency’ which is one of the techniques for converting the analogue
signal to digital using DTMF decoder. The DTMF decoder circuit mostly used in mobile
communications system which recognizes the sequence of DTMF tones from the standard
keypad of the mobile phone.
65
EXPERMENT NO-12
IMPLEMENTATION OF DECIMATION PROCESS.
AIM: -
The objective of this program is To Perform Decimation on the Given Input Sequence.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
In digital signal processing, decimation is the process of reducing the sampling rate of
a signal. Complementary to interpolation, which increases sampling rate, it is a specific case
of sample rate conversion in a multi-rate digital signal processing system. Decimation utilizes
filtering to mitigate aliasing distortion, which can occur when simply down sampling a signal.
A system component that performs decimation is called a decimator.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
%decimation%
clc;
clear all;
close all;
N=input('enter sample value');
n=0:N-1;
m=input('enter down sampling factor');
x=sin(2*pi*0.043*n)+sin(2*pi*0.031*n);
y=decimate(x,m,'fir');
subplot(2,1,1);
stem(n,x(1:N));
xlabel('time');
ylabel('amp');
title('input sequence');
66
t=0:(N/m)-1;
subplot(2,1,2);
stem(t,y(1:N/m));
xlabel('time');
ylabel('amp');
title('output sequence');
OUTPUT:-
enter sample value65

enter down sampling factor3
input sequence
2
1
amp
-1
-2
0 10 20 30 40 50 60 70
time
output sequence
2
1
amp
-1
-2
0 2 4 6 8 10 12 14 16 18 20
time
RESULT:-
This MATLAB program has been written to perform Decimation on the Given
Input Sequence.
67
EXERCISE PROGRAM:-
1. Write a matlab program to illustrate the effect of anti-aliasing filter?

2. Write a matlab program to illustration of upsampling?
3. Write a matlab program to illustration of downsampling?
4. Write a matlab program to illustration of effect of upsampling in frequency
domain?
5. Write a matlab program to illustration of effect of downsampling in frequency
domain?
7. Write a matlab program to plot magnitude response of comb filter?
8. Write a matlab program to plot magnitude response of allpass filter?
9. Write a matlab program to plot magnitude & unwrapped phase response of two
transfer functins?
10. Write a matlab program to design a filter that eliminates high frequency
component in a CT signal?
11. Write a matlab program to illustrate the effect of aliasing filter?
12. Write a matlab program to illustration of upsampling with sampling factor 5 ?
13. Write a matlab program to illustration of downsampling with sampling factor 5 ?
14. Write a matlab program to illustration of effect of upsampling in time domain?
15. Write a matlab program to illustration of effect of downsampling in time domain?
16. Write a matlab program to illustration of upsampling with sampling factor10?
17. Write a matlab program to illustration of downsampling with sampling factor 15 ?
18. Write a matlab program to illustration of upsampling with sampling factor 20?
19. Write a matlab program to illustration of downsampling with sampling factor 20?
20. Write a matlab program to illustration of downsampling with sampling factor 20?
VIVA QUESTIONS:-
1. Define decimation?
2. Define multi rate signal processing?
3. What are the effects of coefficient quantization in FIR filters?
4. What is quantization process?
5. What is transmultiplexer? What is its use?
6. What is the DC gain of a FIR filter?
7. What do you mean by cut-off frequency?
8. Give the difference between analog and digital filter?
9. What is the difference between type 1 and type 2 filter structure?
10. what is the role of delay element in filter design?
11. Explain how the frequency is filter in filters?
12. Differences between Butterworth chebyshev filters?
68
13. Can IIR filters be Linear phase? how to make it linear Phase?
15. What is the special about maximum phase filter?
16. In signal processing, why we are much more interested in orthogonal transform?
REAL TIMEAPPLICATIONS:
It was originally developed for processing broadband, low frequency geophysical data in
real time using a small digital processor. Sine and Cosine transforms are applied to sequences
of data produced by successively applying a low pass digital filter and decimation by two
operator to the original data. The resulting spectra are the average of estimates which are
independent in time and represent constant percentage bandwidths. Thus, the technique is
particularly well suited to applications where it is desireable to reject intervals of bad data or
those which are naturally expressed on a log-frequency scale.
69
EXPERMENT NO-13
IMPLEMENTATION OF INTERPOLATION PROCESS.
AIM: -
The objective of this program is To Perform up sampling on the given Input Sequence.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
Up sampling is interpolation, applied in the context of digital signal processing and

sample rate conversion. When up sampling is performed on a sequence of samples of a
continuous function or signal, it produces an approximation of the sequence that would have
been obtained by sampling the signal at a higher rate (or density, as in the case of a
photograph). For example, if compact disc audio is up sampled by a factor of 5/4, the resulting
sample-rate increases from 44,100 Hz to 55,125 Hz.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
%interpolation%
clc;
clear all;
close all;
N=input('enter sample value');
n=0:N-1;
L=input('enter up sampling factor');
x=sin(2*pi*0.043*n)+sin(2*pi*0.031*n);
y=interp(x,L);
subplot(2,1,1);
stem(n,x(1:N));
xlabel('time');
ylabel('amp');
70
t=0:(N*L)-1;
subplot(2,1,2);
stem(t,y(1:N*L));
xlabel('time');
ylabel('amp');
title('output sequence');
OUTPUT:-
enter sample value50

enter up sampling factor
input sequence
2
1
amp
-1
-2
0 5 10 15 20 25 30 35 40 45 50
time
output sequence
2
1
amp
-1
-2
0 50 100 150 200 250
time
RESULT:-
This MATLAB program has been written to perform interpolation on the Given
Input Sequence.
71
EXERCISE PROGRAM:-

domain?
domain?
transfer functins?
12. Write a matlab program to illustration of up sampling?
13. Write a matlab program to illustration of down sampling?
domain?
domain?
18. Write a matlab program to plot magnitude response of bandpass filter?
19. Write a matlab program to plot magnitude response?
20. Write a matlab program to design a filter that eliminates low frequency component
in a CT signa.
VIVA QUESTIONS:-
1. How aliasing can be avoided?

2. Which type of interpolation is used to reconstruct the signal?
3. What is aliasing?
4. Define interpolation?
5. What is pre-alias filter?
8. How does poly phasefiltering save computations in a decimation filter?
72
12. Can we do decimate in multiple stages?
13. What are "decimation" and "downsampling"?
14. What is the "decimation factor
15. How does poly phase filtering save computations in an interpolation filter?
16. Why do we need I&Q signals?
17. What is Interpolation and decimation filters and why we need it?
18. What are "up sampling" and "interpolation"?
19. Why interpolate, needed for any signal/sequence?
21. What kind of processor is DSP processor?
22. What are the main building
23. blocks of DSP processor?
24. What is the main function of MAC unit?
25. Explain VLIW architecture?
26. Explain the significance of convolution.
27. Define linear convolution.
28. Why linear convolution is called as a periodic convolution?
29. Why zero padding is used in linear convolution?
30. What are the four steps to find linear convolution?
APPLICATIONS:
Interpolation is a technique for obtaining new unknown data points within the range of
discrete known data points and is often used to recover an image from its downsampled
version, or to simply perform image expansion.
73
EXPERMENT NO-14
IMPLEMENTATION OF I/D SAMPLING RATE CONVERTERS
AIM: -
To study sampling rate conversion by a rational form using MATLAB.
SOFTWARE REQURIED:-
1. MATLAB R2010a.
2. Windows XP SP2.
THEORY:-
"Up sampling" is the process of inserting zero-valued samples between original

samples to increase the sampling rate. (This is called "zero-stuffing".) Up sampling adds to the
original signal undesired spectral images which are centered on multiples of the original
sampling rate.
"Interpolation", in the DSP sense, is the process of up sampling followed by filtering.

(The filtering removes the undesired spectral images.) As a linear process, the DSP sense of
interpolation is somewhat different from the "math" sense of interpolation, but the result is
conceptually similar: to create "in-between" samples from the original samples. The result is
as if you had just originally sampled your signal at the higher rate.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
% interpolation/dismation sampling %
clc; clear all;close all;

N=input('enter the sample value');
n=0:N-1;
l=input('enter up sampling factor');
m=input('enter down sampling factor');
74
x=sin(2*pi*0.043*n)+sin(2*pi*0.03*n);
y=resample(x,l,m);
subplot(2,1,1);
stem(n,x(1:N));
xlabel('time');
ylabel('amplitude');
t=0:(N*l/m)-1;
subplot(2,1,2);
stem(t,y(1:N*l/m));
xlabel('time');
title('input sampling sequence');
OUTPUT:-
enter the sample value30

enter up sampling factor7
enter down sampling factor2
input sequence
2
1
amplitude
-1
-2
0 5 10 15 20 25 30
time
input sampling sequence
2
1
amplitude
-1
-2
0 20 40 60 80 100 120
time
75
RESULT:-
Thus sampling rate conversion by a rational form is performed using MATLAB.
EXERCISE PROGRAM:-

domain?
domain?
transfer functins?
12. Write a matlab program to illustration of I/D sampling with sampling factor 5 ?
14. Write a matlab program to illustration of effect of I/D sampling in time domain?
15. Write a matlab program to illustration of effect of I/D sampling in time domain?
16. Write a matlab program to illustration of I/D sampling with sampling factor10?
18. Write a matlab program to illustration of I/D sampling with sampling factor 20?
76
VIVA QUESTIONS:-
1. What is multi rate signal processing?

2. What is the need for anti-imaging filter after up sampling a signal?
3. What is the need for anti-imaging filter prior to down sampling?
4. Define down sampling?
5. What is meant by up sampling?
9. Can we do decimate in multiple stages?
10. What are "decimation" and "downsampling"?
11. What isthe "decimation factor
12. How does polyphase filtering save computations in an interpolation filter?
13. Why do we need I&Q signals?
14. What is Interpolation and decimation filters and why we need it?
15. What are "upsampling" and "interpolation"?
16. Why interpolate,i needed for any signal/sequence?
18. Which signals can be interpolated?
19. Can interpolate will happens in multiple stages? If yes give reason?
21. What is "bit reversal"?
25. What is FIR and IIR filter define, and distinguish between these two?
26. What is window method? How you will design an FIR filter using window method?
27. What are low-pass and band-pass filter and what is the difference between these two?
28. What is the matlab command for Hamming window? Explain.
29. What do you mean by built in function ‘abs’ and where it is used?
30. Explain how the FIR filter are stable?
Interpolation is a technique for obtaining new unknown data points within the range of
discrete known data points and is often used to recover an image from its downsampled
version, or to simply perform image expansion.
It was originally developed for processing broadband, low frequency geophysical data in
real time using a small digital processor. Sine and Cosine transforms are applied to sequences
of data produced by successively applying a low pass digital filter and decimation by two
operator to the original data. The resulting spectra are the average of estimates which are
independent in time and represent constant percentage bandwidths. Thus, the technique is
77
particularly well suited to applications where it is desireable to reject intervals of bad data or
those which are naturally expressed on a log-frequency scale.
EXPERMENT NO-15
IMPULSE RESPONSE OF FIRST ORDER AND SECOND ORDER SYSTEM
AIM: -
To find the impulse response of the given equation

y(n)-y(n-1+0.9y(n-2)=x(n)
SOFTWARE REQURIED:-
1.MATLAB R2010a.
2.Windows XP SP2.
THEORY:-
Second order systems are the systems or networks which contain two or
more storage elements and have describing equations that are second order
differential equations.
The frequency response of second order filters is characterised by three

filter parameters: the gain k, the corner frequency and the quality factor Q.
PROCEDURE:-
 Open MATLAB
 Open new M-file
PROGRAM:-
%To find the impulse response of discrete time system

%y(n)-y(n-1)+0.9y(n-2)=xn)%
clc;
clear all;
close all;
b=input('Enter the coefficeints of x(n):b=');
a=input('Enter the coefficeints of x(n):a=');
N=input('Enter the order of N =');
78
h=impz(b,a,N);
n=0:N-1;
subplot(2,1,1);
stem(n,h);
xlabel('discrete time');
title('Impulse Response');
subplot(2,1,2);
zplane(b,a);
xlabel('real axis');
ylabel('imaginary axis');
title('pole zero in z-plane');
OUTPUT:-
Enter the coefficeints of x(n):b=[1]

Enter the coefficeints of x(n):a=[1 -1 .9]
Enter the order of N =2
Impulse Response
1
amplitude
0.5
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
discrete time
pole zero in z-plane
1
imaginary axis
0.5
2
0
-0.5
-1
-3 -2 -1 0 1 2 3
real axis
79
RESULT:-
Hence the impulse response of the given system is performed.
EXERCISE PROGRAM:-
equation y(n)-7y(n-1)+9y(n-2)=3x(n)-2x(n-1)?
equation 3y(n)+5y(n-1)=9x(n)?
equation 9 y(n)-2y(n-1)+7y(n-2)-3y(n-3)=6x(n)+x(n-1)?
equation 8y(n)+6y(n-1)=4x(n)+2x(n-1)?
equation 3y(n)-8y(n-1)+9y(n-2)=9x(n)+5x(n-1) ?
equation 6y(n)-5y(n-1)=9x(n)+5x(n-1) -7x(n-2)?
equation 9y(n)-8y(n-1)+2y(n-2)=9x(n)-3x(n-1) ?
equation 2y(n)-8y(n-1)=9x(n)+5x(n-1) ?
equation 9y(n)-8y(n-1)+9y(n-2)=9x(n)+5x(n-1)-x(n-2) ?
equation 3y(n)-8y(n-1)=7x(n)-3x(n-1) ?
VIVA QUESTIONS:-
1. What is the commend to find phase angle?

2. What is the commend to find frequency response?
3. What is transition band?
4. What is the formula for Z-transform?
5. What is the relationship b/w impulse response& frequency response?
7. What is the relation between DFT and Z-Transform.
80

20. What is filter?
21. Give mathematical definition of circular convolution.
22. Why circular convolution is called as periodic convolution?
30. Explain modulo N operation
APPLICATIONS:
In practical systems, it is not possible to produce a perfect impulse to serve as input for testing;
therefore, a brief pulse is sometimes used as an approximation of an impulse. Provided that the
pulse is short enough compared to the impulse response, the result will be close to the true,
theoretical, impulse response. In many systems, however, driving with a very short strong pulse
may drive the system into a nonlinear regime, so instead the system is driven with a pseudo-
random sequence, and the impulse response is computed from the input and output signals..
81

DSP Lab Manual 2018-19 PrintCopy MANU

Uploaded by

Copyright:

Available Formats

DSP Lab Manual 2018-19 PrintCopy MANU

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

DSP Lab Manual 2018-19 PrintCopy MANU

Uploaded by

Copyright:

Available Formats

DSP LAB MANUAL Dept of ECE

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

DIGITAL SIGNAL PROCESSING

III B.TECH -II Semester (ECE)

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

VISION AND MISSION OF THE INSTITUTE

1. An atmosphere that facilitates personal commitment to the educational success of

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

VISION AND MISSION OF THE DEPARTMENT

To achieve the Mission the department will

PROGRAMME EDUCATIONAL OBJECTIVES

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

i An ability to function as an individual, and as a member or leader in diverse teams and in

PROGRAM SPECIFIC OUTCOMES

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

COURSE OUTCOMES FOR DIGITAL SIGNAL PROCESSING LAB

1. Generate different types of Signals using MATLAB and DSP

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

DIGITAL SIGNAL PROCESSING LAB

1) To find DFT/IDFT of given DT signal

Dundigal, Quthbullapur (M), Hyderabad-500043.

Department of Electronics & Communication Engineering

INSTRUCTIONS TO THE STUDENT FJHINSTRUCTIOTO

2. Students will work individually in hardware laboratories and in computer laboratories.

SL.NO. EXPERIMENT NAME PAGENO.

1 Generation of sinusoidal waveform /signal based on recursive difference 9

2 Histogram of White Gaussian Noise and Uniformly Distributed Noise. 13

3 To find DFT/IDFT of given DT signal. 18

4 To find frequency response of a given system (Transfer function 25

6 Implementation of FFT of given sequence. 34

8 Implementation of LP FIR filters for a given sequence. 43

9 Implementation of HP IIR filters for a given sequence. 50

10 Generation of Narrow band signal through filtering. 54

11 Generation of DTMF signals. 59

13 Implementation of Interpolation process. 70

14 Implementation of I/D sampling rate converters. 74

Generating of a sinusoidal signal based on recursive difference equation.

Real Life Applications

Trigonometry is especially important in architecture because it allows the architect to

1. Write program to get Discrete time Sinusoidal Signal?

2. Write program to get Fourier Transform of Sinusoidal Signal?

3. Write program to get Inverse Fourier Transform of Sinusoidal Signal?

4. Write Program for the following Function

HISTOGRAM OF GAUSSIAN NOISE

To write a MATLAB program for computation of statistical parameters of Gaussian

8.what is auto correlation?

1.Write a MATLAB program to plot the co-sinusoidal signal.

DFT/IDFT OF A SEQUENCE WITHOUT USING THE INBUILT FUNCTIONS

Given a sequence of N samples f(n), indexed by n = 0..N-1, the Discrete Fourier

F(k) are often called the 'Fourier Coefficients' or 'Harmonics'.

In general, both f(n) and F(k) are complex.

10.0000 -2.0000 + 2.0000i -2.0000 - 0.0000i -2.0000 - 2.0000i

enter the input sequence[10 -2+2j -2 -2-2j]

y= 2pi[1209 1336 1477