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

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

Войти
Регистрация
Восстановить пароль
 
uto
6 / 6 / 0
Регистрация: 07.05.2009
Сообщений: 94
#1

Две матрицы - C++

22.06.2009, 09:58. Просмотров 434. Ответов 9
Метки нет (Все метки)

У меня есть два кода. Нужно код второй матрицы вставить в первый в указаное место (строка 127). Несколько дней уже его мучаю... У меня не получается, помогите, кто-нибудь. :'(

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
144
145
146
147
148
149
#include <vcl.h>
#include <iostream>
#include <string>
#include <sstream> // в некоторых версиях STL называется stringstream
#include <map>
#include <conio.h>
#include <math.h>
using namespace std;
 
//--------------------------------------------------------------
 
const int nn=100;
typedef double arr[nn][nn];
typedef double vect[nn];
 
int _tmain(int argc, _TCHAR* argv[])
{
    // HTTP-заголовок
    cout << "Content-type: text/html" << endl << endl;
 
    // Содержимое HTML-документа
    cout << "<HTML>" << endl;
    cout << "<BODY>" << endl;
 
    //Все переменные окружения (для проверки работоспособности) - extern char **__environ;
    cout << "<h1>CGI Environment</h1>" << endl << "<UL>" << endl;
    char **pe = environ;
    while(*pe != NULL)
    {
        cout << "<LI>" << *pe++ << "<br>" << endl;
    }
    cout << "</UL>" << endl << endl;
 
    // Выделяем строку с GET-параметрами переданными из формы
    printf("<h1>GET parameters:</h1>\n"); 
    string sq;
    char* pq = getenv("QUERY_STRING");
    if (pq!=NULL)
        sq = pq;
    cout << "sq=[" << sq.c_str() << "]";
    cout << endl;
 
    // формирование хеша (масива со строковым индексом) параметров GET
    map<string, string> gp; // результат данной операции - параметры GET
    string name, value;     // имя и зщначение параметра
    bool modename = true;   // режим накопления имени vs режим накопления значения
    sq+='&'; // добавляем ещё один символ '&' в конец
    // этот символ был разделителем пар имя=значение, 
    // а теперь его можно рассматривать как признак конца пары
    for (size_t pos=0; pos<sq.length(); pos++)
    {
        switch (sq[pos])
        {
        case '=' :  // конец имени данной пары
            modename=false;
            break;
        case '&' :  // конец значения данной пары и всей пары
            if (!name.empty())  // если имя определено
            {
                gp[name]=value; // сохраняем пару
            }
            modename=true;
            name.clear();       // важно сбросить накопители имени и значения!
            value.clear();
            break;
        default:    // обычный символ
            if (modename)       // идёт...
            {
                name+=sq[pos];  // либо в имя
            }
            else
            {
                value+=sq[pos]; // либо в значение
            }
        }
    }
 
    // Создан хеш типа строка -> строка
    // формируем из него хеш строка -> число
    map<string, double> gpf; // результат данной операции - параметры GET,со знач-ми float
    for (map<string,string>::const_iterator it=gp.begin(); it!=gp.end(); it++)
    {   // перебор всех пар
        gpf[it->first] = atof(it->second.c_str()); // и преобразование значений
    }
 
    // Вывод в браузер
    cout << "<h1>Pairs of GET-parameters decoded (name=value)</h1>" << endl;
    cout << "<UL>" << endl;
    for (map<string,double>::const_iterator it=gpf.begin(); it!=gpf.end(); it++)
    {
        cout << "<LI>" << it->first << "=" << it->second << "<br>" << endl;
    }
    cout << "</UL>" << endl << endl;
 
    // Приведение в "привычное состояние" для решения СЛАУ
    double a[9], b[3], exp; // параметры уравнения
    int i,j,k=0;
    for (i=1; i<=3; i++)
    {
        for (j=1; j<=3; j++)
        {
            stringstream namea;
            namea << "a" << i << j;
            a[k++]=gpf[namea.str()];
        }
 
        stringstream nameb;
        nameb << "b" << i;
        b[i-1] = gpf[nameb.str()];
    }
    exp=gpf["exp"];
 
    // Контрольный вывод на экран значений перед решением уравнений
    cout << "<h1>Matrix, vector and exp values ready to process:</h1>" << endl;
    for (i=0; i<3; i++) for(j=0; j<3; j++)
    {
        cout << "a[" << i << "," << j << "] = " << a[i*3+j] << "<BR>" << endl;
    }
    for (i=0; i<3; i++) 
    {
        cout << "b[" << i << "] = " << b[i] << "<BR>" << endl;
    }
    cout << "exp = " << exp << "<BR>" << endl << endl;
 
 
    // ... по идее тут надо вставить код для решения СЛАУ (метод Зейделя)...
 
    // **********************************************
 
 
    // Завершение HTML
    cout << "</BODY>" << endl;
    cout << "</HTML>" << endl;
 
        getch();
    return 0;
}
 
    bool exact(const double *px, const double *pp, const int *n)    {
        int i;
        bool f=true;
 
        for (i=0; i<*n; i++)
            if ( fabs( *(px+i) - *(pp+i) ) > exp){
                f=false;
                break;
            }
        return f;
    }
