INSTR F243: Signal & Systems

Instructor: Pawan k.Ajmera

MATLAB ASSIGNMENT

RAJ AARYAN(2017B4A80753P)

Due Friday, May 8, 2020

Contents

1 Question1 2

1.1 CODE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2

1.2 IMAGE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2 Question 2 4

2.1 CODE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2.2 IMAGE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

3 Question3 6

3.1 CODE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

3.2 IMAGE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

4 Question 4 8

4.1 CODE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

4.2 IMAGE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

1
1 Question1

1.1 CODE

1
2 % Specification in problem: Amplitude= 2
3 % Duration= 2 second
4 % frequency = 753 Hertz
5 % Sampling frequency =20kHz
6 clc;
7 close all;
8 clear all;
9 fs = 20000; % Sampling frequency
10 t = 0:1/fs:2; % Duration
11 a = 2; % Amplitude
12 f = 753; % Frequency of signal
13 y = a*cos(2*pi*f*(t)); % Cosine signal
14 % Plot 1−
15 subplot(2,1,1);
16 plot(t,y);
17 axis([0 2/f −2 2]);
18 xlabel('Time (in sec)');
19 ylabel('Amplitude');
20 title('1(a) Cosine Signal − First two cycles');
21 % Delayed signal
22 y2 = a*cos(2*pi*f*(t−0.1));
23 % Plot 2
24 subplot(2,1,2);
25 plot(t,y2);
26 axis([0 2/f −2 2]);
27 xlabel('Time (in sec)');
28 ylabel('Amplitude');
29 title('1(b) Delayed Signal by 100msec − First two cycles');

2
1.2 IMAGE

1(a) Cosine Signal - First two cycles

2

1
Amplitude

-1

-2
0 0.0005 0.001 0.0015 0.002 0.0025
Time (in sec)

1(b) Delayed Signal by 100msec - First two cycles

2

1
Amplitude

-1

-2
0.1 0.1005 0.101 0.1015 0.102 0.1025
Time (in sec)

Figure 1: Caption

3
2 Question 2

2.1 CODE

1 % Specification required : a= 7(Maximum of 7,5 and 3)

2 close all;
3 clear all;
4 clc;
5
6 n1 = 1:1:5;
7 x = [1,1,1,1,1 ];
8 subplot (2, 2, 1);
9 stem (n1,x);
10 xlabel('Samples');
11 ylabel ('Ampltiude');
12 title ('x[n]');
13
14
15 n2 = −7:1:4 ;
16 h1 = [1 1 1 1 1 1 1 1 1 1 1 1 ];
17 subplot (2, 2, 2);
18 stem (n2 ,h1);
19 xlabel('Samples');
20 ylabel ('Ampltiude');
21 title ('h1[n]');
22
23
24 n3 = −6:1:9; %Sum of left sided limit and sum of right Sided Limits
25 y1 = conv(x,h1);
26 subplot (2, 2 ,3);
27 stem (n3,y1)
28 xlabel('Samples');
29 ylabel ('Ampltiude');
30 title ('x[n]*h1[n]');
31
32
33 n3 = −5:1:10; %Sum of left sided limit and sum of right Sided Limits
34 h2 = [1 ];
35 subplot (2, 2, 4);
36 y = conv(h2,y1);
37 stem(n3,y);
38 xlabel('Samples');
39 ylabel ('Ampltiude');
40 title ('Y[n]');

4
2.2 IMAGE

x[n] h1[n]
1 1

0.8 0.8
Ampltiude

Ampltiude
0.6 0.6

0.4 0.4

0.2 0.2

0 0
1 2 3 4 5 -8 -6 -4 -2 0 2 4
Samples Samples

x[n]*h1[n] Y[n]
5 5

4 4
Ampltiude

Ampltiude
3 3

2 2

1 1

0 0
-10 -5 0 5 10 -10 -5 0 5 10
Samples Samples

Figure 2: Caption

5
3 Question3

3.1 CODE

1 %
2 % Specification Mentioned : a=7 b=5 c=3
3
4 close all;
5 clear all;
6 clc;
7
8 num = [7 5 3]; % numerator of transfer function
9 den = [3 5 7]; % denominator of transfer function
10
11 H = tf(num,den) % to find zeros, poles and gain
12
13 pzmap(H); % to plot zeros and poles
14 grid on ;
15 [z,p]= tf2zp(num,den);
16 % stability will be there if all poles lies behind jw axis.

6
3.2 IMAGE

Figure 3: Caption

7
4 Question 4

4.1 CODE

1 close all;
2 clear all;
3 clc;
4
5
6 dur= 2; %Duration of plot is 2 second
7 Fs = 20000; % Sampling frequency is 20kHz
8 dur samples= 0:(1/ Fs):2 ;
9 freq1= 753;
10 freq2= 11753;
11 x = 2*[cos(2*pi*freq1* dur samples)+cos(2*pi*freq2* dur samples)];
12 fft x= fft(x,Fs);
13 frequency = [0:length(fft x)−1]* Fs/length(fft x);
14 plot(frequency,abs(fft x));
15 axis([0 10000 0 25000 ]);
16 xlabel('Frequency(Hz)');
17 ylabel ('Ampltiude');
18 title ('FFT −Magntitude Response');

8
4.2 IMAGE
FFT -Magntitude Response
25000

20000

15000
Ampltiude

10000

5000

0
0 2000 4000 6000 8000 10000
Frequency(Hz)

Figure 4: Caption

COMMENT