Symbolic computation: Mathematica help

  • Context: Mathematica 
  • Thread starter Thread starter confi999
  • Start date Start date
  • Tags Tags
    Computation Mathematica
Click For Summary

Discussion Overview

The discussion revolves around obtaining a simplified equation by setting the determinant of a specific matrix expression to zero using Mathematica. The expression involves an identity matrix and a complex symbolic matrix derived from prior computations. Participants are exploring methods to simplify this expression and the challenges associated with symbolic computation.

Discussion Character

  • Technical explanation
  • Exploratory
  • Mathematical reasoning

Main Points Raised

  • One participant requests assistance with Mathematica code to simplify the determinant of the expression involving an identity matrix and a symbolic matrix G.
  • Another participant inquires about the nature of matrix G, asking whether it is numerical or symbolic, and what types of terms it contains (algebraic, transcendental, etc.).
  • It is noted that simplifying the determinant could be complex if G is symbolic, and using the FullSimplify function with valid assumptions may be necessary.
  • A later reply provides details about G, indicating it is symbolic, algebraic, and complex, with a specific structure that includes many zeroes and long algebraic expressions.
  • Suggestions are made for potential strategies to simplify the expression, such as changing variables, altering the coordinate system, or performing a series expansion around specific points.

Areas of Agreement / Disagreement

Participants express a consensus on the complexity of the expression due to the symbolic nature of G, but there is no agreement on a definitive method for simplification, as various approaches are suggested without a clear resolution.

Contextual Notes

The discussion highlights the challenges of simplifying large symbolic matrices in Mathematica, including potential memory issues and the need for careful handling of assumptions during simplification.

confi999
Messages
18
Reaction score
0
Hi,
I need to get a simplified equation by setting the determinant of the following expression to zero:
exp(j.omega.dt). - [G]

Here, I is an 18x18 identity matrix and G is an 18x18 matrix that I got by doing some symbolic computation in Mathematica (in fact G is huge).

Can anyone help me with Mathematica code to get the simplified equation from
det| exp(j.omega.dt). - [G] | =0

In the literature people used Maple for the above. I have access to Mathematica and would prefer it. Thank you very much
 
Physics news on Phys.org
Is G numerical or symbolic? If it is symbolic are the terms algebraic or transcendtal or irrational or rational or what? Is G symmetric or off diagonal or sparse in any way? Are the terms real or complex?

In general, this will be a hairy expression if G is symbolic. It should evaluate, but simplifying may be a pain. Generally you will just use FullSimplify and explicitly give any assumptions that are valid. If it is particularly complicated then it may take a long time to simplify it or it may run out of memory. Try to do any simplification on the terms of G first and only then calculate the determinant.
 
DaleSpam said:
Is G numerical or symbolic? If it is symbolic are the terms algebraic or transcendtal or irrational or rational or what? Is G symmetric or off diagonal or sparse in any way? Are the terms real or complex?

In general, this will be a hairy expression if G is symbolic. It should evaluate, but simplifying may be a pain. Generally you will just use FullSimplify and explicitly give any assumptions that are valid. If it is particularly complicated then it may take a long time to simplify it or it may run out of memory. Try to do any simplification on the terms of G first and only then calculate the determinant.

Hi DaleSpam,
Thanks a lot for replying.

G is symbolic, algebraic and complex. It was obtained from some symbolic computation in Mathematica (at one stage to get G it involves matrix inversion ). Each column of 18x18 matrix G has roughly 9 zeroes, 1 one and rest of the elements (around 8 / 9 on each column) are long algebraic expression. Example of two such elements are copied below - I showed numerator, denominator by placing a long line

i)

( (v4 + 2 c1 c2 v1 (2 + v4) - 2 c1 c2 v1 (v4 Cos[dx kx] + Cos[dy ky] + Cos[dz kz]))
---------------------------------------------------------------------------------------
-1 - 6 c1 c2 v1 + 2 c1 c2 v1 (Cos[dx kx] + Cos[dy ky] + Cos[dz kz])ii)
-((c1 c2 E^(-I dx kx) (-1 + E^(I dx kx)) (-1 + E^(I dy ky)) v1 (-1 + v4))
--------------------------------------------------------------------------------------------
(-1 - 6 c1 c2 v1 + 2 c1 c2 v1 (Cos[dx kx] + Cos[dy ky] + Cos[dz kz])))

FullSimplify command is already used before getting the above G.
Any advice. Thank you so much
 
Last edited:
Well, your expression is inherently hairy. You can do things like changes of variables or changes of coordinate system, e.g. let dx->Cx/kx and v1->V1/(c1 c2). Also, you can do a Series expansion about some particular point of interest. For instance if you expand to first order about dx=dy=dz=0 then expression i) reduces to -v4.
 

Similar threads

  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 6 ·
Replies
6
Views
7K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
4K