Ask Question, Ask an Expert


Ask MATLAB Expert

1 Errors

a. Run the MATLAB script

sum = single(0); term = single(1);
while sum + term > sum
sum = sum + term;
term = 1/n;

Describe what you think it is trying to do. What would this code do in exact arithmetic? describe why the code produces a finite result and (harder) describe the output. What would happen (don't try it!) if you replaced the  first line with sum=0; term=1;?

b. Use Taylor series to derive the form of the truncation error of the approximation

f'(x) =3f(x) -4f(x - h) + f(x -2h)/2h

assuming f ≡ C3. Given that the round off error has the same form as in the Forward Difference case, find an expression for the step size h that gives the smallest total (absolute) error. Modify ForwardDifference.m to show how the total error in using Equation 1 behaves for various values of h and comment on your results.

2 Operation counts

a. Derive the operation count (multiplications and divisions only) for

(i) Choleski factorization (without pivoting) A detailed algorithm for Choleski factorization can be found via the MAST30028 website on p. 78 of Heath's Lecture Slides --- Chapter 2.

(ii) Gauss elimination (equivalent to LU factorization) without pivoting of an Upper Hessenberg matrix. An Upper Hessenberg matrix is one with zeroes below the subdiagonal. Just modify the derivation for Gauss elimination.

(iii) solution of a tridiagonal system via LU factorization without pivoting (also called the Thomas algorithm) This is the algorithm implemented in the code tridisolve in the Asst1 folder on the server.

b. Verify qualitatively the running times you obtained above, by timing the solution of linear systems of size at least 1000X1000 with matrices having the appropriate structure (symmetric positive definite, Upper Hessenberg or tridiagonal).

Construct the matrices yourself with suitable MATLAB commands such as diag, triu, rather than using particular matrices from the MATLAB gallery.

For all cases, use \ to solve the system. For case c, try a tridiagonal matrix created with diag and one created with spdiags.

describe what you think is going on.

3 Application: Nonlinear systems of equations

You have met Newton's method as a method for solving nonlinear equations in the form f(x) = 0.

xn+1 = xn-f(xn)/f'(xn)

which gives iterates {xn} that converge rapidly to the solution x*  provided the initial guess x0 is sufficiently close to x* . Applying this to the function

f(x) = x2 - a produces the iteration xn+1 = xn - x2n -a/2xn= xn/2 + a/(2xn)

which should look familiar. It can be written in the alternative form: solve in turn

f'(xn)sn = f(xn); xn+1 = xn - sn

For a system of two nonlinear equations f1(x1, x2) = 0; f2(x1,x2) = 0

or  in vector notation

f (x) = 0

Newton's method generalizes in the following way: solve the 2X2 linear system

J(xn)sn = f (xn) then update xn+1 = xn - sn

where J is the Jacobian matrix formed from f

194_Jacobian matrix.jpg

To get an idea how it works, try out the Java applet

a. prepare a MATLAB function to implement Newton's method for a system of 2 nonlinear equations.

It should have the calling sequence [roots, count, resids, history] = ass1Q3(func,x0,tol)


• the first input argument is a function handle to a function that returns the value of the function f and the Jacobian matrix J at the position x e.g. function [ f J] = ex2( xx )

The function ex2 needs to be written by the user, including the expressions for f and J. There is no need for J to be found symbolically from f .

• x0 is the initial guess ( a column vector)

• tol is the absolute tolerance. The iterations should stop once the residual ||f||, measured using the vector ∞-norm, is below tol.

• the  first output argument is a column vector containing the solution x

• count is the number of iterations required

• resids is a vector containing the norms of the residuals at each iteration

• history is a matrix whose columns are the iterates xn i.e. each column is a 2X1 array.

The function should:

• use a default tolerance of 10-10 if none is given

• use a while loop but loop no more than 50 times.

Although I only need it to work for systems of 2 equations, you should be able to prepare it so that it can be used (in principle) on a system of any size. My code is less than 20 lines long.

b. prepare a driver function to test your code on ex 2 of the Newton2D Java applet. The driver should produce output such as
After 5 iterations, the roots are [1.61803 1.61803]

c. Modify your driver to find the solution to

3x31 + 4x22- 145 = 0; 4x21- x32+ 28 = 0

Show results for 2 different initial guesses: x0 = (2, 2),(2,-1).

For each case, plot the convergence history - the residual norms at each iteration - on a suitable plot and (on a different plot) the solution trajectory - the set of iterates xn.
Comment on your results.

MATLAB, Engineering

  • Category:- MATLAB
  • Reference No.:- M91787
  • Price:- $70

Priced at Now at $70, Verified Solution

  • AsyU replied is my first experience and i have to say that you guys are life saver honestly your solution is too awesome which helped me alot thank you so much...

Have any Question? 

Related Questions in MATLAB

Problem 1 the two foremost aerospace companies in the usa

PROBLEM 1: The two foremost aerospace companies in the USA are Lockheed and Boeing. They often compete for government contracts. A major factor in awarding a contract is the claimed reliability of the given system to be ...

Algorithms assignmentthe following 3 tasks are to be

Algorithms Assignment The following 3 tasks are to be completed for Homework 1 - using the following MATLAB syntax for the script file "bike.m" 1) Write additional syntax to create an INPUT variable named time and prompt ...

