'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. EGR 103L Fall Test 2. Michael R. Gustafson II

Size: px

Start display at page:

Download "'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. EGR 103L Fall Test 2. Michael R. Gustafson II"

Esmond O’Neal’
5 years ago
Views:

1 'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ EGR 103L Fall 2017 Test 2 Michael R. Gustafson II Name (please print) NET ID (please print): In keeping with the Community Standard, I have neither provided nor received any assistance on this test. I understand if it is later determined that I gave or received assistance, I will be brought before the Undergraduate Conduct Board and, if found responsible for academic dishonesty or academic contempt, fail the class. I also understand that I am not allowed to speak to anyone except the instructor about any aspect of this test until the instructor announces it is allowed. I understand if it is later determined that I did speak to another person about the test before the instructor said it was allowed, I will be brought before the Undergraduate Conduct Board and, if found responsible for academic dishonesty or academic contempt, fail the class. Signature: Notes You will be turning in each problem in a separate pile. Most of these problems will require working on separate pieces of paper - Make sure that you do not put work for more than any one problem on any one piece of paper. For this test, you will be turning in four different sets of work. Again, Please do not work on multiple problems on the same sheet of paper. Also - please do not put work for one problem on the back of another problem. Be sure your name and NET ID show up on every page of the test. If you are including work on extra sheets of paper, put your name and NET ID on each and be sure to staple them to the appropriate problem. Problems without names will incur at least a 25% penalty for the problem. This first page should have your name, NET ID, and signature on it. It should be stapled on top of and turned in with your submission for Problem I. You will not need and can not use a calculator on this test. For hand calculations you will instead show the set-up of the calculation you need to perform but will not reduce it. You can leave powers, roots, products, sums, differences, and quotients unevaluated - for example, (3)(2) + (8)(11) (6)(24) A = (8) 2 (4) 2 + (2) 2 (7) 2 should be left just like that. You will be asked to write several lines of code on this test. Make sure what you write is MATLAB code and not mathematics. Be very careful with any symbols you use. You do not need to put the honor code statement in your codes. The honor code statement on this page and your NET ID on each problem stands in for that. You do not need to number your lines of code. Please staple in the top left corner of the page. Also, please do not put work too close to the top left corner of the page!

2 Problem I: [25 pts.] This again? (1) Given the following matrices: M = [ ] N = [ ] Oh = Write the results of the following MATLAB commands in the spaces below. (a) A = M*Oh (f) F = M>1 (b) B = N. 2 (g) G = find(m>1) (c) C = N 2 (h) H = M(find(M>1)) (d) D = N(:) (i) I = [N N] (e) E = sum(oh) (j) J = N & N (2) Show what P and R would be after the following command: [P, R] = meshgrid([1 2 3], [4 5]) (3) Show what S, T, U, and V and would be after the following commands: S = [5 2 3; 1 6 8] T = 2; for U=S T=T+1; S(T)=T; V(T)=max(U); end

3 Problem II: [25 pts.] Where is everything? (1) Given the equation: and the two graphs: y(x) = e x x 2 + 8x 5 x y y x x which are of the same function but with two different domains, write MATLAB code to accomplish the tasks below. (a) Write code guaranteed to find the two positive values of x where y(x) = 0. Note that the root at x = 0 is not included in these two. Call these values x1 and x2 (b) Write code guaranteed to find the value and location of the local minimum of the function on the left side of each graph. Assign the value of the minimum to ymin and the location of the minimum to xminloc. (c) Write code guaranteed to find the value and location real global maximum value of the function. Assign the value of the minimum to ymax and the location of the minimum to xmaxloc. (d) From the graphs, you can see that y(x) = 2 at two locations. Write code guaranteed to find both and assign the locations to x2a and x2b from left to right, respectively. (2) Given the function of a surface: f(x,y) = x 2 + x y + y 2 3x + 2y + cos(x) + sin(y) (a) Generate a surface plot of f(x,y) as a function of x and y. x and y should span the domain from -8 to 8 and there should be 50 nodes in each direction. You do not need to label, title, or save this plot. (b) Determine which of the 2500 modes is closest to and find the x, y, f(x,y), and element values at that grid point. Display this information in a manner similar to the following (numerically incorrect) example (which uses %+0.3e for all scientific notation numbers and %0.0f for all integers, as do all in this problem): Max node: f(-7.674e+00,-7.347e+00) = e+02 at node 17 (c) Determine which of the 2500 nodes is closest to and find the x, y, and f(x,y) values at that grid point. Display this information in a manner similar to the following (numerically incorrect) example: Min node: f(+2.976e+00,-2.022e+00) = e+00 at node 1841 (d) Determine the exact location of the minimum of the f(x, y) function. Use the information gained above in (c) as an initial guess. Display this information in a manner similar to the following (numerically incorrect) example: Min val: f(+2.849e+00,-2.015e+00) = e+00 (e) Determine the S t value for the 2500 values of f(x,y) and store it in a variable called St.

