CS310 Team Lab # 1
Engineering Cost Analysis in MATLAB
Numeric Computation

OBJECTIVES

INTRODUCTION

Engineering cost analysis requires decision-making based upon comparison of different financing alternatives.

When you do some work for someone, you expect to get paid. The same is true if you loan money to someone. You expect them to pay you interest for the use of your money. On the other side of the coin (no pun intended), if you need to borrow money, you must pay someone else interest for the use of their money. Rather than doing this transaction individual by individual, society has banks that serve as intermediaries who pay interest for money given to them and who charge interest for money that they loan. Thus in the financial world, money is like labor or goods and services. Businesses that need money must pay for it. This is the way the business, and thus the engineering, world works.

A general principle that applies to all financial transactions is called the time value of money (or future value of money). A dollar that you have today is worth more in the future because you don't hide it under your mattress. You invest it in a savings account, certificate of deposit, stock or bond or some other "financial instrument" and your money earns interest. This works because the institution that takes your money immediately loans it to someone else for their use and charges them interest.

The time value of money is computed using the formula: F = P(1 + a)n where:

P = the present value of your money,
F = the future value of your money after n periods (months or years),
a = the percent interest paid for your money per period (month or year), and
n = the number of periods (months or years) over which the interest is computed (or compounded)

MATLAB is short for Matrix Laboratory. It is a computing environment that makes it easy to perform very complex numeric computations. In numeric computing, we compute initial, intermediate, and final numeric values (numbers) in order to solve problems. This is in contrast with symbolic computing where symbols are used to create expressions and equations and then the equations are manipulated without much thought to the specific underlying value for each symbol.

In numeric computing, the values themselves are needed to find the next step in the solution and the solution is one or more values instead of a symbol or symbolic expression. Note: Variable names are used in MATLAB, but they are not abstract symbolic variables as you might see in a math equation. Rather they are more like nicknames for values that have been computed.

Let's get started!

PROBLEMS

1. Use MATLAB's Command Window to answer the following questions or perform the described calculations

  1. How much money will you have after five years, if you have $10,000 and put it in a bank account that pays 4% interest annually?

    Type the following into the Command Window (don't type the >> -- that's the command prompt):

    >> F = 10000*(1+0.04)^5

  2. For each interest rate from 4% to 7% in 0.5% increments, how much money will you have after five years, if the same amount is invested?

    Hopefully, you are considering how best to do this many calculations without retyping all of the values needed. Hint: First, enter the interest rates into a single vector of interest rates and save it with the name a and then perform an element-wise calculation using all elements in the vector. Type the following to create the vector using the colon operator:

    >> a = 0.04 : 0.005 : 0.07

    Repeat an earlier computation by using the up arrow until the calculation is at the command prompt and then editing as needed for the new calculation. Find your calculation from part a and change the 0.04 to the name of your vector of interest rates, a.

    >> F = 10000*(1+a)^5

  3. What happened?

    Before we explain the error, let's just fix it. Type a period (.) before the ^ symbol in the expression and press the Enter key. The (.^) is an element-wise operator. It indicates that the power operator needs to occur on each value (element) in the matrix, instead of performing matrix algebra.

2. Matrix algebra in MATLAB

The power operator by itself (^) performs the exponential multiplication as defined by matrix algebra. Let's review matrix multiplication first.

  1. Compute [by hand] the values for the blank cells of the result matrix.

    [1 2;3 4;5 6;7 8]*[1 2 3;4 5 6]=[9 12 15;19 26 __;29 __ __; __ __ __]

    SHOW YOUR COMPLETED MATRIX TO A LAB LEADER

    Check your work using MATLAB:

    >> A = [ 1 , 2 ; 3 , 4 ; 5 , 6 ; 7 , 8 ]
    >> B = [ 1, 2, 3 ; 4, 5, 6 ]
    >> C = A * B
    >> D = B * A

    Now what happened?
    The error MATLAB produced is due to the fact that MATLAB tried to multiply two matrices with incompatible dimensions. A 4×2 matrix times a 2×3 matrix is defined, but multiplying a 2×3 matrix with a 4×2 matrix is undefined in matrix algebra, so the B*A operation failed.

  2. Add a scalar value to the matrix named a.

    >> one_plus_a = 1 + a

    Now, one_plus_a contains the quantities of each interest rate plus the value 1. The addition of the scalar value to the matrix of interest rate values (1 + a) works fine. This is because the matrix algebra interpretation of a scalar added to a matrix is to add the scalar to each element of the matrix. Try adding two non-scalar matrices with different dimensions and again an error is the result because that operation is undefined in matrix algebra.

    However, when we tried to raise the quantity (1 + a) to the fifth power, MATLAB produced the error because the matrix dimensions were incompatible for the multiplication operation. Matrix algebra does not allow a 1×N matrix to be multiplied by another 1×N matrix. What we really wanted was each element of the matrix to be raised to the fifth power. This is what element-wise operators do. Thus, replacing the ^ with .^ tells MATLAB to perform the calculation on each element.

  3. How does matrix multiplication work for square matrices in MATLAB?

    Define a square matrix (4×4) and name it x. Then, compute x^2 and x.^2 and store the results in the variables x2_matrix_algebra_mult and x2_elementwise_mult, respectively.

  4. Which operator ^ or .^ will square each interest rate in our computation?

