Survey for Special Math Problem



#6

We are all familiar with the system of linear equations:

a11 x1 + a12 x2 + ... + a1n xn = b1
a21 x1 + a22 x2 + ... + a2n xn = b2

...
an1 x1 + an2 x2 + ... + ann xn = bn

Has anyone seen an application for the variant of the above such that:

a11 x1^2 + a12 x2 + ... + a1n xn = b1
a21 x1 + a22 x2^2 + ... + a2n xn = b2
...
an1 x1 + an2 x2 + ... + ann xn^2 = bn

In other words, where the i'th diagonal element has the square (or even any other power that is not equal to one) of xi instead just xi.

Thanks!

Namir

Edited: 2 June 2012, 5:42 p.m.


#7

Hi, Namir.

I am curious now. You have surely imprinted the question in my mind.

Please, allow me a candid question: do you actually want to know if this sort of 'hybrid' system - not linear, though - is used in some particular application or you already know and would like to check if anyone else has also seen it?

Thanks for bringing the subject up. I confess this is new to me.

Cheers.

Luiz (Brazil)


#8

Luiz,

I have been studying iterative linear methods (both stationary and non-stationary). The problem I am asking about is something that just popped in my mind two days. Interestingly, you can solve for the matrix coefficients a(i,j) for the 'hybrid' system using the iterative Gauss-Seidel method that is also used for linear system.

My question is, does the problem I am asking about model some application?

namir


#9

Quote:
Interestingly, you can solve for the matrix coefficients a(i,j) for the 'hybrid' system using the iterative Gauss-Seidel method that is also used for linear system.
When dealing with electrical distribution networks and load-flow studies, back in the days I was concluding the Electrical Engineering course, I wrote a program for the HP42S to solve a load-flow system by using an alternate Accelerated Gauss-Siedel method. This was done because I decided not to use a computer and FORTRAN back then (1990). The teacher accepted my solution only because I implemented the whole iteration process without using the HP42S SOLVE, only discrete program steps. As you mentioned, Gauss-Siedel is used in this case because the load-flow studies presume a linear system.

Not too much to add, sorry...

Luiz (Brazil)

#10

Can't say I have seen this particular nonlinear system treated as a special case, but on looking at its symmetries it would seem that general solution approaches would be a little simpler to set up.

For example, if you were to take a multidimensional Newton-Raphson approach, the Jacobian matrix that appears in the expression of that problem is simply an nxn matrix with 2*aii*xi along the diagonal and the remaining elements simply constant values of the original matrix. With this sort of redundancy you save a lot of register space when storing elements.

Les


#11

Les,

The iterative solution for these simple nonlinear equations can employ Gauss-Sediel which works for liner system. My question to all is, "Have you seen an application that can be modeled with such system of nonlinear equations."


#12

I just read up on Gauss-Seidel and it is very easy to understand. Can you tell us how the algorithm is modified to accommodate your nonlinear diagonal terms?

Les


#13

To solve the nonlinear/hybrid system, here is the core steps expressed in Excel VBA:

For J = 1 To N
Sum = B(J)
For I = 1 To N
If I <> J Then
Sum = Sum - A(J, I) * X(I)
End If
Next I
X(J) = Sqr(Sum / A(J, J))
Next J

The core steps of the Gauss-Seidel for linear systems is:

For J = 1 To N
Sum = B(J)
For I = 1 To N
If I <> J Then
Sum = Sum - A(J, I) * X(I)
End If
Next I
X(J) = Sum / A(J, J)
Next J

Since we using the squares of the variables n the diagonal elements, calculating them iteratively would require the square root values of Sum/A(I,I).

Using the minor modification of the first code snippet allows Gauss-Seidel method to work with the nonliner/hybrid system of equations.

Namir


Possibly Related Threads…
Thread Author Replies Views Last Post
  Need help understanding math.... cyrille de BrĂ©bisson 9 3,951 12-13-2013, 02:23 AM
Last Post: Didier Lachieze
  OT: a math competition site Pier Aiello 0 1,079 09-16-2013, 06:03 AM
Last Post: Pier Aiello
  Simple Math Question Namir 2 1,431 08-09-2013, 06:13 PM
Last Post: Eddie W. Shore
  Cool math clock Bruce Bergman 28 8,023 04-10-2013, 03:13 AM
Last Post: Siegfried (Austria)
  HP 41CX Survey Pac Keith Burr 7 2,655 03-29-2013, 10:09 AM
Last Post: Keith Burr
  HP-41C Custom Survey Module Identificaiton Help Dan Grelinger 0 953 01-27-2013, 01:51 PM
Last Post: Dan Grelinger
  Math Challenge I could not solve Meindert Kuipers 22 6,339 01-05-2013, 04:43 PM
Last Post: Thomas Klemm
  Error Message/Procedure HP-48SX with TDS Survey SX card bporter101 6 2,074 12-10-2012, 10:18 AM
Last Post: bporter101
  Math Question Namir 0 958 11-06-2012, 07:43 AM
Last Post: Namir
  Survey: Best programming language for next-gen HP? Oliver Unter Ecker 32 8,810 06-06-2012, 03:22 PM
Last Post: David Hayden

Forum Jump: