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

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

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 34, средняя оценка - 4.65
Zurashka
0 / 0 / 0
Регистрация: 06.05.2012
Сообщений: 7
#1

Метод Гаусса решения СЛАУ с полным выбором. C++ - C++

20.04.2013, 00:10. Просмотров 4520. Ответов 1
Метки нет (Все метки)

Написал код реализующий метод Гаусса с полным выбором ведущего элемента, но ответ выдает не правильный. Помогите найти ошибку
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
137
138
139
140
141
142
143
#include "stdafx.h"
#include <conio.h>
#include <stdio.h>
#include <iostream>
#include <iomanip>
#include <math.h>
#define N 6
void PrintArray(double A[][N]){ //вывод матрицы
    for(int i=0; i<N; i++){
        for(int j=0; j<N; j++)
            std::cout<< A[i][j] <<"  ";
    std::cout<<std::endl;
    }
}
void PrintVector(double X[]){ // вывод вектора
    for(int i=0; i<N; i++){
        std::cout << std::fixed << std::setw(22) << std::setprecision(20) << X[i];
        std::cout<<std::endl;
    }
}
void GilbertM(double A[][N]){   //матрица Гильберта
    int i;
    int j;
    for (i = 0;i<N;i++){
        for(j=0;j<N;j++){
            A[i][j]=1./(i+1+j);
        }
    }
}
void VectorB(double B[]){  // заполняем вектор B
    int i;
    B[0] = 1;
    for (i = 1;i<N;i++){
        if (i % 2 == 0){
            B[i]=-(B[i-1]-1);
        }
        else{ B[i] = -(B[i - 1] + 1);}
    }
}
void SwapColumns(double A[][N],int c1,int c2){ // меняем столбцы матрицы
    double temp;
   int j;
   for(j=0; j<N; j++)
    {
     temp=A[j][c1];
     A[j][c1]=A[j][c2];
     A[j][c2]=temp;
    }
}
void SwapLines(double A[][N], double B[],int l1,int l2){ // меняем строки матрицы и вектора
   double temp;
   int j;
   for(j=0; j<N; j++)
    {
     temp=A[l1][j];
     A[l1][j]=A[l2][j];
     A[l2][j]=temp;
     temp=B[l1];
     B[l1]=B[l2];
     B[l2]=temp;
    }
}
void Destroy(double A[][N], double B[], int d){ // зануляем элементы под ведущим элементом
    int i,k;
    double v;
    for(i=d+1; i<N; i++){
        v=(-A[i][d]/A[d][d]);
        for(k=d; k<N; k++){
            A[i][k]=A[i][k]+v*A[d][k];
        }
        B[i]=B[i]+v*B[d];
    }
}
void Solution(double A[][N], double B[], double X[]){ //собственно само решение
    int i,j,s,m,k,l,temp;
    double max;
    int memory[N];
    double Y[N];
    double v;
    for (i=0; i<N; i++)
        memory[i]=i;
    for (i=0; i<N; i++){    
        max=0;
        k=i;
        l=i;
        for (s = i;s<N;s++){//ищем максимум в подматрице
            for (m = i;m<N;m++){
                if (fabs(A[s][m])>max){max=fabs(A[s][m]); k=s; l=m;}
            }
        }
//      if (k!=i)
        SwapLines(A,B,i,k);//меням строки
//      if (l!=i){
        SwapColumns(A,i,l);// меням столбцы и запоминаем какие переменные поменялись
        temp = memory[i];
        memory[i] = memory[l];
        memory[l] = temp;
//      }
        Destroy(A,B,i); // зануление
    }
    for (i=N-1; i>-1; i--){ //обратная подстановку
        v=0;
        for (j=N-1; j>i; j--)
            v=A[i][j]*Y[j]+v;
        Y[i]=(B[i]-v)/A[i][i];
    }
    for (i = 0; i < N; i++)//меняем местами переменные как должно быть
        X[memory[i]] = Y[i];
}
void F(double A[][N], double X[N], double B[N]){// подставляем A,B,X в уравнение AX-B и смотрим чему равен максимальный элемент 
    double L[N];
    double max;
    int i,j;
    max=0;
    for(i=0; i<N; i++){
        L[i]=0;
        for(j=0; j<N; j++){
            L[i]=L[i]+A[i][j]*X[j];
        }
        L[i]=L[i]-B[i];
    }
    for (i=0; i<N; i++){
        max=max+fabs(L[i]);
    }
    std::cout << std::fixed << std::setw(22) << std::setprecision(20) << max<<std::endl<<std::endl;
}
 
 
void main(){
    double B[N];
    double X[N];
    double A[N][N];
    double max;
    GilbertM(A);
    VectorB(B);
    Solution(A,B,X);
    PrintVector(X);
    GilbertM(A);
    VectorB(B);
    std::cout<<std::endl;
    F(A,X,B);
    getch();
}
Добавлено через 2 часа 56 минут
Нашел ошибку. Проблема исчерпана.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.04.2013, 00:10
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод Гаусса решения СЛАУ с полным выбором. C++ (C++):

Метод Гаусса для решения СЛАУ - C++
Здравствуйте! У меня есть алгоритм решения СЛАУ методом Гаусса, но почему-то он не работает((( Я проверил промежуточные матрицы, все ОК! И...

Метод Гаусса (подстановки) решения СЛАУ - C++
Всем доброго времени суток, написал я одну недопрограмму, потому как сделал её наполовину. Не хватает мне опыта и знания математики и опыта...

Метод Гаусса для решения СЛАУ с использованием одномерного массива - C++
Доброго всем времени суток. У меня есть прога для двумерного массива #include &lt;stdio.h&gt; #include &lt;iostream&gt; #include &lt;conio.h&gt; ...

СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя - C++
Помогите ребят. Не могу построить алгоритмы для этих методов Язык C++

СЛАУ методом Гаусса с выбором главного элемента в столбце - C++
СЛАУ методом Гаусса с выбором главного элемента по столбцу.Помогите пожалуйста с кодом программы,ошибки выдает #include &lt;locale&gt; ...

Решения СЛАУ методом Гаусса по шагу - C++
Добрый день! Очередной раз прошу советов! Имеется СЛАУ, которое требуется решить методом Гаусса на с++. Соответственно, на...

1
Lovan
1 / 1 / 0
Регистрация: 10.02.2013
Сообщений: 24
22.09.2013, 01:23 #2
А в чём ошибка-то была?
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
22.09.2013, 01:23
Привет! Вот еще темы с ответами:

Разработать программу для решения СЛАУ методом Гаусса. - C++
Разработать программу для решения СЛАУ методом Гаусса. вот заданный массив 7x1-x2+2x3+x4+x5=7 2x1+5x2-x3-x5=3 ...

Нужен рабочий пример решения СЛАУ методом Гаусса - C++
Нужен рабочий (ну или более менее рабочий) код для решения СЛАУ методом Гауса на С++. Уже замучился искать

Помогите написать программу для решения слау методом Гаусса - C++
Осталось как обычно пару дней до сдачи курсовой, а она еще не готова! Помогите пожалуйста. Вот код, который чуть-чуть работает Главная...

Кто сможет? Написать функцию для решения СЛАУ методом Гаусса - C++
Кто может написать функцию для решения СЛАУ методом Гаусса используя базовые знания.


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

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

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