А вот второй код "Матрица Зейделя", который нужно вставить в первый так, чтобы нормально компилился

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
#include <vcl.h>
#include <math.h>
#include <stdio.h>
#include <iostream.h>
#include <conio.h>
 
//---------------------------------------------------------------------------
 
const int nn=100;
typedef double arr[nn][nn];
typedef double vect[nn];
 
void input(int*, double *);
bool exact(const double *px, const double *pp, const int *n);
 
arr A;
double e;
int main()
{
 
    int n,i,j;
    double sum;
    vect x,p,B;
 
    input(&n,B);             // Ввести n,A,B.
    for (i=0; i<n; i++){     // Начало
        x[i]=1;
        p[i]=1;
    }
 
    do {
        for (i=0; i<n; i++){
            sum = 0;
            for (j=0; j<n; j++)
                if (i!=j) sum+=x[j]*A[i][j];
 
            p[i]=x[i];
            x[i]=(B[i] - sum)/A[i][i];
 
        }
    } while (!exact(x,p, &n)) ;
 
 
    for (j=0; j<n; j++)
        printf("X%d = %.3f\n",j+1,x[j]);
 
    puts("Verification:");
     for (i=0; i<n; i++){
            sum = 0;
            for (j=0; j<n; j++){
                 sum+=x[j]*A[i][j];
                 if ((j==0) && (A[i][j]*x[j]<0)) cout <<"-";
                 printf("%.2f*%.2f",fabs(A[i][j]),fabs(x[j]));
                 if ((j!=n-1)&&(A[i][j+1]*x[j]>=0)) cout <<" + ";
                 else if ((j!=n-1)&&(A[i][j+1]*x[j]<0)) cout <<" - ";
                 else {cout << " = ";
                       printf("%.3f ;\n", sum);
                 }
           getch(); }
     }
 
}
 
void input(int *m, double *PtrB)
{
    int i,j;
 
    cout << "Enter e.\n";
    cin >> e;
    cout << "\nEnter n.\n";
    cin >> *m;
    cout << "\nEnter matrix A.\n";
    for (i=0; i<*m; i++)
        for (j=0; j<*m; j++)
            cin >> A[i][j];
    cout << "\nEnter matrix B.\n";
    for (i=0; i<*m; i++)
        cin >> *(PtrB+i);
    
}
 