3. Future Value of Money (Part 2)

  1. Compute and plot the future value of money vs the interest rates we assigned.

    The up arrow will work, but the command we want now is several commands away. Another way to repeat a previously entered command is to drag it from the Command History window into the Command Window. If the Command History window is not visible, under the HOME tab in the main menu select the Layout drop-down menu. Then under the SHOW heading, select Command History → Docked.

    1. Type format long and press Enter to show more precision in your answers. Note: You won't see any output after entering this command, but you will notice a difference the next time numbers are displayed (as the result of an unsuppressed command). (Typing and executing the command format short will return to standard precision display -- but don't do this now; we want to stay in the long-precision format.)
    2. Find the future value of money calculation command in the Command History.
    3. Click and drag it to the Command Window.
    4. Edit the command so that it reads FV = instead of F =
    5. Press Enter to execute the statement.
    6. Type the command: disp( [ 'The future value is $' , num2str(FV) , '.' ] )
      Note: The disp function is used to provide more control over the display of output. The num2str function is used to convert a numeric value (FV) to its corresponding string of characters. This is necessary for display with other characters. The square brackets, [ ] , are used to put the character strings into a single matrix of characters.
    7. Press Enter to see a more readable output sentence.
    8. Plot the results by typing plot(a,FV) and pressing Enter.
  2. Which of the following investment options is the better investment?

    1. Invest $2,500 in a 3 year CD at 5.7% and then reinvest your money for 3 years at 5.3%.
    2. Invest $2,500 in a CD at 5.5% for 6 years.

    Caution: Do not reenter commands if you can easily drag them from the Command History or repeat them with the up arrow key. Be sure to edit as necessary before pressing Enter.

4. Publishing your work (to HTML)

You have solved some problems, but need to present your work and the results to your boss or instructor.  MATLAB has several tools to help publish your solutions to different forms. In this class, we ask that you publish your homework solutions to HTML or PDF. 

Note: HTML is currently the most common type of web page found on the internet.  If you copy your HTML pages to a computer folder that is available via the internet (on a web server), then you will be able to see your work from anywhere you have internet access.

Here's how to save your work in an m-file:

  1. Scroll to the top of the Command History and highlight (select) all the commands underneath the date and time stamp for today's lab.
  2. Right-click the mouse and select Create Script. This opens a file Editor window with all of your commands in a list.
  3. If the Editor pops up as a separate window (like the Figure window does when you plot something, as opposed to inside the main MATLAB environment like the Workspace or Command Window), dock the Editor window so that it is part of your MATLAB environment. A quick and easy way to do this is to type Ctrl+Shift+D (alternatively, you can find the Dock Editor option by selecting the icon that looks like a point-down triangle in the upper-right corner of the Editor window immediately to the right of the ? icon).
  4. On the L: drive of your computer, find the MATLAB folder within your My Documents folder. Create a new folder within this MATLAB folder named TeamLab1. (Note: the L: drive is also labeled yourCSlogin(\\wfs1\users$\) and shows up as \\wfs1\users$ > yourCSlogin in the location bar in the MATLAB environment).
  5. Save the M-file as TeamLab1.m within the TeamLab1 folder. Caution: Don't use special characters or spaces in your file names.  Digits and the underscore character, _ , are okay as long as they are not the first character in the file name.

Here's how to publish your work (to HTML) using MATLAB:

  1. Download (save) this image file to your TeamLab1 folder.
    You can do this by right-clicking and selecting Save Image As or Save Picture As.
  2. Add this file header as the first line of the TeamLab1.m file and remove the date/time stamp. You will need to type the blank space between %% and Team Lab 1 for the file header to by published correctly.

    %% Team Lab 1
    % <<../tvom_formula.gif>>

  3. Press F5 (or click the green Run icon in the EDITOR tab) to run your script.  Click Change Folder if prompted to change to the current folder.
  4. Fix any errors in your program script (Editor window). Basically, all of the errors we generated earlier in this lab will now have to be removed from your script.
  5. Go to the PUBLISH tab in the menu bar and click the Publish icon.

Now, let's make our work more user-friendly by adding some cell headers and comments. Be sure to edit, remove, or comment out commands that were errors or duplicates from the original command (code) list.

  1. Add the following cell headers before the appropriate lines in your M-file. 
    Don't omit the single blank space after the %% and before the label or the section headers will not be created correctly.

    %% 1.a Numeric Calculations
    %% 2.a Matrix Multiplication
    %% 2.b Add a scalar to a matrix
    %% 2.c Matrix Squared vs Element-wise operator
    %% 3.a Plot the Future Value of Money vs Interest Rate
    %% 3.b Investment Option

  2. Save the file, publish to HTML, and view the results. 
    Notice how each cell header has become a bold heading with a link at the top of the page. If it has not, review the details in the previous steps to trouble-shoot the cause. You will use cell headers to identify each part of your MATLAB homework solutions.
  3. Remove or edit any commands that caused errors.
  4. Save, publish, and view the results.
  5. Delete and edit commands and headers to make your published work easier to read and follow.
  6. Add single line comments (%) immediately after cell headers and in code.
  7. Save, publish, and view the results.
  8. Show a lab leader your work and ask for tips to improve the readability.

5. Publishing multiple plots

When your work requires multiple plots, make sure that each plot is in its own cell division before publishing to HTML.  Otherwise, only the last plot will be displayed on the final html page. Try getting two or more plot figures to display on your HTML page.

Add a new section to your program script and add this command to plot the first row of your squared x matrix. Note: this command will only work if you named matrices as directed in earlier steps.

>> plot(x(1,:),x2_elementwise_mult(1,:));

Save, publish and view the results to confirm that two separate plots are visible.

6. Displaying equations in the published document page

There are lots of things you can do to make your final document a better presentation of your work. One cool thing MATLAB can do is display professionally formatted equations in your page. Even though MATLAB does not allow you to enter equations and solve them, you can display the equations that your work is based on in the final document.

  1. Add these comment lines immediately after the Numeric Calculations cell header in your M-file. Add all three comment lines.

    %
    % $$F = P(1+a)^n$$
    %

  2. Save and publish to view the results.
  3. Create a new cell division and add these lines immediately after its header.

    %
    % $$y = \frac{\sin(x)}{\cos( \frac{\pi}{4} ) }$$
    %
    x = 1:0.1:10;
    y = sin(x)/cos(pi/4);
    plot(x,y)
    title('sin(x)/cos(\pi/4) vs x');
    xlabel('x');
    ylabel('sin(x)/cos(\pi/4)');

  4. Save and publish to view the results.

Note: This syntax is based on LaTex syntax for entering equations to be formatted. While there is no guarantee that all LaTex commands will work the same in MATLAB as in LaTex, many will and it's worth a try.

7. Displaying lists in the published document

It is easy to display a list of items in a cell comment section. MATLAB can produce bulleted lists and numbered lists.

  1. Add these comment lines immediately after the (top) Page Header cell header in your M-file to create a bulleted list. Add all comment lines.

    %
    % * Name 1
    % * Name 2
    % * Name 3
    %

  2. Save and publish to view the results.
  3. Add these comment lines to create a new cell division and numbered list.

    %% Numbered list example
    % # Launch MATLAB
    % # Create script file
    % # Write MATLAB code
    % # Run script
    % # Save and publish
    % # Print or upload as needed

  4. Save and publish to view the results.

8. Publish your work to PDF

  1. Edit the Publish Configuration options.
    1. Select Edit Publishing Options... from the drop-down menu under the Publish icon.
    2. In Edit Configurations window that pops up, find the section for Output settings.
    3. Under Output settings, change the Output file format field from html to pdf (note: it doesn't look like a drop-down menu, but it is).
    4. Click the Close button.
  2. Click the Publish icon to publish your work as a PDF file.

Show your published PDF to a Lab Leader.

9. Find the published files you created

The published files you created are automatically saved for you according to the default configuration settings or those that you edited. If you are asked to print your work, use File→Print while your published document is showing. If we ask you to upload a PDF, then you will need to know where the PDF file is located.

  1. Open your MATLAB script in MATLAB.
  2. Select Edit Publishing Options... to open the Edit Configuration dialog.
  3. Record (or remember) the location of the output folder for each configuration option you configured.
  4. Open up a folder view for your computer's file system.
    In Windows, [Right-click] Start → Open Windows Explorer.
  5. Navigate through your file system to find the html folder and files you created.
  6. This is where you will need to navigate in order to upload your work to the Learn@UW Dropbox.

10. Getting HELP!

We have presented a lot of different things to try in MATLAB, but what if you need help while working in MATLAB?  Help is just a click or a few typed characters away.

If you know the name of the command you need:
Type help plot to see help displayed in the Command Window.
Type doc plot to launch the Help browser.

If you don't know the name of the command:
Type doc or click the blue question mark near the top of your screen to launch the Help browser and then search for a command or topic that will help you solve the problem.

ADDITIONAL PROBLEMS or FURTHER READING

The solution for this lab will be available shortly after 6 pm on Thursday. Please review and ask questions before starting the quiz. Be sure to complete Quiz 1.