Items tagged with nonlinear nonlinear Tagged Items Feed

Hello,

I would like to plot gait diagrams (the lines you can see on the picture belowà from the solutions obtained with a NL oscillator (composed with 8 coupled odes). Here the result that I would like to obtain.

Initial plot:

 

Desired plot

 

 

I would like to obtain 4 lines corresponding to the 4 elliptic trajectories obtained with the NL oscillator. The four lines should be done like this. When the trajectory is above 0, the line should be colored in green. When the trajectory is below 0, the line should be colored in black. 

May you help me to define this kind of graph called gait diagrams from the solution of the NL oscillator ?

Here you can find my maple code:

K:=Matrix([<0, -1, 1, -1>,<-1, 0, -1, 1>,<-1, 1, 0,-1>,<1, -1, -1,0>]);

for i to 4
do
r[i]:=sqrt((u[i](t))^2+(v[i](t))^2):
omega[i]:=omega[sw]/(1+exp(b*v[i](t)))+omega[st]/(1+exp(-b*v[i](t))):
Equ[i]:=diff(u[i](t),t)=Au*(1-r[i]^2)*u[i](t)-omega[i]*v[i](t):
Eqv[i]:=diff(v[i](t),t)=Av*(1-r[i]^2)*v[i](t)+omega[i]*u[i](t)+MatrixVectorMultiply(K,<seq(v[i](t),i=1..4)>)[i]:
EqSys[i]:=[Equ[i],Eqv[i]]:
end do:

paramsCycle:=omega[st]=4*2*Pi,omega[sw]=2*Pi,Au=5,Av=50,b=100;
params:=paramsCycle;

Differential system 
sys:=map(op,eval([seq(EqSys[i],i=1..4)],[params]));
ic:=[u[1](0)=0, v[1](0)=0,u[2](0)=0, v[2](0)=-0.1,u[3](0)=0, v[3](0)=0.1,u[4](0)=0, v[4](0)=0.1];
Résolution1
res:=dsolve([sys[],ic[]],numeric):
Initial boundaries
tcalc:=4;
ic2:=[seq(u[i](0)=eval(u[i](t), res(tcalc)),i=1..4),seq(v[i](0)=eval(v[i](t), res(tcalc)),i=1..4)];
Résolution2
res:=dsolve([sys[],ic2[]],numeric):

tmax:= 40:
numpts:=100*tmax:
plots:-odeplot(res,[t,v[1](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[t,v[2](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[t,v[3](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[t,v[4](t)],0..tmax,thickness=2, view=[0..5, -1.5..1.5],numpoints = numpts);
plots:-odeplot(res,[seq([t,v[i](t)+i*5],i=1..4)],0..tmax,thickness=2,view=[0..5,0..25], numpoints = numpts);

Thanks a lot for your help

I want to obtain polynomial representation of my data:

with(Statistics):

X:=Vector([huge data package-X],datatype=float):

Y:=Vector([huge data package-Y],datatype=float):

NonlinearFit(c*t^2+b*t+a, X, Y, t)

but I can't see any result. What's the problem?

My worksheet: He_p=f(t).mw

I am trying to solve 4 nonlinear equations for four variables using fsolve  and the output that i am getting is basically the same equations repeated after some time.  I even tried reducing one of the equations using assumptions from my side but it results in same behaviour..  Quite new to maple, would like some advice as to this behaviour. Thanks

 Here's the file

fsolve_1.mw

 

PS- using do loop is part of the solving so i cannot remove that

Hi all,

I need to solve det[M]=0 for omega.

M is:

M := Matrix(8, 8, {(1, 1) = BesselJ(n, tp*a), (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = -BesselJ(n, tg*a), (1, 6) = -BesselY(n, tg*a), (1, 7) = 0, (1, 8) = 0, (2, 1) = k*n*BesselJ(n, tp*a)/(tp^2*a), (2, 2) = I*k*n*mu0*omega*(diff(BesselJ(n, tp*a), a))/(tp^2*a), (2, 3) = 0, (2, 4) = 0, (2, 5) = -k*n*BesselJ(n, tg*a)/(tg^2*a), (2, 6) = -k*n*BesselY(n, tg*a)/(tg^2*a), (2, 7) = -I*`&mu;g`*omega*(diff(BesselJ(n, tg*a), a))/tg^2, (2, 8) = -I*`&mu;g`*omega*(diff(BesselY(n, tg*a), a))/tg^2, (3, 1) = 0, (3, 2) = BesselJ(n, tp*a), (3, 3) = 0, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (3, 7) = -BesselJ(n, tg*a), (3, 8) = -BesselY(n, tg*a), (4, 1) = -I*omega*`&varepsilon;p`*(diff(BesselJ(n, tp*a), a))/tp^2, (4, 2) = k*n*BesselJ(n, tp*a)/(tp^2*a), (4, 3) = 0, (4, 4) = 0, (4, 5) = I*`&varepsilon;g`*omega*(diff(BesselJ(n, tg*a), a))/tg^2, (4, 6) = I*`&varepsilon;g`*omega*(diff(BesselY(n, tg*a), a))/tg^2, (4, 7) = -k*n*BesselJ(n, tg*a)/(tg^2*a), (4, 8) = -k*n*BesselY(n, tg*a)/(tg^2*a), (5, 1) = 0, (5, 2) = 0, (5, 3) = k*n*BesselY(n, t0*b)/(t0^2*b), (5, 4) = I*mu0*omega*(diff(BesselY(n, t0*b), b))/t0^2, (5, 5) = -k*n*BesselJ(n, tg*b)/(tg^2*b), (5, 6) = -k*n*BesselY(n, tg*b)/(tg^2*b), (5, 7) = -I*`&mu;g`*omega*(diff(BesselJ(n, tg*b), b))/tg^2, (5, 8) = -I*`&mu;g`*omega*(diff(BesselY(n, tg*b), b))/tg^2, (6, 1) = 0, (6, 2) = 0, (6, 3) = BesselY(n, t0*b), (6, 4) = 0, (6, 5) = -BesselJ(n, tg*b), (6, 6) = -BesselY(n, tg*b), (6, 7) = 0, (6, 8) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = BesselY(n, t0*b), (7, 5) = 0, (7, 6) = 0, (7, 7) = -BesselJ(n, tg*b), (7, 8) = -BesselY(n, tg*b), (8, 1) = 0, (8, 2) = 0, (8, 3) = -I*epsilon0*omega*(diff(BesselY(n, t0*b), b))/t0^2, (8, 4) = -k*n*BesselY(n, t0*b)/(t0^2*b), (8, 5) = I*`&varepsilon;g`*omega*(diff(BesselJ(n, tg*b), b))/tg^2, (8, 6) = I*`&varepsilon;g`*omega*(diff(BesselY(n, tg*b), b))/tg^2, (8, 7) = -k*n*BesselJ(n, tg*b)/(tg^2*b), (8, 8) = -k*n*BesselY(n, tg*b)/(tg^2*b)});

and n=0,1,2.

Except of omega and k ,other parameters is canstant.

After using
with(LinearAlgebra):
detM := Determinant(M):

I used solve(detM=0,omega) and fsolve() but it dosnt work. how can i solve it?

Thanks alot.

hi .please see attached file below and help me

thanks..indices(res_nolist).mw

hi.may  help me for solve this nonlinear equations by numeric solver maple39.d39.pdfocx39.pdf

thanks alot

file format is pdf and word type

 

LE.2a.E.LGM.mwHi, my this programme is executing for linear part but does'nt show the proper results for non linear,plz tell me appropriate code

hi .how i can solve nonlinear equation with unknown prameter omega as below

thanksfrekans.mw

I want to solve numerically the nonlinear pde:

 

u_x+u_t - (u_{xt})^2 = u(x,t)

 

which method do you propose me to use with maple? (I don't mine about which boundary conditions to be used here).

 

hi.I want to dsolve set of nonlinear equations with one unknown parameter ...is this possible with dsolve rule.in matlab this possible with bvp4c rule..please help me for this problem.if we should another rule please attached file reform.Thanks alot12.mw

       Calculation of RSCR mechanism as a  solution to underdetermined system of nonlinear equations.  

 

RSCR.mw 

https://vk.com/doc242471809_376439263
https://vk.com/doc242471809_408704758

RCCC mechanism
https://vk.com/doc242471809_375452868

Hi .please help me for solve this nonlinear equations , that attch below

thanks a lots.....

 

 

dsys3 := {8*(diff(f2(x), x, x, x, x))+9*(diff(f2(x), x, x))+10*f2(x)+11*(diff(f1(x), x, x, x))+12*(diff(f1(x), x))+13*(diff(f3(x), x, x))+14*f3(x)+f3(x)*f3(x)+(diff(f3(x), x))*(diff(f3(x), x))+(diff(f3(x), x, x))*f3(x) = 0, 16*(diff(f3(x), x, x, x, x))+18*(diff(f3(x), x, x))+19*(diff(f3(x), x, x))+22*(diff(f1(x), x))+23*(diff(f1(x), x))+24*(diff(f2(x), x, x))+25*f2(x)+26*f2(x)+27*f3(x)+29*f3(x) = 0, 2*(diff(f1(x), x, x))+3*(diff(f2(x), x, x, x))+4*(diff(f2(x), x))+6*(diff(f3(x), x))+7*f1(x)+(diff(f3(x), x, x))*(diff(f3(x), x))+(diff(f3(x), x))*f3(x)+(diff(f3(x), x))*f3(x) = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0, f3(0) = 0, f3(1) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(1) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(1) = 0}; dsol5 := dsolve(dsys3, 'maxmesh' = 500, abserr = .1, numeric, range = 0 .. 1, output = listprocedure)12.mw

Hi,

I have a system of nonlinear equations in conjuction with boundary conditions, as below:

(all coefficients are constant.)

dsys3 := {a2*(diff(f1(x), x, x))+a3*(diff(f2(x), x, x, x))+a4*(diff(f2(x), x))+a6*(diff(f3(x), x))+a7*f1(x)+a1N*(diff(f3(x), x, x))*(diff(f3(x), x))+a2N*(diff(f3(x), x))*f3(x)+a3N*(diff(f3(x), x))*f3(x) = 0, a8*(diff(f2(x), x, x, x, x))+a9*(diff(f2(x), x, x))+a10*f2(x)+a11*(diff(f1(x), x, x, x))+a12*(diff(f1(x), x))+a13*(diff(f3(x), x, x))+a14*f3(x)+a4N*(f3(x)*f3(x))+a5N*((diff(f3(x), x))*(diff(f3(x), x)))+a6N*(diff(f3(x), x, x))*f3(x) = 0, (diff(f3(x), x, x, x, x))*a16+(diff(f3(x), x, x))*a18+(diff(f3(x), x, x))*a19+(diff(f1(x), x))*a22+(diff(f1(x), x))*a23+(diff(f2(x), x, x))*a24+f2(x)*a25+f2(x)*a26+f3(x)*a27+f3(x)*a29 = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0, f3(0) = 0, f3(1) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(1) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(1) = 0}:
dsol5 := dsolve(dsys3, 'maxmesh' = 500, abserr = .1, numeric, range = 0 .. 1, output = listprocedure)

but i encounter an error as:

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

what is your suggestion to solve this nonlinear equations?

Thank you.

Ok guys this is the problem i have a matrix wich is really the result of a semivariogram so it is like this 
plot([[0, 0], [1, 1], [2, 3], [3, 5], [4, 7], [5, 8], [6, 9], [7, 10], [8, 10.5], [9, 10], [10, 9.8], [11, 9.8], [12, 10.5], [13, 10.89], [14, 11.2], [15, 10.6], [16, 9.85], [17, 9.45], [18, 9.77], [19, 10.15], [20, 10.7], [21, 10.04], [22, 10], [23, 10.85], [24, 10.7], [25, 10.94], [26, 10.81], [27, 10.33]])thi is my semivariogram so i have to find the best model and fit it to it, i can't just throw a polynomial it has to be
spherical= h:= C((3/2)*x/a-(1/2)*x^3/a^3) 
expotencial=Upsilon := C(1-exp(-x/a))
Lineal= y:= p*x
or y:=p*x^a

so i want the maple to show me the variogram plot and the i select wich one of the models i will apply.
the problem is that i don't know how to apply the models so started something that could maybe lead to the values of C and A for the spherical and the xpotencial becouse the others too are easy.
------------------------------------------------------------------.--------------------------------------------------------
restart;
with(plots);
with(CurveFitting);
m3 := Matrix([[0, 0], [1, 2], [2, 5], [3, 7], [4, 10], [5, 14], [6, 15], [7, 15.4], [8, 15.8], [9, 16], [10, 16.4]]);
plot(m3);
Upsilon := C(1-exp(-x/a));
Upsilon := C((3/2)*x/a-(1/2)*x^3/a^3);
for C to 16 do
for a to 100 do
for x to 11 do
H := m3(x, 2);
R[C, a, x] := abs(Upsilon-H)
end do
end do
end do;
for r to 16 do
for j to 100 do
N[r, j] := sum(R[r, j, X], X = 1 .. 11);
if N[r, j] < 55 then
print(r, j)
end if
end do
end do;
k := 14*(3/2*((1/16)*xx)-(1/2)*xx^3/16^3);
multiple(plot, [m3, color = "Green"], [k, xx = 0 .. 11, color = "Blue"])
----------------------------------------------------------...................---------------------------------------
so in this code i rest the values of the model and the values of the matrix for differents C and a so i pick the min of the Sum to see wich C and a give me the min value but i didn't work quite well first is difficult to get the C and a and sometimes the new curve is just not similir to the semivariogram (Matrix) please help i tried 
k := NonlinearFit(C((3/2)*x/a-(1/2)*x^3/a^3), X, Y, x)
and 
k:=LeastSquares`(m3, v, curve = C*(1-exp(-x/a)))
but they don't let me because they aint linear on the parameters
please help 
if you are interested in this topic here is what i am trying to do but explained really well is for geological matters
http://www.kriging.com/PG1979/index.htm#Chapter_2+Part1

if want to help me use this matrix 
m3:= Matrix([[0, 0], [1, 1], [2, 3], [3, 5], [4, 7], [5, 8], [6, 9], [7, 10], [8, 10.5], [9, 10], [10, 9.8], [11, 9.8], [12, 10.5], [13, 10.89], [14, 11.2], [15, 10.6], [16, 9.85], [17, 9.45], [18, 9.77], [19, 10.15], [20, 10.7], [21, 10.04], [22, 10], [23, 10.85], [24, 10.7], [25, 10.94], [26, 10.81], [27, 10.33]])
for that matrix the spherical model is the best i know

goodbye

1 2 3 4 5 6 Page 1 of 6