Calculating Transfer Functions in Simulink: a Step-by-step Approach

by

Transfer functions represent one of the most fundamental concepts in control systems engineering and signal processing. They provide a mathematical representation of the relationship between the input and output of a linear time-invariant (LTI) system in the frequency domain. The Transfer Fcn block models a linear system by a transfer function of the Laplace-domain variable s, making it an essential tool for analyzing system dynamics, designing controllers, and predicting system behavior. Simulink, MATLAB’s graphical programming environment, offers powerful capabilities for working with transfer functions through both direct implementation and analytical extraction methods.

This comprehensive guide explores the complete process of calculating, implementing, and analyzing transfer functions in Simulink. Whether you’re designing a new control system, analyzing an existing model, or extracting transfer functions from complex nonlinear systems, understanding these techniques will significantly enhance your engineering workflow and system analysis capabilities.

Understanding Transfer Functions in Control Systems

Before diving into the practical implementation in Simulink, it’s essential to understand what transfer functions represent and why they’re so valuable in control systems engineering. A transfer function describes the input-output relationship of a system in the Laplace domain, typically expressed as a ratio of polynomials in the complex variable s.

Mathematical Foundation

The transfer function relates system input and outputs through numerator and denominator coefficients in descending powers of s, where the order of the denominator must be greater than or equal to the order of the numerator. This mathematical constraint ensures physical realizability and system causality.

The general form of a transfer function can be written as G(s) = N(s)/D(s), where N(s) represents the numerator polynomial and D(s) represents the denominator polynomial. The coefficients of these polynomials directly correspond to the system’s physical parameters and determine critical characteristics such as poles, zeros, gain, and frequency response.

Types of Transfer Function Representations

Simulink supports multiple approaches for representing transfer functions. The Transfer Fcn block can model single-input single-output (SISO) and single-input multiple-output (SIMO) systems. For more complex scenarios involving time-varying parameters, Simulink provides specialized blocks that accommodate dynamic coefficient changes during simulation.

You can bring in transfer function objects defined in the MATLAB workspace into Simulink by using the LTI System block and specifying the variable name, and a transfer function can also be represented in terms of simple blocks, such as integrators and gains. This flexibility allows engineers to choose the most appropriate representation method based on their specific application requirements and modeling preferences.

Proper environment setup is crucial for efficient transfer function analysis in Simulink. The process begins with launching Simulink and configuring the necessary toolboxes and libraries that provide transfer function functionality.

Start by opening MATLAB and typing “simulink” in the command window or clicking the Simulink icon in the MATLAB toolbar. This opens the Simulink Start Page, where you can create a new blank model. The blank model provides a canvas where you’ll construct your system using various blocks from the Simulink Library Browser.

The Library Browser contains organized categories of blocks, with transfer function-related blocks primarily located in the Continuous library. The Continuous library contains continuous-time system elements including transfer functions, state-space models, and PID controllers. Familiarizing yourself with this library structure will streamline your modeling process.

Essential Blocks for Transfer Function Analysis

A typical transfer function analysis setup requires several key components. The simple model consists of three blocks: Step, Transfer Function, and Scope, where the Step is a Source block from which a step input signal originates, this signal is transferred through the line to the Transfer Function Continuous block, the Transfer Function block modifies its input signal and outputs a new signal to the Scope.

Beyond these basic blocks, you may need additional components depending on your analysis objectives:

  • Signal sources: Step, Ramp, Sine Wave, Chirp Signal, and Random Number blocks for various input scenarios
  • Mathematical operations: Sum, Gain, Product blocks for signal manipulation and system interconnection
  • Visualization tools: Scope, To Workspace, Display blocks for observing and recording simulation results
  • Analysis blocks: Linearization points, measurement blocks for extracting system characteristics

Implementing Transfer Functions Using the Transfer Fcn Block

The Transfer Fcn block serves as the primary method for directly implementing known transfer functions in Simulink. This approach is ideal when you have already derived the mathematical model of your system and want to simulate its behavior.

Configuring Block Parameters

In Simulink, search for the ‘Transfer Fcn’ block and drag and drop the block into your model workspace. Once placed, you need to configure its parameters to match your desired transfer function.

A block can be modified by double-clicking on it, and if you double-click on the Transfer Function block, you will see a dialog box that contains fields for the numerator and the denominator of the block’s transfer function. These fields accept vector inputs representing polynomial coefficients.

