Linearização Carleman 

A rica dinâmica dos sistemas não lineares pode ser apenas parcialmente capturada por linearizações aproximadas. Avalie a eficiência de diferentes ordens de linearização Carleman (Carleman linearization) em vários pontos no espaço de fase de um oscilador de Van der Pol.

O modelo de um oscilador de van der Pol.

In[1]:=

X vanderpol[\[Mu]_] := AffineStateSpaceModel[ {{Subscript[x, 2], -Subscript[x, 1] - \[Mu] (-1 + \!\(\*SubsuperscriptBox[\(x\), \(1\), \(2\)]\)) Subscript[x, 2]}, {{ }, {}}}, {Subscript[x, 1], Subscript[x, 2]}, {}, {Automatic, Automatic}, Automatic , SamplingPeriod -> None]

Sua resposta para 100 unidades de tempo a partir de um ponto inicial p.

In[2]:=

X vanderpolResponse[\[Mu]_, p_] := OutputResponse[{vanderpol[\[Mu]], p}, {}, {t, 0, 100}]

Seu retrato de fase.

In[3]:=

X vanderpolPhasePortrait[\[Mu]_] := With[{v = Normal[vanderpol[\[Mu]]][[1, 1]] /. x_[\[FormalT]] :> x}, StreamPlot[v, {Subscript[x, 1], -3, 3}, {Subscript[x, 2], -3, 3}, StreamStyle -> "PinDart"]]

A resposta de um sistema linearizado por Caleman com ordem de aproximação n.

In[4]:=

X approxResponse[\[Mu]_, p_, n_] := With[{\[ScriptCapitalC] = CarlemanLinearize[vanderpol[\[Mu]], n]}, OutputResponse[{\[ScriptCapitalC]["TransformedSystem"], Last /@ \[ScriptCapitalC]["EmbeddingTransformation"] /. Thread[{Subscript[x, 1], Subscript[x, 2]} -> p]}, {}, {t, 0, 100}]]

Visualize as trajetórias reais e aproximadas a partir de vários pontos iniciais.

In[5]:=

X Manipulate[ Show[vanderpolPhasePortrait[\[Mu]], ParametricPlot[ Evaluate[{vanderpolResponse[\[Mu], p]}~ Join~(approxResponse[\[Mu], p, #] & /@ n)], {t, 0, 100}, PlotLegends -> Join[{"Exact"}, n]], ImageSize -> Medium], {{p, {0.01, 0}}, {-3, -3}, {3, 3}, Locator}, {{\[Mu], 0.6, "Parameter (\[Mu])"}, 0, 1.2, 0.1}, {{n, {1, 4}, "Approx. system(s)"}, CheckboxBar[Dynamic@n, {1, 2, 4, 6, 8, 10}] &}, SaveDefinitions -> True, AutorunSequencing -> {1, 2}]

Reproduzir Parar