Introduction

This document describes how to perform common tasks in 6.302 using Octave and Scilab. These two software packages are free alternatives to MATLAB, and their use is encouraged.

Each software package employs different methods to create, examine, print, and save transfer functions. For the purposes of this tutorial, we will use the following transfer function as an example:

L(s) = (3 * (0.05s + 1)^2) / ((s+1)^3 * (0.01s + 1))

Octave

GNU Octave is a numerical computation software package very similar to MATLAB. In fact, the syntax is almost exactly the same. For help, type help function_name. For an excellent tutorial on control-related functions, type DEMOcontrol.

Creating transfer functions

For the purposes of 6.302, there exist two ways to create transfer functions in Octave. The first is tf2sys, which takes as arguments the coefficients of the numerator and denominator of the transfer function. With the example transfer function, you would type:

L = tf2sys(3 * [0.0025 0.1 1], [0.01 1.03 3.03 3.01 2]);

Another command is zp2sys, which takes as arguments the zeros, poles, and leading coefficient of the transfer function. Consequently, you would type:

L = zp2sys([-20 -20], [-1 -1 -1 -100], 3);

Displaying transfer functions

You will notice that if you merely type L, the output isn't very readable. Use sysout(L, "tf") to display the function in polynomial form and sysout(L, "zp") for zero-pole form.

The graph commands

Root locus

For root locus, use rlocus(L).

L = tf2sys(3 * [0.0025 0.1 1], [0.01 1.03 3.03 3.01 2]);
title("Root locus of L");
rlocus(L);

Bode

For Bode, use bode(L).

L = tf2sys(2*[1, 1], [1, 0, 1]);
bode(L);

Nyquist

For Nyquist, use nyquist(L).

L = tf2sys(2*[1, 1], [1, 0, 1]);
nyquist(L, logspace(-1, 1, 100));
title("Nyquist plot of L");
xlabel("Real axis");
ylabel("Imag axis");
replot;

Nichols

For Nichols, use nichols(L).

L = tf2sys(2*[1, 1], [1, 0, 1]);
nichols(L);
title("Nichols plot of L");
replot;

Auxilary graph commands

As in MATLAB, the commands title, xlabel, and ylabel assign strings to the graphs. Often after these commands, replot is necessary to activate the changes.

However, bode does not work after a replot command; if you find yourself in a situation where a replot is necessary, plot the data from bode manually using plot. (See online help for details.)

Saving and printing graphs

Octave uses Gnuplot to plot its graphs, so the saving and printing process is not as straightforward as clicking a button.

To save a graph in PNG format (which can be read by almost any graphics program), use the following commands:

gset terminal png;
gset output "filename.png";
replot;

Again, with Bode plots, just use bode(L) again instead of replot.

To print a graph, use the following commands:

gset terminal postscript;
gset output "|lpr";
replot;

After saving or printing, reset the environment with the command gset terminal x11, or else any subsequent graphs will be saved/printed without being displayed onscreen.

SciLab

Scilab is a free software package similar to Maple. Its syntax is considerably different from Matlab, but Scilab can perform symbolic manipulation, unlike Matlab. It includes excellent demos, as found by pressing the "Demos" button on the toolbar.

Creating transfer functions

This is perhaps the easiest task of all. Just type the following:


s = poly(0, "s");

L = syslin('c', 2*(s+1)/(s^2+1));

The poly command simply defines a polynomial with the symbol "s." Don't forget it! Simply type L to see the transfer function you've just defined.

The graph commands

Root locus

For root locus, use evans(L).

s = poly(0, "s");
L = syslin('c', 2*(s+1)/(s^2+1));
evans(L, 2.6);