The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Matrix inversion - linear algebra - higher accuracy for some matrix
rows ? (Tikhonov regularization???)

Replies: 6   Last Post: Sep 11, 2012 12:49 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Loren Shure

Posts: 1,067
Registered: 12/7/04
Re: Matrix inversion - linear algebra - higher accuracy for some matrix rows ? (Tikhonov regularization???)
Posted: Sep 11, 2012 11:07 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"someone" <> wrote in message
> Here's the deal:
> A-matrix is approx. 90x90, mostly diagonal but also quite some offdiagonal
> elements here and there.
> Some rows in the matrix are equations that is harder to satisfy than
> others because something is rotating at different speeds - it means that
> those matrix equations (the lower rows in A) that has a physical
> connection to something that rotates really fast, causes some severe
> oscillations (it gives oscillating accelerations and wrong forces) because
> the timestep is very high compared to the rotation speed for the last rows
> in A... Got it?
> Ok, is there any way to make: x = A\b more accurate for the lower
> rows/equations or ???
> Another possibility I've been thinking about - in order to avoid these
> oscillations I talked about before - is to use regularization, maybe? I
> did a project some years ago using Tikhonov regularization, but here
> AFAIR - the idea was that all the elements in the solution-vector x could
> not change too much... This will allow me to use higher timesteps and
> avoid the oscillations... Understand it ?
> Here's an elaboration of why I have a problem with high angular
> velocities: If some part of the problem (top rows in A) rotates very slow,
> i.e. at 1 rad/s and your timestep is 0.01 sec - it means that you have 100
> timesteps per rotation. That's ok... However, the lower rows in A rotates
> much faster - then I only get maybe 5-8 timesteps per rotation. This gives
> some VERY annoying oscillations and forces going in "random" directions
> because the motion is not smoothly captured.
> I want my program to solve x = A\b quickly and efficent - any ideas that
> might help me (and keep the computation time reasonable) ? Any tricks /
> tips ?
> I really hope to hear from some clever experts here... Thanks...

You might try lscov:

It will allow you to weight the observations.


Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.