MATLAB Answers

How to find eigenvalues for a system of lenearized ordinary differential equations?

조회 수: 91(최근 30일)
I have a system of linearized ODEs with corresponding boundary conditions.
%----------------------------system of ODEs--------------------------------------%
%---------------------------boundary conditions----------------------------------%
y(1)=y(2)=y(4)=0 at eta=0
y(2)=y(4)=0 at eta=0;
here Pr phi Ra Da Fr A1 A2 fd1 fd fdd thd are known quantities and 'e' is unknown.
I need to solve the system to find out the eigenvalues (e).
Thanks in advance.

채택된 답변

Pavel Osipov
Pavel Osipov 5 Oct 2019
Tanya, hi.
write so:
dyi/dt =...y1 (t)+...y2 (t)+...+y5(t);
let x (t)=[y1;y2;...;y5]; ->
((V/ve) x=Ah; A - matrix coeff. Your system. Let's formally denote d/dt=p
px-Ax=0; - > (p*E-A) x=0; since x is not 0, then
det(p*E-A)=0. This is the equation for the eigenvalues of p.
  댓글 수: 1
Tanya Sharma
Tanya Sharma 9 Oct 2019
The unkown 'e' is already in the equations y'(3) and y'(4). I want to find all the possible eigenvalues 'e' for this problem. Can you explain again?
Thanks in advance.

댓글을 달려면 로그인하십시오.

추가 답변(2개)

Pavel Osipov
Pavel Osipov 9 Oct 2019
det(p*E-A)=0. This is the equation for the eigenvalues of p. - -> The unkown "p" is solution det(p*E-A)=0. det - is the determinant with dimensions 5x5.
px-Ax=0 ->Ax=px, p is eigenvalues of A MATLAB command [V,D] = eig(A) returns diagonal matrix D of eigenvalues and matrix V whose columns are the corresponding right eigenvectors, so that A*V = V*D. (from MATLAB help).
eigenvalues p is 5x1 vector = liagonal elements D. eigenvectors of A see at columns V.
  댓글 수: 1
Tanya Sharma
Tanya Sharma 21 Oct 2019
Thanks Pavel!
But this will give me only five eigenvalues. As I am solving the eigenvalue of a differential equation and it can have many eigenvalues.
I am attaching the linearized eigenvalue problem. Is there a way I can find the unknown eigenvalues?

댓글을 달려면 로그인하십시오.

Community Treasure Hunt

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

Start Hunting!

Translated by