1
$\begingroup$

I have a square matrix depending on a parameter "a". When "a" has a numeric value (eg 2.1 * 10 ^ 8) through the inverse function Inverse[] I can invert the matrix. If I write the matrix as a function of "a" do not declare its value when the invert Mathematica tells me that it is singular (suspect singular) then not allow me to Invert it. The determinant of the matrix in symbolic form turns out to be a very small value in the order of 6*10 ^ -25000. Do you think I can get the inverse of this matrix in symbolic form? the matrix is generated by:

K0 = Ct.E0.Transpose[Ct];

where Ct is a matrix indipendent from "a" , while "E" is a diagonal matrix which depends on "a" in the following way :

https://i.stack.imgur.com/VXX2A.png

The matrix Ct and E are written manually.

I attach the file containing the matrix mathematica https://drive.google.com/file/d/0B7In1ClJQ-FVZVRuNHMycGotUTg/edit?usp=sharing

I was not able to include them here as a code. Thanks in advance

$\endgroup$
10
  • $\begingroup$ Can you provide the matrix/the code to generate it ? $\endgroup$
    – Sektor
    Mar 15, 2014 at 13:16
  • $\begingroup$ thanks Sektor This is the det of simbolyc matrix : i.imgur.com/OYJjxET.png This is the part of simbolyc matrix (is a large matrix): i.imgur.com/72Sdqrr.png This is the errore reported by mathematica:i.imgur.com/lDqzBrQ.png I don't understand becouse assigning a numerical value to "a" the matrix to be inverted, if I leave everything in a symbolic way, I can't reverse it. $\endgroup$
    – plus91
    Mar 15, 2014 at 13:42
  • 1
    $\begingroup$ No, you should edit your Q and add any relevant code :) BTW Not images, code, because it is tedious to reproduce the problem :) $\endgroup$
    – Sektor
    Mar 15, 2014 at 14:00
  • $\begingroup$ Sektor I edited the Q thx :) $\endgroup$
    – plus91
    Mar 15, 2014 at 14:30
  • 1
    $\begingroup$ Don't post images. Post code instead, please. $\endgroup$ Mar 15, 2014 at 14:46

1 Answer 1

4
$\begingroup$

I think it's just a matter of handling the numerics; you have very small matrix elements which give rows/columns of all zeros when chopped.

You can avoid that by transforming the elements of your matrices to be rationals :

c = Rationalize[Cttt, 10^(-50)];
e = Rationalize[E0, 10^(-50)];
k = Inverse[c.e.Transpose[c]];

Check :

Block[{a = 2.1 10^8}, Det[Cttt.E0.Transpose[Cttt]]]
(* 7.41902*10^100 *)

1/Det[k] /. a -> 2.1 10^8
(* 7.41902084900741*10^100 *)
$\endgroup$
4
  • $\begingroup$ Thank you very much. The procedure that you have shown me, seems to solve the problem. However, I will do some tests on the results obtained from the solution of the system using the inverse matrix of K. I'll let you know soon. $\endgroup$
    – plus91
    Mar 16, 2014 at 0:19
  • $\begingroup$ @ b.gatessucks thank you, if the matrix "k" depends on a single parameter the procedure works. If the matrix depends on several parameters (a, b, c ...) the calculation of the inverse matrix requires a high computational time. Do you think there is another way to overcome the problem? The option ZeroTest associated with the Inverse function [] might be useful? $\endgroup$
    – plus91
    Mar 16, 2014 at 11:14
  • 1
    $\begingroup$ I'd rather spend some time thinking about another way of setting up the problem, so that you don't end up with matrix elements differing by 30+ orders of magnitude. $\endgroup$ Mar 16, 2014 at 16:15
  • $\begingroup$ @ b.gatessucks thanks for your time.The problem concerns the solution of the system to find the nodal displacements of the trusses. The matrices are related to physical parameters that can not vary. At the moment I have no other ideas on how to set up the problem. $\endgroup$
    – plus91
    Mar 16, 2014 at 18:46

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.