bool exact(const double *px, const double *pp, const int *n)
{
    int i;
    bool f=true;
 
    for (i=0; i<*n; i++)
        if ( fabs( *(px+i) - *(pp+i) ) > e){
            f=false;
            break;
        }
    return f;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
22.06.2009, 09:58
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Две матрицы (C++):

Функция: получить две матрицы и вернуть номер матрицы, среднее арифметическое элементов которой максимально - C++
нужна помощь в этой задаче!!! Написать функцию, получающую в качестве параметров два массива целого типа произвольного размера и...

Даны две целочисленные матрицы 4-го порядка. Получить новую * матрицу путем вычитания из элементов каждого столбца первой матрицы - C++
даны две целочисленные матрицы 4-го порядка. Получить новую матрицу путем вычитания из элементов каждого столбца первой матрицы сумму...

Даны две действительные квадратные матрицы порядка n. Получить новую матрицу: а) умножением элементов каждой строки первой матрицы на наибольш - C++
Даны две действительные квадратные матрицы порядка n. Получить новую матрицу: а) умножением элементов каждой строки первой матрицы ...

1)Даны две действительные квадратные матрицы порядка n. Получить новую матрицу : А) умножение элементов каждой строки первой матрицы на наибольшее из - C++
1)Даны две действительные квадратные матрицы порядка n. Получить новую матрицу : А) умножение элементов каждой строки первой матрицы на...

Даны две действительные квадратные матрицы порядка n. Получить новую матрицу: путем умножения элементов каждой строки первой матрицы на наибольшее из - C++
Даны две действительные квадратные матрицы порядка n. Получить новую матрицу: путем умножения элементов каждой строки первой матрицы на...

Перемножить две матрицы - C++
Сформировать два одномерных массива целочисленных из 10 элементов каждый, а потом сформировать матрицу MxN элементами которого являются...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
accept
4821 / 3241 / 165
Регистрация: 10.12.2008
Сообщений: 10,682
22.06.2009, 12:26 #2
у тебя во втором коде используются глобальные переменные, из-за этого теперь сначала надо перевести их в локальные, потом оформить main в виде отдельной функции и ещё отстригать шелуху вроде ввода значений и вывода результатов, ага, и по мере его переписывания, ещё и разбираться в нём, так как переменные названы очень информативно

первый-то нормально подгтовлен (его можно разнести по функциям в любой момент)
uto
6 / 6 / 0
Регистрация: 07.05.2009
Сообщений: 94
22.06.2009, 12:48  [ТС] #3
Так первый то сам писал, а вот второй... ппц, ну не шарю я =/
Словами и я могу... Вы делом помогите я четвертый день с ней мучаюсь, а завтра сдавать
accept
4821 / 3241 / 165
Регистрация: 10.12.2008
Сообщений: 10,682
22.06.2009, 13:59 #4
Цитата Сообщение от uto
Так первый то сам писал
что значит строчка ?
Код
atof(it->second.c_str());
вот именно, что второй надо с нуля писать и по-нормальному
uto
6 / 6 / 0
Регистрация: 07.05.2009
Сообщений: 94
22.06.2009, 23:39  [ТС] #5
Ой ппц, так сказал же первый код не трогайте, это CGI, вставьте в 127 строку второй код, я просто с переменными там чето не могу разобраться, компилятор вечно на чтото жалуется

Добавлено через 5 часов 4 минуты 3 секунды
Так, ну я собсно помудрил еще посидел, и вот до чего в итоге дошел...

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
144
145
146
147
148
149
150
#include <vcl.h>
#include <iostream.h>
#include <string.h>
#include <sstream> // в некоторых версиях STL называется stringstream
#include <map>
 
const int nn=100;
typedef double arr[nn][nn];
typedef double vect[nn];
 
arr A;
 
