non linear eigen value problem

7 次查看（过去 30 天）
rakesh kumar
rakesh kumar 2023-3-19
回答： Christine Tobler 2023-3-27
Ran in:
I want to find the eigen values of
k0=magic(2);
k1=1e-3*k0;
c1=[1 0;0 1]
c1 = 2×2
1 0 0 1
c2=[0 1;1 0]
c2 = 2×2
0 1 1 0
B=kron(c1,k0)+kron(c2,k1)
B = 4×4
1.0000 3.0000 0.0010 0.0030 4.0000 2.0000 0.0040 0.0020 0.0010 0.0030 1.0000 3.0000 0.0040 0.0020 4.0000 2.0000
syms lambda1a lambda1b lambda2a lambda2b
x = sym('x',[size(B,2),1])
x = 
l1=[lambda1a;lambda1b ]
l1 = 
l2=[lambda2a;lambda2b ]
l2 = 
% lambda1 is 2 by 1 matrix and lambda2 is 2 by 1 matrix
% I want to solve the problem
B*x
ans = 
(kron(diag(l1),c1)+kron(diag(l2),c2))*x
ans = 
%I ant to find l1 and l2
  3 个评论
显示 1更早的评论
rakesh kumar
rakesh kumar 2023-3-19
(kron(diag(l1),c1)+kron(diag(l2),c2))
Walter Roberson
Walter Roberson 2023-3-19
l1 and l2 are each 2 x 1 matrices. You cannot take the eigenvalues of non-square matrices.
If you take eig(diag(l1)) and eig(diag(l2)) so that you are making them into 2 x 2 diagonal matrices, then the eigenvalues are just the contents of l1 and l2

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Christine Tobler
Christine Tobler 2023-3-27
This isn't the standard definition of a nonlinear eigenvalue problem, where you would have only one scalar lambda.
Am I understanding correctly that you are looking for four scalar values lambda... for which there exists an non-zero vector x which satisfies
B*x == (kron(diag(l1),c1)+kron(diag(l2),c2))*x
In that case, I would take this to the symbolic toolbox, by trying to solve for values of lambda... for which
det(B == (kron(diag(l1),c1)+kron(diag(l2),c2))) == 0
(using DET is fine for symbolic calculations, we only run into trouble using it when subjected to round-off error).

类别

在 Help Center 和 File Exchange 中查找有关 Linear Algebra 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by