Mathematica Mathematica: ill-conditioned linear system

AI Thread Summary
The discussion centers on challenges faced when solving ill-conditioned linear systems of simultaneous equations using Mathematica. The user successfully solved well-conditioned systems for smaller matrices (m = n = 1 to 7) but encountered issues starting at m = n = 9 due to the matrices becoming ill-conditioned. It is emphasized that when a matrix is ill-conditioned, Mathematica can only report this condition, as the issue lies within the input data rather than the software itself. The user is employing a state space approach to calculate stresses in a plate, which involves summing terms for accuracy. However, problems arise with large numbers leading to unrealistic stress values, highlighting the need for a convergence criterion to prevent adding excessive terms that could cause divergence. The discussion concludes with the recommendation to establish a convergence criterion to manage the summation process effectively and avoid ill-conditioning.
Asad Albostami
Messages
2
Reaction score
0
Hi all, I have problem with regard to ill-conditioned linear system of solving sets of simultaneous equations using Mathematica program. I have tried my best to find a way to solve this but none was successful.
I got results from m =1 and n =1 until m = 7 and n = 7, i,e. the systems are well-conditioned. However, the problem started from when m = 9 and n = 9 due to ill-conditioned systems. Note that m = n = 1,3,5,6,7,9,...99. I am trying to find a treatment so that the matrix would not become a singular.
 
Physics news on Phys.org
If the matrix is ill conditioned there is nothing that Mathematica can do other than to correctly report the ill conditioning. Basically, you are trying to divide by zero, and no math program can make that reasonable, the problem is the input itself not the program.
 
Thank you for your reply.
For my program I used a method it is called state space approach for finding the stresses through the thickness for any plate. From this method the plate has a top and bottom values and from finding the top stresses from an exponential matrix then I can calculate the bottom stresses but here I have a problem for a very big numbers that It gives me (For example if the top stresses is equal to 1 MPa the bottom stress is equal to 10^5). Please see the attached file it shows the formulas that I used. noting that the problem started when the looping numbers is increased (Because this method it is a summation of terms So more terms it means more accurate results).
 

Attachments

The .pdf file doesn't work for me. However, in general you don't want to just blindly add 100 terms. You should specify some convergence criterion and stop adding terms as soon as the criterion is met or the criterion begins to indicate divergence. An ill conditioned matrix should cause any reasonable criterion to diverge, and a good criterion should stop adding terms before it becomes ill conditioned.
 

Similar threads

Replies
1
Views
3K
Replies
5
Views
3K
Replies
1
Views
2K
Replies
1
Views
2K
Replies
2
Views
2K
Replies
6
Views
4K
Back
Top