PROGRAM NO 1(a) Date: 01st Aug,2018

AIM:- To represent basic signals unit impulse.

PROGRAM:-
%program for the generation of unit impulse signal
clc;
clear all;
close all;
t=-2:1:2;
y=[zeros(1,2),ones(1,1),zeros(1,2)];
subplot(2,2,1);
stem(t,y);
ylabel('amplitude-->');
xlabel('(a)n-->');

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 1(b) Date: 01st Aug,2018

AIM:- To represent basic signals unit step.

PROGRAM:-
%program for the generation of unit step signal [u(n)-u(n-N)]
clc;
clear all;
close all;
n=input('enter the N value');
t=0:1:n-1;
y1=ones(1,n);
subplot(2,2,2);
stem(t,y1);
ylabel('amplitude-->');
xlabel('(b)n-->');

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 1(c) Date:01st Aug,2018

AIM:- To represent basic signals unit ramp.

PROGRAM:-
%program for the generation of unit ramp sequence.

clc;

clear all;
close all;
n1=input('enter the length of ramp sequence');
t=0:n1;
subplot(2,2,3);
stem(t,t);
ylabel('amplitude-->');
xlabel('(b)n-->');

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 1(d) Date: 01st Aug,2018
AIM:- To represent basic signals exponential sequence.

PROGRAM:-

%program for the generation of exponential sequence.

clc;
clear all;
close all;
n2=input('enter the length of exponential sequence');
t=0:n2;
a=input('enter the (a)value');
y2=exp(a*t);
subplot(2,2,4);
stem(t,y2);
ylabel('amplitude-->');
xlabel('(b)n-->');

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 1(e) Date: 01st Aug,2018

AIM:- To represent sine sequence

PROGRAM:-
%program for the generation of sin sequence.
t=0:.01:pi;
y=sin(2*pi*t);
figure(2);
subplot(2,1,1);
plot(t,y);
ylabel('AMPLITUDE-->');
xlabel('(a)n-->');

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 1(f) Date: 01st Aug,2018

AIM:- To represent cosine sequence.

PROGRAM:-
%program for the generation of cos sequence.
t=0:.01:pi;
y=cos(2*pi*t);
figure(2);
subplot(2,1,1);
plot(t,y);
ylabel('AMPLITUDE-->');
xlabel('(a)n-->');

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 2 Date: 08th Aug,2018
AIM:- To develop program for discrete convolution.

PROGRAM:-
%program for linear convolution of the sequence X=[1,2] and h=[1,2,4]
clc;
clear all;
close all;
x=input('enter the 1st sequence');
h=input('enter the 2nd sequence');
y=conv(x,h);
figure;
subplot(3,1,1);
stem(x);
ylabel('amplitude-->');
xlabel('(a)n-->');
subplot(3,1,2)
stem(h);
ylabel('amplitude-->');
xlabel('(b)n-->');
subplot(3,1,3)
stem(fliplr(y));
ylabel('amplitude-->');
xlabel('(c)n-->');
disp('the resulant signal is');y

OUTPUT:-

Shivam Pandey
(15ECE27)
Shivam Pandey
(15ECE27)
 PROGRAM NO 3 Date: 22nd Aug,2018

AIM:- To develop program for discrete correlation.

PROGRAM:-
%program for cross correlation of the sequence X=[1,2] and h=[1,2,4]
clc;
clear all;
close all;
x=input('enter the 1st sequence');
h=input('enter the 2nd sequence');
y=xcorr(x,h);
figure;
subplot(3,1,1);
stem(x);
ylabel('amplitude-->');
xlabel('(a)n-->');
subplot(3,1,2)
stem(h);
ylabel('amplitude-->');
xlabel('(b)n-->');
subplot(3,1,3)
stem(fliplr(y));
ylabel('amplitude-->');
xlabel('(c)n-->');
disp('the resulant signal is');y

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 4 Date: 29th Aug,2018

AIM:-To understand stability test.

