DSP LAB MANUAL VITS CONTENTS S.No Experiment Name Page No. 1. Introduction to MATLAB 1 2. Generation of Different Types of Signals 4 3. Generation of Sum of Sinusoidal Signals 8 4. Linear Convolution 10 5. N-POINT FFT 12 6. Auto Correlation & Power Density Spectrum 15 7. To Find Frequency Response of FIR Low Pass /High Pass Filters 17 8. To find frequency response of IIR low pass / High pass filters 29 9. Architecture and Instruction Set of DSPCHIP-TMS320C6713 32 10. CODE COMPOSER STUDIO 54 11. Linear Convolution using CC Studio 58 12. Circular Convolution using CC Studio 63 0
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
DSP LAB MANUAL VITS
CONTENTS
S.No Experiment Name Page No.
1. Introduction to MATLAB 1
2. Generation of Different Types of Signals 4
3. Generation of Sum of Sinusoidal Signals 8
4. Linear Convolution 10
5. N-POINT FFT 12
6. Auto Correlation & Power Density Spectrum 15
7. To Find Frequency Response of FIR Low Pass /High Pass Filters 17
8. To find frequency response of IIR low pass / High pass filters 29
9. Architecture and Instruction Set of DSPCHIP-TMS320C6713 32
10. CODE COMPOSER STUDIO 54
11. Linear Convolution using CC Studio 58
12. Circular Convolution using CC Studio 63
13. FIR FILTER DESIGN USING TMS320C6713 DSP PROCESSOR 67
[RECTANGULER/TRIANGULAR/KAISER WINDOW]
14. IIR FILTER DESIGN USING TMS320C6713 DSP PROCESSOR 76
15. N-POINT Fast Fourier Transform (FFT) 85
16. To Compute Power Density Spectrum of A Sequence 93
[USING TMS320C6713 DSP PROCESSOR]
ADDITIONAL EXPERIEMENTS
1. To Find Frequency Response of Band pass filter using FIR
2. To Find Frequency Response of Band pass filter using IIR
3. To Find Frequency Response of Band reject filter using FIR
4. To Find Frequency Response of Band reject filter using IIR
Annexure –I Viva Questions 100
0
DSP LAB MANUAL VITS
1
DSP LAB MANUAL VITS
EXP.NO: 1INTRODUCTION TO MATLAB
MATLAB (MATrix LABoratory):
MATLAB is a software package for high-performance language for technical
computing. It integrates computation, visualization, and programming in an easy-to-use
environment where problems and solutions are expressed in familiar mathematical
notation. Typical uses include the following
Math and computation
Algorithm development
Data acquisition
Modeling, simulation, and prototyping
Data analysis, exploration, and visualization
Scientific and engineering graphics
Application development, including graphical user interface building
The name MATLAB stands for matrix laboratory. MATLAB was originally written
to provide easy access to matrix software developed by the LINPACK and EISPACK
projects. Today, MATLAB engines incorporate the LAPACK and BLAS libraries,
embedding the state of the art in software for matrix computation.
MATLAB has evolved over a period of years with input from many users. In
university environments, it is the standard instructional tool for introductory and
advanced courses in mathematics, engineering, and science. In industry, MATLAB is the
tool of choice for high-productivity research, development, and analysis.
MATLAB features a family of add-on application-specific solutions called
toolboxes. Very important to most users of MATLAB, toolboxes allow learning and
applying specialized technology. Toolboxes are comprehensive collections of MATLAB
2
DSP LAB MANUAL VITS
functions (M-files) that extend the MATLAB environment to solve particular classes of
problems. Areas in which toolboxes are available include Image processing, signal
processing, control systems, neural networks, fuzzy logic, wavelets, simulation, and many
others.
The main features of MATLAB
1. Advance algorithm for high performance numerical computation ,especially in the
Field matrix algebra
2. A large collection of predefined mathematical functions and the ability to define
one’s own functions.
3. Two-and three dimensional graphics for plotting and displaying data
4. A complete online help system
5. Powerful, matrix or vector oriented high level programming language for individual
applications.
6. Toolboxes available for solving advanced problems in several application areas
% BASIC OPERATIONS ON MATRICES
clc; % clear the command windowclear all; %clear the workspaceclose all; %clear the figure window
%generaion of exponentialn=0:1:20;y=exp(-.3*n);figure;subplot(2,2,1)stem(n,y);title('discrete exponential');xlabel('time');ylabel('Amplitude');subplot(2,2,2);plot(n,y);title('continuous exponential');xlabel('time');ylabel('Amplitude');
%generaion of sinusoidal signaln=0:1:50;y=sin(0.1*pi*n);subplot(2,2,3)stem(n,y);title('discrete sinusoidal signal');xlabel('time');ylabel('Amplitude');subplot(2,2,4);plot(n,y);title('continuous sinusoidal signal');xlabel('time');ylabel('Amplitude');
6
DSP LAB MANUAL VITS
Generation of signals
output:
7
DSP LAB MANUAL VITS
8
DSP LAB MANUAL VITS
EXP.NO: 3GENERATION OF SUM OF SINUSOIDAL SIGNALS
Aim: To generate sum of sinusoidal signals Using MATLAB Software.
EQUIPMENTS:
MATLAB Software
Program:clc;close all;clear all;
t=0:0.05:3*pi;x1=sin(t*5);%sine wave with period 5x2=sin(t*9);%sine wave with period 9x3=x1+x2;%sum of x1 and x2subplot(3,1,1);plot(t,x1)xlabel('time'),ylabel('amplitude')title('sin signal-1');subplot(3,1,2);plot(t,x2)xlabel('time'),ylabel('amplitude')title('sin signal-2');subplot(3,1,3);plot(t,x3)xlabel('time'),ylabel('amplitude')title('sum of sin signals');
9
DSP LAB MANUAL VITS
sum of sinusoidal signals output:
10
DSP LAB MANUAL VITS
EXP.NO: 4 LINEAR CONVOLUTION
Aim: To find the out put with linear convolution operation Using MATLAB
Software.
EQUIPMENTS:
MATLAB Software
Theory:
Linear Convolution involves the following operations.1. Folding2. Multiplication3. Addition4. Shifting
These operations can be represented by a Mathematical Expression as follows:
x[ ]= Input signal Samplesh[ ]= Impulse response co-efficient.y[ ]= Convolution output. n = No. of Input samples h = No. of Impulse response co-efficient.
title('linear convolution')disp('The resultant signal is');disp(y)
linear convolution output:
enter input sequence[1 4 3 2]enter impulse response[1 0 2 1]The resultant signal is 1 4 5 11 10 7 2
12
DSP LAB MANUAL VITS
EXP.NO: 5 N-POINT FFT
Aim: To compute N-point FFTEQUIPMENTS:
MATLAB Software
Theory:DFT of a sequence
Where N= Length of sequence.
K= Frequency Coefficient.
n = Samples in time domain.
FFT : -Fast Fourier transform .
There are Two methods.
1.Decimation in time (DIT FFT).
2. Decimation in Frequency (DIF FFT).
Why we need FFT ?
The no of multiplications in DFT = N2.
The no of Additions in DFT = N(N-1).
For FFT.
The no of multiplication = N/2 log 2N.
The no of additions = N log2 N.
13
DSP LAB MANUAL VITS
Program:clc;close all;clear all;x=input('enter the sequence');N=length(x);n=0:1:N-1;y=fft(x,N)subplot(2,1,1);stem(n,x);title('input sequence');xlabel('time index n----->');ylabel('amplitude x[n]----> ');subplot(2,1,2);stem(n,y);title('output sequence');xlabel(' Frequency index K---->');ylabel('amplitude X[k]------>');
14
DSP LAB MANUAL VITS
N-point FFT Output:enter the sequence[2 3 1 4]
y =
10.0000 1.0000 + 1.0000i -4.0000 1.0000 - 1.0000i
15
DSP LAB MANUAL VITS
EXP.NO: 6 AUTO CORRELATION & POWER DENSITY SPECTRUM
Aim: To compute auto correlation and then find the power density spectrum.EQUIPMENTS:MATLAB SoftwareSignal Processing Toolbox
Programclc;close all;clear all;x=input('enter the sequence');N=length(x);n=0:1:N-1;y=xcorr(x,x);subplot(3,1,1);stem(n,x);xlabel(' n----->');ylabel('Amplitude--->');title('input seq');subplot(3,1,2);N=length(y);n=0:1:N-1;stem(n,y);xlabel('n---->');ylabel('Amplitude----.');title('autocorr seq for input');disp('autocorr seq for input');disp(y)p=fft(y,N);subplot(3,1,3);stem(n,p);xlabel('K----->');ylabel('Amplitude--->');title('psd of input');disp('the psd fun:');disp(p)
16
DSP LAB MANUAL VITS
auto correlation & power density spectrum enter the sequence[1 0 1 2]autocorr seq for input 2 1 2 6 2 1 2
the psd fun:Columns 1 through 5 16.0000 -4.0048 - 1.9286i 3.6174 + 4.5361i -0.6126 - 2.6840i -0.6126 + 2.6840iColumns 6 through 7 3.6174 - 4.5361i -4.0048 + 1.9286i
17
DSP LAB MANUAL VITS
EXP.NO: 7TO FIND FREQEUNCY RESPONSE OF FIR LOW PASS / HIGH
PASS FILTERS
Aim: To find frequency response of FIR low pass/ high pass filter.
EXP.NO: 8TO FIND FREQEUNCY RESPONSE OF IIR LOW PASS / HIGH
PASS FILTERS
Aim: To find frequency response of IIR low pass/ high pass filter.
EQUIPMENTS:MATLAB Software
PROGRAM:% IIR filters LPF & HPFclc;clear all;close all;disp('enter the IIR filter design specifications');rp=input('enter the passband ripple');rs=input('enter the stopband ripple');wp=input('enter the passband freq');ws=input('enter the stopband freq');fs=input('enter the sampling freq');w1=2*wp/fs;w2=2*ws/fs;[n,wn]=buttord(w1,w2,rp,rs,'s');c=input('enter choice of filter 1. LPF 2. HPF \n ');if(c==1)disp('Frequency response of IIR LPF is:');[b,a]=butter(n,wn,'low','s');endif(c==2)disp('Frequency response of IIR HPF is:');[b,a]=butter(n,wn,'high','s');endw=0:.01:pi;[h,om]=freqs(b,a,w);m=20*log10(abs(h));an=angle(h);figure,subplot(2,1,1);plot(om/pi,m);title('magnitude response of IIR filter is:');xlabel('(a) Normalized freq. -->');ylabel('Gain in dB-->');subplot(2,1,2);plot(om/pi,an);title('phase response of IIR filter is:');xlabel('(b) Normalized freq. -->');
30
DSP LAB MANUAL VITS
ylabel('Phase in radians-->');
OUTPUT:enter the IIR filter design specificationsenter the passband ripple.15enter the stopband ripple60enter the passband freq1500enter the stopband freq3000enter the sampling freq7000enter choice of filter 1. LPF 2. HPF 1Frequency response of IIR LPF is:
31
DSP LAB MANUAL VITS
enter the IIR filter design specificationsenter the passband ripple.15enter the stopband ripple60enter the passband freq1500enter the stopband freq3000enter the sampling freq7000enter choice of filter 1. LPF 2. HPF 2Frequency response of IIR HPF is:
32
DSP LAB MANUAL VITS
EXP.No:9
ARCHITECTURE AND INSTRUCTION SET OFDSPCHIP-TMS320C6713
Aim: To Study the architecture and Instruction set of DSP chips
Features of Highest-Performance Floating-Point Digital Signal Processor Tms320c6713
Enhanced Harvard Architecture
VLIW Parallel Architecture
Rich Addressing modes
Two general purpose Register files (A0-A15 & B0-B15)
DPINT (Convert Double-Precision Floating-Point Value to Integer)DPINT .L1 A1:A0,A4
MVD (Move From Register to Register, Delayed)
MVK/MVKL (Move a Signed Constant Into a Register and Sign-Extend)
48
DSP LAB MANUAL VITS
PACKHL2 (Pack 16 MSB, 16 LSB Into Packed 16-Bit):
Floating Point Addition
.text
NUM1 .FLOAT -2.5 ;
NUM1 .FLOAT 8.6;
Start:
mvkl NUM1,a7
mvkh NUM1,a7
mvkl NUM2,a8
mvkh NUM2,a8
ldw *a7,a4
ldw *a8,a5
mvkl 90030000h,a3
mvkh 90030000h,a3
addsp a4,a5,a2
stw a2,*a3++
.end Result: 6.1 ; 40C3 3334h
Floating Point Multiplication
49
DSP LAB MANUAL VITS
.text
NUM1 .FLOAT -2.5 ;
NUM1 .FLOAT 8.6;
Start:
mvkl NUM1,a7
mvkh NUM1,a7
mvkl NUM2,a8
mvkh NUM2,a8
ldw *a7,a4
ldw *a8,a5
mvkl 90030000h,a3
mvkh 90030000h,a3
addsp a4,a5,a2
stw a2,*a3++
.end Result: 6.1 ; 40C3 3334h
50
DSP LAB MANUAL VITS
TMS C6713DSK
The 6713 DSK is a low-cost standalone development platform that enables customers to evaluate and develop applications for the TI C67XX DSP family. The DSK also serves as a hardware reference design for the TMS320C6713 DSP. Schematics, logic equations and application notes are available to ease hardware development and reduce time to market.
The DSK uses the 32-bit EMIF for the SDRAM (CE0) and daughtercard expansion interface (CE2 and CE3). The Flash is attached to CE1 of the EMIF in 8-bit mode.
An on-board AIC23 codec allows the DSP to transmit and receive analog signals. McBSP0 is used for the codec control interface and McBSP1 is used for data. Analog audio I/O is done through four 3.5mm audio jacks that correspond to microphone input, line input, line output and headphone output. The codec can select the microphone or the line input as the active input. The analog output is driven to both the line out (fixed gain) and headphone (adjustable gain) connectors. McBSP1 can be re-routed to the expansion connectors in software.
A programmable logic device called a CPLD is used to implement glue logic that ties the board components together. The CPLD has a register based user interface that lets the user configure the board by reading and writing to the CPLD registers. The registers reside at the midpoint of CE1.
The DSK includes 4 LEDs and 4 DIPswitches as a simple way to provide the user with interactive feedback. Both are accessed by reading and writing to the CPLD registers.
An included 5V external power supply is used to power the board. On-board voltage regulators provide the 1.26V DSP core voltage, 3.3V digital and 3.3V analog voltages. A voltage supervisor monitors the internally generated voltage, and will hold the board in reset until the supplies are within operating specifications and the reset button is released. If desired, JP1 and JP2 can be used as power test points for the core and I/O power supplies.
Code Composer communicates with the DSK through an embedded JTAG emulator with a USB host interface. The DSK can also be used with an external emulator through the external JTAG connector.
TMS320C6713 DSP FEATURES Highest-Performance Floating-Point Digital Signal Processor (DSP):
Eight 32-Bit Instructions/Cycle 32/64-Bit Data Word 300-, 225-, 200-MHz (GDP), and 225-, 200-, 167-MHz (PYP) Clock Rates 3.3-, 4.4-, 5-, 6-Instruction Cycle Times 2400/1800, 1800/1350, 1600/1200, and 1336/1000 MIPS /MFLOPS
51
DSP LAB MANUAL VITS
Rich Peripheral Set, Optimized for Audio Highly Optimized C/C++ Compiler Extended Temperature Devices Available
Advanced Very Long Instruction Word (VLIW) TMS320C67x™ DSP Core Eight Independent Functional Units:
Two ALUs (Fixed-Point) Four ALUs (Floating- and Fixed-Point) Two Multipliers (Floating- and Fixed-Point)
Load-Store Architecture With 32 32-Bit General-Purpose Registers Instruction Packing Reduces Code Size All Instructions Conditional
Instruction Set Features Native Instructions for IEEE 754
Boot Mode: HPI, 8-, 16-, 32-Bit ROM Boot Endianness: Little Endian, Big Endian
32-Bit External Memory Interface (EMIF) Glueless Interface to SRAM, EPROM, Flash, SBSRAM, and SDRAM 512M-Byte Total Addressable External Memory Space
Enhanced Direct-Memory-Access (EDMA) Controller (16 Independent Channels) 16-Bit Host-Port Interface (HPI) Two Multichannel Audio Serial Ports (McASPs)
Two Independent Clock Zones Each (1 TX and 1 RX) Eight Serial Data Pins Per Port:
Individually Assignable to any of the Clock Zones Each Clock Zone Includes:
Programmable Clock Generator Programmable Frame Sync Generator TDM Streams From 2-32 Time Slots Support for Slot Size:
8, 12, 16, 20, 24, 28, 32 Bits Data Formatter for Bit Manipulation
Wide Variety of I2S and Similar Bit Stream Formats Integrated Digital Audio Interface Transmitter (DIT) Supports:
S/PDIF, IEC60958-1, AES-3, CP-430 Formats
52
DSP LAB MANUAL VITS
Up to 16 transmit pins Enhanced Channel Status/User Data
Extensive Error Checking and Recovery Two Inter-Integrated Circuit Bus (I2C Bus™) Multi-Master and Slave Interfaces Two Multichannel Buffered Serial Ports:
Shut down and Power off the PC Connect the supplied USB port cable to the board Connect the other end of the cable to the USB port of PC
Note: If you plan to install a Microphone, speaker, or Signal generator/CRO these must be plugged in properly
Before you connect power to the DSK Plug the power cable into the board Plug the other end of the power cable into a power outlet The user LEDs should flash several times to indicate board is operational When you connect your DSK through USB for the first time on windows loaded
PC the new hardware found wizard will come up. So, Install the drivers (The CCS CD contains the require drivers for C5416 DSK).
Install the CCS software for C5416 DSK
54
DSP LAB MANUAL VITS
EXP.No 10: CODE COMPOSER STUDIO
INTRODUCTION TO CODE COMPOSER STUDIO
Code Composer is the DSP industry's first fully integrated development environment
(IDE) with DSP-specific functionality. With a familiar environment liked MS-based C+
+TM, Code Composer lets you edit, build, debug, profile and manage projects from a
single unified environment. Other unique features include graphical signal analysis,
injection/extraction of data signals via file I/O, multi-processor debugging, automated
testing and customization via a C-interpretive scripting language and much more.
CODE COMPOSER FEATURES INCLUDE: IDE
Debug IDE
Advanced watch windows
Integrated editor
File I/O, Probe Points, and graphical algorithm scope probes
Advanced graphical signal analysis
Interactive profiling
Automated testing and customization via scripting
Visual project management system
Compile in the background while editing and debugging
Multi-processor debugging
Help on the target DSP
To create a system configuration using a standard configuration file:
Step 1: Start CCS Setup by double clicking on the Setup CCS desktop icon.
Step 2: select Family c67xx
Platform simulator
Endianness little
.
Step 3: Click the Import button (File import) to import our selection (c67xx_sim.ccs)
to the system configuration currently being created in the CCS Setup window.
55
DSP LAB MANUAL VITS
Step 4: Click the Save and Quit button to save the configuration in the System Registry.
Step 5: Click the Yes button to start the CCS IDE when we exit CCS Setup. The CCS
Setup closes and the CCS IDE automatically opens using the configuration we just
created.
56
DSP LAB MANUAL VITS
PROCEDURE TO WORK ON CODE COMPOSER STUDIO
Step 1: Creating a New Project From the Project menu, choose New.
In the Project Name field, type the name we want for our project. Each project we
create must have a unique name, and Click Finish. The CCS IDE creates a project
file called projectname.pjt. This file stores our project settings and references the
various files used by our project.
The Project Creation wizard window displays.
Step 2: Creating a source fileCreate a new source file using ‘File new source file ‘ pull down menu and
save the source file with .c extension in the current project name directory.
Run: Run the program using the ‘Debug-Run’ pull down menu or by clicking the
shortcut icon on the left side of program window.
Step 5: observe output using graph Choose View GraphTime/Frequency.
In the Graph Property Dialog, change the Graph Title, Start Address, and
Acquisition Buffer Size, Display Data Size, DSP Data Type, Auto scale, and
Maximum Y- Value properties to the values.
58
DSP LAB MANUAL VITS
EXP.No:11LINEAR CONVOLUTION
AIM: Verify the linear convolution operation Using DSK Code composer studioEQUIPMENTS:
TMS 320C6713 Kit.RS232 Serial CablePower CordOperating System – Windows XPSoftware – CCStudio_v3.1
THEORY:
Convolution is a formal mathematical operation, just as multiplication, addition, and integration. Addition takes two numbers and produces a third number, while convolution takes two signals and produces a third signal. Convolution is used in the mathematics of many fields, such as probability and statistics. In linear systems, convolution is used to describe the relationship between three signals of interest: the input signal, the impulse response, and the output signal.
In this equation, x(k), h(n-k) and y(n) represent the input to and output from the system at time n. Here we could see that one of the inputs is shifted in time by a value every time it is multiplied with the other input signal. Linear Convolution is quite often used as a method of implementing filters of various types.
ALGORITHM
Step 1 Declare three buffers namely Input buffer, Temporary Buffer, Output Buffer.Step 2 Get the input from the CODEC, store it in Input buffer and transfer it to the first location of the Temporary buffer.Step 3 Make the Temporary buffer to point to the last location.Step 4 Multiply the temporary buffer with the coefficients in the data memory and accumulate it with the previous output.Step 5 Store the output in the output buffer.Step 6 Repeat the steps from 2 to 5.
ASSEMBLY PROGRAM TO IMPLEMENT LINEAR CONVOLUTIONconv.asm:
//Linear convolution program in c language using CCStudio
#include<stdio.h>
60
DSP LAB MANUAL VITS
int x[15],h[15],y[15];main(){int i,j,m,n;printf("\n enter first sequence length m:");scanf("%d",&m);printf("\n enter second sequence length n:");scanf("%d",&n);printf("Enter i/p sequence for x(n):\n");for(i=0;i<m;i++)scanf("%d",&x[i]);printf("Enter i/p sequence for h(n): \n");for(i=0;i<n; i++)scanf("%d",&h[i]);// padding of zeorsfor(i=m;i<=m+n-1;i++)x[i]=0;for(i=n;i<=m+n-1;i++)h[i]=0;/* convolution operation */for(i=0;i<m+n-1;i++){y[i]=0;for(j=0;j<=i;j++){y[i]=y[i]+(x[j]*h[i-j]);}}//displaying the o/pprintf("Output (Linear Convolution) sequence is:\n ");for(i=0;i<m+n-1;i++)printf("y[%d]=%d\t",i,y[i]);}
PROCEDURE:
Open Code Composer Studio, make sure the DSP kit is turned on.
Start a new project using ‘Project-new ‘ pull down menu, save it in a
separate directory(c:\ti\myprojects) with name lconv.pjt.
Add the source files conv.asm.
to the project using ‘Projectadd files to project’ pull down menu.
Add the linker command file hello.cmd.
61
DSP LAB MANUAL VITS
(Path: c:\ti\tutorial\dsk6713\hello1\hello.cmd)
Add the run time support library file rts6700.lib.
(Path: c:\ti\c6000\cgtools\lib\rts6700.lib)
Compile the program using the ‘Project-compile’ pull down menu or by
clicking the shortcut icon on the left side of program window.
Build the program using the ‘Project-Build’ pull down menu or by
clicking the shortcut icon on the left side of program window.
Load the program (lconv.out) in program memory of DSP chip using the
‘File-load program’ pull down menu.
To View output graphically
Select view graph time and frequency.
62
DSP LAB MANUAL VITS
OUTPUT FOR LINEAR CONVOLUTION:
enter first sequence length m:4enter second sequence length n:3Enter i/p sequence for x(n):1 2 3 4Enter i/p sequence for h(n): 2 3 1Output (Linear Convolution) sequence is: y[0]=2 y[1]=7 y[2]=13 y[3]=19 y[4]=15 y[5]=4
Graph Property Graph
LINEAR CONVOLUTION
63
DSP LAB MANUAL VITS
EXP.No:12CIRCULAR CONVOLUTION
AIM: To verify the circular convolution operation Using DSK Code composer studio EQUIPMENTS:
TMS 320C6713 Kit.RS232 Serial CablePower CordOperating System – Windows XPSoftware – CCStudio_v3.1
THEORY:Circular convolution is another way of finding the convolution sum of two input
signals. It resembles the linear convolution, except that the sample values of one of the input signals is folded and right shifted before the convolution sum is found. Also note that circular convolution could also be found by taking the DFT of the two input signals and finding the product of the two frequency domain signals. The Inverse DFT of the product would give the output of the signal in the time domain which is the circular convolution output. The two input signals could have been of varying sample lengths. But we take the DFT of higher point, which ever signals levels to.. This process is called circular convolution.
0
64
DSP LAB MANUAL VITS
/*program to implement circular convolution */
#include<stdio.h> int m,n,x[30],h[30],y[30],i,j, k,x2[30],a[30];void main(){ printf(" enter the length of the first sequence\n"); scanf("%d",&m); printf(" enter the length of the second sequence\n"); scanf("%d",&n); printf(" enter the first sequence\n"); for(i=0;i<m;i++) scanf("%d",&x[i]); printf(" enter the second sequence\n"); for(j=0;j<n;j++) scanf("%d",&h[j]);
if(m-n!=0) /*If length of both sequences are not equal*/ {
if(m>n) /* Pad the smaller sequence with zero*/ { for(i=n;i<m;i++) h[i]=0; n=m; } for(i=m;i<n;i++) x[i]=0; m=n;
y[k]+=x[i]*x2[i]; } } /*displaying the result*/ printf(" the circular convolution is\n"); for(i=0;i<n;i++) printf("%d \t",y[i]); }
PROCEDURE:
Open Code Composer Studio; make sure the DSP kit is turned on.
Start a new project using ‘Project-new ‘ pull down menu, save it in a
separate directory(c:\ti\myprojects) with name cir conv.pjt.
Add the source files Circular Convolution.C.
to the project using ‘Projectadd files to project’ pull down menu.
Add the linker command file hello.cmd .
(Path: c:\ti\tutorial\dsk6713\hello1\hello.cmd)
Add the run time support library file rts6700.lib
(Path: c:\ti\c6000\cgtools\lib\rts6700.lib)
Compile the program using the ‘Project-compile’ pull down menu or by
clicking the shortcut icon on the left side of program window.
Build the program using the ‘Project-Build’ pull down menu or by
clicking the shortcut icon on the left side of program window.
Load the program(lconv.out) in program memory of DSP chip using the
‘File-load program’ pull down menu.
66
DSP LAB MANUAL VITS
OUTPUT FOR CIRCULAR CONVOLUTION:
enter the length of the first sequence4enter the length of the second sequence3enter the first sequence1 2 3 4enter the second sequence1 2 3the circular convolution is18 16 10 16
USING TMS320C6713 DSP PROCESSOR[RECTANGULER/TRIANGULAR/KAISER WINDOW]
AIM:
The aim of this laboratory exercise is to design and implement a Digital FIR Filter & observe its frequency response. In this experiment we design a simple FIR filter so as to stop or attenuate required band of frequencies components and pass the frequency components, which are outside the required band.
EQUIPMENTS:
TMS 320C6713 Kit.Oscilloscope & Function GeneratorRS232 Serial CablePower CordOperating System – Windows XPSoftware – CCStudio_v3.1
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.
(1)
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).
DSK6713_AIC23_Config config = {\ 0x0017, /* 0 DSK6713_AIC23_LEFTINVOL Leftline input channel volume */\ 0x0017, /* 1 DSK6713_AIC23_RIGHTINVOL Right line input channel volume*/\ 0x00d8, /* 2 DSK6713_AIC23_LEFTHPVOL Left channel headphone volume */\ 0x00d8, /* 3 DSK6713_AIC23_RIGHTHPVOL Right channel headphone volume */\ 0x0011, /* 4 DSK6713_AIC23_ANAPATH Analog audio path control */\ 0x0000, /* 5 DSK6713_AIC23_DIGPATH Digital audio path control */\ 0x0000, /* 6 DSK6713_AIC23_POWERDOWN Power down control */\ 0x0043, /* 7 DSK6713_AIC23_DIGIF Digital audio interface format */\ 0x0081, /* 8 DSK6713_AIC23_SAMPLERATE Sample rate control */\ 0x0001 /* 9 DSK6713_AIC23_DIGACT Digital interface activation */ \};
/* * main() - Main code routine, initializes BSL and generates tone */
void main(){ DSK6713_AIC23_CodecHandle hCodec; Uint32 l_input, r_input,l_output, r_output; /* Initialize the board support library, must be called first */ DSK6713_init(); /* Start the codec */ hCodec = DSK6713_AIC23_openCodec(0, &config); DSK6713_AIC23_setFreq(hCodec, 1);
70
DSP LAB MANUAL VITS
while(1) { /* Read a sample to the left channel */
while (!DSK6713_AIC23_read(hCodec, &l_input));
/* Read a sample to the right channel */while (!DSK6713_AIC23_read(hCodec, &r_input));
STEPS TO IMPLEMENT FILTERS Do diagnostic test Open DSKC6713 studio Open new project and give it a name. see whether the project is in C drive or D
drive Go to file. Say new and select DSPBIOS configuration Select DSK6713 or DSK6711 whatever is available (double click) Click on configuration window and save
Path: C:/ CCStudio_v3.1/myprojects/ project nameProjectAdd files to project Filter type as: configuration file(*.cdb)Configl—open
Observe under your project .cdb file under DSP and two configure files under generated files
Close configl window Go to file say new/source file/copy paste the filter program onto untitled window
The filter program is on desktop in programs for 6713 File/save/your project/file type is c source file and give a name to the file and say
save After saving the file add to the project/the file name i.e the saved file and say save. Add command file to the project
Project Add files to projectPath:CCStudio_v3.1/myprojects/ our project nameFiles of type: Linker command Files (*.cmd*,*.lcf*)
Add library file to the projectProject Add files to projectPath: CCStudio_v3.1\C6000\dsk6713\lib\dsk6713bsl.libFiles of type: Object and Library Files (*.o*,*.l*)
Add header file to the projectOpen configlcfg-c.c(under generated files)/Copy #include configlcfg.h (close) and paste it to our programRemove filter.h file
Go to project/ add files to projectPath: CCStudio_v3.1\C6000\dsk6713Files of type: All filesSelect include/copy the first two header files(dsk 6713,dsk6713_aic23) and paste in our project i.ePath: CCStudio_v3.1/myprojects/ our project name Save
Building and Running the Program (compile\ Build\ Load Program\ Run)Compile: Compile the program using the ‘Project-compile’ pull down menu or
by clicking the shortcut icon on the left side of program window.
72
DSP LAB MANUAL VITS
Build: Build the program using the ‘Project-Build’ pull down menu or by
clicking the shortcut icon on the left side of program window.
Load Program: Load the program in program memory of DSP chip using the
Run: Run the program using the ‘Debug-Run’ pull down menu or by clicking the
shortcut icon on the left side of program window.
73
DSP LAB MANUAL VITS
RESULT:FREQUENCY RESPONSE
High Pass FIR filter(Fc= 800Hz).
Low Pass FIR filter
(Fc=1000Hz)
74
DSP LAB MANUAL VITS
Model 2: ’c’ program to implement FIR Filter
#include<stdio.h>#include<math.h>#define pi 3.1415int n,N,c;float wr[64],wt[64];void main(){printf("\n enter no. of samples,N= :");scanf("%d",&N);printf("\n enter choice of window function: 1.rect 2. triang \n c= :");scanf("%d",&c);printf("\n elements of window function are:");switch(c){case 1:for(n=0;n<=N-1;n++){wr[n]=1;printf(" \n wr[%d]=%f",n,wr[n]);}break;case 2:for(n=0;n<=N-1;n++){wt[n]=1-(2*(float)n/(N-1));printf("\n wt[%d]=%f",n,wt[n]);}break;}}
75
DSP LAB MANUAL VITS
Output for FIR Filter Rectangular Windowenter no. of samples,N= :64 enter choice of window function:1.rect 2. triang c= 1
Graph Property GraphFIR FILTER RECTANGULAR WINDOW
Output for FIR Filter Triangular Window enter no. of samples,N= :64enter choice of window function:1.rect 2. triang c=2
Graph Property GraphFIR FILTER TRIANGULAR WINDOW
76
DSP LAB MANUAL VITS
EXP. No:14IIR FILTER DESIGN USING TMS320C6713 DSP PROCESSOR
AIM:
The aim of this laboratory exercise is to design and implement a Digital IIR Filter & observe its frequency response. In this experiment we design a simple IIR filter so as to stop or attenuate required band of frequencies components and pass the frequency components which are outside the required band
EQUIPMENTS: TMS 320C6713 Kit.Oscilloscope & Function GeneratorRS232 Serial CablePower CordOperating System – Windows XPSoftware – CCStudio_v3.1
INTRODUCTIONGENERAL CONSIDERATIONS:
In the design of frequency – selective filters, the desired filter characteristics are specified in the frequency domain in terms of the desired magnitude and phase response of the filter. In the filter design process, we determine the coefficients of a causal IIR filter that closely approximates the desired frequency response specifications.
IMPLEMENTATION OF DISCRETE-TIME SYSTEMS:
Discrete time Linear Time-Invariant (LTI) systems can be described completely by constant coefficient linear difference equations. Representing a system in terms of constant coefficient linear difference equation is it’s time domain characterization. In the design of a simple frequency–selective filter, we would take help of some basic implementation methods for realizations of LTI systems described by linear constant coefficient difference equation.
BACKGROUND CONCEPTS:
An Infinite impulse response (IIR) filter possesses an output response to an impulse which is of an infinite duration. The impulse response is "infinite" since there is feedback in the filter, that is if you put in an impulse ,then its output must produced for infinite duration of time. The IIR filter can realize both the poles and zeroes of a system because it has a rational transfer function, described by polynomials in z in both the numerator and the denominator:
77
DSP LAB MANUAL VITS
(1)
The difference equation for such a system is described by the following:
(2)
M and N are order of the two polynomials bk and ak are the filter coefficients. These filter coefficients are generated using FDS (Filter Design software or Digital Filter design package).
ALGORITHM TO IMPLEMENT:
We need to realize the Butter worth band pass IIR filter by implementing the difference
equation y[n] = b0x[n] + b1x[n-1]+b2x[n-2]-a1y[n-1]-a2y[n-2] where b0 – b2, a0-a2 are feed
forward and feedback word coefficients respectively [Assume 2nd order of filter].These
coefficients are calculated using MATLAB.A direct form I implementation approach is
taken.
Step 1 - Initialize the McBSP, the DSP board and the on board codec.
“Kindly refer the Topic Configuration of 6713Codec using BSL”
Step 2 - Initialize the discrete time system , that is , specify the initial conditions.
Generally zero initial conditions are assumed.
Step 3 - Take sampled data from codec while input is fed to DSP kit from the signal
generator. Since Codec is stereo , take average of input data read from left and right
channel . Store sampled data at a memory location.
Step 4 - Perform filter operation using above said difference equation and store
filter Output at a memory location .
Step 5 - Output the value to codec (left channel and right channel) and view the
/* Codec configuration settings */DSK6713_AIC23_Config config = { \ 0x0017, /* 0 DSK6713_AIC23_LEFTINVOL Left line input channel volume */ \ 0x0017, /* 1 DSK6713_AIC23_RIGHTINVOL Right line input channel volume */\ 0x00d8, /* 2 DSK6713_AIC23_LEFTHPVOL Left channel headphone volume */ \ 0x00d8, /* 3 DSK6713_AIC23_RIGHTHPVOL Right channel headphone volume */\ 0x0011, /* 4 DSK6713_AIC23_ANAPATH Analog audio path control */ \ 0x0000, /* 5 DSK6713_AIC23_DIGPATH Digital audio path control */ \ 0x0000, /* 6 DSK6713_AIC23_POWERDOWN Power down control */ \ 0x0043, /* 7 DSK6713_AIC23_DIGIF Digital audio interface format */ \ 0x0081, /* 8 DSK6713_AIC23_SAMPLERATE Sample rate control */ \ 0x0001 /* 9 DSK6713_AIC23_DIGACT Digital interface activation */ \};
/* * main() - Main code routine, initializes BSL and generates tone */void main(){ DSK6713_AIC23_CodecHandle hCodec; int l_input, r_input, l_output, r_output; /* Initialize the board support library, must be called first */ DSK6713_init(); /* Start the codec */ hCodec = DSK6713_AIC23_openCodec(0, &config); DSK6713_AIC23_setFreq(hCodec, 3);
80
DSP LAB MANUAL VITS
while(1) { /* Read a sample to the left channel */
while (!DSK6713_AIC23_read(hCodec, &l_input));
/* Read a sample to the right channel */while (!DSK6713_AIC23_read(hCodec, &r_input));
/* Send a sample to the left channel */ while (!DSK6713_AIC23_write(hCodec, l_output));
/* Send a sample to the right channel */ while (!DSK6713_AIC23_write(hCodec, r_output)); } /* Close the codec */ DSK6713_AIC23_closeCodec(hCodec);}
signed int IIR_FILTER(const signed int * h, signed int x1){
static signed int x[6] = { 0, 0, 0, 0, 0, 0 }; /* x(n), x(n-1), x(n-2). Must be static */ static signed int y[6] = { 0, 0, 0, 0, 0, 0 }; /* y(n), y(n-1), y(n-2). Must be static */ int temp=0;
/* Shuffle values along one place for next time */ y[2] = y[1]; /* y(n-2) = y(n-1) */ y[1] = y[0]; /* y(n-1) = y(n) */ x[2] = x[1]; /* x(n-2) = x(n-1) */ x[1] = x[0]; /* x(n-1) = x(n) */
/* temp is used as input next time through */ return (temp<<2); }PROCEDURE:
Switch on the DSP board.
Open the Code Composer Studio.
Create a new project
Project New (File Name. pjt , Eg: FIR.pjt)
Initialize on board codec.
“Kindly refer the Topic Configuration of 6713 Codec using BSL”
Add the given above ‘C’ source file to the current project (remove codec.c source
file from the project if you have already added).
Connect the speaker jack to the input of the CRO.
Build the program.
Load the generated object file(*.out) on to Target board.
Run the program using F5.
Observe the waveform that appears on the CRO screen.
82
DSP LAB MANUAL VITS
RESULT:
83
DSP LAB MANUAL VITS
MODEL 2:‘C’ PROGRAM TO IMPLEMENT IIR FILTER
#include<stdio.h>#include<math.h>int i,w,wc,c,N;float H[100];float mul(float, int);void main(){printf("\n enter order of filter ");scanf("%d",&N);printf("\n enter the cutoff freq ");scanf("%d",&wc);printf("\n enter the choice for IIR filter 1. LPF 2.HPF ");scanf("%d",&c);switch(c){case 1:for(w=0;w<100;w++){H[w]=1/sqrt(1+mul((w/(float)wc),2*N));printf("H[%d]=%f\n",w,H[w]);}break;case 2:for(w=0;w<=100;w++){H[w]=1/sqrt(1+mul((float)wc/w,2*N));printf("H[%d]=%f\n",w,H[w]);}break;}}float mul(float a,int x){for(i=0;i<x-1;i++)a*=a;return(a);}
84
DSP LAB MANUAL VITS
OUTPUT FOR IIR FILTER LPFenter order of filter 2 enter the cutoff freq 50 enter the choice for IIR filter 1. LPF 2.HPF: 1
Graph Property Graph
IIR FILTER LOW PASS FILTER
OUTPUT FOR IIR FILTER HPFenter order of filter 2 enter the cutoff freq 50 enter the choice for IIR filter 1. LPF 2.HPF: 2
Graph Property Graph
IIR FILTER HIGH PASS FILTER
85
DSP LAB MANUAL VITS
EXP.No:16N-POINT FAST FOURIER TRANSFORM (FFT)
AIM: To find the DFT of a sequence using FFT algorithm.
EQUIPMENTS:
TMS 320C6713 Kit.Oscilloscope & Function GeneratorRS232 Serial CablePower CordOperating System – Windows XPSoftware – CCStudio_v3.1
THEORY:The Fast Fourier Transform is useful to map the time-domain sequence into a continuous function of a frequency variable. The FFT of a sequence {x(n)} of length N is given by a Complex-valued sequence X (k).
The above equation is the mathematical representation of the DFT. As the number of computations involved in transforming a N point time domain signal into its corresponding frequency domain signal was found to be N2 complex multiplications, an alternative algorithm involving lesser number of computations is opted. When the sequence x(n) is divided into 2 sequences and the DFT performed separately, the resulting number of computations would be N2/2. (i.e.)
Consider x(2n) be the even sample sequences and x(2n+1) be the odd sample sequence derived form x(n).
(N/2)2multiplication’s
an other (N/2)2 multiplication's finally resulting in (N/2)2 + (N/2)2
=
86
DSP LAB MANUAL VITS
Further solving Eg. (2)
Dividing the sequence x(2n) into further 2 odd and even sequences would reduce the computations.
WN is the twiddle factor
Employing this equation, we deduce
(13)
87
DSP LAB MANUAL VITS
(14)
The time burden created by this large number of computations limits the usefulness of DFT in many applications. Tremendous efforts devoted to develop more efficient ways of computing DFT resulted in the above explained Fast Fourier Transform algorithm. This mathematical shortcut reduces the number of calculations the DFT requires drastically. The above mentioned radix-2 decimation in time FFT is employed for domain transformation.
Dividing the DFT into smaller DFTs is the basis of the FFT. A radix-2 FFT divides the DFT into two smaller DFTs, each of which is divided into smaller DFTs and so on, resulting in a combination of two-point DFTs. The Decimation -In-Time (DIT) FFT divides the input (time) sequence into two groups, one of even samples and the other of odd samples. N/2 point DFT are performed on the these sub-sequences and their outputs are combined to form the N point DFT.
The above shown mathematical representation forms the basis of N point FFT and is called the Butterfly Structure.
88
DSP LAB MANUAL VITS
ALGORITHM
Step 1 sample the input (N) of any desired frequency. Convert it to fixed-point format and scale the input to avoid overflow during manipulation.
Step 2 Declare four buffers namely real input, real exponent, imaginary exponent and imaginary input.
Step 3 Declare three counters for stage, group and butterfly.
Step 4 Implement the Fast Fourier Transform for the input signal.
Step 5 Store the output (Real and Imaginary) in the output buffer.
Step 6 Decrement the counter of butterfly. Repeat from the Step 4 until the counter reaches zero.
Step 7 If the butterfly counter is zero, modify the exponent value.
Step 8 Repeat from the Step 4 until the group counter reaches zero.
Step 9 If the group counter is zero, multiply the butterfly value by two and divide the group value by two.
Step 10 Repeat from the Step 4 until the stage counter reaches zero.
Step 11 Transmit the FFT output through line out port.
PROGRAMfft256.c#include <math.h> #define PTS 64 //# of points for FFT #define PI 3.14159265358979
typedef struct {float real,imag;} COMPLEX;
void FFT(COMPLEX *Y, int n); //FFT prototypefloat iobuffer[PTS]; //as input and output bufferfloat x1[PTS]; //intermediate buffer short i; //general purpose index variable short buffercount = 0; //number of new samples in iobuffer short flag = 0; //set to 1 by ISR when iobuffer full COMPLEX w[PTS]; //twiddle constants stored in w COMPLEX samples[PTS]; //primary working buffer
main(){ for (i = 0 ; i<PTS ; i++) // set up twiddle constants in w { w[i].real = cos(2*PI*i/(PTS*2.0)); //Re component of twiddle constants w[i].imag =-sin(2*PI*i/(PTS*2.0)); //Im component of twiddle constants } for (i = 0 ; i < PTS ; i++) //swap buffers {
89
DSP LAB MANUAL VITS
iobuffer[i] = sin(2*PI*10*i/64.0);/*10- > freq, 64 -> sampling freq*/ samples[i].real=0.0; samples[i].imag=0.0; } for (i = 0 ; i < PTS ; i++) //swap buffers { samples[i].real=iobuffer[i]; //buffer with new data } for (i = 0 ; i < PTS ; i++) samples[i].imag = 0.0; //imag components = 0
FFT(samples,PTS); //call function FFT.c
for (i = 0 ; i < PTS ; i++) //compute magnitude { x1[i] = sqrt(samples[i].real*samples[i].real
+ samples[i].imag*samples[i].imag); } } //end of main
fft.c:
#define PTS 64 //# of points for FFTtypedef struct {float real,imag;} COMPLEX;extern COMPLEX w[PTS]; //twiddle constants stored in w
void FFT(COMPLEX *Y, int N) //input sample array, # of points { COMPLEX temp1,temp2; //temporary storage variables int i,j,k; //loop counter variables int upper_leg, lower_leg; //index of upper/lower butterfly leg int leg_diff; //difference between upper/lower leg int num_stages = 0; //number of FFT stages (iterations) int index, step; //index/step through twiddle constant i = 1; //log(base2) of N points= # of stages do { num_stages +=1; i = i*2; }while (i!=N); leg_diff = N/2; //difference between upper&lower legs step = (PTS*2)/N; //step between values in twiddle.h for (i = 0;i < num_stages; i++) //for N-point FFT {
90
DSP LAB MANUAL VITS
index = 0; for (j = 0; j < leg_diff; j++) { for (upper_leg = j; upper_leg < N; upper_leg += (2*leg_diff)) { lower_leg = upper_leg+leg_diff;
+temp2.imag*(w[index]).real; (Y[upper_leg]).real = temp1.real; (Y[upper_leg]).imag = temp1.imag; } index += step; } leg_diff = leg_diff/2; step *= 2; } j = 0; for (i = 1; i < (N-1); i++) //bit reversal for resequencing data { k = N/2; while (k <= j) { j = j - k; k = k/2; } j = j + k; if (i<j) { temp1.real = (Y[j]).real; temp1.imag = (Y[j]).imag; (Y[j]).real = (Y[i]).real; (Y[j]).imag = (Y[i]).imag; (Y[i]).real = temp1.real; (Y[i]).imag = temp1.imag; } } return;}
91
DSP LAB MANUAL VITS
PROCEDURE:
Open Code Composer Studio, make sure the DSP kit is turned on.
Start a new project using ‘Project-new ‘ pull down menu, save it in a
separate directory(c:\ti\myprojects) with name “FFT.pjt”.
Add the source files “fft256.c“ and “fft.C” in the project using
‘Projectadd files to project’ pull down menu.
Add the linker command file “hello.cmd”.
Add the rts file “rts6700.lib” .
Compile the program using the ‘Project-compile’ pull down menu or by
clicking the shortcut icon on the left side of program window.
Load the program in program memory of DSP chip using the ‘File-load program’
pull down menu.
Run the program and observe output using graph utility.
92
DSP LAB MANUAL VITS
GRAPHS
Graph Property GraphFFT INPUT
Graph Property GraphFFT OUTPUT
93
DSP LAB MANUAL VITS
EXP.No:16TO COMPUTE POWER DENSITY SPECTRUM OF
A SEQUENCE [USING TMS320C6713 DSP PROCESSOR]
AIM: To find the PSD of a sequence.
EQUIPMENTS:
TMS 320C6713 Kit.Oscilloscope & Function GeneratorRS232 Serial CablePower CordOperating System: Windows XPSoftware: CCStudio_v3.1
INTRODUCTION:The total or the average power in a signal is often not of as great an interest. We are most
often interested in the PSD or the Power Spectrum. We often want to see is how the input
power has been redistributed by the channel and in this frequency-based redistribution of
power is where most of the interesting information lies. The total area under the Power
Spectrum or PSD is equal to the total avg. power of the signal. The PSD is an even
function of frequency or in other words.
To compute PSD:
The value of the auto-correlation function at zero-time equals the total power in the
signal. To compute PSD we compute the auto-correlation of the signal and then take its
FFT. The auto-correlation function and PSD are a Fourier transform pair. (Another
estimation method called “period gram” uses sampled FFT to compute the PSD.).
E.g.: For a process x(n) correlation is defined as:
Power Spectral Density is a Fourier transform of the auto correlation.
94
DSP LAB MANUAL VITS
ALGORITHM TO IMPLEMENT PSD:
Step 1 - Select no. of points for FFT(Eg: 64).
Step 2 – Generate a sine wave of frequency ‘f ‘ (eg: 10 Hz with a sampling rate =
No. of Points of FFT(eg. 64)) using math library function.
Step 3 - Compute the Auto Correlation of Sine wave .
Step4 - Take output of auto correlation, apply FFT algorithm .
Step 4 - Use Graph option to view the PSD.
Step 5 - Repeat Step-1 to 4 for different no. of points & frequencies.
‘C’ PROGRAM TO IMPLEMENT PSD: PSD.c:/************************************************************* * FILENAME * Non_real_time_PSD.c * DESCRIPTION * Program to Compute Non real time PSD * using the TMS320C6711 DSK.
*************************************************************** * DESCRIPTION * Number of points for FFT (PTS) * x --> Sine Wave Co-Efficients * iobuffer --> Out put of Auto Correlation. * x1 --> use in graph window to view PSD /*===========================================================*/
#include <math.h> #define PTS 128 //# of points for FFT #define PI 3.14159265358979
typedef struct {float real,imag;} COMPLEX;
void FFT(COMPLEX *Y, int n); //FFT prototypefloat iobuffer[PTS]; //as input and output bufferfloat x1[PTS],x[PTS]; //intermediate buffer short i; //general purpose index variable short buffercount = 0; //number of new samples in iobuffer short flag = 0; //set to 1 by ISR when iobuffer full float y[128];COMPLEX w[PTS]; //twiddle constants stored in w COMPLEX samples[PTS]; //primary working buffer
main(){
95
DSP LAB MANUAL VITS
float j,sum=0.0 ; int n,k,i,a; for (i = 0 ; i<PTS ; i++) // set up twiddle constants in w { w[i].real = cos(2*PI*i/(PTS*2.0));
/*Re component of twiddle constants*/ w[i].imag =-sin(2*PI*i/(PTS*2.0)); /*Im component of twiddle constants*/ } /****************Input Signal X(n) *************************/ for(i=0,j=0;i<PTS;i++) { x[i] = sin(2*PI*5*i/PTS); // Signal x(Fs)=sin(2*pi*f*i/Fs); samples[i].real=0.0; samples[i].imag=0.0; } /********************Auto Correlation of X(n)=R(t) ***********/ for(n=0;n<PTS;n++) { sum=0; for(k=0;k<PTS-n;k++) { sum=sum+(x[k]*x[n+k]); // Auto Correlation R(t) } iobuffer[n] = sum; } /********************** FFT of R(t) ***********************/
for (i = 0 ; i < PTS ; i++) //swap buffers { samples[i].real=iobuffer[i]; //buffer with new data } for (i = 0 ; i < PTS ; i++) samples[i].imag = 0.0; //imag components = 0
FFT(samples,PTS); //call function FFT.c
/******************** PSD ********************/ for (i = 0 ; i < PTS ; i++) //compute magnitude
96
DSP LAB MANUAL VITS
{ x1[i] = sqrt(samples[i].real*samples[i].real
+ samples[i].imag*samples[i].imag); }
} //end of main
FFT.c:
#define PTS 128 //# of points for FFT
typedef struct {float real,imag;} COMPLEX;
extern COMPLEX w[PTS]; //twiddle constants stored in w
void FFT(COMPLEX *Y, int N) //input sample array, # of points { COMPLEX temp1,temp2; //temporary storage variables int i,j,k; //loop counter variables int upper_leg, lower_leg; //indexof upper/lower butterfly leg int leg_diff; //difference between upper/lower leg int num_stages = 0; //number of FFT stages (iterations) int index, step; //index/step through twiddle constant i = 1; //log(base2) of N points= # of stages do { num_stages +=1; i = i*2; }while (i!=N); leg_diff = N/2; //difference between upper&lower legs step = (PTS*2)/N; //step between values in twiddle.h// 512 for (i = 0;i < num_stages; i++) //for N-point FFT { index = 0; for (j = 0; j < leg_diff; j++) {
Figure : Magnitude response of the FIR bandpass filter designed by the window method.
104
DSP LAB MANUAL VITS
EXP.NO: 2TO FIND FREQEUNCY RESPONSE OF BAND PASS FILTER
UISNG IIR
Aim: To find frequency response of Band pass filter using IIR.
EQUIPMENTS:MATLAB Software, CCS Studio
PROGRAM:MATLAB CODE%steps of designing a digital filter implied in the called functions:% step 1: estimate the minimum order of the filter from specifications% step 2: get analog, pre-warped frequencies% step 3: convert to low-pass prototype estimate% step 4: Get n-th order analog lowpass prototype with desired filter characters% step 5: Transform to lowpass, bandpass, highpass, or bandstop of desired Wn% step 6: Use Bilinear transformation to find discrete equivalent:
clear all; format;Fs = 100000; %sampling frequency(Hz). Wp = [20500 23500]/(Fs/2); %passband edge frequency normalised by Fs/2. Ws = [19000 25000]/(Fs/2); %stopband edge frewquency normalised by Fs/2.Rp = 0.25; %passband attenuation in dB. Rs = 45; %stopband attenuation in Dbctype = 4; %character of filterftype = 'bandpass'; %type of filterif ctype==1 [n,Wn]=buttord(Wp,Ws,Rp,Rs); [b,a]=butter(n,Wn,ftype); elseif ctype==2 [n,Wn]=cheb1ord(Wp,Ws,Rp,Rs); [b,a]=cheby1(n,Rp,Wn,ftype); elseif ctype==3 [n,Wn]=cheb2ord(Wp,Ws,Rp,Rs); [b,a]=cheby2(n,Rs,Wn,ftype); elseif ctype==4
105
DSP LAB MANUAL VITS
[n,Wn]=ellipord(Wp,Ws,Rp,Rs); [b,a]=ellip(n,Rp,Rs,Wn,ftype); end
%Output the resultdisp('Numerator coefficients (in descending powers of z):'); disp(b);disp('Denominator coefficients (in descending powers of z):'); disp(a);freqz(b,a,1024,Fs);
=======================================================function Impinv %Impulse invariance method of anolog-to-digital filter conversion %a,b -- s-plane coefficients%az,bz -- digital filter coefficients
clear all;b = 1; a = [1.84496 1.920675 1]; [bz,az]=impinvar(b,a) %get z-plane coefficients using impulse Inv.freqz(bz,az,1024);
C CODE#include <math.h>#include <string.h>#include "st_i.h"
/* Private data for Bandpass effect */typedef struct bandstuff {
float center;float width;double A, B, C;double out1, out2;short noise;/* 50 bytes of data, 52 bytes long for allocation purposes. */
} *band_t;
/* * Process options */int st_band_getopts(eff_t effp, int n, char **argv) {
/* * Do anything required when you stop reading samples. * Don't close input file! */int st_band_stop(eff_t effp){
return (ST_SUCCESS); /* nothing to do */}
108
DSP LAB MANUAL VITS
RESULT
109
DSP LAB MANUAL VITS
EXP.NO: 3TO FIND FREQEUNCY RESPONSE OF BAND REJECT FILTER
UISNG FIR
Aim: To find frequency response of Band reject filter using FIR.
EQUIPMENTS:MATLAB Software
PROGRAM:MATLAB CODE%IIR Chebyshev type-I bandstop filterclc;clear all;close all;rp=input('enter the passband ripple ');rs=input('enter the stopband ripple ');wp=input('enter the passband frequency ');ws=input('enter the stopband frequency ');fs=input('enter the sampling frequency ');w1=2*wp/fs;w2=2*ws/fs;[n]=cheb1ord(w1,w2,rp,rs);wn=[w1 w2];[b,a]=cheby1(n,rp,wn,'stop');w=0:0.01:pi;[h,om]=freqz(b,a,w);m=20*log10(abs(h));an=angle(h);subplot(2,1,1);plot(om/pi,m);ylabel('Gain in dB--->');xlabel('Normalised frequency--->');title('CHEBYSHEV TYPE-I BAND STOP MAGNITUDE RESPONSE');subplot(2,1,2);plot(om/pi,an);ylabel('Phase in radians--->');xlabel('Normalised frequency--->');title('CHEBYSHEV TYPE-I BAND STOP PHASE RESPONSE');SAMPLE INPUTSenter the passband ripple 0.25enter the stopband ripple 40enter the passband frequency 2500
110
DSP LAB MANUAL VITS
enter the stopband frequency 2750enter the sampling frequency 7000
RESULT
111
DSP LAB MANUAL VITS
EXP.NO: 4TO FIND FREQEUNCY RESPONSE OF BAND REJECT FILTER
UISNG IIR
Aim: To find frequency response of Band reject filter using IIR.
EQUIPMENTS:MATLAB Software
PROGRAM:MATLAB CODE
%IIR Butterworth bandstop filter clc;clear all;close all; rp=input('enter the passband ripple '); rs=input('enter the stopband ripple '); wp=input('enter the passband frequency '); ws=input('enter the stopband frequency '); fs=input('enter the sampling frequency '); w1=2*wp/fs; w2=2*ws/fs; [n,wn]=buttord(w1,w2,rp,rs); wn=[w1,w2]; [b,a]=butter(n,wn,'stop'); w=0:0.01:pi; [h,om]=freqz(b,a,w); m=20*log10(abs(h)); an=angle(h); subplot(2,1,1);
112
DSP LAB MANUAL VITS
plot(om/pi,m); ylabel('Gain in dB--->'); xlabel('Normalised frequency--->'); title('BUTTERWORTH BANDSTOP MAGNITUDE RESPONSE'); subplot(2,1,2); plot(om/pi,an); ylabel('Phase in radians--->'); xlabel('Normalised frequency--->'); title('BUTTERWORTH BANDSTOP PHASE RESPONSE'); SAMPLE INPUTS enter the passband ripple .4 enter the stopband ripple 46 enter the passband frequency 1100 enter the stopband frequency 2200 enter the sampling frequency 6000 RESULT
113
DSP LAB MANUAL VITS
VIVA QUESTIONS
GENERATION OF SINUSOIDAL SIGNAL QUESTIONS
1. what is the difference between sin & cos signals?
2. What is meant by signal?
3. What is the difference between time domain & frequency domain signal?
4. What is the difference between periodic & a periodic signal.
5. What is the difference between orthogonal and orthonormal signals?
6. What is the need for Fourier series & Fourier transform?
7. What is the difference between discrete & digital signals?
8. What is the difference between even signal & odd signal?
9. What is the difference between power signal & energy signal?
10. What is the difference between amplitude scaling & time scaling of a signal?
11. What is the difference between deterministic & random signal?
LINEAR CONVOLUTION QUESTIONS
1. What is the requirement for convolution?.
2. What is the difference between convolution & correlation?
3. What is meant by impulse response?
4. Is it possible to represent any discrete time signal in terms of impulses? If yes,
represent by using example.
5. Draw the h(2n-k) & h(n-2k) for the following sequence
h(n) = { 4 3 2 1} assume (i) k= 3 (ii) k =5.
6. Write the expressions for LTI system convolution formula & causal LTI system
convolution formula.
7. What us the length of linear convolution if length of input & impulse responses
are N1 & N2 respectively?
8. What is the difference between continuous and discrete convolution?
114
DSP LAB MANUAL VITS
CIRCULAR CONVOLUTION QUESTIONS
1. Why we need circular convolution?
2. What is the difference between circular & linear convolution?
3. What is the length of output sequence after circular convolution if the lengths of
input & impulse responses are M1 & M2 respectively?
4. State the circular convolution property of DFT?
5. Where we required convolution property?
6. What does zero padding mean? Where we required this concept?
7. What is difference between linear shifting & circular shifting of signal? Show
with example.
8. What is difference between linear & circular folding of signal? Show with
example.
9. What is the advantage with sectioned convolution?.
FAST FOURIER TRANSFORMS QUESTION
1. What is the difference between continuous time & discrete time Fourier transform?
2. What is the condition for convergence of Fourier transform?
3. What is the difference between discrete Time Fourier Transform (DTFT)& DFT?
4. What is the difference between Z transform & DFT?
5. State convolution property of the DFT? Where we could use the convolution
property?
6. State Parseval’s theorem.?
7. State correlation property of the DFT.?
8. What is the difference between radix 2 & radix4 FFT algorithms?
9. Why we need FFT?.
10. What is the difference between decimation in time (DIT FFT) & Decimation in
frequency(DIFFFT) algorithms?
11. What is meant by ‘in-place’ computation in DIF & DIF algorithms?
12. Which properties are used in FFT to reduce no of computations?
115
DSP LAB MANUAL VITS
FIR FILTER QUESTIONS
1. What are the advantages of FIR as compared to IIR?
2. How many types of FIR design methods are used in real time?.
3. What is meant by Gibbs Phenomenon? Where we found such type of effect in FIR
filters?
4. What are the advantages& disadvantages of Rectangular window FIR filter as
compared to remaining window techniques?
5. Which window technique having less peak amplitude of side lobe as compared to
all?
6. What do you understand by linear phase responce?
7. To design all types of filters what would be the expected impulse response?
8. What are the properties of FIR filter?.
9. How the zeros in FIR filter is located?
10. What are the desirable characteristics of the window?
11. What are the specifications required to design filter
IIR FILTER QUESTIONS
1. What is meant by IIR filter?
2. What is the difference between recursive & non-recursive systems?
3. Write the difference equation for IIR system.
4. What are the mapping techniques in IIR filter design? Discuss the advantage &
disadvantages of them.
5. What are IIR analog filters? What are the advantages & disadvantages of them?
6. What is the disadvantage in impulse invariance method?
7. What does warping effect mean? Where we found this effect? How can we
eliminate warping effect
8. Explain the pole mapping procedure of Impulse invariant & bilinear
transformation method.
9. For given same specification which difference we found in Butter worth &
Tchebyshev filter.
10. What is the difference between type I & type II Tchebyshev filters?.
116
DSP LAB MANUAL VITS
11. Where the poles are located for Butter worth & Tchedbyshev filters?
12. What is meant by spectral transformation?
13. Why we need spectral transformation in IIR filter?
POWER SPECTRUM DENSITY QUESTION
1. What is the difference between correlation & auto correlation function?
2. What is the difference between PSD & ESD?
3. What is the unit for energy density spectrum?
4. What is the formula for PSD of a function?
5. “Same power density spectrum signals always have same magnitude & phase
spectrums” Is above statement true (or) False: Justify your answer.
6. If we know the impulse response of the system, the How can you find output
signal power density from the input signal?
7. What is the unit for power density spectrum?
8. What is the relation between auto correlation & PSD of a function?
DSP PROCESSORS QUESTIONS
1. How many types of DSP processors are available in the market?
2. TMS 320C6X, ‘C’ stands for what?
3. What are the features of TMS 320C6X processor?
4. What is meant by VLIW architecture? Why we required in DSP processor?
5. How many functional units are in TMS 320C6X DSP processor?
6. What is meant by Circular addressing mode how is it useful for DSP?
7. Which instruction is used to move 16 bit constant in to the upper bits of a register?
8. What is the difference between Von Neumann architecture & Harvard
architecture?
9. Which architecture is used in DSP processor?
10. How many instructions can we execute per cycle in TMS320C6X DSP processor?
11. What are the applications for the TMS320 DSP’s?
12. Which soft ware tool is required to compile and run the DSP assembly program ?
13. What is the difference between full version Code composer studio &DSK CCS?