Unable to evaluate Eigenvalues and Eigenvectors for a matrix (2)
$begingroup$
I have posted a similar question last year pertaining to this issue. Here's a link to my post together with the solution given: Unable to evaluate Eigenvalues and Eigenvectors for a matrix
I have tried the methods in my previous posts but to no avail. Here's the problem: I have the following 3x3 matrix
m = {{-γ/2, -I*g1, -I*Exp[-I*α]*g3}, {-I*g1, -(κ1)/2, -I*g2}, {-I*Exp[I*α]*g3, -I*g2, -(κ2)/2]}}
where I represents the complex identity Sqrt[-1]. I wish to find the eigenvectors for the matrix for two different alpha values. For α = π/2, simply doing (after manually replacing α with π/2)
Eigenvectors[m, Cubics->True]
Returns the appropriate (albeit long) eigenvectors. Now however, if I change my α to α = π and run
Eigenvectors[m, Cubics->True]
I am returned with
...Eigenvectors: Unable to find all eigenvectors
Which is the similar issue encountered in the link that I provided above a while ago. I proceed to perform the same fix detailed in that question. Namely
Simplify[Eigenvectors[mchiral /. Complex[0, -1] -> mi, Cubics -> True] /. mi -> -I];
and I am still returned with the same error. Namely
...Eigenvectors: Unable to find all eigenvectors
What is the problem here?
matrix eigenvalues
edited 4 hours ago
corey979
20.9k64282
asked 4 hours ago
kowalskikowalski
1559
$endgroup$
add a comment |
$begingroup$
I have posted a similar question last year pertaining to this issue. Here's a link to my post together with the solution given: Unable to evaluate Eigenvalues and Eigenvectors for a matrix
I have tried the methods in my previous posts but to no avail. Here's the problem: I have the following 3x3 matrix
m = {{-γ/2, -I*g1, -I*Exp[-I*α]*g3}, {-I*g1, -(κ1)/2, -I*g2}, {-I*Exp[I*α]*g3, -I*g2, -(κ2)/2]}}
where I represents the complex identity Sqrt[-1]. I wish to find the eigenvectors for the matrix for two different alpha values. For α = π/2, simply doing (after manually replacing α with π/2)
Eigenvectors[m, Cubics->True]
Returns the appropriate (albeit long) eigenvectors. Now however, if I change my α to α = π and run
Eigenvectors[m, Cubics->True]
I am returned with
...Eigenvectors: Unable to find all eigenvectors
Which is the similar issue encountered in the link that I provided above a while ago. I proceed to perform the same fix detailed in that question. Namely
Simplify[Eigenvectors[mchiral /. Complex[0, -1] -> mi, Cubics -> True] /. mi -> -I];
and I am still returned with the same error. Namely
...Eigenvectors: Unable to find all eigenvectors
What is the problem here?
matrix eigenvalues
edited 4 hours ago
corey979
20.9k64282
asked 4 hours ago
kowalskikowalski
1559
$endgroup$
$begingroup$
There is no problem: imgur.com/a/ALdYCou Except that you have some brackets misplaced in the definition ofmthat I fixed – but check if the form is the desired one.
$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Interesting. I am on version 11.3 for macOS and I do get the error message (even after fixing the brackets).
$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
How many times have I advocated for providing the$Versionone is using... I'm on 10.4 and there is no problem, as showed. Indeed, there is an error in 11.3. No idea what version the OP is using. Unless he clarifies there is virtually no problem.
$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Apologies but I'm using version 11.3. The error still persists even after the bracket fix. I don't know what the problem is
$endgroup$
– kowalski
3 hours ago
add a comment |
$begingroup$
I have posted a similar question last year pertaining to this issue. Here's a link to my post together with the solution given: Unable to evaluate Eigenvalues and Eigenvectors for a matrix
I have tried the methods in my previous posts but to no avail. Here's the problem: I have the following 3x3 matrix
m = {{-γ/2, -I*g1, -I*Exp[-I*α]*g3}, {-I*g1, -(κ1)/2, -I*g2}, {-I*Exp[I*α]*g3, -I*g2, -(κ2)/2]}}
where I represents the complex identity Sqrt[-1]. I wish to find the eigenvectors for the matrix for two different alpha values. For α = π/2, simply doing (after manually replacing α with π/2)
Eigenvectors[m, Cubics->True]
Returns the appropriate (albeit long) eigenvectors. Now however, if I change my α to α = π and run
Eigenvectors[m, Cubics->True]
I am returned with
...Eigenvectors: Unable to find all eigenvectors
Which is the similar issue encountered in the link that I provided above a while ago. I proceed to perform the same fix detailed in that question. Namely
Simplify[Eigenvectors[mchiral /. Complex[0, -1] -> mi, Cubics -> True] /. mi -> -I];
and I am still returned with the same error. Namely
...Eigenvectors: Unable to find all eigenvectors
What is the problem here?
matrix eigenvalues
edited 4 hours ago
corey979
20.9k64282
asked 4 hours ago
kowalskikowalski
1559
$endgroup$
I have posted a similar question last year pertaining to this issue. Here's a link to my post together with the solution given: Unable to evaluate Eigenvalues and Eigenvectors for a matrix
I have tried the methods in my previous posts but to no avail. Here's the problem: I have the following 3x3 matrix
m = {{-γ/2, -I*g1, -I*Exp[-I*α]*g3}, {-I*g1, -(κ1)/2, -I*g2}, {-I*Exp[I*α]*g3, -I*g2, -(κ2)/2]}}
where I represents the complex identity Sqrt[-1]. I wish to find the eigenvectors for the matrix for two different alpha values. For α = π/2, simply doing (after manually replacing α with π/2)
Eigenvectors[m, Cubics->True]
Returns the appropriate (albeit long) eigenvectors. Now however, if I change my α to α = π and run
Eigenvectors[m, Cubics->True]
I am returned with
...Eigenvectors: Unable to find all eigenvectors
Which is the similar issue encountered in the link that I provided above a while ago. I proceed to perform the same fix detailed in that question. Namely
Simplify[Eigenvectors[mchiral /. Complex[0, -1] -> mi, Cubics -> True] /. mi -> -I];
and I am still returned with the same error. Namely
...Eigenvectors: Unable to find all eigenvectors
What is the problem here?
matrix eigenvalues
matrix eigenvalues
edited 4 hours ago
corey979
20.9k64282
asked 4 hours ago
kowalskikowalski
1559
edited 4 hours ago
corey979
20.9k64282
asked 4 hours ago
kowalskikowalski
1559
edited 4 hours ago
corey979
20.9k64282
edited 4 hours ago
corey979
20.9k64282
edited 4 hours ago
corey979
20.9k64282
20.9k64282
asked 4 hours ago
kowalskikowalski
1559
asked 4 hours ago
kowalskikowalski
1559
asked 4 hours ago
kowalskikowalski
1559
1559
$begingroup$
There is no problem: imgur.com/a/ALdYCou Except that you have some brackets misplaced in the definition ofmthat I fixed – but check if the form is the desired one.
$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Interesting. I am on version 11.3 for macOS and I do get the error message (even after fixing the brackets).
$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
How many times have I advocated for providing the$Versionone is using... I'm on 10.4 and there is no problem, as showed. Indeed, there is an error in 11.3. No idea what version the OP is using. Unless he clarifies there is virtually no problem.
$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Apologies but I'm using version 11.3. The error still persists even after the bracket fix. I don't know what the problem is
$endgroup$
– kowalski
3 hours ago
add a comment |
$begingroup$
There is no problem: imgur.com/a/ALdYCou Except that you have some brackets misplaced in the definition ofmthat I fixed – but check if the form is the desired one.
$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Interesting. I am on version 11.3 for macOS and I do get the error message (even after fixing the brackets).
$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
How many times have I advocated for providing the$Versionone is using... I'm on 10.4 and there is no problem, as showed. Indeed, there is an error in 11.3. No idea what version the OP is using. Unless he clarifies there is virtually no problem.
$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Apologies but I'm using version 11.3. The error still persists even after the bracket fix. I don't know what the problem is
$endgroup$
– kowalski
3 hours ago
$begingroup$
There is no problem: imgur.com/a/ALdYCou Except that you have some brackets misplaced in the definition of
m that I fixed – but check if the form is the desired one.$endgroup$
– corey979
3 hours ago
$begingroup$
There is no problem: imgur.com/a/ALdYCou Except that you have some brackets misplaced in the definition of
m that I fixed – but check if the form is the desired one.$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Interesting. I am on version 11.3 for macOS and I do get the error message (even after fixing the brackets).
$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
@corey979 Interesting. I am on version 11.3 for macOS and I do get the error message (even after fixing the brackets).
$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
How many times have I advocated for providing the
$Version one is using... I'm on 10.4 and there is no problem, as showed. Indeed, there is an error in 11.3. No idea what version the OP is using. Unless he clarifies there is virtually no problem.$endgroup$
– corey979
3 hours ago
$begingroup$
How many times have I advocated for providing the
$Version one is using... I'm on 10.4 and there is no problem, as showed. Indeed, there is an error in 11.3. No idea what version the OP is using. Unless he clarifies there is virtually no problem.$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Apologies but I'm using version 11.3. The error still persists even after the bracket fix. I don't know what the problem is
$endgroup$
– kowalski
3 hours ago
$begingroup$
@corey979 Apologies but I'm using version 11.3. The error still persists even after the bracket fix. I don't know what the problem is
$endgroup$
– kowalski
3 hours ago
add a comment |
1 Answer
1
active
oldest
votes
$begingroup$
I have o clue why this did not work. However, this old-fashioned method seems to work
m = {
{-γ/2, -I g1, -I Exp[-I α] g3},
{-I g1, -(κ1)/2, -I g2},
{-I Exp[I α] g3, -I g2, -(κ2)/2}
};
a = m /. α -> π;
λ = Eigenvalues[a] // ToRadicals;
U = Flatten[NullSpace[a - # IdentityMatrix[3]] & /@ λ, 1];
Simplify[a.Transpose[U] - Transpose[U].DiagonalMatrix[λ]]
{{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
$endgroup$
$begingroup$
The eigenvalues seems to be fine. It's long but it prints. The eigenvectors on the other hand, are still unobtainable. I'm starting to think if this is alpha dependent since it works for pi/2 but not pi. But that's just silly and shouldn't happen
$endgroup$
– kowalski
3 hours ago
$begingroup$
@kowalskiUproduced by the code above contains the eigenvectors.
$endgroup$
– Henrik Schumacher
3 hours ago
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "387"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f193391%2funable-to-evaluate-eigenvalues-and-eigenvectors-for-a-matrix-2%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
$begingroup$
I have o clue why this did not work. However, this old-fashioned method seems to work
m = {
{-γ/2, -I g1, -I Exp[-I α] g3},
{-I g1, -(κ1)/2, -I g2},
{-I Exp[I α] g3, -I g2, -(κ2)/2}
};
a = m /. α -> π;
λ = Eigenvalues[a] // ToRadicals;
U = Flatten[NullSpace[a - # IdentityMatrix[3]] & /@ λ, 1];
Simplify[a.Transpose[U] - Transpose[U].DiagonalMatrix[λ]]
{{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
$endgroup$
$begingroup$
The eigenvalues seems to be fine. It's long but it prints. The eigenvectors on the other hand, are still unobtainable. I'm starting to think if this is alpha dependent since it works for pi/2 but not pi. But that's just silly and shouldn't happen
$endgroup$
– kowalski
3 hours ago
$begingroup$
@kowalskiUproduced by the code above contains the eigenvectors.
$endgroup$
– Henrik Schumacher
3 hours ago
add a comment |
$begingroup$
I have o clue why this did not work. However, this old-fashioned method seems to work
m = {
{-γ/2, -I g1, -I Exp[-I α] g3},
{-I g1, -(κ1)/2, -I g2},
{-I Exp[I α] g3, -I g2, -(κ2)/2}
};
a = m /. α -> π;
λ = Eigenvalues[a] // ToRadicals;
U = Flatten[NullSpace[a - # IdentityMatrix[3]] & /@ λ, 1];
Simplify[a.Transpose[U] - Transpose[U].DiagonalMatrix[λ]]
{{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
$endgroup$
$begingroup$
The eigenvalues seems to be fine. It's long but it prints. The eigenvectors on the other hand, are still unobtainable. I'm starting to think if this is alpha dependent since it works for pi/2 but not pi. But that's just silly and shouldn't happen
$endgroup$
– kowalski
3 hours ago
$begingroup$
@kowalskiUproduced by the code above contains the eigenvectors.
$endgroup$
– Henrik Schumacher
3 hours ago
add a comment |
$begingroup$
I have o clue why this did not work. However, this old-fashioned method seems to work
m = {
{-γ/2, -I g1, -I Exp[-I α] g3},
{-I g1, -(κ1)/2, -I g2},
{-I Exp[I α] g3, -I g2, -(κ2)/2}
};
a = m /. α -> π;
λ = Eigenvalues[a] // ToRadicals;
U = Flatten[NullSpace[a - # IdentityMatrix[3]] & /@ λ, 1];
Simplify[a.Transpose[U] - Transpose[U].DiagonalMatrix[λ]]
{{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
$endgroup$
I have o clue why this did not work. However, this old-fashioned method seems to work
m = {
{-γ/2, -I g1, -I Exp[-I α] g3},
{-I g1, -(κ1)/2, -I g2},
{-I Exp[I α] g3, -I g2, -(κ2)/2}
};
a = m /. α -> π;
λ = Eigenvalues[a] // ToRadicals;
U = Flatten[NullSpace[a - # IdentityMatrix[3]] & /@ λ, 1];
Simplify[a.Transpose[U] - Transpose[U].DiagonalMatrix[λ]]
{{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
edited 3 hours ago
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
answered 4 hours ago
Henrik SchumacherHenrik Schumacher
56.7k577157
56.7k577157
$begingroup$
The eigenvalues seems to be fine. It's long but it prints. The eigenvectors on the other hand, are still unobtainable. I'm starting to think if this is alpha dependent since it works for pi/2 but not pi. But that's just silly and shouldn't happen
$endgroup$
– kowalski
3 hours ago
$begingroup$
@kowalskiUproduced by the code above contains the eigenvectors.
$endgroup$
– Henrik Schumacher
3 hours ago
add a comment |
$begingroup$
The eigenvalues seems to be fine. It's long but it prints. The eigenvectors on the other hand, are still unobtainable. I'm starting to think if this is alpha dependent since it works for pi/2 but not pi. But that's just silly and shouldn't happen
$endgroup$
– kowalski
3 hours ago
$begingroup$
@kowalskiUproduced by the code above contains the eigenvectors.
$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
The eigenvalues seems to be fine. It's long but it prints. The eigenvectors on the other hand, are still unobtainable. I'm starting to think if this is alpha dependent since it works for pi/2 but not pi. But that's just silly and shouldn't happen
$endgroup$
– kowalski
3 hours ago
$begingroup$
The eigenvalues seems to be fine. It's long but it prints. The eigenvectors on the other hand, are still unobtainable. I'm starting to think if this is alpha dependent since it works for pi/2 but not pi. But that's just silly and shouldn't happen
$endgroup$
– kowalski
3 hours ago
$begingroup$
@kowalski
U produced by the code above contains the eigenvectors.$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
@kowalski
U produced by the code above contains the eigenvectors.$endgroup$
– Henrik Schumacher
3 hours ago
add a comment |
Thanks for contributing an answer to Mathematica Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f193391%2funable-to-evaluate-eigenvalues-and-eigenvectors-for-a-matrix-2%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
$begingroup$
There is no problem: imgur.com/a/ALdYCou Except that you have some brackets misplaced in the definition of
mthat I fixed – but check if the form is the desired one.$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Interesting. I am on version 11.3 for macOS and I do get the error message (even after fixing the brackets).
$endgroup$
– Henrik Schumacher
3 hours ago
$begingroup$
How many times have I advocated for providing the
$Versionone is using... I'm on 10.4 and there is no problem, as showed. Indeed, there is an error in 11.3. No idea what version the OP is using. Unless he clarifies there is virtually no problem.$endgroup$
– corey979
3 hours ago
$begingroup$
@corey979 Apologies but I'm using version 11.3. The error still persists even after the bracket fix. I don't know what the problem is
$endgroup$
– kowalski
3 hours ago