PROGRAM:-
%program for stability test
clc;
clear all;
close all;
b=input('enter the denominator coefficients of the filter');
k=poly(b);
knew=fliplr(k);
s=all(abs(knew));
if(s==1)
disp('"stable system"');
else
disp("'non-stable system"');
end

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO 5 Date: 12th Sep,2018

AIM:- To understand sampling theorem.

PROGRAM:-
%program to understand sapling theorem
f1=1/128;
f2=5/128;
n=0:255;
fc=50/128;
x=cos(2*pi*f1*n)+cos(2*pi*f2*n);
xa=cos(2*pi*fc*n);
xamp=x.*xa;
subplot(2,2,1);
plot(n,x);
title('x(n)');
xlabel('Sn-->');
ylabel('amplitude');
subplot(2,2,2);
plot(n,xa);
title('xa(n)');
xlabel('n-->');
ylabel ('amplitude');
subplot(2,2,3);
plot(n,xamp);
xlabel('n-->');
ylabel('amplitude');

OUTPUT:-

Shivam Pandey
(15ECE27)
 PROGRAM NO-6 Date: 26th Sep,2018

AIM:-To design analog butterworth low pass filter.

PROGRAM:-
%program for the design of butterworth low pass filter
clc;
format long
rs=input('enter the stopband ripple...');
rp=input('enter the passband ripple...');
ws=input('enter the stopband freq...');
wp=input ('enter the passband freq...');
fs=input('enter the sampling freq...');
w1=2*wp/fs;
w2=2*ws/fs;
[n,wn]=buttord(w1,w2,rp,rs);
[z,p,k]=butter(n,wn);
[b,a]=zp2tf(z,p,k);
[b,a]=butter(n,wn);
w=0:0.01:pi;
[h,om]=freqs(b,a,w);
m=20*log10(abs(h));
an=angle(h);
subplot(2,1,1);
plot(om/pi,m);
ylabel('gainindB-->');
xlabel('(a)normalised frequency-->');
subplot(2,1,2);
plot(om/pi,an);
xlabel('(b)normalised frequency-->');
ylabel('phase in radians-->');

OUTPUT:-
enter the stopband ripple...60

enter the stopband freq...3000

enter the passband freq...1500

enter the sampling freq…7000

Shivam Pandey
(15ECE27)
Shivam Pandey
(15ECE27)
 PROGRAM NO-7 Date: 03rd Oct,2018

AIM:-Program for design of butterworth high pass analog filter.

PROGRAM:-
%program for the design of butterworth high pass filter
clc;
format long
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');
[b,a]=butter(n,wn,'high','s');
w=0:0.01:pi;
[h,om]=freqs(b,a,w);
m=20*log10(abs(h));
an=angle(h);
subplot(2,1,1);
plot(om/pi,m);
ylabel('gainindB-->');
xlabel('(a)normalised frequency-->');
subplot(2,1,2);
plot(om/pi,an);
xlabel('(b)normalised frequency-->');
ylabel('phase in radians-->');

OUTPUT:-
enter the passband ripple...0.2

enter the stopband ripple...40

enter the passband freq...2000

enter the stopband freq...3500

enter the sampling freq...8000

Shivam Pandey
(15ECE27)
Shivam Pandey
(15ECE27)
 PROGRAM NO :-8 Date: 03rd Oct,2018

Aim:-Program for the design of butter worth band pass analog filter.

PROGRAM:-
%program for the design of butterworth band pass filter
clc;
format long
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]=buttord(w1,w2,rp,rs,'s');
wn=[w1 w2];
[b,a]=butter(n,wn,'bandpass','s');
w=0:0.01:pi;
[h,om]=freqs(b,a,w);
m=20*log10(abs(h));
an=angle(h);
subplot(2,1,1);
plot(om/pi,m);
ylabel('gainindB-->');
xlabel('(a)normalised frequency-->');
subplot(2,1,2);
plot(om/pi,an);
xlabel('(b)normalised frequency-->');
ylabel('phase in radians-->');

OUTPUT:-
enter the passband ripple...0.36
enter the stopband ripple...36
enter the passband freq...1500
enter the stopband freq...2000
enter the sampling freq...6000

Shivam Pandey
(15ECE27)
Shivam Pandey
(15ECE27)
 PROGRAM NO:9 Date: 03rd Oct,2018

AIM:-Program for the design of butter worth band stop analog filter.

PROGRAM:-
%program for the design of butterworth band stop filter
clc;
format long
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]=buttord(w1,w2,rp,rs,'s');
wn=[w1 w2];
[b,a]=butter(n,wn,'stop','s');
w=0:0.01:pi;
[h,om]=freqs(b,a,w);
m=20*log10(abs(h));
an=angle(h);
subplot(2,1,1);
plot(om/pi,m);
ylabel('gain in dB-->');
xlabel('(a)normalised frequency-->');
subplot(2,1,2);
plot(om/pi,an);
xlabel('(b)normalised frequency-->');
ylabel('phase in radians-->');

OUTPUT:-
enter the passband ripple...0.28

enter the stopband ripple...28

enter the passband freq...1000

enter the stopband freq...14000

enter the sampling freq...5000

Shivam Pandey
(15ECE27)
Shivam Pandey
(15ECE27)