Advanced Search

Question: solving homogeneous algebraic equations for non-trivial solution

Question:solving homogeneous algebraic equations for non-trivial solution

sunit 105 Maple
solve
September 03 2016
0

Hi,

I am trying to solve a set of homogeneous equations for the non-trivial solutions. Mathematically it is possible to get it. But is there any way to get it in Maple. Please find the attached maple sheet for the question. Please help me regarding this.

Regards

Sunit

> 

restart
> 

with(plots):
> 

with(LinearAlgebra):
> 

eq[1] := diff(x[1](t), t)-x[2](t)

diff(x[1](t), t)-x[2](t)

 (1)
> 

eq[2] := diff(x[2](t), t)+2*zeta*beta*x[2](t)+beta^2*x[1](t)+n*psi*(-v*(phi[1](t)-phi[1](t-2*Pi/(n*omega0)))+x[1](t)-x[1](t-2*Pi/(n*omega0)))

diff(x[2](t), t)+2*zeta*beta*x[2](t)+beta^2*x[1](t)+n*psi*(-v*(phi[1](t)-phi[1](t-2*Pi/(n*omega0)))+x[1](t)-x[1](t-2*Pi/(n*omega0)))

 (2)
> 

eq[3] := diff(phi[1](t), t)-phi[2](t)

diff(phi[1](t), t)-phi[2](t)

 (3)
> 

eq[4] := diff(phi[2](t), t)+2*kappa*phi[2](t)+phi[1](t)+n*(-v*(phi[1](t)-phi[1](t-2*Pi/(n*omega0)))+x[1](t)-x[1](t-2*Pi/(n*omega0)))

diff(phi[2](t), t)+2*kappa*phi[2](t)+phi[1](t)+n*(-v*(phi[1](t)-phi[1](t-2*Pi/(n*omega0)))+x[1](t)-x[1](t-2*Pi/(n*omega0)))

 (4)
> 

for k to 4 do eqn[k] := simplify(coeff(map(expand, eval(eq[k], [x[1] = (proc (t) options operator, arrow; x[1]*exp(lambda*t) end proc), x[2] = (proc (t) options operator, arrow; x[2]*exp(lambda*t) end proc), phi[1] = (proc (t) options operator, arrow; phi[1]*exp(lambda*t) end proc), phi[2] = (proc (t) options operator, arrow; phi[2]*exp(lambda*t) end proc)])), exp(lambda*t))) end do

x[1]*lambda-x[2]

  

x[2]*lambda+2*zeta*beta*x[2]+beta^2*x[1]-n*psi*v*phi[1]+n*psi*v*phi[1]*exp(-2*lambda*Pi/(n*omega0))+n*psi*x[1]-n*psi*x[1]*exp(-2*lambda*Pi/(n*omega0))

  

phi[1]*lambda-phi[2]

  

phi[2]*lambda+2*kappa*phi[2]+phi[1]-n*v*phi[1]+n*v*phi[1]*exp(-2*lambda*Pi/(n*omega0))+n*x[1]-n*x[1]*exp(-2*lambda*Pi/(n*omega0))

 (5)
> 

A, b := GenerateMatrix([seq(eqn[k], k = 1 .. 4)], [x[1], x[2], phi[1], phi[2]])

A, b := Matrix(4, 4, {(1, 1) = lambda, (1, 2) = -1, (1, 3) = 0, (1, 4) = 0, (2, 1) = beta^2+n*psi-n*psi*exp(-2*lambda*Pi/(n*omega0)), (2, 2) = 2*Zeta*beta+lambda, (2, 3) = n*psi*v*exp(-2*lambda*Pi/(n*omega0))-n*psi*v, (2, 4) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = lambda, (3, 4) = -1, (4, 1) = n-n*exp(-2*lambda*Pi/(n*omega0)), (4, 2) = 0, (4, 3) = -n*v+1+n*v*exp(-2*lambda*Pi/(n*omega0)), (4, 4) = 2*kappa+lambda}), Vector(4, {(1) = 0, (2) = 0, (3) = 0, (4) = 0})

 (6)
> 

> 

right_vector := Matrix(4, 1, [r[1], r[2], r[3], r[4]])

right_vector := Matrix(4, 1, {(1, 1) = r[1], (2, 1) = r[2], (3, 1) = r[3], (4, 1) = r[4]})

 (7)
> 

junk := MatrixVectorMultiply(subs(lambda = I*omega, A), right_vector)

junk := Matrix(4, 1, {(1, 1) = I*omega*r[1]-r[2], (2, 1) = (beta^2+n*psi-n*psi*exp(-(2*I)*omega*Pi/(n*omega0)))*r[1]+(2*Zeta*beta+I*omega)*r[2]+(n*psi*v*exp(-(2*I)*omega*Pi/(n*omega0))-n*psi*v)*r[3], (3, 1) = I*omega*r[3]-r[4], (4, 1) = (n-n*exp(-(2*I)*omega*Pi/(n*omega0)))*r[1]+(-n*v+1+n*v*exp(-(2*I)*omega*Pi/(n*omega0)))*r[3]+(2*kappa+I*omega)*r[4]})

 (8)
> 

junk(1)

I*omega*r[1]-r[2]

 (9)
> 

for k to 4 do eqnn[k] := junk(k) end do

I*omega*r[1]-r[2]

  

(beta^2+n*psi-n*psi*exp(-(2*I)*omega*Pi/(n*omega0)))*r[1]+(2*zeta*beta+I*omega)*r[2]+(n*psi*v*exp(-(2*I)*omega*Pi/(n*omega0))-n*psi*v)*r[3]

  

I*omega*r[3]-r[4]

  

(n-n*exp(-(2*I)*omega*Pi/(n*omega0)))*r[1]+(1-n*v+n*v*exp(-(2*I)*omega*Pi/(n*omega0)))*r[3]+(2*kappa+I*omega)*r[4]

 (10)
> 

solve({seq(eqnn[k], k = 1 .. 4)}, {seq(r[k], k = 1 .. 4)})

{r[1] = 0, r[2] = 0, r[3] = 0, r[4] = 0}

 (11)
> 

``
> 

``
> 

``

 

Download question4.mw
Please Wait...