How to find eigenvalue of a matrix (attached) having three variables?

조회 수: 3 (최근 30일)
Ismaeel
Ismaeel 2019년 6월 18일
댓글: Ismaeel 2019년 7월 26일
I have a 33X33 matrix with three variables (KA1, KA2, KA3). I want to take its eigenvalues. The eigenvalue should be a function of these three variables. Any idea? Thanks.
A =
[ -0.1754, -0.01261, 0.005666, 0.03756, 0.02875, 0.01006, -0.04169, -0.002683, 0.000833, -0.02971, -0.02183, -0.01673, 0.004094, -0.002681, -0.001413, 0, 0, 0, 0.1116, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -0.2057, -13.31, -0.0657, 0.03443, 13.41, -0.1167, 0.5673, 0.00206, -0.000514, 0.4043, 0.01676, 0.01032, -0.1771, 0.2988, -0.1217, 0, 0, 0, 0, 0.1667, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0.1378, -0.09722, -1.388, -0.02306, 0.2217, 1.044, -0.3221, -0.004694, -1.429e-5, -0.2296, -0.03819, 0.0002871, 0.1011, 0.06514, -0.1663, 0, 0, 0, 0, 0, 0.1698, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 10.12, 0.5683, -0.2554, -8.933, -1.296, -0.4534, 1.879, 0.121, -0.03755, 1.339, 0.984, 0.7541, -0.1846, 0.1209, 0.06369, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 1.076, 30.76, 0.3436, -0.1801, -32.16, 0.6101, -2.967, -0.01077, 0.002688, -2.114, -0.08764, -0.05399, 0.9263, -1.563, 0.6364, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 5.039, -3.556, 21.52, -0.8434, 8.11, -27.88, -11.78, -0.1717, -0.0005228, -8.397, -1.397, 0.0105, 3.699, 2.383, -6.082, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, -3.226, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -0.3163, 0.07269, -0.02432, 0.05294, -0.1658, -0.04318, -0.3983, -16.94, -0.01262, -0.2839, -15.62, 0.2534, 0.112, -0.2546, 0.1426, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0.1245, -0.001107, 0.01228, -0.02083, 0.002524, 0.0218, 0.2527, -0.01684, -10.82, 0.1801, -0.137, -8.9, -0.0731, -0.08933, 0.1624, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -0.3569, 0.9731, -0.3902, 0.05974, -2.219, -0.6928, -38.3, -0.03862, 0.02101, -36.0, -0.3142, -0.4219, 2.197, -1.111, -1.086, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -2.2, 0.5057, -0.1692, 0.3683, -1.153, -0.3004, -2.771, -13.36, -0.08779, -1.975, -17.17, 1.763, 0.7792, -1.771, 0.992, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -1.458, 0.01296, -0.1438, 0.244, -0.02957, -0.2554, -2.961, 0.1973, -5.191, -2.11, 1.605, -8.184, 0.8565, 1.047, -1.903, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -2.86, -11.95, 4.807, 0.4787, 27.26, 8.535, 86.55, 0.4028, -0.2332, 61.68, 3.277, 4.683, -26.26, 13.27, 12.99, -0.1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7.974, 0, 0, 0, 0, 0]
[ 3.478, 74.73, 10.82, -0.5822, -170.4, 19.21, -170.5, -4.851, -1.489, -121.5, -39.46, 29.91, 51.65, -95.45, 43.8, 0, -0.2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 29.45, 0, 0, 0, 0]
[ 18.45, -61.83, -61.2, -3.087, 141.0, -108.7, -349.2, 6.465, 5.145, -248.9, 52.59, -103.3, 106.1, 92.13, -198.2, 0, 0, -0.3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 62.62, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -3.365, 0, 0, 0, 0, 0, 3.185, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -3.947, 0, 0, 0, 0, 0, 3.185, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -3.536, 0, 0, 0, 0, 0, 3.185, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5143, 0, 0, -2.857, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5143, 0, 0, -2.857, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.5143, 0, 0, -2.857, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -5.322*KA1, 0.2824*KA1, -0.1263*KA1, 0.8907*KA1, -0.6439*KA1, -0.2243*KA1, -0.2677*KA1, 0.05718*KA1, -0.01764*KA1, -0.1908*KA1, 0.4652*KA1, 0.3543*KA1, -0.06187*KA1, 0.04441*KA1, 0.01746*KA1, 0, 0, 0, -0.9*KA1, 0, 0, 5.0*KA1, 0, 0, -5.0, 0, 0, 0, 0, 0, 0, 0, 0]
[ -1.192*KA2, 1.977*KA2, -0.2315*KA2, 0.1995*KA2, -4.509*KA2, -0.411*KA2, 0.816*KA2, 0.4063*KA2, -0.02598*KA2, 0.5816*KA2, 3.305*KA2, 0.5217*KA2, -0.2788*KA2, 0.3427*KA2, -0.06391*KA2, 0, 0, 0, 0, -0.9*KA2, 0, 0, 5.0*KA2, 0, 0, -5.0, 0, 0, 0, 0, 0, 0, 0]
[ -1.237*KA3, 0.541*KA3, -0.8733*KA3, 0.207*KA3, -1.234*KA3, -1.55*KA3, 0.8115*KA3, 0.09036*KA3, -0.1336*KA3, 0.5783*KA3, 0.7351*KA3, 2.683*KA3, -0.2786*KA3, -0.01786*KA3, 0.2965*KA3, 0, 0, 0, 0, 0, -0.9*KA3, 0, 0, 5.0*KA3, 0, 0, -5.0, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -10.0, 0, 0, 10.0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -10.0, 0, 0, 10.0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -10.0, 0, 0, 10.0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1.061, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -20.0, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1.061, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -20.0, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1.061, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -20.0]

답변 (1개)

johnson wul
johnson wul 2019년 7월 26일
If your matrix is not symmetric, you should change it into a symmetric one by using the formular:
symmetric_Matrix(i,j) = ur_(matrix(i,j)+ur_matrix(i,j)')/2
After doing this just use eig() function to obtain eigenvalues

카테고리

Help CenterFile Exchange에서 Linear Algebra에 대해 자세히 알아보기

태그

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by