Entering Coefficient Vectors

By entering a vector containing the coefficients of the desired numerator or denominator polynomial, the desired transfer function can be entered. The coefficients must be arranged in descending powers of s, starting with the highest order term.

For example, consider implementing the transfer function G(s) = (s + 3)/(s² + 2s). Enter the numerator coefficients [1 3] (representing s + 3) and enter the denominator coefficients [1 2 0] (representing s² + 2s). Note that you must include zeros for any missing terms in the polynomial to maintain proper coefficient alignment.

Block Display Options

If each is specified as an expression, a vector, or a variable enclosed in parentheses, the icon shows the transfer function with the specified coefficients and powers of s, and if you specify a variable in parentheses, the variable is evaluated. This feature allows for flexible display options that can make your model more readable and maintainable.

You can specify coefficients in three ways: as explicit numerical vectors (e.g., [1 3]), as MATLAB workspace variables (e.g., num and den), or as expressions that evaluate to vectors. Each method has advantages depending on whether you prioritize model clarity, parameter flexibility, or computational efficiency.

Handling Initial Conditions

Simulink presets the initial conditions of the Transfer Fcn block to zero, and to specify initial conditions for a given transfer function, convert the transfer function to its controllable, canonical state-space realization using tf2ss. This limitation exists because transfer functions have infinitely many time-domain realizations, and most state representations don’t correspond to physically meaningful initial conditions.

When non-zero initial conditions are required, the recommended approach involves converting to state-space representation, setting the initial state vector, and using the State-Space block instead. This provides explicit control over initial conditions while maintaining mathematical equivalence to the transfer function representation.

One of Simulink’s most powerful capabilities is extracting transfer functions from existing models, including complex nonlinear systems. This process, called linearization, approximates system behavior around an operating point using a linear transfer function model.

Understanding Model Linearization

A more direct way to find a transfer function from a Simulink model is to linearize it to a state-space model, then convert it to a transfer function using sys = linearize(“mymodel”,io_points); sys_tf = tf(sys). This two-step process first extracts a linear state-space representation and then converts it to transfer function form.

Linearization works by computing the Jacobian matrices of the system equations at a specified operating point. For nonlinear systems, this produces a linear approximation that’s valid in a neighborhood around that point. For systems that are already linear, linearization extracts the exact transfer function regardless of the operating point chosen.

Specifying Analysis Points

Before linearizing a model, you must identify which signals represent inputs and outputs for the transfer function you want to extract. To specify the portion of the model to linearize, first open the Linearization tab by clicking Linearization Manager in the Apps gallery in the Simulink window.

To specify an analysis point for a signal, click the signal in the model, then on the Linearization tab, in the Insert Analysis Points gallery, select the type of analysis point. Simulink supports several analysis point types:

  • Input Perturbation: Adds a small signal at the specified location to measure system response
  • Output Measurement: Records the signal value without affecting system dynamics
  • Loop Break: Opens a feedback loop for open-loop analysis
  • Open-Loop Output: Combines output measurement with loop opening

To specify the portion of the model to linearize, create an array of linearization I/O objects using the linio command, such as creating an input perturbation analysis point at the output of the PID Controller block. This programmatic approach offers more flexibility than interactive point placement, especially for batch processing or automated workflows.

Using the Model Linearizer App

Model Linearizer lets you perform linear analysis of nonlinear Simulink models. This powerful application provides a comprehensive interface for linearization tasks, operating point specification, and result analysis.

You can open the Model Linearizer App by going to the Simulink Toolstrip: On the Apps tab, under Control Systems, click Model Linearizer. Once opened, the app displays your model’s analysis points and allows you to configure linearization settings.

The Model Linearizer interface includes several key sections:

  • Analysis I/Os: Displays and manages input/output point configurations
  • Operating Point: Specifies the state and input values around which linearization occurs
  • Linearization Results: Shows extracted models and their characteristics
  • Plot Options: Generates frequency and time domain response plots

Performing the Linearization

For this example, use the model operating point for linearization by leaving Model Initial Condition selected in the Operating Point drop-down list, and to linearize the system and generate a response plot for analysis, in the Linearize section, click a response.

