gibson:teaching:spring-2018:math445:lab2
====== Differences ====== This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
gibson:teaching:spring-2018:math445:lab2 [2018/01/30 02:59] gibson |
gibson:teaching:spring-2018:math445:lab2 [2018/01/31 16:07] (current) gibson [Part 3: Plots] |
||
|---|---|---|---|
| Line 1: | Line 1: | ||
| - | ====== Math 445 lab 2: plotting and scripts====== | + | ====== Math 445 lab 2: vectors, plots, and scripts====== |
| - | Do the following in Matlab, saving your work to a file with the ''diary'' command, and labeling each problem with a Matlab comment (e.g. ''% problem 1''). Reduce unnecessary whitespace with the ''format compact'' command. When you're done, edit the diary file to eliminate errors, so that your diary shows just the correct work. You do not need to turn in the plots of 12, 13, and 14. Just describe them verbally using Matlab comments. For 15, just provide your estimates of the zeroes of the polynomial. Turn in your diary file on Canvas. | + | Do the problems for parts 1 and 2 at the Matlab prompt, saving your work to a file with the ''diary'' command. |
| - | 1. Make a plot of a unit circle, i.e. a curve that satisfies $x^2 + y^2 = 1$. Draw the circle with a thick red line, label the axes, and give the plot a title. Make sure the circle looks like a circle and not an oval. Hint: don't try to draw the plot using the equation $x^2 + y^2 = 1$. Instead parameterize the curve in terms of the angle $\theta$, i.e. calculate $x$ and $y$ in terms of a parameter $\theta$. | + | Be kind to your teaching assistant! |
| + | * Label each problem with a Matlab comment (e.g. <nowiki>% ==== Problem 1 ====</nowiki>). | ||
| + | * For questions with free-form response, write your answer as a Matlab comment. Write in complete sentences! | ||
| + | * Reduce unnecessary whitespace with the ''format compact'' command. | ||
| + | * When you're done, edit the diary file to eliminate errors, so that your diary shows just the correct work. | ||
| + | * Add and subtract blank lines as necessary to make your work organized and easiy readable. | ||
| - | 2. Make a plot of the polynomial $f(x) = x^3 -5x^2 + 2x + 3$ as a blue line. In the same plot, draw a line along the $x$ axis. (Hint: in Matlab create a zero vector of the same length as your $x$ vector with ''y=0*x'', then plot $y$ versus $x$ along with $f$ versus $x$.) Using this plot, estimate the zeros of the polynomial, i.e. the values of $x$ for which $f(x) = 0$. Make sure to find all the zeros of $f$ by adjusting limits of the plot until all intersectionf of $f$ with the $x$-axis are visible. | + | ==== Part 1: Vector manipulation ==== |
| + | Relevant Matlab | ||
| + | <code Matlab> | ||
| + | [] % square brackets | ||
| + | , % row separator | ||
| + | ; % col separator | ||
| + | : % colon operator | ||
| + | ' % apostrophe | ||
| + | linspace | ||
| + | </code> | ||
| - | 3. | + | **1.** Assign a row vector with elements 3, 4, 5, 9 to the variable u. |
| - | 1. Assign a row vector with elements 3, 4, 5, 9 to the variable u. | + | |
| - | 2. Assign the transpose of u to the variable w. | + | **2.** Assign the transpose of u to the variable w. |
| - | 3. Change the third element of w to 10. | + | **3.** Change the third element of w to 10. |
| - | 4. Assign a column vector with the elements 7, 1, -2, 3 to the variable z. | + | **4.** Assign a column vector with the elements 7, 1, -2, 3 to the variable z. |
| - | 5. Add w and z. Does the result make sense? | + | **5.** Add w and z. Does the result make sense? |
| - | 6. Add u and z. Does the result make sense? | + | **6.** Add u and z. Does the result make sense? |
| - | 7. Create a vector of the even integers between 2 and 14, inclusive. | + | **7.** Create a vector of the even integers between 2 and 14, inclusive. |
| - | 8. Create a vector of the odd integers between 7 and 19, inclusive. | + | **8.** Create a vector of the odd integers between 7 and 19, inclusive. |
| - | 9. What is the sum of every third number between 3 and 27, inclusive? | + | **9.** Make x a vector from 0 to $2\pi$ in increments of 0.02. |
| - | 10. What is the sum of the integers between 1 and 100, inclusive? | + | ==== Part 2: Computing sums ==== |
| - | 11. Make x a vector from 0 to $2\pi$ in increments of 0.02. | + | Relevant Matlab syntax |
| + | <code matlab> | ||
| + | : % colon operator | ||
| + | sum % sum function | ||
| + | .* ./ .^ % dot syntax | ||
| + | </code> | ||
| - | 12. Run ''plot(x, sin(x))''. What happens? | ||
| - | 13. Run ''plot(x, sin(x), 'r--')''. What happens? | ||
| - | 14. Run | + | **10.** What is the sum of every third number between 3 and 27, inclusive? |
| + | **11.** What is the sum of the integers between 1 and 100, inclusive? | ||
| + | |||
| + | **12.** The infinite series | ||
| + | |||
| + | \begin{equation*} | ||
| + | \sum_{n=1}^{N} \frac{1}{n} | ||
| + | \end{equation*} | ||
| + | |||
| + | diverges as $ | ||
| + | N \rightarrow \infty$. Demonstrate this in Matlab by computing the sum for several values of $N$, e.g. $N=10$ to $N=10^5$ by powers of ten. | ||
| + | |||
| + | **13.** The $\sin$ function can be calculated from the infinite series | ||
| + | |||
| + | \begin{equation*} | ||
| + | \sin x = \sum_{n=0}^{\infty} (-1)^n \frac{x^{2n+1}}{(2n+1)!} | ||
| + | \end{equation*} | ||
| + | |||
| + | Of course we have to truncate this infinite sum to a finite number of terms in order to calculate it on a computer. How many terms do you need to keep in order to compute $\sin \: \pi/2 = 1$ to sixteen digits accuracy? | ||
| + | Hint: use ''format long'' to see all sixteen digits of the computation. | ||
| + | |||
| + | Are you amazed or what? | ||
| + | |||
| + | |||
| + | **14.** Now use the same infinite series for $\sin x$ to calculate $\sin \: \pi = 0$. How many terms do you need to keep in order get the correct answer to sixteen digits accuracy? | ||
| + | |||
| + | Are you perplexed or what? What is strange about this calculation? Can you explain what happened? | ||
| + | |||
| + | |||
| + | ==== Part 3: Plots ==== | ||
| + | |||
| + | Relevant Matlab | ||
| <code matlab> | <code matlab> | ||
| - | plot(x, sin(x), 'r--', x, cos(x), 'm-') | + | .* ./ .^ % dot syntax |
| - | xlabel('x') | + | plot % plot function |
| - | legend('sin x', 'cos x') | + | axis % set limits on plot or aspect ratio |
| + | xlim % set x limits | ||
| + | ylim % set y limits | ||
| + | grid % turn grid on/off | ||
| + | help % help function | ||
| </code> | </code> | ||
| + | **15.** Make a plot of the polynomial $f(x) = x^3 -5x^2 + 2x + 3$ as a blue line. In the same plot, draw a line along the $x$ axis. (Hint: in Matlab create a zero vector of the same length as your $x$ vector with ''y=0*x'', then plot $y$ versus $x$ along with $f$ versus $x$.) Using this plot, estimate the zeros of the polynomial, i.e. the values of $x$ for which $f(x) = 0$. Make sure to find all the zeros of $f$ by adjusting limits of the plot until all intersection of $f$ with the $x$-axis are visible. | ||
| + | |||
| + | |||
| + | **16.** Make a plot of a unit circle, i.e. a curve that satisfies $x^2 + y^2 = 1$. Draw the circle with a thick red line, label the axes, and give the plot a title. Make sure the circle looks like a circle and not an oval. Hint: don't try to draw the plot using the equation $x^2 + y^2 = 1$. Instead parameterize the curve in terms of the angle $\theta$, i.e. calculate $x$ and $y$ from $\theta$. | ||
| + | |||
| + | (I previously had a problem 17, revisiting the sum of $1/n$...to be written. But on second thought, this problem belongs in lab 3. So expect it there.) | ||
gibson/teaching/spring-2018/math445/lab2.1517309948.txt.gz · Last modified: 2018/01/30 02:59 by gibson
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 3.0 Unported
