12/15/2023 0 Comments Eigenvalues matlab dont match mathcad![]() TimWescott at 17:12 I have already did this. Then instead of going by some dubious authority on the Internet, you're actually giving a mathematical proof. For A not exactly symmetric it probably takes route (1) and gives you a basis of each subspace, but not necessarily with orthogonal vectors. 1 I would suggest taking the FFT of sinusoids of various frequencies, and showing that the resulting FFT has a single bin with the same amplitude as the input sinusoid. Matlab is probably taking route (2) (thus forcing V_a to be orthogonal) only if A is symmetric. (2) However, since every subspace has an orthonormal basis,you can find orthonormal bases for each eigenspace, so you can find an orthonormal basis of eigenvectors. ![]() Eigenvectors corresponding to the same eigenvalue need not be orthogonal to each other. This is my bet as to what's happening: as noted by Eigenvectors corresponding to distinct eigenvalues of a symmetric matrix must be orthogonal to each other. Still, it's surprising that so wildly different results are obtained for V_A when A is symmetric and when A is nearly symmetric. If you force A to be exactly symmetric you'll get an orthogonal V_A, up to numerical errrors on the order of eps: > A = (A+A.')/2 The differences are on the order of eps, as expected from numerical errors. But your input matrix A is not exactly symmetric. The eigenvectors of a real symmetric matrix are orthogonal. This seems to be a numerical precision issue. I was under the impression that if my matrix was real and symmetric then Matlab would return an orthogonal matrix of eigenvectors, so what is the issue here? V_A*V_A' is not equalling the identity matrix (taking into account rounding errors). The values of that satisfy the equation are the eigenvalues. The eigenvalue problem is to determine the solution to the equation Av v, where A is an n-by-n matrix, v is a column vector of length n, and is a scalar. ![]() ![]() I am finding that my matrix of eigenvectors V_A is not orthogonal i.e. V,D,W eig(A) also returns full matrix W whose columns are the corresponding left eigenvectors, so that WA DW. % less than the given number of repeats if we ever input the same % Put one random value into another (note this sometimes will result in % Now generate a vector of eigenvalues with the given number of repeats % First generate a random matrix of eignevectors that is orthonormal The variable-precision counterparts are E eig (vpa (A)) and V,E eig (vpa (A)). When I try to find the eigen-decomposition of a matrix in Matlab that has a repeated eigenvalue but is NOT defective, it is not returning an orthonormal matrix of eignevectors. Answers Trial Software Product Updates Eigenvalues The symbolic eigenvalues of a square matrix A or the symbolic eigenvalues and eigenvectors of A are computed, respectively, using the commands E eig (A) and V,E eig (A). ![]()
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |