Advanced Search

Question: How to solve this equation and get time history plot? There is something wrong with my code. Could someone help me check it? Thx.

Question:How to solve this equation and get time history plot? There is something wrong with my code. Could someone help me check it? Thx.

ssllys 35 Maple 2018
differential-equation dsolve numeric
May 30 2018
1


 

> 

 

> 

C_t := -A(t)*omega[0]*(sigma-(diff(C(t), t)))-(1/4)*A(t)^5*beta[2]*omega[0]^2-(1/4)*A(t)^3*beta[1]*omega[0]^2-(1/2)*F[0]*cos(C(t))+5*alpha[2]*A(t)^5*(1/16)+3*alpha[1]*A(t)^3*(1/8) = 0;

-A(t)*omega[0]*(sigma-(diff(C(t), t)))-0.3297500000e-1*A(t)^5*omega[0]^2-0.8345000000e-1*A(t)^3*omega[0]^2-cos(C(t))+0.4059375000e-1*A(t)^5+.1249125000*A(t)^3 = 0

 (1)

 

> 

A_t := (1/2)*mu*A(t)*omega[0]+(diff(A(t), t))*omega[0]-(1/2)*F[0]*sin(C(t)) = 0;

0.5000000000e-2*A(t)*omega[0]+(diff(A(t), t))*omega[0]-sin(C(t)) = 0

 (2)
> 

 

> 

F[0] :=2: w0 := 1: alpha[1]:=0.3331:alpha[2]:=0.1299:beta[1]:=0.3338:beta[2]:=0.1319:mu:=0.01:sigma=0：

sigma = `0：`

 (3)
> 

 

> 

C_t;A_t;

-A(t)*omega[0]*(sigma-(diff(C(t), t)))-0.3297500000e-1*A(t)^5*omega[0]^2-0.8345000000e-1*A(t)^3*omega[0]^2-cos(C(t))+0.4059375000e-1*A(t)^5+.1249125000*A(t)^3 = 0

  

0.5000000000e-2*A(t)*omega[0]+(diff(A(t), t))*omega[0]-sin(C(t)) = 0

 (4)
> 

ICS:=[C(0)=0,D(C)=0,A(0)=1,D(A)=0]:
> 

sol:=dsolve({C_t,A_t,ICS},[C(t2),D(C),A(t2),D(A)],type=numeric):

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

  
> 

pRange:=0..20:
> 

 Atime:=odeplot( sol, [t, A],t=pRange, numpoints=10000 ):
> 

Atime;

odeplot(sol, [t, A], t = 0 .. 20, numpoints = 10000)

 (5)
> 

 


 

Download phase_plot.mw
Please Wait...