Linearizing dynamics - derivative of matrix wrt vector

Click For Summary
SUMMARY

The discussion centers on the linearization of non-linear equations of motion represented by the equation x_dot(x,u) = B(x)*u + A(x). The user seeks to compute the derivative of the matrix B with respect to the vector x, specifically focusing on the Jacobian. The user clarifies that dA/du = 0 and d(Bu)/du = B, leading to the conclusion that d(Bu)/dx is indeed the Jacobian of the vector with respect to the vector. This realization indicates a solid understanding of the mathematical principles involved in the linearization process.

PREREQUISITES
  • Understanding of non-linear dynamics and equations of motion
  • Familiarity with matrix calculus and Jacobians
  • Knowledge of state-space representation in control theory
  • Proficiency in linearization techniques for differential equations
NEXT STEPS
  • Study the derivation and applications of the Jacobian matrix in dynamic systems
  • Explore linearization methods for non-linear control systems
  • Learn about state-space models and their implications in control theory
  • Investigate the role of perturbation analysis in system dynamics
USEFUL FOR

Engineers, control theorists, and researchers involved in dynamic systems modeling and analysis will benefit from this discussion, particularly those focusing on linearization techniques in non-linear dynamics.

alexjbuck
Messages
8
Reaction score
0
Hi,

I'm working on a project (not homework) where I have to linearize a set of non-linear equations of motion

They take the form:

x_dot(x,u) = B(x)*u + A(x)

where A and B are functions of state vector x. x_dot is a time derivative of x.
A - n x 1
B - n x m
u - m x 1

(wish I could get the latex to work... Its probably not broken, but rather I am...)

I also know that: (linearizing about some arbitrary x and u)

delta(x_dot) = (d(Bu)/dx + dA/dx)*delta x + (d(Bu)/du + dA/du)*delta u

I know dA/du = 0 and also that d(Bu)/du = B, but then the question is, how do I evaluate the first half? If I'm not mistaken it involves the derivative of a matrix with respect to a vector and I do not know how to go about this.

Thanks!
 
Physics news on Phys.org
I just realized that Bu is an n x m times a m x 1 thus the product is an n x 1 vector.

d(Bu)/dx is then just the derivative of a vector wrt a vector, which is the Jacobian if I'm not mistaken...

looks like I just answered my own question (I think)
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 13 ·
Replies
13
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
4K
Replies
1
Views
1K
  • · Replies 4 ·
Replies
4
Views
4K