Eigensystem
- To: mathgroup at smc.vnet.net
- Subject: [mg39671] Eigensystem
- From: Angela Birk <angela.birk at hwwa.de>
- Date: Fri, 28 Feb 2003 04:46:32 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
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? Third, after having the vector, I'd like to generate the 12x12 matrix with the following command Table[f[a,b[[k]]], {a, 1, 4}, {k, 1, 3}] // MatrixForm But actually this gives a 4x3x4x3 matrix (see attached file). Fourth, if I transpose the above matrix and try to calculate the Eigensystem, I get the failure message (see attached file). I think, the mistake has something to do with the Table command. Does somebody know, how I have to change the commands to get the desired results? Thank's in advance for the help. Angela -- Angela Birk University of Hamburg HWWA 20347 Hamburg, Germany +49 40 42834 - 478