MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Asymmetric quadratic coefficient problem

  • To: mathgroup at
  • Subject: [mg49909] Re: Asymmetric quadratic coefficient problem
  • From: "Dr. Wolfgang Hintze" <weh at>
  • Date: Thu, 5 Aug 2004 09:22:36 -0400 (EDT)
  • References: <cequfb$k3m$>
  • Sender: owner-wri-mathgroup at

As to your first question:
You can easily convice yourself that the asymmetric part of the matrix A 
just drops out of the term z'Az.


S=(A+Transpose[A])/2  (symmetric part of A)
N=(A-Transpose[A])/2  (asymmetric part of A)

then, obviously

A = S+N

and we have

z'Az = z'Sz + z'Nz

but the term x=z'Nz vanishes because (not in strict Mathematica style)

x=z'Nz = Sum[ z[i] N[i,k] z[k] ]
= Sum[ z[k] N[k,i] z[i] ] (just renaming summation indices)
= - Sum [ z[i] N[i,k] z[k] ] (because of N[k,i] = - N[i,k] )
= - x

Hences x=0.

Hope this helps.


Xiao Huang wrote:

> Dear All,
> Does any one know how to deal with asymmetric quadratic coefficient in 
> QuadraticFormDistribution[ {A,b,c}, {mu, sigma}] for z'Az+b'z+c?
> Mathematica assumes A is a symmetric matrix, but I now have an asymmetric 
> and singular matrix A.
> Also, I'd really appreciate if any expert here can give me some 
> clue/reference about the following problem: how can I calculate different 
> quantiles of a mixture distribution? One way is to do the simulation, but is 
> it possible to use the function Quantile[.] in Mathematica? if yes, how?
> Thank you!
> Xiao
> _________________________________________________________________
> Overwhelmed by debt? Find out how to 'Dig Yourself Out of Debt' from MSN 
> Money.

  • Prev by Date: Re: populate a list with random numbers from normaldistribution?
  • Next by Date: Re: Matrix + Simplify
  • Previous by thread: Re: Asymmetric quadratic coefficient problem
  • Next by thread: 'NonlinearFit` confusion