int main()  {
    // HTTP-заголовок
    cout << "Content-type: text/html\n\n";
    
    // Содержимое HTML-документа
    cout << "<HTML>" << endl;
    cout << "<BODY>" << endl;
    
    //Все переменные окружения (для проверки работоспособности) - extern char **__environ;
    cout << "<h1>CGI Environment</h1>" << endl << "<UL>" << endl;
    char **pe = environ;
    while(*pe != NULL)
    {
        cout << "<LI>" << *pe++ << "<br>" << endl;
    }
    cout << "</UL>" << endl << endl;
    
    // Выделяем строку с GET-параметрами переданными из формы
    cout << "<h1>GET parameters:</h1>\n";
    string sq;
    char* pq = getenv("QUERY_STRING");
    if (pq!=NULL)
        sq = pq;
    cout << "sq=[" << sq.c_str() << "]";
    cout << endl;
    
    // формирование хеша (масива со строковым индексом) параметров GET
    map <string, string> gp;    // результат данной операции - параметры GET
    string name, value;     // имя и зщначение параметра
    bool modename = true;   // режим накопления имени vs режим накопления значения
    sq+='&'; // добавляем ещё один символ '&' в конец
    // этот символ был разделителем пар имя=значение, 
    // а теперь его можно рассматривать как признак конца пары
    for (size_t pos=0; pos<sq.length(); pos++)
    {
        switch (sq[pos])
        {
        case '=' :  // конец имени данной пары
            modename=false;
            break;
        case '&' :  // конец значения данной пары и всей пары
            if (!name.empty())  // если имя определено
            {
                gp[name]=value; // сохраняем пару
            }
            modename=true;
            name.clear();       // важно сбросить накопители имени и значения!
            value.clear();
            break;
        default:    // обычный символ
            if (modename)       // идёт...
            {
                name+=sq[pos];  // либо в имя
            }
            else
            {
                value+=sq[pos]; // либо в значение
            }
        }
    }
    
    // Создан хеш типа строка -> строка
    // формируем из него хеш строка -> число
    map <string, double> gpf; // результат данной операции - параметры GET,со знач-ми float
    for (map <string,string>::const_iterator it=gp.begin(); it!=gp.end(); it++)
    {   // перебор всех пар
        gpf[it->first] = atof(it->second.c_str()); // и преобразование значений
    }
    
    // Вывод в браузер
    cout << "<h1>Pairs of GET-parameters decoded (name=value)</h1>" << endl;
    cout << "<UL>" << endl;
    for (map <string,double>::const_iterator it=gpf.begin(); it!=gpf.end(); it++)
    {
        cout << "<LI>" << it->first << "=" << it->second << "<br>" << endl;
    }
    cout << "</UL>" << endl << endl;
    
    // Приведение в "привычное состояние" для решения СЛАУ
    double a[9], b[3], exp; // параметры уравнения
    int i,j,k=0;
    for (i=1; i<=3; i++)
    {
        for (j=1; j<=3; j++)
        {
            stringstream namea;
            namea << "a" << i << j;
            a[k++]=gpf[namea.str()];
        }
 
        stringstream nameb;
        nameb << "b" << i;
        b[i-1] = gpf[nameb.str()];
    }
    exp=gpf["exp"];
    
    // Контрольный вывод на экран значений перед решением уравнений
    cout << "<h1>Matrix, vector and exp values ready to process:</h1>" << endl;
    for (i=0; i<3; i++) for(j=0; j<3; j++)
    {
        cout << "a[" << i << "," << j << "] = " << a[i*3+j] << "<BR>" << endl;
    }
    for (i=0; i<3; i++) 
    {
        cout << "b[" << i << "] = " << b[i] << "<BR>" << endl;
    }
    cout << "exp = " << exp << "<BR>" << endl << endl;
    
    
    // ... по идее тут надо вставить код для решения СЛАУ (метод Зейделя)...
    vect x, p, B;
    double sum;
for (i=0; i<3; i++){     // Начало
        x[i]=1;
        p[i]=1;
    }
 
    do {
        for (i=0; i<3; i++){
            sum = 0;
            for (j=0; j<3; j++)
                if (i!=j) sum+=x[j]*A[i][j];
 
            p[i]=x[i];
            x[i]=(B[i] - sum)/A[i][i];
 
        }
    } 
 
for (j=0; j<3; j++)
printf("X%d = %.3f\n",j+1,x[j]);
    // **********************************************
    
    
    // Завершение HTML
    cout << "</BODY>" << endl;
    cout << "</HTML>" << endl;
}
P.s. Хоть это и не компилируемый вариант, но уже почти всё...
Компилятор ругается, говорит что в 141 строке нужна функция while, но я не знаю, что мне писать (хотя мне и нужно где-то вписать точность по эпсиланту).
Что делать??? Я в ступоре

