Crthuky

What is a ^ b and (a & b) << 1?

Is honey really a supersaturated solution? Does heating to un-crystalize redissolve it or melt it?

As a new Ubuntu desktop 18.04 LTS user, do I need to use ufw for a firewall or is iptables sufficient?

Is it good practice to use Linear Least-Squares with SMA?

Meme-controlled people

Can I use USB data pins as power source

Why does overlay work only on the first tcolorbox?

This word with a lot of past tenses

How to deal with taxi scam when on vacation?

Describing a chess game in a novel

Are ETF trackers fundamentally better than individual stocks?

Is "upgrade" the right word to use in this context?

Aluminum electrolytic or ceramic capacitors for linear regulator input and output?

Is there a hypothetical scenario that would make Earth uninhabitable for humans, but not for (the majority of) other animals?

Why do passenger jet manufacturers design their planes with stall prevention systems?

Is there a symmetric-key algorithm which we can use for creating a signature?

Is it true that good novels will automatically sell themselves on Amazon (and so on) and there is no need for one to waste time promoting?

What is the purpose or proof behind chain rule?

Why does energy conservation give me the wrong answer in this inelastic collision problem?

Why do tuner card drivers fail to build after kernel update to 4.4.0-143-generic?

World War I as a war of liberals against authoritarians?

Recruiter wants very extensive technical details about all of my previous work

Is it normal that my co-workers at a fitness company criticize my food choices?

How difficult is it to simply disable/disengage the MCAS on Boeing 737 Max 8 & 9 Aircraft?

Unable to evaluate Eigenvalues and Eigenvectors for a matrix (2)

Unable to evaluate Eigenvalues and Eigenvectors for a matrixProblem with Eigenvectors when given a matrix containing approximate numbers and symbolsIs Eigensystem::eivin message a bug?Seemingly wrong eigenvectors for numerical matrix whose elements differ in scale by orders of magnitudeInverse of a 3x3 MatrixConvert, using the Pauli matrices, an $n times m$ matrix of quaternions into a $2 n times 2 m$ matrix with complex entries, and vice versaEigenvalues and eigenvectors of tensorsNo eigenvectors coming for a very simple* matrixUnable to evaluate Eigenvalues and Eigenvectors for a matrixBasis for unstable manifold of a matrixDensity map for complex and imaginary parts of eigenvalues on one graph

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

share|improve this question

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 of m that I fixed – but check if the form is the desired one.
  $endgroup$
  – corey979
  4 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
  4 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
  4 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
  4 hours ago
  

  
  
  
  

add a comment | 

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

share|improve this question

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 of m that I fixed – but check if the form is the desired one.
  $endgroup$
  – corey979
  4 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
  4 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
  4 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
  4 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

share|improve this question

edited 4 hours ago

corey979

20.9k64282

asked 4 hours ago

kowalskikowalski

1559

$endgroup$

m = -γ/2, -I*g1, -I*Exp[-I*α]*g3, -I*g1, -(κ1)/2, -I*g2, -I*Exp[I*α]*g3, -I*g2, -(κ2)/2]

Eigenvectors[m, Cubics->True]

 Eigenvectors[m, Cubics->True]

...Eigenvectors: Unable to find all eigenvectors

Simplify[Eigenvectors[mchiral /. Complex[0, -1] -> mi, Cubics -> True] /. mi -> -I];

...Eigenvectors: Unable to find all eigenvectors

share|improve this question

edited 4 hours ago

corey979

20.9k64282

asked 4 hours ago

kowalskikowalski

1559

share|improve this question

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

asked 4 hours ago

kowalskikowalski

1559

asked 4 hours ago

kowalskikowalski

1559

1 Answer
1

active

oldest

votes

2

$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

share|improve this answer

edited 4 hours ago

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
  4 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @kowalski U produced by the code above contains the eigenvectors.
  $endgroup$
  – Henrik Schumacher
  4 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
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');

var $window = $(window),
onScroll = function(e)
var $elem = $('.new-login-left'),
docViewTop = $window.scrollTop(),
docViewBottom = docViewTop + $window.height(),
elemTop = $elem.offset().top,
elemBottom = elemTop + $elem.height();
if ((docViewTop elemBottom))
StackExchange.using('gps', function() StackExchange.gps.track('embedded_signup_form.view', location: 'question_page' ); );
$window.unbind('scroll', onScroll);

;
$window.on('scroll', onScroll);

);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

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

Name

Email

Required, but never shown

2

$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

share|improve this answer

edited 4 hours ago

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
  4 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @kowalski U produced by the code above contains the eigenvectors.
  $endgroup$
  – Henrik Schumacher
  4 hours ago
  

  
  
  
  

add a comment | 

2

$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

share|improve this answer

edited 4 hours ago

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
  4 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @kowalski U produced by the code above contains the eigenvectors.
  $endgroup$
  – Henrik Schumacher
  4 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

share|improve this answer

edited 4 hours ago

answered 4 hours ago

Henrik SchumacherHenrik Schumacher

56.7k577157

$endgroup$

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[λ]]

share|improve this answer

edited 4 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