ENGINEER · PILOT · BUILDER

Research & Theory · 2024

Model-Based Learning & Control Studies

Built three MATLAB studies of nonlinear dynamics: a data-driven Van der Pol predictor, operator-based analysis, and iterative trajectory steering.

Project type
Individual graduate nonlinear-dynamics studies
Role
Individual MATLAB implementation, numerical experiments, and analysis
Maturity
Nonlinear dynamics studies · MATLAB
Collaboration
Individual study
Local flow predictor
Delaunay
Operator analysis
FP / EDMD
Iterative steering
LQ
Primary nonlinear system
Van der Pol
Read the triangulation predictor report
01

The studies

This page connects three MATLAB studies of nonlinear systems. The first learned local Van der Pol motion from sampled transitions. The second used Frobenius–Perron, Koopman, and clustering views to reveal structure in the same oscillator. The third repeatedly linearized, steered, and rolled out nonlinear trajectories on a nonholonomic integrator and a controlled Van der Pol system.

SYSTEM ARCHITECTURE

Learn local motion

1,000 sampled transitions → Delaunay mesh → local affine maps → rollout

Each triangle approximated one step of the Van der Pol flow. Repeated local predictions were then compared with direct ode45 integration.

Expose nonlinear structure

snapshot data → density transport + EDMD observables + clustered local flow

Operator and clustering methods showed where trajectories concentrate, which lifted quantities may remain structured, and where the local flow behaves similarly.

Steer the nonlinear plant

nominal trajectory → local linearization → steering update → nonlinear rollout → repeat

The iterative controller first satisfied the endpoint task and then reduced input energy, checking every local update against the nonlinear dynamics.

02

What I implemented

  1. STAGE 01

    Build the sampled flow map

    Sampled 1,000 Van der Pol states, advanced each by 0.1 seconds with ode45, and triangulated the paired snapshots into a piecewise-affine flow map.

  2. STAGE 02

    Test multi-step prediction

    Iterated the learned map as a multi-step predictor and compared its trajectory with direct integration from the same initial condition.

  3. STAGE 03

    Analyze the nonlinear structure

    Estimated density transport with a Monte Carlo Frobenius–Perron matrix, applied EDMD with nine polynomial observables, and clustered related local flow maps.

  4. STAGE 04

    Iterate nonlinear steering

    Implemented a 200-step linearize–steer–rollout loop on a nonholonomic integrator, then applied the same structure to controlled Van der Pol dynamics.

03

Interpretation and scope

The learned map remained local

Prediction quality depended on sample density, triangle geometry, and staying inside the sampled region. A close overlay on one trajectory verified the implementation but did not establish accuracy across the full state space.

Numerical modes were candidates, not proofs

The dominant density mode and eigenvalue-one EDMD observable depended on the grid, sample count, and chosen polynomial dictionary. They reveal computational structure but still require analytical verification before being called exact invariants.

Every steering update needed nonlinear validation

A successful linearized update could still miss when applied to the nonlinear plant. Rolling out both trajectories at every iteration exposed that gap and kept endpoint steering separate from later energy reduction.

04

What worked

  • The triangulation predictor remained visually close to ode45 on the documented multi-step trajectory toward the Van der Pol limit cycle.
  • Frobenius–Perron evolution moved a uniform density toward the oscillator's limit cycle, while the clustering curve showed diminishing improvement beyond roughly eight groups.
  • Nonholonomic and Van der Pol animations compare the nominal, linearized, and nonlinear trajectories across successive steering updates.
  • Input-energy plots decreased after the endpoint task was satisfied in the demonstrated iterative-control runs.
05

Takeaway

A learned dynamics model is only trustworthy inside the region supported by its data, so multi-step rollout against the known plant mattered more than a clean one-step fit.

Local linearization was useful for proposing control updates, but the nonlinear rollout remained the deciding test for whether each update actually improved the trajectory.

06

Study reports

NEXT PROJECT / 08

Manipulator Kinematics & Dynamics Studies