Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
2 views

Expt-3-Transfer-Function

This document outlines an experiment focused on determining the transfer function of a mechanical rotational system using MatLab. It details the steps to create symbolic expressions, compute transfer functions, and analyze system responses, including finding poles, zeros, and DC gain. The document also provides instructions for plotting the system's output over time and deriving the output signal equation.

Uploaded by

s2302810
Copyright
© © All Rights Reserved
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
2 views

Expt-3-Transfer-Function

This document outlines an experiment focused on determining the transfer function of a mechanical rotational system using MatLab. It details the steps to create symbolic expressions, compute transfer functions, and analyze system responses, including finding poles, zeros, and DC gain. The document also provides instructions for plotting the system's output over time and deriving the output signal equation.

Uploaded by

s2302810
Copyright
© © All Rights Reserved
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 4

EXPERIMENT NO.

3
TRANSFER FUNCTION

I. INTRODUCTION

1. In MatLab polynomials can be manipulated using symbolic expressions. The constructor is


“syms”
Example: To declare a polynomial Y = s2 + 4s + 3 (variable is “s”)
>> syms s
>> Y = s^2 + 4*s + 3

The first derivative of Y is:


>> diff(Y)

2. System transfer function can be created using the tf() function. First, create the numerator
and denominator polynomials expressing them as row vectors.

Example: s+1
H = -------------------------- can be created as follows:
s3 + 2s2 + 3s + 2

>>numerator = [1 1] → s+1
>>denominator = [1 2 3 2] → s3 + 2s3 + 3s + 2

3. Zeros are values of ‘s’ that would make the system transfer function zero. Or the roots of the
numerator of the transfer function that are not common with the poles.

4. Poles are values of ‘s’ that would make the system transfer function infinite. Or the roots of
the denominator of the transfer function that are not common with the zeros.
5. System response of the system can be view using the plot() function.

II. OBJECTIVES

1. Determine the transfer function of the rational system


2. Use MatLab to derive the output signal equation in differential equation
3. Ploy and analyze system response
4. Determine the pole, zero locations and dc gain of the mechanical system

III. MATERIALS/EQUIPMENT

Matlab Software
Personal Computer
Objective A and B. Determine the transfer function by computation given the mechanical rotational
system:

Transfer function: θ(s)/τ(s) = _____________________________

Concept
Objective C.

1. Run MatLab Software.


2. Using symbolic expression constructor “syms” construct “s” as a symbolic expression >>syms s
3. Write the transfer function, save the result to variable TF

>>TF = 1 / (s^2 + 2*s + 2)

output: TF = ____________________________________

4. Use function pretty() to display a more understandable equation >>pretty(TF)

output: _________________________________________

5. If the input torque is 10 u(t) what is τ(s) ? __________________


6. Store the value of input torque to MatLab variable TORQUE

7. The output of the system is TF x Input. Store the output signal equation to
OUTPUT
>>OUTPUT = TF * TORQUE

output: OUTPUT = _______________________

8. The inverse laplace of an equation can be solve using ilaplace() function. Store the inverse laplace of
OUTPUT OUTPUT_TIME variable

>>OUTPUT_TIME = ilaplace(OUTPUT)

output: OUTPUT_TIME = ___________________________

9. To plot the equation in time is to substitute values of “t”. Create a series of time values from 0 to 15
seconds incremented by 0.1 using row vector declaration.

>>t = [0 : 0.1 : 15]

10. Substitute the values of “t” to OUTPUT_TIME equation using subs() function. This will create 150
values or points. Store the answer to variable RESULT

>>RESULT = subs(OUTPUT_TIME, t)

Tip: So as not to display the 150 values add a semicolon (;) at the end of the statement.

11. Plot the values of OUTPUT_TIME using plot(x,y) function.

>>plot(t, RESULT)

Draw the response curve:

What is the maximum amplitude? _______________

What is the approximate transient time? ________________

What is the steady state value/amplitude? _____________________


Objective D. Determine poles, zeros and dc gain

1. Declare the numerator and denominator as a row vector:

>>numerator = ________________
>>denominator = ________________

2. Create the transfer function using tf()


>>SYS = tf(numerator, denominator)

3. Solve for the pole location of the transfer function


>>pole(SYS)

ans:

4. Solve for the zero location of the transfer function


>>zero(SYS)

ans:

5. Using tf2zp() function pole, zero and dc gain can also be found.
>> [z,p,g] = tf2zp(num,den)

ans: z=

p=

g=

V. CONCLUSION

You might also like