Скачать или смотреть Understanding the Eigenvalue Issue with Scipy: Solving for Correct Frequencies ω

This guide explains a common mistake when calculating eigenvalues in Python using Scipy, particularly in vibration analysis. We take a closer look at a specific example to clarify how to derive the correct frequency values from the computed eigenvalues.
---
This video is based on the question https://stackoverflow.com/q/68231509/ asked by the user 'buğra caner' ( https://stackoverflow.com/u/14021838/ ) and on the answer https://stackoverflow.com/a/68231684/ provided by the user 'Warren Weckesser' ( https://stackoverflow.com/u/1217358/ ) at 'Stack Overflow' website. Thanks to these great users and Stackexchange community for their contributions.

Visit these links for original content and any more details, such as alternate solutions, latest updates/developments on topic, comments, revision history etc. For example, the original title of the Question was: Scipy is not showing correct eigenvalues?

Also, Content (except music) licensed under CC BY-SA https://meta.stackexchange.com/help/l...
The original Question post is licensed under the 'CC BY-SA 4.0' ( https://creativecommons.org/licenses/... ) license, and the original Answer post is licensed under the 'CC BY-SA 4.0' ( https://creativecommons.org/licenses/... ) license.

If anything seems off to you, please feel free to write me at vlogize [AT] gmail [DOT] com.
---
Understanding the Eigenvalue Issue with Scipy

When working with numerical computations in Python, particularly in scientific contexts such as vibration analysis, you might run into some confusing results. A frequent query from users revolves around the eigenvalues computed using the Scipy library, specifically when dealing with systems modeled by mass and stiffness matrices. One such situation arises when users expect to see specific eigenvalues and instead observe different results. This post explores a typical scenario to understand why this occurs and how to solve it.

The Problem at Hand

Suppose you have a mass matrix M and a stiffness matrix K, which together describe a dynamic system. The goal is to compute the eigenvalues of the system represented by the equation M^-1 * K. Upon running the provided code snippet, results do not match the expected eigenvalues of 3.17, 9.13, 13.98, and 17.15. Instead, an unexpected array is returned.

Example Code

Here's the code provided for context:

[[See Video to Reveal this Text or Code Snippet]]

Understanding the Output

Upon examining the output from the eigenvalue calculation (la), you would get:

[[See Video to Reveal this Text or Code Snippet]]

The confusion arises from interpreting these eigenvalues. The key point to notice is that these values represent the squares of the actual frequencies (ω²).

Extracting the Frequencies

To derive the expected frequencies, you need to take the square root of your computed eigenvalues:

[[See Video to Reveal this Text or Code Snippet]]

The output from this operation will yield:

[[See Video to Reveal this Text or Code Snippet]]

As you can see, these values match the expected frequencies that were initially stated!

The Root of the Confusion

The main takeaway from this analysis is recognizing the relationship between eigenvalues and the physical quantities they represent:

Eigenvalues computed through matrix formulations related to mechanical systems often correspond to frequencies squared.

The codes provided and the mathematical framework utilized correctly implement the algorithm to extract the eigenvalues, but the interpretation error arose in not realizing that a square root was needed to determine the actual frequencies.

Conclusion

In summary, when working with eigenvalues in Python using Scipy, particularly within the context of systems involving mass and stiffness matrices, keep in mind the mathematical principles at play. If the computed eigenvalues seem off, check if they are in a squared format and remember to apply the square root to extract the expected physical quantities, like vibration frequencies.

This guide aimed to clarify why incorrect eigenvalues might be shown and how to ensure that you are calculating and interpreting them correctly. If you're embarking on your analytical journey with eigenvalues, follow these guidelines, and you'll be well on your way to accurate results.