Re: Eigensystem
- To: mathgroup at smc.vnet.net
- Subject: [mg39664] Re: Eigensystem
- From: Sseziwa Mukasa <mukasa at jeol.com>
- Date: Sat, 1 Mar 2003 02:47:15 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
On Friday, February 28, 2003, at 04:46 AM, Angela Birk wrote:
> Dear MathGroup Members,
>
> I would like to generate a 12x12 matrix for that I can caluclate the
> Eigensystem by transposing a 12x12 matrix. This 12x12 matrix is the
> result of a model that includes a bunch of functions and variables.
> (Note the attached file doesn't work, because the model is missing. But
> one can see that the mistake occurs.)
> [contact the author to get the attachment - moderator]
>
> I try to derive the 12x12 matrix by using different commands:
>
> First, I define the first row of the matrix (see attached file):
>
> f[a_,b_] := f[a,b]=Table[(1-CDF[h[a,bValue,W]])*
> SETP[bValue,b[[i]]]*ATP[a, j] + CDF[h[a,bValue,W]]*
> SETP[bValue,b[[i]]]*ATP[a, j], {j, 1, 4}, {i, 1,3}]
>
> where CDF[h[a,bValue,W]] is a probability that depends on the function
> h
> and SETP[bValue,b[[i]]] as well as ATP[a, j] are given matrices having
> the dimensions 3x3 and 4x4, respectively.
>
> Second, if I execute the command
>
> f[1, b[[1]]] // MatrixForm
>
> this should give a 1x12 vector. But actually I get a 4x3 matrix (see
> attached file). This 1x12 vector is the first row of the 12x12 matrix,
> and the sum of row is one. (Each row has to sum to one.)
>
> Could someone suggest how the above command(s) should be changed to get
> the 1x12 vector that sums to one?
>
>
f[a_,b_]:=f[a,b]=Flatten[Table[(1-CDF[h[a,bValue,W]])*
SETP[bValue,b[[i]]]*ATP[a, j] + CDF[h[a,bValue,W]]*
SETP[bValue,b[[i]]]*ATP[a, j], {j, 1, 4}, {i, 1,3}]]
Regards,
Sseziwa