<< Chapter < Page Chapter >> Page >
Introduction to MATLAB and script files. This development of these labs was supported by the National Science Foundation under Grant No. DUE-0511635. Any opinions, conclusions orrecommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation.

Introduction

The goal of this lab is to provide exercises that will help you get started learning MATLAB. You will learn about the help function, vectors, complex numbers, simple math operations, and 2-D plots. You may find it useful to try some of the built-in demos in Matlab. Type demo to see the choices. In particular, look at the demo on "Basic matrix operations" (under "Mathematics") and on "2-D" plots (under "Graphics"). We will also look at script files in MATLAB, which we will refer to as M-files and have the file extension *.m .

Getting started

Start MATLAB by clicking on it in the start menu. Once MATLAB is running you will see a screen similar to Figure 1. The command window, (A), is where you can enter commands. The current working directory, (B), displays the directory that MATLAB will look first for any commands. For example, if you made a new MATLAB function called myfunc.m , then this will need to be placed in the current working directory. You can change the working directory by typing it in the box, clicking the "..." button to browse to a folder, or typing cd [directory name] (i.e. cd 'H:\ee235\lab0\' ), which is similar to the DOS/Linux cd command).

MATLAB supports tab completion. This means that you can type part of a command and then press tab and it will fill in the rest of the command automatically.
The workspace displays information about all the variables currently active and is shown in (C). The files in the current directory can also be displayed in (C) by clicking on the tab labeled Current Directory . A history of your commands is shown in (D). If you find that you do not need some of these windows open you can close them by clicking on the small x in that section of the window.

The matlab gui

(a) Command window, (b) working directory, (c) workspace, (d) command history.

There are a number of different ways to use MATLAB on Linux. Typing matlab at the command prompt will run MATLAB in X-Windows (warning, MATLAB in X-Windows can be slow when connecting off campus). To run MATLAB without X-Windows type matlab -nodisplay . You can also run MATLAB using the current terminal for commands and use X-Windows for everything else (like figures) by typing matlab -nodesktop .

Matlab commands

MATLAB works with matrices and therefore treats all variables as a matrix. To enter the matrix x = ( 3 1 5 6 4 1 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIhacqGH9aqpdaqadaqaauaabeqacmaaaeaacaaIZaaabaGaaGymaaqaaiaaiwdaaeaacaaI2aaabaGaaGinaaqaaiaaigdaaaaacaGLOaGaayzkaaaaaa@3DF1@ type the command x = [3 1 5; 6 4 1] . We can represent an array with a vector, which is a special case of a matrix. A vector can be entered in by typing y = [1 2 3] . Now try entering z = [1 2 3]' . Is the output what you expect?

  • Familiarize yourself with the help command. Typing help gives you a list of all help topics. Typing help<topicname> gives help on a specific MATLAB function. For example, use help plot to learn about the plot command.

    More useful commands

    • whos lists all variables
    • clear clears all variables
  • Perform the following operations in MATLAB:
    • Generate the following column vectors as MATLAB variables: x = ( 2 4 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIhacqGH9aqpdaqadaqaauaabeqaceaaaeaacaaIYaaabaGaaGinaaaaaiaawIcacaGLPaaaaaa@3AF5@ and y = ( 3 1 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadMhacqGH9aqpdaqadaqaauaabeqaceaaaeaacaaIZaaabaGaaGymaaaaaiaawIcacaGLPaaaaaa@3AF4@
    • Using the computer, issue the following MATLAB commands x * y' x' * yx .* y Be sure you understand the differences between each of these and you know what the ' , * , and .* operators do.
    • Convince yourself that the answer makes sense by checking the matrix dimension and computing each result by hand.
  • Plot and subplot

    The Matlab command plot allows you to graphically display vector data (in our case here, the two signals). For example, if you had the variables t for time, and y for the signal, typing the command plot(t, y); will display a plot of t vs. y . See help plot for more information.

    Annotating your plots is very IMPORTANT! Here are a few annotation commands.

    • title('Here is a title'); - Adds the text "Here is a title" to the top of the plot.
    • xlabel('Control Voltage (mV)'); - Adds text to the X-axis.
    • ylabel('Current (mA)'); - Adds text to the Y-axis.
    • grid on; - Adds a grid to the plot.

    In order to display multiple plots in one window you must use the subplot command. This command takes three arguments, subplot(m, n, p) . The first two breaks the window into a m by n matrix of smaller plot windows. The third argument, p, selects which plot is active.

    For example, if we have three signals x, y, and z, that we want to plot against time, t, then we can use the subplot command to produce a single window with three plots stacked on top of each other. subplot(3,1,1); plot(t,x);subplot(3,1,2); plot(t,y);subplot(3,1,3); plot(t,z);

    See help subplot and help plot for much more information.

    Create and plot a signal x 0 ( t ) = t e | t | MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIhadaWgaaWcbaGaaGimaaqabaGccaGGOaGaamiDaiaacMcacqGH9aqpcaWG0bGaamyzamaaCaaaleqabaGaeyOeI0IaaGymamaaemaabaGaamiDaaGaay5bSlaawIa7aaaaaaa@42FA@ over the time range [-10,10] using the following MATLAB commands: t = -10:0.1:10; xo = t .* exp(-abs(t));plot(t, xo); grid; The first command defines an array with time values having an 0.1 increment. The ";" is used to suppress printout of the arrays (which are large), and the "grid" command makes the plot easier to read. Now create the signals:
    x e t t e t
    x ( t ) = 0.5 * [ x 0 ( t ) + x e ( t ) ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIhacaGGOaGaamiDaiaacMcacqGH9aqpcaaIWaGaaiOlaiaaiwdacaGGQaWaamWaaeaacaWG4bWaaSbaaSqaaiaaicdaaeqaaOGaaiikaiaadshacaGGPaGaey4kaSIaamiEamaaBaaaleaacaWGLbaabeaakiaacIcacaWG0bGaaiykaaGaay5waiaaw2faaaaa@4892@
    Plot all signals together using 3 plots stacked on top of each other with the subplotcommand. subplot(3,1,1); plot(t,xo);subplot(3,1,2); plot(t,xe);subplot(3,1,3); plot(t,x); Note that x o ( t ) and x e ( t ) are the odd and even components, respectively, of x t t e t u t
  • Complex Numbers: One of the strengths of MATLAB is that most of its commands work with complex numbers. Perform the following computations in MATLAB.
    • MATLAB recognizes i as an imaginary number. Try entering sqrt(-1) into MATLAB, does the result make sense?
    • MATLAB uses the letter i instead of j by default. Electrical Engineers prefer using j however, and MATLAB will recognize that as well. Try entering i+j , does this make sense.
      If you are using complex numbers in your code, it's a good idea to avoid using i and j as variables to prevent confusion.
    • Define z 1 = 1 + j . Find the magnitude, phase, real and imaginary parts of z (using abs() , angle() , real() , imag() , respectively). Is the phase in radians or degrees?
    • Find the magnitude of z 1 + z 2 where z 2 2 e 3
    • Compute the value of j j MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadQgadaahaaWcbeqaaiaadQgaaaaaaa@37ED@ . Is the result what you expect?
  • Complex Functions: MATLAB also handles complex time functions in the same way (again, implemented as vectors) . Create a signal x 1 t t e jt over the range [-10,10], as in part 3. Next plot the real and imaginary parts of the signal in two plots, one over the other using the subplot command. Notice that oneplot is odd and one is even. Try proving to your self analytically that this is what you would expect.
  • Playing and Plotting a Sound Load the built-in data named "handel" and play it: load handel; plot(linspace(0,9,73113),y);sound(y);
You can use the clear command in MATLAB to clear all of the varibles

Script files

Scripts are m-files files that contain a sequence of commands that are executed exactly as if they were manually typed into the MATLAB console. Script files are useful for writing things in MATLAB that you want to save and run later. You can enter all the commands into a script file that you can run at a later time, and have the ability to go back and edit it later.

You need to use a text editor to create script files, e.g. Notepad on the PC's ( pico , emacs , or vi on Linux machines). MATLAB also has an internal editor that you can start by clicking on a .m file within MATLAB's file browser. All are easy to learn and will produce text files that MATLAB can read.

Click here to download the dampedCosine.m script and be sure to save it with that name to follow the instructions here exactly. It is very important that script filenames end in .m . Be sure that MATLAB's working directory is set to the location of where you saved the script file. Type dampedCosine at the MATLAB prompt. Look at the m-file in a text editor and verify that you get the plot predicted in the comment field of the script.

The % character marks the rest of the line as a comment.

Scripts

Now we are going to edit the dampedCosine.m file to create our own script file. At the top of the file, you will see the following commands

diary 'your_name_Lab1.txt' disp('NAME: your name')disp('SECTION:your section')

  • Edit the dampedCosine.m (download from link above) script and enter your name and section where indicated. Save this new version of the script as yourName_dampedCosine.m
  • Edit the script to create a second signal where the cosine with twice the period (which gives half the frequency) of the first.
  • Add to the script the commands to plot these together with the first signal on top and the second on the bottom. In other words, you should have a single figure with two different plots, one on top and one on bottom. You will need to use subplot and plot . Save this plot as yourName_dampedCosine.fig .
  • Show the TA your dampedCosine plot. What is the period of the cosine?

Complex exponentials

Download and run compexp.m , which includes a 3-D plot of a complex exponential, y t , as well as 2-D magnitude/phase and real/imaginary plots. You need 2 2-D plots to have the same information as the 3-D plot. How would you change the script to make the oscillation frequency lower by half? How would you change the script to make the decay faster? Show the TA your plots.

Questions & Answers

A golfer on a fairway is 70 m away from the green, which sits below the level of the fairway by 20 m. If the golfer hits the ball at an angle of 40° with an initial speed of 20 m/s, how close to the green does she come?
Aislinn Reply
cm
tijani
what is titration
John Reply
what is physics
Siyaka Reply
A mouse of mass 200 g falls 100 m down a vertical mine shaft and lands at the bottom with a speed of 8.0 m/s. During its fall, how much work is done on the mouse by air resistance
Jude Reply
Can you compute that for me. Ty
Jude
what is the dimension formula of energy?
David Reply
what is viscosity?
David
what is inorganic
emma Reply
what is chemistry
Youesf Reply
what is inorganic
emma
Chemistry is a branch of science that deals with the study of matter,it composition,it structure and the changes it undergoes
Adjei
please, I'm a physics student and I need help in physics
Adjanou
chemistry could also be understood like the sexual attraction/repulsion of the male and female elements. the reaction varies depending on the energy differences of each given gender. + masculine -female.
Pedro
A ball is thrown straight up.it passes a 2.0m high window 7.50 m off the ground on it path up and takes 1.30 s to go past the window.what was the ball initial velocity
Krampah Reply
2. A sled plus passenger with total mass 50 kg is pulled 20 m across the snow (0.20) at constant velocity by a force directed 25° above the horizontal. Calculate (a) the work of the applied force, (b) the work of friction, and (c) the total work.
Sahid Reply
you have been hired as an espert witness in a court case involving an automobile accident. the accident involved car A of mass 1500kg which crashed into stationary car B of mass 1100kg. the driver of car A applied his brakes 15 m before he skidded and crashed into car B. after the collision, car A s
Samuel Reply
can someone explain to me, an ignorant high school student, why the trend of the graph doesn't follow the fact that the higher frequency a sound wave is, the more power it is, hence, making me think the phons output would follow this general trend?
Joseph Reply
Nevermind i just realied that the graph is the phons output for a person with normal hearing and not just the phons output of the sound waves power, I should read the entire thing next time
Joseph
Follow up question, does anyone know where I can find a graph that accuretly depicts the actual relative "power" output of sound over its frequency instead of just humans hearing
Joseph
"Generation of electrical energy from sound energy | IEEE Conference Publication | IEEE Xplore" ***ieeexplore.ieee.org/document/7150687?reload=true
Ryan
what's motion
Maurice Reply
what are the types of wave
Maurice
answer
Magreth
progressive wave
Magreth
hello friend how are you
Muhammad Reply
fine, how about you?
Mohammed
hi
Mujahid
A string is 3.00 m long with a mass of 5.00 g. The string is held taut with a tension of 500.00 N applied to the string. A pulse is sent down the string. How long does it take the pulse to travel the 3.00 m of the string?
yasuo Reply
Who can show me the full solution in this problem?
Reofrir Reply
Got questions? Join the online conversation and get instant answers!
Jobilize.com Reply

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Continuous time linear systems laboratory (ee 235). OpenStax CNX. Sep 28, 2007 Download for free at http://cnx.org/content/col10374/1.8
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Continuous time linear systems laboratory (ee 235)' conversation and receive update notifications?

Ask