Transforming association matrix to a functional matrix

  • Context: Graduate 
  • Thread starter Thread starter Pythagorean
  • Start date Start date
  • Tags Tags
    Functional Matrix
Click For Summary
SUMMARY

The discussion focuses on transforming an association matrix into a functional matrix for a system of n elements, specifically using MATLAB for operations like circular shifts and transposes. The initial association matrix is defined for three elements, leading to a functional matrix that captures unidirectional diffusion coupling between electrophysiological neurons. A general formula is proposed: F = transpose(A) - diag(sum(rows(A))), which allows for the computation of the functional matrix from any association matrix, accommodating both unidirectional and two-way diffusion scenarios.

PREREQUISITES
  • Understanding of association and functional matrices in mathematical modeling
  • Familiarity with MATLAB commands, specifically circshift and transpose
  • Knowledge of diffusion processes in biological systems
  • Basic linear algebra concepts, including matrix operations and diagonal matrices
NEXT STEPS
  • Explore advanced MATLAB techniques for matrix manipulation and optimization
  • Research the mathematical principles behind diffusion processes in biological systems
  • Learn about the implications of two-way diffusion in neural networks
  • Investigate numerical methods for solving differential equations related to diffusion
USEFUL FOR

Researchers, mathematicians, and engineers working on mathematical modeling of biological systems, particularly those focused on neural dynamics and diffusion processes.

Pythagorean
Science Advisor
Messages
4,430
Reaction score
327
I have n elements. Say n = 3.

Suppose I have an association matrix that gives the relationship between each element

[tex] \begin{array}{cc} <br /> 0 & 0 & D3\\ <br /> D1 & 0 & 0\\<br /> 0 & D2 & 0 <br /> \end{array}[/tex]

I have a function in mind now, I want to operate and the physical variables representing my three elements, [tex]\vec{y} = [y_1 y_2 y_3][/tex]. My functional matrix would look like:

[tex] \begin{array}{cc} <br /> -D1 & D1 & 0\\ <br /> 0 & -D2 & D2\\<br /> D3 & 0 & -D3 <br /> \end{array}[/tex]

so that I get [tex]\vec{y} = [D1(y_2 - y_1) D2(y_3-y_2) D3(y_1-y_3)[/tex]

If you're interested in the physical/biological motivation, we basically have a unidirectional diffusion coupling between electrophysiological neurons here so you're seeing a numerical second derivative. Now, using MATLAB commands (circshift and transpose) I do circular shifts and transposes on the association matrix to nudge it into my functional shape.But I'm having trouble with the more general case. What if I have two-way diffusion, but the diffusion is stronger in one direciton than in the other? Now the association matrix is:

[tex] \begin{array}{cc} <br /> 0 & D4 & D3\\ <br /> D1 & 0 & D5\\<br /> D6 & D2 & 0 <br /> \end{array}[/tex]

and what we for functional is:[tex] \begin{array}{cc} <br /> -(D1+D4) & D1 & D4\\ <br /> D5 & -(D2+D5) & D2\\<br /> D3 & D6 & -(D3 + D6) <br /> \end{array}[/tex]

so that [tex]\vec{y} = [D1(y_3-y_1) + D4(y2_y1); D2(y_1- y_2) + D5(y_3-y_2) + ...][/tex]
Now, I can design another series of circshifts and tranposes, but it won't work for the case above. I can't find a general set of operations that works for both

This should work in general, for an nxn matrix.

Thank you for your help.
 
Physics news on Phys.org
A friend and co-worker came up with this (F = functional matrix, A = association matix)

F = transpose(A) - diag(sum(rows(A)))

that is, if you take the sum of the rows of A and place it on a diagonal of an otherwise 0 matrix, then subtract that from the transpose of A, you have the functional matrix of diffusive coupling for that association matrix!
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 23 ·
Replies
23
Views
2K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
8K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 1 ·
Replies
1
Views
5K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K