Consider the equation Δυ = = Uxx + Uyy = f on D = [0, 1] × [0,1] u = 0 on aD. Approximate values of u at the set of grid points (xi, Yj), i, j = 0,...,n+1 will be denoted by Ui, j ≈ u(xi, yj). The grid points cover the domain D, so that, calling h the mesh-size, we have (n + 1)h = 1. A common second-order approximation of the Laplacian is (1) Ui+1,j + Ui−1,j + Ui, j+1 +Ui, j-1 - 4ui,j Au≈ (Au)i,j == h² Setting ui,o the form = Ui,n+1 = uo,j = Un+1,0 = where uT = (u1,1, U2,1, , un, 1, u1,2, An is an n² x n² matrix. (2) = 0, the discrete analog of equation (1) can be expressed in Anu = f, Un,n), ƒT = (f1,1, f2,1, · · ·, fn,1, f1,2, . . ., fn,n) and where For this assignment take n = 10, 20, 40 and some other values you consider interesting. For each value form the matrix An and use appropriate MATLAB commands to evaluate the eigenvalues of An. Consider the following questions: (1) What are the eigenvalues of the continuous Laplacian with vanishing Dirichlet boundary conditions? (The eigenvectors are of the form sin(kлx) sin(lπy).) (2) What are the minimum and maximum eigenvalues of An. How do they change with n? Do they seem to be converging to the eigenvalues of the continuous Laplacian? If so, at what rate? Is this expected?
Consider the equation Δυ = = Uxx + Uyy = f on D = [0, 1] × [0,1] u = 0 on aD. Approximate values of u at the set of grid points (xi, Yj), i, j = 0,...,n+1 will be denoted by Ui, j ≈ u(xi, yj). The grid points cover the domain D, so that, calling h the mesh-size, we have (n + 1)h = 1. A common second-order approximation of the Laplacian is (1) Ui+1,j + Ui−1,j + Ui, j+1 +Ui, j-1 - 4ui,j Au≈ (Au)i,j == h² Setting ui,o the form = Ui,n+1 = uo,j = Un+1,0 = where uT = (u1,1, U2,1, , un, 1, u1,2, An is an n² x n² matrix. (2) = 0, the discrete analog of equation (1) can be expressed in Anu = f, Un,n), ƒT = (f1,1, f2,1, · · ·, fn,1, f1,2, . . ., fn,n) and where For this assignment take n = 10, 20, 40 and some other values you consider interesting. For each value form the matrix An and use appropriate MATLAB commands to evaluate the eigenvalues of An. Consider the following questions: (1) What are the eigenvalues of the continuous Laplacian with vanishing Dirichlet boundary conditions? (The eigenvectors are of the form sin(kлx) sin(lπy).) (2) What are the minimum and maximum eigenvalues of An. How do they change with n? Do they seem to be converging to the eigenvalues of the continuous Laplacian? If so, at what rate? Is this expected?
Linear Algebra: A Modern Introduction
4th Edition
ISBN:9781285463247
Author:David Poole
Publisher:David Poole
Chapter4: Eigenvalues And Eigenvectors
Section4.2: Determinants
Problem 13AEXP
Related questions
Question
help me with part (1) and (2) please. please also provide the matlab code
Expert Solution
This question has been solved!
Explore an expertly crafted, step-by-step solution for a thorough understanding of key concepts.
Step by step
Solved in 2 steps with 5 images
Recommended textbooks for you
Linear Algebra: A Modern Introduction
Algebra
ISBN:
9781285463247
Author:
David Poole
Publisher:
Cengage Learning
Elementary Geometry for College Students
Geometry
ISBN:
9781285195698
Author:
Daniel C. Alexander, Geralyn M. Koeberlein
Publisher:
Cengage Learning
Algebra & Trigonometry with Analytic Geometry
Algebra
ISBN:
9781133382119
Author:
Swokowski
Publisher:
Cengage
Linear Algebra: A Modern Introduction
Algebra
ISBN:
9781285463247
Author:
David Poole
Publisher:
Cengage Learning
Elementary Geometry for College Students
Geometry
ISBN:
9781285195698
Author:
Daniel C. Alexander, Geralyn M. Koeberlein
Publisher:
Cengage Learning
Algebra & Trigonometry with Analytic Geometry
Algebra
ISBN:
9781133382119
Author:
Swokowski
Publisher:
Cengage
Elementary Geometry For College Students, 7e
Geometry
ISBN:
9781337614085
Author:
Alexander, Daniel C.; Koeberlein, Geralyn M.
Publisher:
Cengage,
Elementary Linear Algebra (MindTap Course List)
Algebra
ISBN:
9781305658004
Author:
Ron Larson
Publisher:
Cengage Learning