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

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

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

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

18.11.2012, 18:00. Просмотров 7372. Ответов 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}

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

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

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 104
Записей в блоге: 1
Завершенные тесты: 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();}
Кстати, этот алгоритм является общим. При исключении проверки на симметричность он будет правильно искать собственные вектора и собственные значения любой действительной матрицы.
2
williamvolkov
0 / 0 / 0
Регистрация: 17.11.2010
Сообщений: 23
28.11.2012, 00:51 #3
Цитата Сообщение от Nuclear_Razor Посмотреть сообщение
stdafx.h
а где данный файл?
0
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 104
Записей в блоге: 1
Завершенные тесты: 1
28.11.2012, 01:15 #4
Что? В смысле? Это библиотека, необходимая для компиляции в MVS.
0
williamvolkov
0 / 0 / 0
Регистрация: 17.11.2010
Сообщений: 23
28.11.2012, 01:49 #5
Цитата Сообщение от Nuclear_Razor Посмотреть сообщение
Что? В смысле? Это библиотека, необходимая для компиляции в MVS.
работаю в elipse и он жалуется что нет файлика) если можешь выложи библиотечку
0
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 104
Записей в блоге: 1
Завершенные тесты: 1
28.11.2012, 01:56 #6
Чувак, меня просили код на С++. А проблемы заголовочных файлов - это не ко мне. Скачай с сайта Майкрософт MVS 2012 Express убого вставь код и запусти, в чем проблема?)
0
MrGluck
Модератор
Эксперт CЭксперт С++
7239 / 4407 / 642
Регистрация: 29.11.2010
Сообщений: 11,927
28.11.2012, 01:58 #7
C++
1
2
3
4
#include <iostream>
#include <cmath>
 
int main()
И не забудьте память вернуть delete-ом.
0
Nuclear_Razor
49 / 2 / 0
Регистрация: 22.07.2012
Сообщений: 104
Записей в блоге: 1
Завершенные тесты: 1
28.11.2012, 02:01 #8
Ну. На счет этого спорить не буду =)
0
williamvolkov
0 / 0 / 0
Регистрация: 17.11.2010
Сообщений: 23
01.12.2012, 14:11 #9
Цитата Сообщение от Nuclear_Razor Посмотреть сообщение
Чувак, меня просили код на С++. А проблемы заголовочных файлов - это не ко мне. Скачай с сайта Майкрософт MVS 2012 Express убого вставь код и запусти, в чем проблема?)
ок, нашел)все работает..а можно ли как нибудь сделать ввод из файлов? если знаешь, покажи пожалуйста...

Добавлено через 15 часов 0 минут
с файлами уже тож понял..ток единственно тут вектора невязки не правильно находит. а так в целом норм алгоритм)
0
Katya_dmtr
Сообщений: n/a
26.09.2014, 16:26 #10
тогда это программа только для симметричных матриц?
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
26.09.2014, 16:26
Привет! Вот еще темы с ответами:

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

Итерационные методы. Метод Якоби - 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 // число...


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

Или воспользуйтесь поиском по форуму:
Yandex
Объявления
26.09.2014, 16:26
Ответ Создать тему
Опции темы

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