Scribd
Upload
48 views

Tutorial 1

(1) This document provides instructions for 4 exercises introducing numerical methods and MATLAB/Octave programming. Students are asked to practice representing real numbers in different formats, performing basic operations, plotting functions, finding intersections, and creating a text file with a value table. (2) The exercises involve building matrices using loops and vector operations, computing matrix-vector products, plotting sin and cos curves with different step sizes, finding the intersections of 3 functions, and making a 3D surface plot. (3) Students are instructed to complete the exercises and submit their MATLAB code by the specified deadline.

Uploaded by

anon_545085293
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
48 views

Tutorial 1

(1) This document provides instructions for 4 exercises introducing numerical methods and MATLAB/Octave programming. Students are asked to practice representing real numbers in different formats, performing basic operations, plotting functions, finding intersections, and creating a text file with a value table. (2) The exercises involve building matrices using loops and vector operations, computing matrix-vector products, plotting sin and cos curves with different step sizes, finding the intersections of 3 functions, and making a 3D surface plot. (3) Students are instructed to complete the exercises and submit their MATLAB code by the specified deadline.

Uploaded by

anon_545085293
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Universität Duisburg-Essen Fall Term 2010/11

Computational Mechanics
Campus Essen

Dr. S. Vanis 14. October 2010

Introduction to Numerical Methods


Tutorial 1
This first tutorial wants to introduce you to the matlab/octave programming
language. These exercises are not counted as programming exercises. The pro-
gramming exercises will be explicitly tagged out as such.

Exercise 1:
(i) Give the representation of the real numbers 43 , 0.0000012345 and
1.2345e−06 in the following formats
(1) short ; short e ; short g
(2) long ; long e ; long g
(3) rat
and explain the results. What is the meaning of ”e” in this context?
(ii) What happens if you use the following commands in matlab/octave
(1) ”f 1 = 3 ∗ 3 ; ”
(2) ”f 2 = 3ˆ3 ; ” ”g1 = [3 , 1; 1 , 3]ˆ3”
(3) ”f 3 = 3.ˆ2 ” ”g2 = [3 , 1; 1 , 3] .ˆ3”
(4) ”3 + 3 ”
(5) ” disp(’Hallo world!’); disp(’Introduction to Numerical Methods’); ”
What is the difference between g1 and g2, cf., (2) and (3) ?

Exercise 2:
Let A = (aij )i,j=1...n be the matrix with entries aij = i/j and x the vector xj = j
for j = 1 . . . n, i.e.,
1 21 1
· · · n1
   
3 1
2 .. 2
 2 1 . n   2 
   
3
A=
 .. ..   3 
 x= .
3 . . n3 
 3 2  .. 
 . .
 .. . . . . . . . . ...   . 

n n2 · · · · · · 1 n
(i) Try three different ways to build the matrix A, i.e.,

(1) using 2 nested for-loops


(for k=1... for l=1... A(k,l)=k/l; end end)
(2) using one vector and one for-loop
(k=1:n; for l=1... A(:,l) = k’/l; end)
(3) using the product of two vectors
(k=(1:n)’; l=ones(1,n)./k’; A = k * l;)

With the commands tic and toc the time needed for the operations in bet-
ween these commands can be determined (tic; operations; toc;). Use these
commands to determine the time matlab/octave needs to build the matrix
A in the ways described above for n = 100, 500, 1000, 1500, 2500(, 5000).
What do you recognize?
Hint: Write an m-file which takes the size of the matrix n as input and in
which all three ways are included one after another. Determine the separate
times by enclosing each part in tic ... toc.

(ii) Build the matrix A and the vector x in matlab/octave and compute the
matrix-vector product y = Ax for n = 3, 5, 10, 100, 1000.
Verify the solution by your theoretical knowledge on matrix-vector multi-
plication.

(∗)-Exercise 3: (2 + 5 + 3 = 10 points)

(i) Build a vector x which starts with 0 and ends at 2π and use it to plot with
matlab sin(x) and cos(x) on the interval [0, 2π] ⊂ R. Plot the functions
once into two separate diagrams and once into only one diagram.
Use different step sizes for the vector x, i.e., step sizes 1; 0.5; 0.1, and com-
pare the results. Just plot the sin-curves for all three step sizes into a single
diagram.

(ii) Use the possibility to plot different graphs into one diagram to obtain the
intersection points of the following three functions on the interval [−2, 2] ⊂
R

• exp(−x2 ) ,
• sin(x2 ) ,
1
• x3 +9
.

Just use the ”Data Cursor” to select the intersection points and read off
the values.
Try different step sizes in this case as well. Use the command legend (in
”matlab”, when you use ”octave” you have to find an alternative to this
command), to label the different functions in the diagram.
(iii) Use the command mesh(x,y,z) to make a 3 dimensional diagramm for x ∈
[−2, 2], y ∈ [−2, 2] and z(x, y) = x2 − y 2 . Choose a step size of 0.1 in both
directions, i.e., in x-direction as well as in y-direction.

Exercise 4:
Create a file ”Name.txt” (where you replace name by your first name). Write an
arbitrary sentence in the file. Leave two lines free and write a value table for the
function exp(x3 ) on the interval [−1.5, 3.5] ⊂ R using the step size 0.25.

Delivery: The exercises have to be handed in before the tutorial on Thursday,


21 October 2010. Exercises which are handed in later may be corrected but
will give no points for the (∗)- and programming exercises. The source codes have
to arrive by E-Mail not later than 12 o’clock on the date of delivery.

You might also like