Discrete communicationsobjectiverevision of fundamental

Discrete Communications Objective: Revision of fundamental concepts and demonstration of necessary foundation skills. Assessment: The assignment will comprise 20% of your final mark and all ten problems will be of equal ...

1 you are testing a theory which says that the output

1. You are testing a theory which says that the output displacement of a measurement device should vary with time according to the equation: D(x) = {x 2 cos (x)} / (x +1) So far, you have measured the following set of ex ...

Part -1introduction to programming with matlab1 the

Part -1: Introduction to Programming with MATLAB 1. The function move_me is defined like this: function w = move_me(v,a). The first input argument v is a row-vector, while a is a scalar. The function moves every element ...

Pre-labobjectiveslearn how to use various matlab built-in

Pre-Lab Objectives ? Learn how to use various MATLAB built-in functions ? Become familiar with vector and matrix multiplication in MATLAB The input function is a MATLAB built-in function that allows a program to accept u ...

Project 1complete parts 1a 1b and 1cprepare a report

Project 1 Complete parts 1A, 1B, and 1C Prepare a report describing your steps and results Send by email a zip or rar file that contains the report and all MATLAB files you have used. Project 1A: In this project, we cons ...

Qestion 11 introductionmeasurements have been made of a

Question 1: 1 Introduction: Measurements have been made of a cantilever beam's deflection (stored in the file ass2q1in.csv). You must determine what type of loading has been applied and the magnitude of the load. There a ...

Assignment - instructions for the overall material

Assignment - Instructions for the Overall Material Balance Work with your learning pod to perform an overall material balance for fuel-grade ethanol production.  This balance is much simplified, but will help you to revi ...

1 write a matlab script that will solve the system of

1. Write a Matlab script that will solve the system of linear equations shown below. This can be done by solving for the vector of unknowns (X) in the equation AX=B where A is the coefficient matrix and B is the vector o ...

  • 4,153,160 Questions Asked
  • 13,132 Experts
  • 2,558,936 Questions Answered

Ask Experts for help!!

Looking for Assignment Help?

Start excelling in your Courses, Get help with Assignment

Write us your full requirement for evaluation and you will receive response within 20 minutes turnaround time.

Ask Now Help with Problems, Get a Best Answer

A cola-dispensing machine is set to dispense 9 ounces of

A cola-dispensing machine is set to dispense 9 ounces of cola per cup, with a standard deviation of 1.0 ounce. The manuf

What is marketingbullwhat is marketing think back to your

What is Marketing? • "What is marketing"? Think back to your impressions before you started this class versus how you

Question -your client david smith runs a small it

QUESTION - Your client, David Smith runs a small IT consulting business specialising in computer software and techno

Inspection of a random sample of 22 aircraft showed that 15

Inspection of a random sample of 22 aircraft showed that 15 needed repairs to fix a wiring problem that might compromise

Effective hrmquestionhow can an effective hrm system help

Effective HRM Question How can an effective HRM system help facilitate the achievement of an organization's strate