Добавлено через 3 часа 51 минуту 43 секунды
Аааааа, я сделал, помогите откомпилировать программу!!! Плиииз, там одно предупреждение, а из-за него программа не работает =(

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
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
#include <vcl.h>
#include <iostream.h>
#include <string.h>
#include <sstream> // в некоторых версиях STL называется stringstream
#include <map>
#include <math.h>
#include <conio.h>
 
const int nn=100;
typedef double arr[nn][nn];
typedef double vect[nn];
 
bool exact(const double *px, const double *pp);
 
arr A;
double eps;
 
int main()  {
    // HTTP-заголовок
    cout << "Content-type: text/html\n\n";
    
    // Содержимое HTML-документа
    cout << "<HTML>" << endl;
    cout << "<BODY>" << endl;
    
    //Все переменные окружения (для проверки работоспособности) - extern char **__environ;
    cout << "<h1>CGI Environment</h1>" << endl << "<UL>" << endl;
    char **pe = environ;
    while(*pe != NULL)
    {
        cout << "<LI>" << *pe++ << "<br>" << endl;
    }
    cout << "</UL>" << endl << endl;
    
    // Выделяем строку с GET-параметрами переданными из формы
    cout << "<h1>GET parameters:</h1>\n";
    string sq;
    char* pq = getenv("QUERY_STRING");
    if (pq!=NULL)
        sq = pq;
    cout << "sq=[" << sq.c_str() << "]";
    cout << endl;
    
    // формирование хеша (масива со строковым индексом) параметров GET
    map <string, string> gp;    // результат данной операции - параметры GET
    string name, value;     // имя и зщначение параметра
    bool modename = true;   // режим накопления имени vs режим накопления значения
    sq+='&'; // добавляем ещё один символ '&' в конец
    // этот символ был разделителем пар имя=значение, 
    // а теперь его можно рассматривать как признак конца пары
    for (size_t pos=0; pos<sq.length(); pos++)
    {
        switch (sq[pos])
        {
        case '=' :  // конец имени данной пары
            modename=false;
            break;
        case '&' :  // конец значения данной пары и всей пары
            if (!name.empty())  // если имя определено
            {
                gp[name]=value; // сохраняем пару
            }
            modename=true;
            name.clear();       // важно сбросить накопители имени и значения!
            value.clear();
            break;
        default:    // обычный символ
            if (modename)       // идёт...
            {
                name+=sq[pos];  // либо в имя
            }
            else
            {
                value+=sq[pos]; // либо в значение
            }
        }
    }
    
    // Создан хеш типа строка -> строка
    // формируем из него хеш строка -> число
    map <string, double> gpf; // результат данной операции - параметры GET,со знач-ми float
    for (map <string,string>::const_iterator it=gp.begin(); it!=gp.end(); it++)
    {   // перебор всех пар
        gpf[it->first] = atof(it->second.c_str()); // и преобразование значений
    }
    
    // Вывод в браузер
    cout << "<h1>Pairs of GET-parameters decoded (name=value)</h1>" << endl;
    cout << "<UL>" << endl;
    for (map <string,double>::const_iterator it=gpf.begin(); it!=gpf.end(); it++)
    {
        cout << "<LI>" << it->first << "=" << it->second << "<br>" << endl;
    }
    cout << "</UL>" << endl << endl;
    
    // Приведение в "привычное состояние" для решения СЛАУ
    double a[9], b[3]; // параметры уравнения
    int i,j,k=0;
    for (i=1; i<=3; i++)
    {
        for (j=1; j<=3; j++)
        {
            stringstream namea;
            namea << "a" << i << j;
            a[k++]=gpf[namea.str()];
        }
 
        stringstream nameb;
        nameb << "b" << i;
        b[i-1] = gpf[nameb.str()];
    }
    eps=gpf["eps"];
    
    // Контрольный вывод на экран значений перед решением уравнений
    cout << "<h1>Matrix, vector and eps values ready to process:</h1>" << endl;
    for (i=0; i<3; i++) for(j=0; j<3; j++)
    {
        cout << "a[" << i << "," << j << "] = " << a[i*3+j] << "<BR>" << endl;
    }
    for (i=0; i<3; i++) 
    {
        cout << "b[" << i << "] = " << b[i] << "<BR>" << endl;
    }
    cout << "eps = " << eps << "<BR>" << endl << endl;
    
    
    // ... по идее тут надо вставить код для решения СЛАУ (метод Зейделя)...
    vect x, p, B;
    double sum;
for (i=0; i<3; i++){     // Начало
        x[i]=1;
        p[i]=1;
    }
 
    do {
        for (i=0; i<3; i++){
            sum = 0;
            
            for (j=0; j<3; j++)
                if (i!=j) sum+=x[j]*A[i][j];
            p[i]=x[i];
            x[i]=(B[i] - sum)/A[i][i];
 
        }
    }  while (!exact(x,p)) ;
 
for (j=0; j<3; j++)
printf("X%d = %.3f\n",j+1,x[j]);
    // **********************************************
     getch();
    
    // Завершение HTML
    cout << "</BODY>" << endl;
    cout << "</HTML>" << endl;
}
 
bool exact(const double *px, const double *pp)
{
    int i;
    bool done=true;
 
    for (i=0; i<3; i++)
        if ( fabs( *(px+i) - *(pp+i) ) > eps){
            done=false;
            break;
        }
    return done;
}
142ая строка...
Пишет мол A[i][i] не может быть равна нулю.... Я и так е и сяк, ну не хочет она корректно работать, помогите, плиииз. То она у меня выдает еденицы, то тупо подвисает
insolent
23.06.2009, 00:32
  #6

Не по теме:

Цитата Сообщение от uto Посмотреть сообщение
Пишет мол A[i][i] не может быть равна нулю
Правильно пишет - на ноль делить нельзя, ещё в младших классах учат

uto
6 / 6 / 0
Регистрация: 07.05.2009
Сообщений: 94
23.06.2009, 00:39  [ТС] #7
Да ёк макарек, это я и так вижу. Вы подскажите как исправить. У меня уже температура поднялась, голова болит, ничего не соображаю. Ужас... Помогите, пожаааалуйста!!!! :'(
accept
4821 / 3241 / 165
Регистрация: 10.12.2008
Сообщений: 10,682
23.06.2009, 03:54 #8
так он даже не работает, оказывается, в таком виде (пространства имён нет), я запустил его, может быть он даже правильно считает

вывод
Код
[guest@station tmp]$ ./test
Enter e.
0.1

Enter n.
3

Enter matrix A.
1 2 3
1 2 3
1 2 3

Enter matrix B.
5
5
5
X1 = 0.000
X2 = 1.000
X3 = 1.000
Verification:
1.00*0.00 + 2.00*1.00 +
3.00*1.00 = 5.000 ;

1.00*0.00 +
2.00*1.00 +
3.00*1.00 = 5.000 ;

1.00*0.00 +
2.00*1.00 +
3.00*1.00 = 5.000 ;

[guest@station tmp]$
Добавлено через 58 минут 18 секунд
ещё неправильная verification (показывает минус, хотя там плюс)

Код
[guest@station tmp]$ ./test
Enter exponent.
0.1

Enter n.
2

Enter matrix A.
1 2
3 4

Enter matrix B.
1
1
X1 = -1.000
X2 = 1.000
Verification:
-1.00 * 1.00 - 2.00 * 1.00 = 1.000 ;
-3.00 * 1.00 - 4.00 * 1.00 = 1.000 ;
[guest@station tmp]$
Добавлено через 59 минут 58 секунд
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
#include <cmath>
#include <iostream>
 
using namespace std; 
 
void MethodZeidelya(
         double *px,
         double (*pa)[100],
         double *pb,
         const int n,
         const double expon);
bool exact(
         const double *px,
         const double *pp,
         const int n,
         const double e);
 
int main()
{
    const double expon = 0.1;
    double a[100][100],
           b[100],
           x[100],
           sum;
    int n, i, j;
    
    cout << "Enter dimension: ";
    cin >> n;
    
    cout << "Enter matrix A: " << endl;
    for (i = 0; i < n; i++)
        for (j = 0; j < n; j++)
            cin >> a[i][j];
    
    cout << "Enter comumn B: " << endl;
        for (i = 0; i < n; i++)
            cin >> b[i];
            
    cout << "Exponent = "
         << expon
         << ","
         << " " "Dimension = "
         << n
         << endl;
            
    MethodZeidelya(x, a, b, n, expon);
 
    for (i = 0; i < n; i++)
        cout << "x"
             << i+1
             << " = "
             << x[i]
             << endl;
        
    return 0;    
}
 
void MethodZeidelya(
         double *px,
         double (*pa)[100],
         double *pb,
         const int n,
         const double expon)
{
    int i, j;
    double sum, p[100];
    
    for (i = 0; i < n; i++)
        px[i] = p[i] = 1.0;
    
    do
        for (i = 0; i < n; i++) {
            sum = 0.0;
            
            for (j = 0; j < n; j++)
                if (i != j)
                    sum += px[j] * pa[i][j];
            
            p[i] = px[i];
            px[i] = (pb[i] - sum) / pa[i][i];
 
        }
    while (!exact(px, p, n, expon));
}
 
bool exact(
         const double *px,
         const double *pp,
         const int n,
         const double e)
{
    int i;
    bool flag = true;
 
    for (i = 0; i < n; i++)
        if (fabs(*(px+i)-*(pp+i)) > e) {
            flag = false;
            break;
        }
    return flag;
}
на данном этапе код уже автономный

Добавлено через 1 минуту 57 секунд
и вывод к нему
Код
[guest@station tmp]$ ./test
Enter dimension: 3
Enter matrix A:
1 2 3
1 2 3
1 2 3
Enter comumn B:
5
5
5
Exponent = 0.1, Dimension = 3
x1 = 0
x2 = 1
x3 = 1
[guest@station tmp]$
uto
6 / 6 / 0
Регистрация: 07.05.2009
Сообщений: 94
23.06.2009, 08:08  [ТС] #9
Проверка мне не нужна... верификэйшн (я ее удалил), мне бы перехватить как-нибудь ошибку, чтобы прога спокойно компилировалась и делила на ноль...
insolent
23.06.2009, 13:24     Две матрицы
  #10

Не по теме:

Цитата Сообщение от uto Посмотреть сообщение
делила на ноль...
делить на ноль нельзя...то-есть можно, но получится бесконечность и далее...

MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
23.06.2009, 13:24
Привет! Вот еще темы с ответами:

Как сложить две матрицы ? - C++
Есть такой код #include &lt;iostream&gt; #include &lt;ctime&gt; #include &lt;cstdlib&gt; #include &lt;iomanip&gt; using namespace std; int main() { ...

Функция разбиения матрицы на две - C++
Доброго времени суток уважаемые форумчане. Столкнулся с проблемой. Необходимо написать функцию разбиения матрицы на две матрицы, в качестве...

Даны две матрицы произвольной размерности - C++
При выполнении этого задания необходимо написать две функции. Одна из этих функций должна получать и обрабатывать матрицу целиком. Другая...

Считать с файла две матрицы и перемножить их - C++
Считать с файла две матрицы и перемножить их


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

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

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