The linearization process compiles the model, evaluates it at the specified operating point, and computes the linear approximation. For large or complex models, this may take several seconds. The resulting linear model appears in the Linear Analysis Workspace, where you can examine its properties, generate plots, and export it to MATLAB for further analysis.

You can also export the linearized model to the MATLAB workspace by right-clicking linsys1 in the data browser and selecting Export to MATLAB Workspace. This enables additional analysis using Control System Toolbox functions like bode, step, margin, and others.

Working with Time-Varying Transfer Functions

Many real-world systems have parameters that change over time or operating conditions. Simulink provides specialized blocks for handling transfer functions with variable coefficients, enabling accurate simulation of gain-scheduled controllers, adaptive systems, and parameter-varying plants.

The Varying Transfer Function Block

Use this block and the other blocks in the Linear Parameter Varying library to implement common control elements with variable parameters or coefficients, for more information, see Model Gain-Scheduled Control Systems in Simulink.

Transfer function order N is specified as a positive integer, and this parameter determines the number of coefficient input ports on the block. Each coefficient can be supplied as a time-varying signal, allowing the transfer function characteristics to adapt during simulation.

Configuring Variable Coefficients

Enable the b0 input port for a transfer function with direct feedthrough, and for a zero-feedthrough transfer function, clear this checkbox, as disabling the port for zero-feedthrough models is numerically more reliable than feeding a zero-constant into the port.

When implementing time-varying transfer functions, careful attention to numerical stability is essential. Rapid coefficient changes can cause numerical issues, especially if the denominator polynomial approaches zero or if poles move significantly between time steps. Consider implementing rate limiters or smoothing filters on coefficient signals to maintain simulation stability.

Avoiding Algebraic Loops

Avoid making the transfer-function coefficients depend on the block output y, because if you have such dependence, the resulting transfer function causes an algebraic loop, since computing the block output value requires knowing the block output value, and this algebraic loop is prone to instability and divergence, so instead of the output, try expressing the coefficients in terms of the time t and the block input u.

Algebraic loops occur when Simulink cannot determine the order of block execution because outputs depend on inputs in a circular fashion. While Simulink can sometimes solve algebraic loops iteratively, they significantly slow simulation and may fail to converge, especially with nonlinear dependencies or discontinuous coefficient changes.

Advanced Linearization Techniques

Beyond basic linearization, Simulink offers advanced techniques for extracting transfer functions from complex systems, handling multiple operating points, and analyzing frequency-dependent behavior.

Batch Linearization

If you want to obtain multiple open-loop or closed-loop transfer functions from the linearized system without recompiling the model, you can specify linear analysis points using an slLinearizer interface, for more information, see Mark Signals of Interest for Batch Linearization.

The slLinearizer interface provides a programmatic approach to linearization that’s particularly valuable for parametric studies, Monte Carlo analysis, or automated design optimization. You can define multiple input-output pairs, sweep through parameter ranges, and collect linearization results without manual intervention.

Operating Point Specification

The operating conditions are the values of the state and input vectors about which the linearization of the analysis model occurs, and initial conditions are the values of the state and input vectors at time t=0. For nonlinear systems, the choice of operating point significantly affects the resulting transfer function.

Simulink provides several methods for specifying operating points:

  • Model Initial Conditions: Uses the initial state and input values defined in the model
  • Steady-State Operating Point: Computes equilibrium conditions where derivatives equal zero
  • Simulation Snapshot: Captures system state at a specific simulation time
  • Custom Operating Point: Manually specifies state and input values

For control system design, steady-state operating points often provide the most meaningful linearization. These represent the conditions around which the controller will regulate the system, making the linear model directly applicable to controller synthesis and stability analysis.

Frequency Response Estimation

To estimate the transfer function of a system in Simulink, use the Discrete Transfer Function Estimator block, which implements the Welch’s average modified periodogram method and uses the measured input and output data for estimation.

This approach differs from linearization by using actual simulation data rather than analytical derivatives. It’s particularly useful when:

  • The system contains components that don’t support analytical linearization
  • You want to validate a linearization against simulation results
  • The system includes stochastic elements or measurement noise
  • You need frequency response data at specific frequencies

The relationship between the input x and output y is modeled by the linear, time-invariant transfer function Txy, and the transfer function is the ratio of the cross power spectral density of x and y, Pyx, to the power spectral density of x, Pxx. This frequency-domain approach provides robust estimates even in the presence of noise and disturbances.

Simulation and Analysis Workflow

Once you’ve implemented or extracted a transfer function, the next step involves simulation and analysis to understand system behavior and validate performance against requirements.

Configuring Simulation Parameters

In the model window, select Model Configuration Parameters from the Simulation menu, and there are many simulation parameter options including start and stop times, which tell Simulink over what time period to perform the simulation.

Key simulation parameters to consider include:

  • Solver Type: Variable-step solvers automatically adjust time steps for accuracy, while fixed-step solvers maintain constant intervals suitable for real-time applications
  • Relative Tolerance: Controls the acceptable error in state calculations, with smaller values increasing accuracy but slowing simulation
  • Absolute Tolerance: Sets the threshold below which state values are considered negligible
  • Maximum Step Size: Limits how large time steps can grow, preventing the solver from missing fast dynamics

Variable-step solvers use absolute and relative tolerances when choosing the step size to determine whether the error in state calculations is acceptable, and to inherit the absolute tolerance from the Absolute tolerance configuration parameter, specify this parameter value as auto or -1.

Visualizing Results with Scopes

The most complicated of these three blocks is the Scope block, and double-clicking on this brings up a blank oscilloscope screen, and when a simulation is performed, the signal which feeds into the scope will be displayed in this window.

Modern Simulink versions offer enhanced scope functionality including:

  • Multiple input channels with independent scaling
  • Cursor measurements for precise value reading
  • Signal statistics (mean, RMS, peak values)
  • Data logging to workspace or file
  • Trigger conditions for capturing specific events
  • FFT analysis for frequency content examination

Exporting Data to MATLAB

For detailed analysis beyond what Simulink’s built-in tools provide, export simulation data to the MATLAB workspace. Use To Workspace blocks to capture signals of interest, or configure scopes to log data automatically. Once in MATLAB, you can apply the full range of Control System Toolbox functions for analysis.

Common post-processing tasks include:

  • Computing performance metrics (rise time, settling time, overshoot, steady-state error)
  • Generating publication-quality plots with custom formatting
  • Comparing multiple simulation runs or parameter variations
  • Performing statistical analysis on Monte Carlo simulation results
  • Fitting transfer function models to simulation data

Practical Applications and Examples

Understanding transfer function calculation in Simulink becomes more concrete through practical examples that demonstrate common applications in control systems engineering.

Controller Design and Analysis

The controllers transfer function is implemented using the transfer function block, which is what we use to represent the engine and the actuator as well, and we can see that the model is able to follow step inputs with some overshoot and zero steady-state error.

A typical controller design workflow involves:

  1. Modeling the plant (system to be controlled) as a transfer function
  2. Specifying performance requirements (bandwidth, phase margin, disturbance rejection)
  3. Designing a controller transfer function to meet these requirements
  4. Implementing the closed-loop system in Simulink
  5. Simulating with various inputs and disturbances
  6. Iterating the design based on simulation results

The Simulink Control Design toolbox offers the functionality to extract a model from Simulink into the MATLAB workspace, which is especially useful for complicated, or nonlinear simulation models, and is also useful for generating discrete-time (sampled) models.

System Identification

When you have experimental data or simulation results but don’t know the underlying transfer function, system identification techniques can estimate a transfer function model that matches the observed behavior. This process involves:

  • Collecting input-output data from the system
  • Choosing an appropriate model structure (order, delays)
  • Estimating parameters that minimize prediction error
  • Validating the model against independent data sets

Simulink’s Discrete Transfer Function Estimator block automates much of this process for frequency-domain identification, while the System Identification Toolbox provides comprehensive capabilities for time-domain and frequency-domain parameter estimation.

Multi-Domain System Modeling

Transfer functions excel at representing systems from various engineering domains—mechanical, electrical, thermal, hydraulic—using a unified mathematical framework. Simulink enables you to combine transfer functions from different domains into integrated system models.

For example, an electromechanical actuator system might include:

  • Electrical transfer function relating voltage to current
  • Electromechanical transfer function relating current to torque
  • Mechanical transfer function relating torque to position
  • Sensor transfer function relating position to measurement signal

By connecting these transfer function blocks in series and feedback configurations, you create a complete system model that captures the essential dynamics while remaining computationally efficient.

Troubleshooting Common Issues

Working with transfer functions in Simulink occasionally presents challenges. Understanding common issues and their solutions w