login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit

(the maths inside is not meant to be taken seriously; 'tis a silly idea that can't work)

Note that dropping comutativity below we get quantum probablility. Franz Lehner added packages supporting computations for several popular versions. In principle one can treat commutative case as a special case of quantum probablility...

from a recent email by Peter Broadbery

Random variables are assumed to have the following properties:

  1. complex constants are random variables;
  2. the sum of two random variables is a random variable;
  3. the product of two random variables is a random variable;
  4. addition and multiplication of random variables are both commutative; and
  5. there is a notion of conjugation of random variables, satisfying:

    (ab)<em> = b</em> a<em> and a** = a

    for all random variables a, b, and coinciding with complex conjugation if a is a constant.

This means that random variables form complex abelian *-algebras. If a = a*, the random variable a is called "real".

An expectation E on an algebra A of random variables is a normalized, positive linear functional. What this means is that

  1. E(1) = 1;
  2. E(a <em> a) >= 0 for all random variables a;
  3. E(a + b) = E(a) + E(b) for all random variables a and b; and
  4. E(za) = zE(a) if z is a constant.

*-algebra

From Wikipedia, the free encyclopedia

In mathematics, a *-algebra is an associative algebra over the field of complex numbers with an antilinear, antiautomorphism <em> : A \rightarrow A which is an involution. More precisely, * is required to satisfy the following properties:

  • (x + y)^<em> = x^</em> + y^<em> \quad
  • (z x)^<em> = \overline{z} x^</em>
  • (x y)^<em> = y^</em> x^<em> \quad
  • (x^<em>)^</em> = x \quad

for all x,y in A, and all z in C.

The most obvious example of a *-algebra is the field of complex numbers C where * is just complex conjugation. Another example is the algebra of nn matrices over C with * given by the conjugate transpose.

An algebra homomorphism f : A \rightarrow B is a *-homomorphism if it is compatible with the involutions of A and B, i.e.

  • f(a <em> ) = f(a) </em> for all a in A.

An element a in A is called self-adjoint if a* = a.

aldor
#include "axiom"
RandomAlgebra(F: Field): Category == with { Algebra F; E: % -> F; sample: % -> F; }
local PolyHelper(F: Field): with { expand: SparseUnivariatePolynomial F -> Generator Cross(F, NonNegativeInteger); } == add { expand(p: SparseUnivariatePolynomial F): Generator Cross(F, NonNegativeInteger) == generate { default m: SparseUnivariatePolynomial F; import from SparseUnivariatePolynomial F; import from List SparseUnivariatePolynomial F; for m in monomials p repeat { yield (leadingCoefficient m, degree m); } } }
UnivariateNormalRandomAlgebra: RandomAlgebra Float with { X: () -> %; variance: % -> Float; } == add { Rep ==> SparseUnivariatePolynomial Float; import from Rep;
0: % == per 0; 1: % == per 1; X(): % == per(monomial(1$Float,1$NonNegativeInteger)$Rep);
characteristic(): NonNegativeInteger == 0;
-(x: %): % == per(-rep x); (a: %) = (b: %): Boolean == rep(a) = rep(b);
(a: %) + (b: %): % == per(rep(a) + rep(b)); (a: %) * (b: %): % == per(rep(a) * rep(b)); (a: Float) * (b: %): % == per(a * rep(b));
coerce(x: Integer): % == per(x::Rep); coerce(x: Float): % == per(x::Rep);
coerce(x: %): OutputForm == coerce rep(x);
E(X: %): Float == { import from PolyHelper Float; z: Float := 0; for p in expand rep(X) repeat { (a, b) := p; z := z + a * E(b); } z }
-- should be a random sampling of x. sample(X: %): Float == { import from PolyHelper Float; import from Float; u := uniform01()$RandomFloatDistributions; x: Float := 0; for p in expand rep(X) repeat { (a, b) := p; x := x + a * u^b; } return x; }
variance(X: %): Float == { A := (X-E(X)*1); E(A*A); }
-- return expected value of X^n local E(n: NonNegativeInteger): Float == { p: Rep := 1; -- yuck. There must be a nicer way than this.. for i in 1..n repeat p := differentiate(p) + monomial(1,1)*p; coefficient(p,0); } }
aldor
   Compiling FriCAS source code from file 
      /var/lib/zope2.10/instance/axiom-wiki/var/LatexWiki/2618023199496535153-25px001.as
      using Aldor compiler and options 
-O -Fasy -Fao -Flsp -lfricas -Mno-ALDOR_W_WillObsolete -DFriCAS -Y $FRICAS/algebra -I $FRICAS/algebra
      Use the system command )set compiler args to change these 
      options.
   The )library system command was not called after compilation.

fricas
a := X()$UnivariateNormalRandomAlgebra
UnivariateNormalRandomAlgebra is not a valid type.




  Subject:   Be Bold !!
  ( 15 subscribers )  
Please rate this page: