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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 49, средняя оценка - 4.98
victory19933
2 / 2 / 0
Регистрация: 24.04.2011
Сообщений: 70
#1

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

18.11.2012, 18:00. Просмотров 7298. Ответов 9
Метки нет (Все метки)

Используя метод Якоби найти с точностью 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++
помогите, может у когото имеется приложение (на языке Pascal или C++), Решение СЛАУ методом вращений (Якоби). дали задание, уже...

метод вращений - C++
нужно найти собственные значения и векторы: вот прога тока она кажется путает индесы элемента a_{ij} из-за этого не правильно считает,...

Метод вращений с построением КЮЭР-разложения - C++
Я написал только Метод вращения #include<math.h> #include<stdio.h> #include<stdlib.h> #include<time.h> void f(double *mass,...

Метод Якоби - C++
здраствуйте. помогите пожалуйста... нужно прокоментировать метод Якоби.(если можно каждую строчку) заранее спасибо) вот мой код: ...

Метод Якоби - C++
Ребят помогите реализовать метод Якоби по алгоритму (находится во вложениях) проблема с циклами (не стал выкладывать)доп требование это...

Классический метод якоби - C++
Дана матрица размера 5*5: a=e^(abs(i-j)) и столбец свободных членов:{1 0 0 0 0} Помогите пожалуйста реализовать алгоритмы решения...

Метод итерации( Якоби) - C++
Доброго времени суток!) может у кого есть исходник или поможет решит поставленную задачу. Нужно написать программу которая решит...

После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 104
Записей в блоге: 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
Сообщений: 104
Записей в блоге: 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
Сообщений: 104
Записей в блоге: 1
28.11.2012, 01:56     Метод вращений Якоби с++ #6
Чувак, меня просили код на С++. А проблемы заголовочных файлов - это не ко мне. Скачай с сайта Майкрософт MVS 2012 Express убого вставь код и запусти, в чем проблема?)
MrGluck
Модератор
Эксперт CЭксперт С++
7178 / 4344 / 634
Регистрация: 29.11.2010
Сообщений: 11,820
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
Сообщений: 104
Записей в блоге: 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++
Не могу найти ошибку, скорее всего что-то с циклами. Формула самих итераций: Xk+1=AXk+Y Ребяят, помогите пожалуйста, уже несколько дней...

Метод Якоби. Выводит результат -1.INF и -1.IND - C++
Вот сделал прогу, но почему-то вместо нормального результата выводит это: result|x: -1.#INF -1.#IND -1.#INF -1.#IND -1.#IND ...

КУЭР-разложение методом вращений - C++
Я программу создал через двухмерный массив работает #include &lt;stdio.h&gt; #include &lt;math.h&gt; #include &lt;time.h&gt; #define _N 10 // число...

Вычисление символа якоби - C++
Нужно вычислить символ якоби. пользовался этим алгоритмом. Вот код: int a,b,t,c; cin &gt;&gt; a; cin &gt;&gt; b; r=1; if...

Решение СЛАУ методом Якоби - C++
Решить СЛАУ методом Якоби. Вывести значения решения, график зависимости нормы невязки от номера итерации и его значение, при котором...


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

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

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