Compute approximations of the largest eigenvalue

Assignment Help MATLAB Programming
Reference no: EM131068140

Problem - We consider the approximations to the eigenvalues and eigenfunctions of the one-dimensional Laplace Operator L[u] = -(d2u/dx2) on the unit interval [0, 1] with boundary u(0) = u(1) = 0. A scalar λ is an eigenvalue of L if there exists a twice -differentiable function u s. t. -u''(x) = λu(x) on [0, 1] with u(0) = u(1) = 0. In this case u is said to be an eigenfunction of L corresponding to the eigenvalue λ. This continuous problem can be approximated by the discrete eigenvalue problem

h-2TNu = λu

where we have set-

518_Figure.png

with ui = u(xi). It can be shown that the N x N matrix TN has eigenvalues vj = 2(1 - cosxj/N+1) for j = 1 : N, corresponding to the eigenvectors uj(k) = √(2/N+1)sin(jkπ/N+1) is the kth entry in uj. Notice that the eigenvectors are normalized w.r.t. the 2-norm. Also notice that eigenvalues of Tn lie in the interval (0, 4). Hence, the eigenvalues of h-2TN lie in the interval (0, 4(N + 1)2).

(a) Express the spectral condition number κ = λN1 of TN (which is the same as the condition number of h-2TN) as a simple function of N for N → ∞.

(b) Use Matlab to plot the eigenvalues of TN for N = 21 (from smallest to largest).

(c) Again using N = 21, use Matlab to plot the eigenvectors uj of TN for j = 1, 2, 3, 5, 11, 21. The plot should be of the form (k, uj(k)) where k = 1 : 21. Note that as j increases, the behavior of the eigenvectors (approximate eigenfunctions) becomes increasingly oscillatory; these eigenfunctions are often called "high energy modes".

(d) Write a Matlab code implementing the power method to compute approximations of the largest eigenvalue λN of L and corresponding eigenfunctions. Use N = 500 and apply the algorithm to h-2TN. Start with a random initial guess, use normalization in the 2-norm, and stope when the difference of two subsequent approximate eigenvectors has 2-norm less than some prescribed tolerance r. (e. g. τ = 10-6).

(i) Report the number of iterations performed

(ii) Compared approximations with exact eigenvalue/eigenvector pair.

(iii) Use matalb to plot the computed eigenvector.

Reference no: EM131068140

Questions Cloud

Component cost of preferred stock : The firm's marginal tax rate (combined federal and state) is 40 percent, and the firm plans to maintain its current capital structure relationship into the future. The component cost of preferred stock to Lei-Feng, Inc. would be closest to.
Million principal outstanding : Now assuming that at the beginning of the 5th year, Tranche A has $20 million principal outstanding, and the CMO trustee receives total cash flows of $36 million, how are they distributed?
Evaluate kanymoon universitys benchmarked position : Evaluate, as far as possible, Kanymoon University's benchmarked position. Define what you understand by the term, 'Strategic management accounting'; critically evaluate any academic literature and refer to recent developments on the subject.
Bankruptcy and agency costs : What is the value of the tax shield if the value of the firm is $5 million, its value if unlevered would be $4.78 million, and the present value of bankruptcy and agency costs is $360,000?
Compute approximations of the largest eigenvalue : Write a Matlab code implementing the power method to compute approximations of the largest eigenvalue λN of L and corresponding eigenfunctions. Use N = 500 and apply the algorithm to h-2TN
Did the ads address issues or candidate perceptions or both : Please briefly summarize one newspaper article related to topics discussed in the chapter and discuss your opinion on the issue. You should copy and paste the article below your response so that other group members can read your article if they wa..
Generate posterior predictive datasets y : Posterior predictive checks: Let's investigate the adequacy of the Poisson model for the tumor count data. Following the example in Section
Brought up provocative ideas in our discussions : Since you have probably already written a number of essays on assigned topics in other English courses, and since many of you have brought up provocative ideas in our discussions, I would like each of you to develop your own topic. I have only one..
Determining the appropriate discount rate : The preferred stock of Rolling Sand Company pays a $4.50 annual dividend. What is the value of the stock if your required rate of return is 10%? Preferred stock of Leaping Dolphin Inc. pays 8% on the $100 par value. What is the value of the stock..

Reviews

Write a Review

MATLAB Programming Questions & Answers

  Finite difference method

Use the finite difference method to calculate the temperature at the point specified since it is easier.

  Determine the necessary shell temperature

In a shell-and-tube heat exchanger, one fluid passes through a central tube while another fluid flows through an outer shell in the opposite direction. The purpose is to heat the fluid passing through the central tube.

  Find the integral of a function at an arbitrary location

Write a Matlab function to perform numerical integration of a set of evenly spaced data points using the trapezoidal rule

  Compute the speed of single-stage planetary gear train

Write a MATLAB function [speed] = planetary (N, emesh, first, last, arm) that computes the speed of a given link in a single-stage planetary gear train.

  Calculate and plot the error in the numerical derivative

Write a program to calculate and plot the error in the numerical estimate of the derivative.

  Create the graph using matlab functions

Create the graph, which contains a piecewise function where a line exists in the first interval, a parabola in the second interval, and the sine function in the third interval.

  Develop a simulation program

Develop a simulation program

  Create a vector in matlab

Create a three dimensional diagram of function.

  Open a named pipe and to read data from the pipe

Open a named pipe and to read data from the pipe in matlab

  Write the commands that will create the matrix

Write the commands that will create the matrix.

  Lagrange interpolating polynomial of degree

Lagrange interpolating polynomial of degree

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