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

  Write the commands that will create the matrix

Write the commands that will create the matrix.

  Deterministic dimensions and rebar area

Performance function for bending of a rectangular reinforced concrete beam under moments due to dead and live loads

  Is the interest rate expressed as a fractional amount

Create a MATLAB® function called future_value with three inputs: the investment (present value), the interest rate expressed as a fraction, and the number of compounding periods.

  Newtons and lagrange interpolation polynomials

What conclusion can you draw from the comparison between Newton's and Lagrange interpolation polynomials? Explain your findings.

  Write a script to integrate a function

Write a script to integrate a function in three different ways: by eyeballing it, by calling trapz with a list of points, and by calling quad with a function.

  Create a project to detect the iris and pupil

Create a Project to do the following: Detect the Iris and pupil, Measure the location "Centroid" of iris/pupil and Measure the diameter of the iris and pupil.

  What is the average growth rates of the economy

What is the average growth rates of the economy along the path (kH, kL, kH kL, kH, ...) and how it compares with the average growth rate along the the balanced growth path?

  Does the threshold increase or decrease with axon diameter

Write down your threshold amplitude and duration values in a table and enter into an mfile. Plot the current versus duration. This is called strength-duration curve.

  Design the digital receiver system

Objective of the project is to design the digital receiver system of the 100-Base TX Ethernet provided in given that the transmitter is provided.

  Write a matlab program that plots x1(t)

Plot (by hand) the signal. Write a MatLab program that plots x1(t) , |t|

  Create a function for calculating british pound

What is the name of the dependant variable and what is the name of the independent variable - Create a function for calculating british pound to danish kroner.

  Prompts the user to enter grades in a loop

Define a function that takes two input parameters, a vector of test scores and a scalar passing grade. The function returns the fraction of passing grades as a number in the range

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