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

Поиск обратной матрицы методом Гаусса - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 105, средняя оценка - 4.69
прасковья
0 / 0 / 0
Регистрация: 13.09.2009
Сообщений: 3
28.09.2009, 00:41     Поиск обратной матрицы методом Гаусса #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
#include<iostream>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
using namespace std;
static int n;
void obrat(double **a)
{double **d,vrem,max,tryam;
int i,j,*p,jmax,temp;
d=new double *[2*n];//videlenie dlinnoi stroki
for (i=0;i<n;i++) d[i]=new double[n];
for(i=0;i<n;i++)
{for(j=0;j<n;j++)d[i][j]=a[i][j];
  for(j=n;j<2*n;j++)
   {if(j==n+i)d[i][j]=1;
     else d[i][j]=0;
   }
}
p=new int[n];
for(j=0;j<n;j++)p[j]=j;
//k-nomer pervogo uravn
for(int k=0;k<n-1;k++)
 {max=fabs(d[k][k]);
  jmax=j;
  for(j=k;j<n;j++)
   if(fabs(d[k][j])>max)
    {max=fabs(d[k][j]);
     jmax=j;
      temp=p[k];
      p[k]=p[jmax];
      p[jmax]=temp;
      }
   for(i=k;i<n;i++)//menyaem mestami stolbiki
       {tryam=d[i][k];//tryam - vrem peremennaya
        d[i][k]=d[i][jmax];
        d[i][jmax]=tryam;
        }
 
for(i=k+1;i<n;i++)
 {vrem=d[i][k]/d[k][k];
  for(j=k;j<2*n;j++) d[i][j]-=vrem*d[i][j];
   }
 }
for(int k=n-1;k>0;k--)
{for(i=k-1;i>=0;i--)
 vrem=d[i][k]/d[k][k];
 for(j=k;j<2*n;j++)d[i][j]-=vrem*d[i][j];
}
 
 
for(i=0;i<n;i++)
 {for(j=n;j<2*n;j++)
  d[i][j]=d[i][j]/d[i][i];
    d[i][i]=1;
  }
for(i=0;i<n;i++)
 {for(j=n;j<2*n;j++)
  cout<<d[p[i]][j]<<" ";
   cout<<endl;
 }
}
 
 
int main ()
{int i,j;
 cout<<"vvedite razmer matrici"<<endl;
 cin>>n;
 double **a;
 a=new double*[n];
for (i=0;i<n;i++) a[i]=new double[n];//i-stroki,j-stolbci
for(i=0;i<n;i++) 
 {for(j=0;j<n;j++)
   {cout<<" vvedite elementi matrici a["<<i<<"]["<<j<<"]"<<endl;
            cin>>a[i][j];
   }
  
 }
     for (i=0;i<n;i++)//vivod matrici
    {
        for(j=0;j<n;j++)
            cout<<a[i][j]<<" ";
    }
obrat(a);
getchar ();
getchar ();
return 0;
}
Добавлено через 1 минуту
помогите пожалуйста разобраться с кодом, проблема заключается в этом цикле:
C++
1
2
3
4
5
for(int k=n-1;k>0;k--)
{for(i=k-1;i>=0;i--)
vrem=d[i][k]/d[k][k];
for(j=k;j<2*n;j++)d[i][j]-=vrem*d[i][j];
}
Лучшие ответы (1)
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
-=ЮрА=-
Заблокирован
Автор FAQ
17.11.2009, 11:20     Поиск обратной матрицы методом Гаусса #2
Лови мой проект в вижуале, разберёшся...
Вложения
Тип файла: rar MetodGayssa_visual.rar (21.3 Кб, 1908 просмотров)
прасковья
0 / 0 / 0
Регистрация: 13.09.2009
Сообщений: 3
17.11.2009, 22:15  [ТС]     Поиск обратной матрицы методом Гаусса #3
Спасибо большое!
Quart
0 / 0 / 0
Регистрация: 04.06.2011
Сообщений: 3
10.06.2011, 22:31     Поиск обратной матрицы методом Гаусса #4
Цитата Сообщение от прасковья Посмотреть сообщение
Спасибо большое!
Если вы разобрались, скиньте пожалуйста готовый код программы. А то сам метод у меня получается некорректным ((
Malagos
0 / 0 / 1
Регистрация: 23.02.2012
Сообщений: 4
23.02.2012, 17:27     Поиск обратной матрицы методом Гаусса #5
если можно и мне кинь прогу.
vanomux
3 / 3 / 0
Регистрация: 08.02.2012
Сообщений: 17
27.02.2012, 18:04     Поиск обратной матрицы методом Гаусса #6
можно отсюда скачать http://algorithmlib.org/gauss_inverse или http://algorithmlib.org/lu_inverse смотря какой метод нужен
Malagos
0 / 0 / 1
Регистрация: 23.02.2012
Сообщений: 4
28.02.2012, 18:04     Поиск обратной матрицы методом Гаусса #7
Она у меня что-то не пашет.
-=ЮрА=-
Заблокирован
Автор FAQ
28.02.2012, 18:10     Поиск обратной матрицы методом Гаусса #8
Вот готовый код консольного проекта
Методом Гаусса решить систему n линейных алгебраических уравнений
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
14.03.2012, 20:50     Поиск обратной матрицы методом Гаусса
Еще ссылки по теме:

Определитель матрицы методом Гаусса C++
C++ Решения матрицы методом Гаусса
Решение СЛАУ методом обратной матрицы C++

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

Или воспользуйтесь поиском по форуму:
Malagos
0 / 0 / 1
Регистрация: 23.02.2012
Сообщений: 4
14.03.2012, 20:50     Поиск обратной матрицы методом Гаусса #9
Сообщение было отмечено автором темы, экспертом или модератором как ответ
вот в принцепи, но она у меня зацыкливается. В чем проблема подскажите кто увидел.

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
#include<iostream.h>
#include<stdio.h>
#include<conio.h>
#include<stdlib.h>
int n,i,j, A[2][4]; int k, j1; float S, b[20][20], x[20][20];
 
void vvod()
{
cout<<"vvedite rozmirnist matricy A";
cin>>n;
for(int i=0;i<=n-1;i++)
for(int j=0;j<=n-1;j++)
{cout<<"\n A["<<i+1<<", "<<j+1<<"]";
cin>>A[i][j];
}
}
 
void obern ()
{
for(k=0;k<n;k++)
{if(A[k][k]==0)
{cout<<"Oshibka vvedenih dannih";
exit(1);
}
for(j=0;j<2*n;j++)
b[k][j]=A[k][j]/A[k][k];
for(i=k+1;j<n;i++)
for(j=k+1;j<2*n;j++)
A[i][j]=A[i][j]-A[i][k]*b[k][j];
}
for(j=0;j<n;j++)
x[n-1][j]=b[n-1][j+n];
for(k=n-2;k>=0;k--)
for(j1=0;j1<n;j1++)
{S=0;
for(i=k+1;i<n;i++)
for(j=k+1;j<n;j++)
S+=b[k][j]*x[i][j-n];
x[k][j1]=b[k][j1+n]-S;
}
}
void vuvod()
{for(i=0;n-1;1)
for(j=0;n-1;1)
cout<<x[i][j]<<"\t";
cout<<endl;
}
void main()
{int r;
do
{
clrscr();
cout<<"\n\t\t menu\n\t 1-vvedenna\n\t 2-obchuslenna\n\t 3-vuvod\n\t 4-exit\n\tvkagit regum\n";
cin>>r;
switch(r)
{case 1:{vvod(); break;}
 case 2:{obern(); break;}
 case 3:{vuvod();break;}
 case 4:exit(1);
 default:cout<<"\n ne pravelno, povtor vvod\n";
 }
 }
 while(1);
 }
Yandex
Объявления
14.03.2012, 20:50     Поиск обратной матрицы методом Гаусса
Ответ Создать тему
Опции темы

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