Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
ZigangPan`RobustBackSteppingCancellation`
RobustBackSteppingCancellationG  |
|
| | ||||
Examples
(1)
Basic Examples
(1)
In[1]:=
Needs["ZigangPan`Examples`"]
In[2]:=
Needs["ZigangPan`DifferentialEquationSolver`"]
In[3]:=
xc={x1,x2,x3,x4};$Assumptions=assumeReal[xc]
Out[3]=
x1∈&&x2∈&&x3∈&&x4∈
In[4]:=
fe[x1_,x2_,x3_,x4_]:={x3+x1^2,x1*x2+x4,x1*x3,x2*x4}
In[5]:=
ge[x1_,x2_,x3_,x4_]:={{0,0},{0,0},{1,0},{0,-1}}
In[6]:=
he[x1_,x2_,x3_,x4_]:={{1,0},{-1,x1},{0,0},{1,1}}
In[7]:=
V0=0;α0={0,0};γ=4;f=Apply[fe,xc];g=Apply[ge,xc];h=Apply[he,xc];
In[8]:=
{V1,α1,x,σ}=
[V0,α0,{},{},f〚1;;2〛/.{x30,x40},D[f〚1;;2〛,{{x3,x4}}],h〚1;;2,All〛,{},{x1,x2},{},γ^2*{{1,1/2},{1/2,2}},{x1,x2},IdentityMatrix[2]/2]
 | RobustBackSteppingCancellationG |
Out[8]=
+,-x1-+----x2,-x2-x1x2+---x1-+++x1x2,{x1,x2},+--x2,-++x2
2
x1
2
2
x2
2
2
x1
1
4
x1
14
1
14
x1
56
1
4
1
14
x1
56
1
14
x1
56
1
56
x1
28
1
2
x1
14
1
14
x1
56
1
2
x1
56
1
56
x1
28
In[9]:=
{V,α,x,σ}=
[V1,α1,f〚1;;2〛,h〚1;;2,All〛,f〚3;;4〛,g〚3;;4,All〛,h〚3;;4,All〛,{x1,x2},{x3,x4},{},γ^2*{{1,1/2},{1/2,2}},{x1,x2}-α1,1/2*IdentityMatrix[2]];
 | RobustBackSteppingCancellationG |
In[10]:=
closedloopdynamics=FormulaToFunction[xc,f+g.α+h.{0,0}];
In[11]:=
sol=RungeKutta45[closedloopdynamics,0,{-2,2.,1,0},20,1/10];
10% completed.
20% completed.
30% completed.
40% completed.
50% completed.
60% completed.
70% completed.
80% completed.
90% completed.
100% completed.
In[12]:=
Plot[sol[t],{t,0,20},PlotRangeAll]
Out[12]=
In[13]:=
Plot[Apply[FormulaToFunction[xc,α],sol[t]],{t,0,20},PlotRangeAll]
Out[13]=
""