4 Problem III: [25 pts.] Fall & Spring(s) (1) Assume you have the following equations for a system: Hr = Is + J Kr + Ls + M = 0 and assuming H through M are known constants while r and s are your unknowns, (a) Fill in the following framework for a linear algebra system (there are six total blanks): [A]{x} = {b} [ ]{ } { } r = s (b) By hand, determine a formula for the determinant of the coefficient matrix, A. (c) By hand, determine a formula for the determinant of the inverse of the coefficient matrix, A 1. (d) By hand, determine solutions to the vector of unknowns r and s in terms of known values H through M. (2) Three masses M i are held apart at (unknown) positions x i between two walls by a series of large springs with known spring constants K 1 through K 4. There are known external forces F i acting on each of the masses. x 1 x 2 x 3 F 1 F 2 F 3 M 1 M 2 M 3 K 1 K 2 K 3 K 4 The equilibrium equations for the positions of the masses as functions of the applied forces are: (K 1 + K 2 )x 1 = F 1 + K 2 x 2 (K 2 + K 3 )x 2 = F 2 + K 2 x 1 + K 3 x 3 (K 3 + K 4 )x 3 = F 3 + K 3 x 2 (a) Fill in the following framework for a linear algebra system that could be used to solve for the positions (x 1, x 2, and x 3 ); there are 12 total blanks: x 1 x 2 x 3 = (b) Assuming you know that K 1 =K 2 =K 3 =1000 and F 1 =F 2 =F 3 =2000, write code that will allow you to solve for and plot the values of x 2 for 200 linearly spaced values of K 4 between 500 and Give your graph proper axis labels and use a dashed blue line. You may assume the code clear; K1 = 1000; K2=1000; K3=1000; F1=2000; F2=2000; F3=2000; is already written. (c) Assume your system measurements have six significant figures. It turns out the condition number of your coefficient matrix, based on the 2-norm, is 5.05 when K 4 = What does this mean for the solutions to your linear algebra problem when K 4 = 2000?

5 Problem IV: [25 pts.] Under Pressure This problem is adapted from Chapra The dependent variable has been transformed to something called LP versus the original p. The Antoine equation describes the relation between vapor pressure and temperature for pure components as LP = A B C + T where LP is related to the pressure, T is temperature, and A, B, and C are component-specific constants. Use MATLAB to determine the best values for the constants for carbon monoxide based on the following measurements where you may assume the following lines of code are given: T = [ ]; LP = [ ]; Write the rest of the code to: (1) Determine the coefficients A, B, and C that give the mathematically best version of the Antoine equation for this data set. At the end of this part of the code, you should have variables named A, B, and C. Of the three coefficients, in looking at how the formula changes with T, you should be able to come up with a guess for A based on the data. Describe your process for making a guess of A. Use 50 as an initial guess for B and -2 as an initial guess for C. (2) Determine the r 2 value for this fit and call it r2a for Antoine. (3) Ask the user for a single integer called Ord between 1 and 8. Keep asking the user for an input until the user correctly gives a single integer with a value between 1 and 8 (inclusive). (4) Determine the coefficients for a polynomial fit or order Ord. Store them in a variable called Coefs2. (5) Determine the r 2 value for this polynomial fit and call it r2p for Polynomial. (6) Make a graph that has the original data points as black circles, the Antoine equation model as a solid blue line that uses 200 values of T that span the domain of the original data set, and the polynomial fit as a dashed red line that also uses 200 values of T that span the domain of the original data set. Put proper axis labels on this graph along with a meaningful legend. (7) Print a statement about which fit was mathematically better. Your program should either say: Antoine fit is better: r2=9.993e-01 vs e-01 or Order 4 polynomial fit is better: r2=9.997e-01 vs e-01 where you need to indicate the order of the polynomial, or Fits have same r2=9.993e-01 depending on which is better or if they are the same. Use %0.0f for integers and %0.3e for other numbers.

'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. EGR 103L Fall Test 2 Solutions. Michael R. Gustafson II

$'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. EGR 103L Fall Test 2 Solutions. Michael R. Gustafson II$ 'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ EGR 103L Fall 2017 Test 2 Solutions Michael R. Gustafson II Name (please print) NET ID (please print): In keeping with the Community Standard, I have neither