Форум программистов, компьютерный форум, киберфорум
Наши страницы
Matlab
Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 10, средняя оценка - 4.90
HMB
0 / 0 / 0
Регистрация: 25.02.2013
Сообщений: 17
#1

Какая функция в matlab является аналогичной функции eigenvec(M, z) в mathcad? - Matlab

25.02.2013, 17:57. Просмотров 1465. Ответов 7
Метки нет (Все метки)

В mathcad я использую функцию eigenvec(M, z) для вывода из матрицы M, нормированный собственный вектор, соответствующий собственному значению z. Но мне нужно то же самое выполнить в Matlab. Как это сделать? Я попробовал функции eig(M) и eigs(M) но результат не тот. Пожалуйста кто знает помогите!!! Очень нужно.
http://www.cyberforum.ru/matlab/thread1955126.html
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
25.02.2013, 17:57
Я подобрал для вас темы с готовыми решениями и ответами на вопрос Какая функция в matlab является аналогичной функции eigenvec(M, z) в mathcad? (Matlab):

Из MathCad в MatLab
Здравствуйте! Пытаюсь функцию MathCad переделать в Matlab,но что-то не...

Из mathcad в Matlab
Здравствуйте, ооочень нужна помощь, есть программка в маткаде с функцией...

Из mathCAD в matlab
Здравствуйте! Есть вот такое выражение в маткаде: Пытаюсь сделать в...

Портирование из Mathcad в Matlab
Ребят, кто-нибудь может, пожалуйста, конвертировать программу(во вложении) в...

Интерполяция (MathCad to Matlab)
Здравствуйте, Уважаемые форумчане! Подскажите, пожалуйста, как данная запись...

7
Зосима
4915 / 3284 / 308
Регистрация: 02.04.2012
Сообщений: 6,191
Записей в блоге: 15
Завершенные тесты: 1
25.02.2013, 18:25 #2
А можешь показать в числах, что должно получиться?
Ф-ция eig может возвращать и собственные значения и собственный вектор:
Matlab M
1
2
M = magic(3) % некая матрица
[V,D] = eig(M)
Результат:
Код
M =
     8     1     6
     3     5     7
     4     9     2

V =
   -0.5774   -0.8131   -0.3416
   -0.5774    0.4714   -0.4714
   -0.5774    0.3416    0.8131

D =
   15.0000         0         0
         0    4.8990         0
         0         0   -4.8990
Единственное, вектор собственных значений не нормированный но это нетрудно исправить:
Matlab M
1
Z = diag(D)/norm(diag(D))
Также можно подшаманить и собственные вектора
0
HMB
0 / 0 / 0
Регистрация: 25.02.2013
Сообщений: 17
26.02.2013, 09:56  [ТС] #3
Спасибо за ответ но это я пробовал результат не тот. В matlab функция eig(KM) используется для вывода собственных значений и собственных векторов матрицы KM. А функция eigenvec(KM,s) в mathcad выводит собственный вектор матрицы KM по соответствующему собственному значению s.
Исходные данные и результаты полученные в Mathcad представлены в файле Пример.
0
Вложения
Тип файла: doc Пример.doc (57.5 Кб, 5 просмотров)
Зосима
4915 / 3284 / 308
Регистрация: 02.04.2012
Сообщений: 6,191
Записей в блоге: 15
Завершенные тесты: 1
26.02.2013, 13:01 #4
Вообщем гляди
Вначале вводим матрицу KM
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
clear all
clc
 
