Форум программистов, компьютерный форум CyberForum.ru

Метод вращений Якоби с++ - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 49, средняя оценка - 4.98
victory19933
2 / 2 / 0
Регистрация: 24.04.2011
Сообщений: 70
18.11.2012, 18:00     Метод вращений Якоби с++ #1
Используя метод Якоби найти с точностью 0.0001 все собственные значения и векторы матрицы А:
{1.953, 1.535, -3.220 , 0.826, 1.592 , 0.042, 0.570},
{1.535, -0.893, 0.783, 2.583, -0.117, -1.178, 1.240},
{-3.220, 0.783, 1.007 , 3.216, -1.908 , 0.305 , 0.342},
{0.826, 2.583, 3.216, 0.421, 2.132, 3.130, -2.218},
{1.592, -0.117, -1.908, 2.132, 2.530, -0.316, 0.368},
{0.042, -1.178, 0.305, 3.130, -0.316, -2.476, -0.374},
{0.570, 1.240, 0.342 , -2.218, 0.368, -0.374, 1.916}

(перечитала уже кучу информации, начиная с Википедии, заканчивая форумами. никак не могу собрать всё и сообразить.просмотрела несколько программ из интернета, но все они никак не хотят работать...)Помогите плиз!без вашей помощи не справиться
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.11.2012, 18:00     Метод вращений Якоби с++
Посмотрите здесь:

C++ метод вращений
C++ метод якоби
C++ Метод итерации( Якоби)
Метод вращений с построением КЮЭР-разложения C++
C++ КУЭР-разложение методом вращений
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 103
Записей в блоге: 1
19.11.2012, 21:55     Метод вращений Якоби с++ #2
Вот рабочий код:

C++
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
#include "stdafx.h"
#include "stdio.h"
#include <iostream>
#include "math.h"
 
const double PI = 3.1415926536;
 
bool isSimmetrial( double **coefficients, int numberOfEquation ) {
bool result = true;
int i, j;
for ( i = 0; i < numberOfEquation; i++ ) {
for ( j = i + 1; j < numberOfEquation; j ++ ) {
if ( coefficients[i][j] != coefficients[j][i] ) {
result = false;
break;}}
if ( !result ) {break;}}
return result;}
int wrachenie( double **coefficients, int numberOfEquation,
double **solution, double precision ) {
int result = 1;
int i, j, k;
int maxI, maxJ;
double max, fi;
double** matricaPoworota;
matricaPoworota = new double*[numberOfEquation];
for ( i = 0; i < numberOfEquation; i ++ ) {
matricaPoworota[i] = new double[numberOfEquation];}
double** temp;
temp = new double*[numberOfEquation];
for ( i = 0; i < numberOfEquation; i ++ ) {
temp[i] = new double[numberOfEquation];}
double fault = 0.0;
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = i + 1; j < numberOfEquation; j ++ ) {
fault = fault + coefficients[i][j]*coefficients[i][j];}}
fault = sqrt( 2*fault );
while ( fault > precision ) {
max = 0.0;
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = i + 1; j < numberOfEquation; j ++ ) {
if ( coefficients[i][j] > 0 && coefficients[i][j] > max ) {
max = coefficients[i][j];
maxI = i;
maxJ = j;}
else if ( coefficients[i][j] < 0 && - coefficients[i][j] > max ) {
max = - coefficients[i][j];
maxI = i;
maxJ = j;}}}
for ( i = 0; i < numberOfEquation; i ++ ){
for ( j = 0; j < numberOfEquation; j ++ ){
matricaPoworota[i][j] = 0;}
matricaPoworota[i][i] = 1;}
if ( coefficients[maxI][maxI] == coefficients[maxJ][maxJ] ) {
matricaPoworota[maxI][maxI] = matricaPoworota[maxJ][maxJ] =
matricaPoworota[maxJ][maxI] = sqrt( 2.0 ) / 2.0;
matricaPoworota[maxI][maxJ] = - sqrt( 2.0 ) / 2.0;}
else {
fi = 0.5 * atan( ( 2.0 * coefficients[maxI][maxJ] ) /
( coefficients[maxI][maxI] - coefficients[maxJ][maxJ] ) );
matricaPoworota[maxI][maxI] = matricaPoworota[maxJ][maxJ] = cos( fi );
matricaPoworota[maxI][maxJ] = - sin( fi );
matricaPoworota[maxJ][maxI] = sin( fi );}
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = 0; j < numberOfEquation; j ++ ) {
temp[i][j] = 0.0;}}
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = 0; j < numberOfEquation; j ++ ) {
for ( k = 0; k < numberOfEquation; k ++ ) {
temp[i][j] = temp[i][j] + matricaPoworota[k][i] * coefficients[k][j];}}}
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = 0; j < numberOfEquation; j ++ ) {
coefficients[i][j] = 0.0;}}
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = 0; j < numberOfEquation; j ++ ) {
for ( k = 0; k < numberOfEquation; k ++ ) {
coefficients[i][j] = coefficients[i][j] +
temp[i][k] * matricaPoworota[k][j];}}}
fault = 0.0;
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = i + 1; j < numberOfEquation; j ++ ) {
fault = fault + coefficients[i][j]*coefficients[i][j];}}
fault = sqrt( 2*fault );
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = 0; j < numberOfEquation; j ++ ) {
temp[i][j] = 0.0;}}
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = 0; j < numberOfEquation; j ++ ) {
for ( k = 0; k < numberOfEquation; k ++ ) {
temp[i][j] = temp[i][j] + solution[i][k] * matricaPoworota[k][j];}}}
for ( i = 0; i < numberOfEquation; i ++ ) {
for ( j = 0; j < numberOfEquation; j ++ ) {
solution[i][j] = temp[i][j];}}
result++;}
return result;}
 
using namespace std;
 
void main() 
{
setlocale(LC_ALL,"rus");
int i, j;
int size;
double **coefficients, **solution, precision;
cout << "Метод вращения Якоби.\nВведите размерность матрицы: ";
cin >> size;
cout << "\nВведите элементы матрицы: \n";
coefficients = new double*[size];
solution = new double*[size];
for ( i = 0; i < size; i++ ) {
coefficients[i] = new double[size];
solution[i] = new double[size];}
for ( i = 0; i < size; i ++ ){
for ( j = 0; j < size; j ++ ){
solution[i][j] = 0;}
solution[i][i] = 1;}
for ( i = 0; i < size; i ++ ){
cout << "Enter " << i + 1 << " row: ";
for ( j = 0; j < size; j ++ ){
cin >> coefficients[i][j];}}
cout << "Введите точность расчета: ";
cin >> precision;
if ( !isSimmetrial( coefficients, size ) ) {
cout << "Матрица не симметричная";}
else {
int steps = wrachenie( coefficients, size, solution, precision );
cout << "Reshenie:\n";
for ( i = 0; i < size; i++ ) {
cout << "Собственный вектор k " << i + 1 << ":\n";
for ( j = 0; j < size; j ++ ){
cout << solution[j][i] << "\n";}}
cout << "Собственные значения:\n";
for ( i = 0; i < size; i++ ) {
cout << coefficients[i][i] << "\n";}
cout << "Общее число шагов: " << steps;}
cout << "\nPress \"Enter\" to continue..." << endl;
getchar();}
Кстати, этот алгоритм является общим. При исключении проверки на симметричность он будет правильно искать собственные вектора и собственные значения любой действительной матрицы.
williamvolkov
0 / 0 / 0
Регистрация: 17.11.2010
Сообщений: 23
28.11.2012, 00:51     Метод вращений Якоби с++ #3
Цитата Сообщение от Nuclear_Razor Посмотреть сообщение
stdafx.h
а где данный файл?
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 103
Записей в блоге: 1
28.11.2012, 01:15     Метод вращений Якоби с++ #4
Что? В смысле? Это библиотека, необходимая для компиляции в MVS.
williamvolkov
0 / 0 / 0
Регистрация: 17.11.2010
Сообщений: 23
28.11.2012, 01:49     Метод вращений Якоби с++ #5
Цитата Сообщение от Nuclear_Razor Посмотреть сообщение
Что? В смысле? Это библиотека, необходимая для компиляции в MVS.
работаю в elipse и он жалуется что нет файлика) если можешь выложи библиотечку
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 103
Записей в блоге: 1
28.11.2012, 01:56     Метод вращений Якоби с++ #6
Чувак, меня просили код на С++. А проблемы заголовочных файлов - это не ко мне. Скачай с сайта Майкрософт MVS 2012 Express убого вставь код и запусти, в чем проблема?)
MrGluck
Ворчун
Эксперт С++
 Аватар для MrGluck
4919 / 2662 / 243
Регистрация: 29.11.2010
Сообщений: 7,399
28.11.2012, 01:58     Метод вращений Якоби с++ #7
C++
1
2
3
4
#include <iostream>
#include <cmath>
 
int main()
И не забудьте память вернуть delete-ом.
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 103
Записей в блоге: 1
28.11.2012, 02:01     Метод вращений Якоби с++ #8
Ну. На счет этого спорить не буду =)
williamvolkov
0 / 0 / 0
Регистрация: 17.11.2010
Сообщений: 23
01.12.2012, 14:11     Метод вращений Якоби с++ #9
Цитата Сообщение от Nuclear_Razor Посмотреть сообщение
Чувак, меня просили код на С++. А проблемы заголовочных файлов - это не ко мне. Скачай с сайта Майкрософт MVS 2012 Express убого вставь код и запусти, в чем проблема?)
ок, нашел)все работает..а можно ли как нибудь сделать ввод из файлов? если знаешь, покажи пожалуйста...

Добавлено через 15 часов 0 минут
с файлами уже тож понял..ток единственно тут вектора невязки не правильно находит. а так в целом норм алгоритм)
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
26.09.2014, 16:26     Метод вращений Якоби с++
Еще ссылки по теме:

C++ СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя
Метод Якоби. Выводит результат -1.INF и -1.IND C++
C++ Итерационные методы. Метод Якоби

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

Или воспользуйтесь поиском по форуму:
Katya_dmtr
Сообщений: n/a
26.09.2014, 16:26     Метод вращений Якоби с++ #10
тогда это программа только для симметричных матриц?
Yandex
Объявления
26.09.2014, 16:26     Метод вращений Якоби с++
Ответ Создать тему
Опции темы

Текущее время: 23:19. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru