Scribd
Upload
7 views4 pages

Transfer Function Analysis and Stability

The document contains MATLAB code for analyzing different linear time-invariant systems. It calculates the poles and zeros of transfer functions, plots their pole-zero maps, and converts transfer functions between different representations. Several examples are provided to demonstrate these analyses and conversions between transfer function parameter forms.

Uploaded by

3126328310
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
7 views4 pages

Transfer Function Analysis and Stability

The document contains MATLAB code for analyzing different linear time-invariant systems. It calculates the poles and zeros of transfer functions, plots their pole-zero maps, and converts transfer functions between different representations. Several examples are provided to demonstrate these analyses and conversions between transfer function parameter forms.

Uploaded by

3126328310
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd

LAB TASK# 3

Q# 1.
b=[0 -1.618 1]
a= [1 -1.5161 0.878]
roots(a)
roots(b)
zplane(b,a);

1
0.8
0.6

Imaginary Part

0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
-1

-0.5

0
Real Part

0.5

Yes this system is stable because poles are inside the unit circle.

Q # 2.
b=[1 2 1];
a= [1 -1 0.3561];
freqz(b,a);

Name: SHAHBAZ IQBAL

Roll No: BEES/S11/0150

Page 1

Magnitude (dB)

50

-50

-100
0

0.1

0.2

0.3
0.4
0.5
0.6
0.7
0.8
Normalized Frequency ( rad/sample)

0.9

0.1

0.2

0.3
0.4
0.5
0.6
0.7
0.8
Normalized Frequency ( rad/sample)

0.9

Phase (degrees)

0
-50
-100
-150
-200
0

Q# 3 (a)
clc; close all; clear all;
num = input('2 5 9 5 3');
den = input('5 -45 2 1 1');
%Conversion from rational to Factored form
[z,p,k]=tf2zp(num,den);
%Convert transfer function parameters to zero-pole-gain form
disp('Zeros are at');
disp(z);
disp('Poles are at');
disp(p);
disp('Gain Constant');
disp(k); %Determination of radius of the poles
radius=abs(p);
disp('Radius of the poles');
disp(radius); %Pole Zero Map using numerator and denominator coefficients
zplane(num,den) %Conversion from factored to second ordered factored
sos=zp2sos(z,p,k);
disp('Second Order Sections');
disp(sos);
figure;

RESPONSE:
2 5 9 5 3
5 -45 2 1 1
Zeros are at
Poles are at
Gain Constant
Radius of the poles
Second Order Sections
1
0
0

Name: SHAHBAZ IQBAL

Roll No: BEES/S11/0150

Page 2

1
0.8
0.6

Imaginary Part

0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
-1

-0.5

0
Real Part

0.5

Q# 3.(b)
clc; close all; clear all;
num = input('6 -9');
den = input('1 -2.5 1');
%Conversion from rational to Factored form
[z,p,k]=tf2zp(num,den);
%Convert transfer function parameters to zero-pole-gain form
disp('Zeros are at');
disp(z);
disp('Poles are at');
disp(p);
disp('Gain Constant');
disp(k); %Determination of radius of the poles
radius=abs(p);
disp('Radius of the poles');
disp(radius); %Pole Zero Map using numerator and denominator coefficients
zplane(num,den) %Conversion from factored to second ordered factored
sos=zp2sos(z,p,k);
disp('Second Order Sections');
disp(sos);
figure;

Name: SHAHBAZ IQBAL

Roll No: BEES/S11/0150

Page 3

1
0.8
0.6

Imaginary Part

0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
-1

-0.5

0
Real Part

0.5

RESPONSE:
6 -9
1 -2.5 1
Zeros are at
Poles are at
Gain Constant
Radius of the poles
Second Order Sections
1

Name: SHAHBAZ IQBAL

Roll No: BEES/S11/0150

Page 4

You might also like