M1 = [1.2286    1.1867    1.1556    1.1196    1.0819    1.0500    1.0471
    1.1867    1.2039    1.1699    1.1363    1.0990    1.0666    1.0425
    1.1556    1.1699    1.1941    1.1570    1.1219    1.0891    1.0632
    1.1196    1.1363    1.1570    1.1789    1.1408    1.1100    1.0835
    1.0819    1.0990    1.1219    1.1408    1.1619    1.1279    1.1033
    1.0500    1.0666    1.0891    1.1100    1.1279    1.1525    1.1238
    1.0471    1.0425    1.0632    1.0835    1.1033    1.1238    1.1523
    1.0522    1.0410    1.0409    1.0589    1.0778    1.1003    1.1248
    1.0655    1.0377    1.0313    1.0295    1.0467    1.0689    1.0961
    1.0944    1.0632    1.0400    1.0302    1.0266    1.0465    1.0727
    1.1297    1.0906    1.0634    1.0375    1.0263    1.0251    1.0486
    1.1494    1.1247    1.0904    1.0601    1.0326    1.0242    1.0272
    1.1592    1.1444    1.1260    1.0874    1.0548    1.0310    1.0281
    1.1522    1.1499    1.1414    1.1193    1.0788    1.0501    1.0320
    1.1380    1.1526    1.1563    1.1428    1.1182    1.0810    1.0573
    1.1123    1.1375    1.1568    1.1568    1.1414    1.1193    1.0860
    1.0843    1.1053    1.1342    1.1515    1.1506    1.1375    1.1188
    1.0633    1.0817    1.1065    1.1327    1.1487    1.1499    1.1401
    1.0528    1.0521    1.0739    1.0976    1.1233    1.1417    1.1464
    1.0556    1.0478    1.0504    1.0704    1.0931    1.1208    1.1421
    1.0766    1.0490    1.0445    1.0455    1.0646    1.0894    1.1201
    1.1137    1.0815    1.0578    1.0493    1.0482    1.0694    1.0971
    1.1610    1.1226    1.0955    1.0662    1.0548    1.0562    1.0811
    1.2053    1.1618    1.1291    1.0971    1.0654    1.0571    1.0633
    1.2159    1.1986    1.1614    1.1245    1.0905    1.0625    1.0600
    1.2013    1.2070    1.1960    1.1550    1.1162    1.0860    1.0638
    1.1850    1.1942    1.2061    1.1911    1.1482    1.1130    1.0883
    1.1368    1.1678    1.1830    1.1927    1.1767    1.1377    1.1083
    1.0974    1.1265    1.1617    1.1750    1.1838    1.1706    1.1355
    1.0728    1.0896    1.1227    1.1558    1.1681    1.1795    1.1699
    1.0740    1.0709    1.0904    1.1215    1.1537    1.1675    1.1813
    1.0875    1.0689    1.0691    1.0867    1.1169    1.1510    1.1680
    1.1128    1.0826    1.0676    1.0657    1.0821    1.1144    1.1519
    1.1568    1.1172    1.0913    1.0720    1.0679    1.0866    1.1223
    1.1993    1.1591    1.1241    1.0941    1.0727    1.0710    1.0934
    1.2290    1.1975    1.1629    1.1235    1.0914    1.0731    1.0762];
M2 = [1.0522    1.0655    1.0944    1.1297    1.1494    1.1592    1.1522
    1.0410    1.0377    1.0632    1.0906    1.1247    1.1444    1.1499
    1.0409    1.0313    1.0400    1.0634    1.0904    1.1260    1.1414
    1.0589    1.0295    1.0302    1.0375    1.0601    1.0874    1.1193
    1.0778    1.0467    1.0266    1.0263    1.0326    1.0548    1.0788
    1.1003    1.0689    1.0465    1.0251    1.0242    1.0310    1.0501
    1.1248    1.0961    1.0727    1.0486    1.0272    1.0281    1.0320
    1.1548    1.1215    1.1018    1.0763    1.0525    1.0336    1.0312
    1.1215    1.1462    1.1206    1.0999    1.0737    1.0501    1.0285
    1.1018    1.1206    1.1559    1.1274    1.1075    1.0851    1.0577
    1.0763    1.0999    1.1274    1.1604    1.1324    1.1153    1.0897
    1.0525    1.0737    1.1075    1.1324    1.1660    1.1411    1.1203
    1.0336    1.0501    1.0851    1.1153    1.1411    1.1792    1.1494
    1.0312    1.0285    1.0577    1.0897    1.1203    1.1494    1.1830
    1.0420    1.0322    1.0442    1.0690    1.1025    1.1390    1.1628
    1.0645    1.0423    1.0436    1.0523    1.0782    1.1161    1.1486
    1.0867    1.0605    1.0453    1.0452    1.0539    1.0812    1.1167
    1.1229    1.0855    1.0678    1.0503    1.0507    1.0623    1.0866
    1.1371    1.1163    1.0842    1.0658    1.0478    1.0482    1.0578
    1.1478    1.1344    1.1200    1.0864    1.0683    1.0518    1.0499
    1.1422    1.1441    1.1366    1.1210    1.0874    1.0703    1.0517
    1.1301    1.1458    1.1583    1.1472    1.1332    1.1048    1.0839
    1.1121    1.1363    1.1666    1.1739    1.1653    1.1588    1.1251
    1.0910    1.1132    1.1515    1.1775    1.1864    1.1836    1.1720
    1.0685    1.0875    1.1233    1.1582    1.1849    1.1980    1.1902
    1.0634    1.0637    1.0955    1.1283    1.1636    1.1937    1.2022
    1.0683    1.0597    1.0728    1.1015    1.1349    1.1739    1.1994
    1.0848    1.0582    1.0592    1.0709    1.0989    1.1327    1.1684
    1.1072    1.0788    1.0595    1.0592    1.0707    1.0994    1.1307
    1.1361    1.1027    1.0819    1.0610    1.0608    1.0737    1.0998
    1.1728    1.1352    1.1078    1.0854    1.0650    1.0667    1.0775
    1.1828    1.1700    1.1391    1.1102    1.0880    1.0692    1.0685
    1.1700    1.1802    1.1746    1.1420    1.1134    1.0931    1.0717
    1.1622    1.1733    1.1951    1.1857    1.1547    1.1316    1.1071
    1.1315    1.1642    1.1869    1.2050    1.1971    1.1712    1.1439
    1.1010    1.1310    1.1766    1.1956    1.2150    1.2118    1.1813];
 
 
M3 =[1.1380    1.1123    1.0843    1.0633    1.0528    1.0556    1.0766
    1.1526    1.1375    1.1053    1.0817    1.0521    1.0478    1.0490
    1.1563    1.1568    1.1342    1.1065    1.0739    1.0504    1.0445
    1.1428    1.1568    1.1515    1.1327    1.0976    1.0704    1.0455
    1.1182    1.1414    1.1506    1.1487    1.1233    1.0931    1.0646
    1.0810    1.1193    1.1375    1.1499    1.1417    1.1208    1.0894
    1.0573    1.0860    1.1188    1.1401    1.1464    1.1421    1.1201
    1.0420    1.0645    1.0867    1.1229    1.1371    1.1478    1.1422
    1.0322    1.0423    1.0605    1.0855    1.1163    1.1344    1.1441
    1.0442    1.0436    1.0453    1.0678    1.0842    1.1200    1.1366
    1.0690    1.0523    1.0452    1.0503    1.0658    1.0864    1.1210
    1.1025    1.0782    1.0539    1.0507    1.0478    1.0683    1.0874
    1.1390    1.1161    1.0812    1.0623    1.0482    1.0518    1.0703
    1.1628    1.1486    1.1167    1.0866    1.0578    1.0499    1.0517
    1.2075    1.1809    1.1545    1.1284    1.0864    1.0645    1.0544
    1.1809    1.2207    1.1850    1.1631    1.1279    1.0911    1.0676
    1.1545    1.1850    1.2201    1.1870    1.1598    1.1281    1.0903
    1.1284    1.1631    1.1870    1.2254    1.1856    1.1625    1.1295
    1.0864    1.1279    1.1598    1.1856    1.2203    1.1833    1.1594
    1.0645    1.0911    1.1281    1.1625    1.1833    1.2211    1.1832
    1.0544    1.0676    1.0903    1.1295    1.1594    1.1832    1.2201
    1.0732    1.0703    1.0741    1.1012    1.1314    1.1662    1.1883
    1.1150    1.0964    1.0803    1.0900    1.1047    1.1416    1.1742
    1.1489    1.1325    1.1022    1.0917    1.0900    1.1112    1.1460
    1.1891    1.1612    1.1346    1.1096    1.0884    1.0931    1.1125
    1.2043    1.1991    1.1620    1.1403    1.1053    1.0903    1.0933
    1.2177    1.2155    1.2007    1.1686    1.1368    1.1080    1.0912
    1.2014    1.2186    1.2110    1.1998    1.1608    1.1339    1.1038
    1.1720    1.2038    1.2165    1.2116    1.1950    1.1597    1.1318
    1.1370    1.1764    1.2030    1.2186    1.2079    1.1951    1.1587
    1.1082    1.1431    1.1778    1.2068    1.2176    1.2098    1.1960
    1.0845    1.1132    1.1432    1.1806    1.2044    1.2184    1.2097
    1.0766    1.0903    1.1137    1.1466    1.1782    1.2055    1.2186
    1.0944    1.0934    1.0969    1.1250    1.1482    1.1852    1.2108
    1.1281    1.1099    1.0990    1.1072    1.1257    1.1543    1.1896
    1.1637    1.1426    1.1141    1.1084    1.1061    1.1307    1.1573];
