Azzera filtri
Azzera filtri

Eigen value eigen vectors in matlab

6 visualizzazioni (ultimi 30 giorni)
Aykut Albayrak
Aykut Albayrak il 21 Gen 2018
Risposto: mercy charles il 19 Feb 2022
I want to solve this question.I'm going to atteched my code but I don't know this code is true.
if true
clc;
clear all;
A=[4 -5;-2 1];
Eig=eig(A)
F=poly(A)
[V,D]=eig(A)
V1=V(:,1)
V11=V1/V(1,1)
end
  1 Commento
Walter Roberson
Walter Roberson il 21 Gen 2018
What are you doing with the Eig and F values you calculate?
I do not see any initial value problem in the question presented: it gives the initial value. Initial value problems require that information about the original condition is missing and needs to be calculated, which is not the case here.

Accedi per commentare.

Risposta accettata

Birdman
Birdman il 21 Gen 2018
Your approach is nice, but not enough of course. First, you need to construct a differential equation for your second question. Third question comes with the second one. The answer for your first question is included in this line:
[V,D]=eig(A)
Full solution may look like as follows:
%%1st
A=[4 -5;-2 1];
[V,D]=eig(A)
%%2nd
syms x1(t) x2(t)
x0=[2.9;2.6];
eq1=diff(x1,t)==A(1,:)*[x1(t);x2(t)];
eq2=diff(x2,t)==A(2,:)*[x1(t);x2(t)];
solx=dsolve([eq1 eq2],[x1(0)==2.9 x2(0)==2.6]);
x1=solx.x1
x2=solx.x2
%%3rd
t=subs(t,0:0.1:1);
x1=subs(x1,t);
x2=subs(x2,t);
plot(t,x1,'o',t,x2,'-*')
  2 Commenti
Aykut Albayrak
Aykut Albayrak il 21 Gen 2018
Thanks for helping.I really appreciate for helping
Birdman
Birdman il 21 Gen 2018
You are welcome. You can officially thank me by accepting my answer.

Accedi per commentare.

Più risposte (2)

johnson wul
johnson wul il 23 Lug 2019
Nice that. is this formular for eigenvector:
plot(t,x1,'o',t,x2,..,t,x10, '-*')
have been used to for square matrix 10*10 for example if some one want to plot differents eigenvectors?
And how to plot in 3 dimensions the previous eigenvectors? just a step

mercy charles
mercy charles il 19 Feb 2022
%%1st
A=[4 -5;-2 1];
[V,D]=eig(A)
%%2nd
syms x1(t) x2(t)
x0=[2.9;2.6];
eq1=diff(x1,t)==A(1,:)*[x1(t);x2(t)];
eq2=diff(x2,t)==A(2,:)*[x1(t);x2(t)];
solx=dsolve([eq1 eq2],[x1(0)==2.9 x2(0)==2.6]);
x1=solx.x1
x2=solx.x2
%%3rd
t=subs(t,0:0.1:1);
x1=subs(x1,t);
x2=subs(x2,t);
plot(t,x1,'o',t,x2,'-*')

Community Treasure Hunt

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

Start Hunting!

Translated by