Compute a matrix vector multiplication

Assignment Help Engineering Mathematics
Reference no: EM131299227

Math Programming Exercise Assignment

 Use the finite difference method to approximate the solution of the Poisson problem:

-?u := -∂xxu - ∂yyu = f, in ? = (0, 1) × (0, 1),

  u = 0 ,  on ∂?,

where f is a given function. Specifically, solve the approximation problem

1612_Figure.png

where n ≥ 2 is an integer, h := 1/n, fi,j := f (ih, jh), and wi,j ≈ u(ih, jh). The object wi,j is a called a grid function. It is a discrete approximation of the continuous variable u. This approximation problem can be written as a linear system Aw = f, where A is an (n - 1)2 × (n - 1)2 matrix.

Write a code in Matlab to solve the approximation problem Aw = f using five different iterative schemes: (a) Jacobi, (b) Gauss-Seidel, (c) relaxation, (d) Steepest Descent, (e) Conjugate Gradient, and, finally, (f) preconditioned Steepest Descent using the incomplete Choleski factorization of A.

You should not explicitly form the stiffness matrix A. Instead, if you need to compute a matrix vector multiplication Ax, where x is an arbitrary vector, you should instead use the stencil form as above. In other words, translate x into a grid function with zero boundary values.

The methods above are iterative. They only stop iterating when you supply a stopping tolerance. For this use

649_Figure1.png

where wi,j(k) is the kth iterate and τ is the stopping tolerance, which you should make an input parameter for your code. Compute and plot Rk for each iteration.

In your tests, use n = 64 and n = 128, wi,j(0) = 0, and

f(x, y) = sin3(3πx)cos(πy).

Plot your final iterate for each method.

Reference no: EM131299227

Questions Cloud

Why then does the circuit of figure oscillate : In our study of oscillators, we concluded that a loop containing only two poles cannot oscillate (unless both are at the origin). Why then does the circuit of Figure oscillate?
Discuss the campus newspaper’s arithmetic : The answer given was that "the figures for the fall semester are approximately 12 percent for the girls and 15-18 percent for the guys, which produces a ‘Greek' figure of approximately 27-30 percent of the undergraduates." Discuss the campus newsp..
What would the equilibrium price be : Assume that the market demand for milkshakes is as follows:
Write a research report that cover a description of each erp : Write a minimum of three page research report, using current APA formatting. The report should cover the previously listed MS Dynamics ERP systems, a description of each ERP and list and description of each associated models.
Compute a matrix vector multiplication : Math 571 Programming Exercise Assignment. You should not explicitly form the stiffness matrix A. Instead, if you need to compute a matrix vector multiplication Ax, where x is an arbitrary vector, you should instead use the stencil form as above. In..
Critical slides of your presentation and placement : Create a 12-slide presentation. Follow the outline on pp. 343-346 for the critical slides of your presentation and their placement - Did not submit or incompletely followed the outline on pp. 343-346 for the critical slides of your presentation and..
Must the clock transition be abrupt for the d latch : Must the clock transition be abrupt for the D latch of Figure to operate properly? Consider a clock transition time on the order of the time constant at X and Y.
How about ethernet sata firewire usb media cards : How about Ethernet, SATA, FireWire, USB (2.0 or 3.0), media cards? Think of the data transfer/exchange requirements and what kind of speeds are necessary to make them work effectively.
What is the population for this sample survey : Find the margin of error. Should the commissioner of baseball punish Major League Baseball players named in the Mitchell report on steroid use in baseball for having used steroids? A poll of 413 professional baseball fans found that 248 said "Yes...

Reviews

Write a Review

Engineering Mathematics Questions & Answers

  Prime number theorem

Dirichlet series

  Proof of bolzano-weierstrass to prove the intermediate value

Every convergent sequence contains either an increasing, or a decreasing subsequence.

  Antisymmetric relations

How many relations on A are both symmetric and antisymmetric?

  Distributed random variables

Daily Airlines fies from Amsterdam to London every day. The price of a ticket for this extremely popular flight route is $75. The aircraft has a passenger capacity of 150.

  Prepare a system of equations

How much money will Dave and Jane raise for charity

  Managing ashland multicomm services

This question is asking you to compare the likelihood of your getting 4 or more subscribers in a sample of 50 when the probability of a subscription has risen from 0.02 to 0.06.]  Talk about the comparison of probabilities in your explanation.

  Skew-symmetric matrices

Skew-symmetric matrices

  Type of taxes and rates in spokane wa

Describe the different type of taxes and their rates in Spokane WA.

  Stratified random sample

Suppose that in the four player game, the person who rolls the smallest number pays $5.00 to the person who rolls the largest number. Calculate each player's expected gain after one round.

  Find the probability density function

Find the probability density function.

  Develop a new linear programming for an aggregate production

Linear programming applied to Aggregate Production Planning of Flat Screen Monitor

  Discrete-time model for an economy

Discrete-time model for an economy

Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd