Question: eigenvector equation discrepancies

restart

with(LinearAlgebra)

A := `<,>`(`<|>`(-2*`&omega;__0`^2, `&omega;__0`^2), `<|>`(`&omega;__0`^2, -`&omega;__0`^2))

Matrix(%id = 36893488151921005020)

(1)

AA := subs(`&omega;__0` = 3, A)

Matrix(%id = 36893488151921002244)

(2)

evalf([Eigenvectors(AA)])

[Vector[column](%id = 36893488151920990196), Matrix(%id = 36893488151920990316)]

(3)

NULL

Below is an eigenvector equation for AA for the eigenvalue -3.4376...

 

S := AA+3.43769410*IdentityMatrix(2)

Matrix(%id = 36893488151920987180)

(4)

When I solve the system below, I expect to get the vector `<,>`(.6180, 1)

Vector[column](%id = 36893488151920981284)

(5)

, just like the result of the Eigenvectors command used above.

 

LinearSolve(S, `<,>`(0, 0))

Vector[column](%id = 36893488151920975988)

(6)

Instead I get the zero vector.

 

In addition, when I compute the reduced row echelon form, I expect to have one of the rows be zero since matrix S has rank 1.

 

Rank(S)

1

(7)

ReducedRowEchelonForm(S)

Matrix(%id = 36893488151959135044)

(8)

NULL

Therefore, I guess my question is maybe about calculations and not about Maple. That being said, I did the calculations manually, then I asked chatgpt, which agreed with me.

 

Why are the calculations not coinciding with Maple?

 

I can do the calculations "manually" with Maple:

 

solve(-14.5623059000000*x+9*y = 0)

{x = x, y = 1.618033989*x}

(9)

solve(9*x-5.56230590000000*y = 0)

{x = x, y = 1.618033988*x}

(10)

 

Now, if I try to solve them simultaneously, it seems I have the same issue as with LinearSolve

NULL

solve({9*x-5.56230590000000*y = 0, -14.5623059000000*x+9*y = 0})

{x = 0., y = 0.}

(11)

I tried something else that worked, but the questions above still remain.

 

We start with A again

 

A = Matrix(%id = 36893488151921005020)NULL

In my problem I have the variable

 

`&omega;__a` := (1/2)*`&omega;__0`*(sqrt(5)-1) = (1/2)*omega__0*(5^(1/2)-1)NULL

NULL

and it turns out that -`#msub(mi("omega",fontstyle = "normal"),mi("a"))`^2 is an eigenvalue of A.

In particular, for `&omega;__0` = 3 we get the values I showed previously.

 

For example, here is the eigenvalue we were looking at beforeevalf(subs(`&omega;__0` = 3, -`&omega;__a`^2))

-3.437694099

(12)

Okay so the eigenvector equation without subbing in a value for `&omega;__0` is (A+I*`&omega;__a`^2)*v = 0.

 

Here is what the matrix on the left looks like

NULL

A+`&omega;__a`^2*IdentityMatrix(2)

Matrix(%id = 36893488151964015780)

(13)

 

And here we solve for the eigenvectors

sol := LinearSolve(A+`&omega;__a`^2*IdentityMatrix(2), `<,>`(0, 0))

Vector[column](%id = 36893488151964001572)

(14)

evalf(subs(_t2[1] = 1, sol))

Vector[column](%id = 36893488151964010244)

(15)

Which is the expected result.

 

So, Maple can solve the eigenvector equation in this more general case but not if I specify `&omega;__0` before asking it to solve.NULL

Download evcalc.mw

Please Wait...