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

Перемножения матриц (обратную на вектор). Найти ошибку в коде - C++

Восстановить пароль Регистрация
 
cj3687
2 / 2 / 0
Регистрация: 10.10.2013
Сообщений: 43
10.02.2014, 22:30     Перемножения матриц (обратную на вектор). Найти ошибку в коде #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
#include <iostream>
#include <locale.h>
 
using namespace std;
 
int gaus_obr(int cnt_str,double **mass,double **&A_obr)
{
   int i,j,k;
    //ñîçäГ*Г*ГЁГҐ åäèГ*ГЁГ·Г*îé Г¬Г*òðèöû
   A_obr=new double* [cnt_str];
   for(i=0;i<cnt_str;i++)
   {
   A_obr[i]=new double [cnt_str];
    for(j=0;j<cnt_str;j++)A_obr[i][j]=0;
    A_obr[i][i]=1;
   }
    //ïðÿìîé õîä ìåòîäîì ГѓГ*ГіГ±Г±Г*
    double a,b;
    for(i=0;i<cnt_str;i++)
   {
      a=mass[i][i];
      for(j=i+1;j<cnt_str;j++)
      {
        b=mass[j][i];
        for(k=0;k<cnt_str;k++)
         {
           mass[j][k]=mass[i][k]*b-mass[j][k]*a;
          A_obr[j][k]=A_obr[i][k]*b-A_obr[j][k]*a;
         }
      }
   }
     //îáðГ*ГІГ*ûé õîä âû÷èñëåГ*ГЁГї ýëåìåГ*òîâ îáðГ*ГІГ*îé Г¬Г*òðèöû
    double sum;
    for(i=0;i<cnt_str;i++)
    {
      for(j=cnt_str-1;j>=0;j--)
     {
      sum=0;
      for(k=cnt_str-1;k>j;k--)
      sum+=mass[j][k]*A_obr[k][i];
      if(mass[j][j]==0)
       {
        for(i=0;i<cnt_str;i++)
          delete []A_obr[i];
          delete []A_obr;
       return 0;
       }
       A_obr[j][i]=(A_obr[j][i]-sum)/mass[j][j];
     }
  }
return 1;
}
 
int main()
{
    setlocale(LC_CTYPE,"Russian");
//ГђГ*çìåð Г¬Г*òðèö:
int N;
cout<<"Size=";
cin>>N;
//Г€Г*äåêñГ*ûå ïåðåìåГ*Г*ûå:
int i,j,k;
//ÏåðâГ*Гї Г¬Г*òðèöГ*:
cout<<"Matrix A:\n";
double** A1 =new double*[N];
     for(i=0;i<N;i++)
     {
        A1[i]=new double [N];
        for(j=0;j<N;j++)
        {
            cout<<"A["<<i<<"]["<<j<<"]= ";
            cin>>A1[i][j];
        }
     }
//ÂòîðГ*Гї Г¬Г*òðèöГ*:
double B[N][N];
//Òðåòüÿ Г¬Г*òðèöГ* (ðåçóëüòГ*ГІ):
double C[N][N];
//Ââîä (ïîñòðî÷Г*ûé) ýëåìåГ*òîâ ïåðâîé Г¬Г*òðèöû:
 
//Ââîä (ïîñòðî÷Г*ûé) ýëåìåГ*òîâ ïåðâîé Г¬Г*òðèöû:
cout << "Matrix B:\n";
for(i=0; i < N; i++)
{
 for(j=0; j < 1; j++)
 {
 cout<<"B["<<i+1<<"]["<<j+1<<"]= ";
 cin>>B[i][j];
 }
}
cout<<"\n";cout<<"\nГЋГЎГ°Г*ГІГ*Г*Гї Г¬Г*òðèöГ* ГЂ:\n";
double **A_OBR;
     if(gaus_obr(N,A1,A_OBR))
     {
          for(i=0;i<N;i++)
          {
            for(j=0;j<N;j++)
            {
            cout<<A_OBR[i][j]<<'\t';
            }
            cout<<endl;
          }
          for(i=0;i<N;i++)
               delete[]A_OBR[i];
          delete[] A_OBR;
 
     }else
        cout<<"Error\n";
     for(i=0;i<N;i++)
           delete[]A1[i];
     delete[] A1;
 
cout<<"\n";cout<<"\n";
//Âû÷èñëåГ*ГЁГҐ ïðîèçâåäåГ*ГЁГї Г¬Г*òðèö:
cout << "Matrix C=AB:\n";
for(i=0; i < N; i++)
{
 for(j=0; j < 1; j++)
 {
 C[i][j]=0;
 for(k=0; k < N; k++)
 C[i][j]+=((A_OBR[i][k])*(B[k][j]));
 //Âûâîä Г§Г*Г*Г·ГҐГ*ГЁГї ýëåìåГ*ГІГ* Г*Г* ГЅГЄГ°Г*Г*:
 cout << C[i][j] << '\t';
 }
 cout << endl;
}
return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
10.02.2014, 22:30     Перемножения матриц (обратную на вектор). Найти ошибку в коде
Посмотрите здесь:

Сложение матриц,не могу найти ошибку в коде! C++
Найти обратную матрицу и умножить ее на вектор C++
C++ Функция перемножения матриц.
C++ Не могу написать программу перемножения матриц
Массив: перемножения двух матриц A,B размерности n*n C++
Написать процедуру перемножения матриц. C++
Алгоритм Штрассена для быстрого перемножения матриц C++
Умножение матриц. Не могу найти ошибку C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

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