M4 = [1.1137    1.1610    1.2053    1.2159    1.2013    1.1850    1.1368
    1.0815    1.1226    1.1618    1.1986    1.2070    1.1942    1.1678
    1.0578    1.0955    1.1291    1.1614    1.1960    1.2061    1.1830
    1.0493    1.0662    1.0971    1.1245    1.1550    1.1911    1.1927
    1.0482    1.0548    1.0654    1.0905    1.1162    1.1482    1.1767
    1.0694    1.0562    1.0571    1.0625    1.0860    1.1130    1.1377
    1.0971    1.0811    1.0633    1.0600    1.0638    1.0883    1.1083
    1.1301    1.1121    1.0910    1.0685    1.0634    1.0683    1.0848
    1.1458    1.1363    1.1132    1.0875    1.0637    1.0597    1.0582
    1.1583    1.1666    1.1515    1.1233    1.0955    1.0728    1.0592
    1.1472    1.1739    1.1775    1.1582    1.1283    1.1015    1.0709
    1.1332    1.1653    1.1864    1.1849    1.1636    1.1349    1.0989
    1.1048    1.1588    1.1836    1.1980    1.1937    1.1739    1.1327
    1.0839    1.1251    1.1720    1.1902    1.2022    1.1994    1.1684
    1.0732    1.1150    1.1489    1.1891    1.2043    1.2177    1.2014
    1.0703    1.0964    1.1325    1.1612    1.1991    1.2155    1.2186
    1.0741    1.0803    1.1022    1.1346    1.1620    1.2007    1.2110
    1.1012    1.0900    1.0917    1.1096    1.1403    1.1686    1.1998
    1.1314    1.1047    1.0900    1.0884    1.1053    1.1368    1.1608
    1.1662    1.1416    1.1112    1.0931    1.0903    1.1080    1.1339
    1.1883    1.1742    1.1460    1.1125    1.0933    1.0912    1.1038
    1.2381    1.2141    1.1950    1.1623    1.1265    1.1083    1.0961
    1.2141    1.2749    1.2443    1.2190    1.1833    1.1488    1.1166
    1.1950    1.2443    1.2980    1.2609    1.2328    1.1985    1.1510
    1.1623    1.2190    1.2609    1.3079    1.2681    1.2415    1.1949
    1.1265    1.1833    1.2328    1.2681    1.3125    1.2743    1.2363
    1.1083    1.1488    1.1985    1.2415    1.2743    1.3202    1.2704
    1.0961    1.1166    1.1510    1.1949    1.2363    1.2704    1.3084
    1.1089    1.1043    1.1203    1.1505    1.1932    1.2355    1.2635
    1.1387    1.1195    1.1105    1.1224    1.1512    1.1948    1.2305
    1.1664    1.1500    1.1277    1.1157    1.1264    1.1558    1.1940
    1.2032    1.1771    1.1569    1.1311    1.1178    1.1293    1.1528
    1.2177    1.2151    1.1850    1.1611    1.1338    1.1214    1.1265
    1.2377    1.2452    1.2371    1.2019    1.1756    1.1494    1.1260
    1.2288    1.2638    1.2655    1.2523    1.2147    1.1895    1.1525
    1.2074    1.2548    1.2831    1.2787    1.2629    1.2267    1.1896];
 
M5 = [1.0974    1.0728    1.0740    1.0875    1.1128    1.1568    1.1993
    1.1265    1.0896    1.0709    1.0689    1.0826    1.1172    1.1591
    1.1617    1.1227    1.0904    1.0691    1.0676    1.0913    1.1241
    1.1750    1.1558    1.1215    1.0867    1.0657    1.0720    1.0941
    1.1838    1.1681    1.1537    1.1169    1.0821    1.0679    1.0727
    1.1706    1.1795    1.1675    1.1510    1.1144    1.0866    1.0710
    1.1355    1.1699    1.1813    1.1680    1.1519    1.1223    1.0934
    1.1072    1.1361    1.1728    1.1828    1.1700    1.1622    1.1315
    1.0788    1.1027    1.1352    1.1700    1.1802    1.1733    1.1642
    1.0595    1.0819    1.1078    1.1391    1.1746    1.1951    1.1869
    1.0592    1.0610    1.0854    1.1102    1.1420    1.1857    1.2050
    1.0707    1.0608    1.0650    1.0880    1.1134    1.1547    1.1971
    1.0994    1.0737    1.0667    1.0692    1.0931    1.1316    1.1712
    1.1307    1.0998    1.0775    1.0685    1.0717    1.1071    1.1439
    1.1720    1.1370    1.1082    1.0845    1.0766    1.0944    1.1281
    1.2038    1.1764    1.1431    1.1132    1.0903    1.0934    1.1099
    1.2165    1.2030    1.1778    1.1432    1.1137    1.0969    1.0990
    1.2116    1.2186    1.2068    1.1806    1.1466    1.1250    1.1072
    1.1950    1.2079    1.2176    1.2044    1.1782    1.1482    1.1257
    1.1597    1.1951    1.2098    1.2184    1.2055    1.1852    1.1543
    1.1318    1.1587    1.1960    1.2097    1.2186    1.2108    1.1896
    1.1089    1.1387    1.1664    1.2032    1.2177    1.2377    1.2288
    1.1043    1.1195    1.1500    1.1771    1.2151    1.2452    1.2638
    1.1203    1.1105    1.1277    1.1569    1.1850    1.2371    1.2655
    1.1505    1.1224    1.1157    1.1311    1.1611    1.2019    1.2523
    1.1932    1.1512    1.1264    1.1178    1.1338    1.1756    1.2147
    1.2355    1.1948    1.1558    1.1293    1.1214    1.1494    1.1895
    1.2635    1.2305    1.1940    1.1528    1.1265    1.1260    1.1525
    1.3045    1.2611    1.2309    1.1928    1.1518    1.1316    1.1300
    1.2611    1.3036    1.2625    1.2310    1.1933    1.1589    1.1375
    1.2309    1.2625    1.3063    1.2644    1.2333    1.2013    1.1662
    1.1928    1.2310    1.2644    1.3072    1.2657    1.2406    1.2078
    1.1518    1.1933    1.2333    1.2657    1.3088    1.2740    1.2480
    1.1316    1.1589    1.2013    1.2406    1.2740    1.3292    1.2931
    1.1300    1.1375    1.1662    1.2078    1.2480    1.2931    1.3470
    1.1546    1.1344    1.1441    1.1715    1.2139    1.2667    1.3103];
M6 = [1.2290
    1.1975
    1.1629
    1.1235
    1.0914
    1.0731
    1.0762
    1.1010
    1.1310
    1.1766
    1.1956
    1.2150
    1.2118
    1.1813
    1.1637
    1.1426
    1.1141
    1.1084
    1.1061
    1.1307
    1.1573
    1.2074
    1.2548
    1.2831
    1.2787
    1.2629
    1.2267
    1.1896
    1.1546
    1.1344
    1.1441
    1.1715
    1.2139
    1.2667
    1.3103
    1.3629];
 
KM = 1.0e+12*[M1 M2 M3 M4 M5 M6];

Ну а дальше считаем собственные вектора и значения, и находим подходящие:
Matlab M
1
2
3
4
5
6
7
8
9
10
% заданное собственное значение
s = 40822611056652.8; 
 % находим собственные вектора и значения
[V, D] = eig(KM);
 % вектор собственных значений
Z = diag(D);
% находим собственный вектор, которы сообтветствует 
% собственному значению, наиболее близкому к заданному
q = abs(Z-s);
U = V(find(q==min(q)),:)
Единственное, тут я нахожу строку в матрице векторов V, но я не силен в матричной теории, поэтому если нужен столбец, то двоеточие переместиться:
Matlab M
1
U = V(:,find(q==min(q)))
*норма полученного вектора norm(U) получается равной 1. Для подстраховки можно добавить строку:
Matlab M
1
U = U/norm(U)
тогда U гарантированно будет нормированным
0
HMB
0 / 0 / 0
Регистрация: 25.02.2013
Сообщений: 17
26.02.2013, 14:12  [ТС] #5
Спасибо Зосима за ответ. Это моя не внимательность ведь eig выдает сразу собственные числа и собственные векторы и из этих СЧ я мог сопоставить и выделить нужный СВ. Кстати указанное в Примере СЧ является 36-м. Но тут уже возникла другая проблема почему знаки разные. Например для (36 СЧ матрицы КМ) при использовании функции mathcad результаты с минусом а при eig они с плюсом. А для (35 СЧ) по обеим функциям совпадают. Я проверил по 6 последним СЧ (31-36) в двух случаях (32,35) знаки совпадают а в остальных случаях значения одинаковые а знаки разные. Интересно в чем же причина. Я попробовал найти ответ в нете, на других форумах но толком ничего не нашел.
0
vital792
1997 / 1269 / 60
Регистрация: 05.06.2010
Сообщений: 2,213
04.03.2013, 17:57 #6
Цитата Сообщение от HMB Посмотреть сообщение
Интересно в чем же причина.
очевидно в разных алгоритмах. Вообще, в задачах сингулярного анализа, главных компонент и других задачах где используют разложение на собственные вектора, знак не важен, важно чтобы эти вектора образовывали полный ортонормированный базис. А что вектор указывает в противоположную сторону - это обычно не принципиально, если совпадают прямые на которых лежат эти вектора
0
HMB
0 / 0 / 0
Регистрация: 25.02.2013
Сообщений: 17
05.03.2013, 10:41  [ТС] #7
Спасибо за ответ. Теперь понятно.
0
vital792
1997 / 1269 / 60
Регистрация: 05.06.2010
Сообщений: 2,213
05.03.2013, 13:27 #8
это кстати можно показать на простеньком примере. Возьмем симметричную матрицу 2х2:

Matlab M
1
2
a = [5.5 -7;
     -7  3.2];
Для разложения на собственные числа, ее надо привести к диагональной форме. Для этого проще всего выполнить преобразование подобия T*A*T' (разложение Шура), где Т - ортогональная матрица. В качестве этой матрицы можно использовать плоское вращение
Код
[cos(phi) -sin(phi);
 sin(phi)   cos(phi)];
Угол phi надо выбрать такой, чтобы обнулить внедиагональные элементы, т.е. как бы посмотреть на пространство, натянутое на вектора матрицы А под другим углом - под которым концы векторов лежат на осях. Сам угол напрямую вычислять не обязательно, можно сразу найти синус и косинус. Для этого можно используя простые тригонометрические тождества решить квадратное уравнение и выбрать меньший угол(метод Рутисхаузера):
Matlab M
1
2
3
4
5
6
7
8
ctg = 0.5*(a(2,2)-a(1,1))/a(1,2);
t = sign(ctg) / (abs(ctg) + sqrt(1 + ctg^2));
c = 1/sqrt(1+t^2);
s = t*c;
 
rot = [c -s;
       s c];
rot * a * rot'
Полученные вектора [sin cos] образуют полный ортонормированный базис в пространстве R(2x2) и они точно являются собственными векторами матрицы А. Но сравним результаты с функцией eig():

Matlab M
1
2
3
4
[v,d] = eig(a)
 
[a*v(:,1) a*[s; c]]
[a*v(:,2) a*[-c; s]]
Имеем
Matlab M
1
2
3
4
5
6
7
8
9
ans =
 
   1.775971052399038  -1.775971052399040
   2.091544706533572  -2.091544706533570
 
ans =
 
  -8.723297932607865  -8.723297932607865
   7.407120948153297   7.407120948153297
Видно, что первый выбранный вектор противоположного знака, но это абсолютно не принципиально...
2
05.03.2013, 13:27
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.03.2013, 13:27
Привет! Вот еще темы с решениями:

Импорт данных из Mathcad в Matlab
Здравствуйте. Помогите разобраться. Из Mathcad экспортированы данные в формате...

Дифференциальные уравнения в matlab и mathcad
Напишите пожалуйста решение дифференциального уравнения различными методами...

В чем отличие Mathcad и matlab
Добрый день! В чем отличие?

Решение задачи из Mathcad в Matlab
Данную задачу в Mathcad я решил, а вот в Matlab - не знаю как решать, так как...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
8
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru