VPF::Матрица - Форум программистов

Agentx86

Дата 4.5.2006, 19:26 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 110
Регистрация: 27.4.2006

Репутация: нет
Всего: нет

Помогите пожалуйста. Нашел код в книге который находит обратную матрицу:

Код


# include <iostream.h>
# include <math.h>
# include <conio.h>


const int dim=20;
float a [dim][dim];
void InvMatrG (int n, float e, float d)
{
int i,j,k;
float y,w;
float b[dim], c[dim];
int z[dim];
d=1;
for (j=0; j<n; j++)
z[j]=j;
for (i=0; i<n; i++)
{
k=i;
y=a [i][j];
for (j=i+1; j<n; j++)
{
w=a[i][j];
if (fabs (w)>fabs (y))
{
k=j;
y=w;
}
}
d*=y;
if (fabs(y)<e)
{
cout<<" Matrica virozdena\n";
getch ();
}
for (j=0; j<n; j++)
{
c[j]=a[j][k];
a[j][k]=a[j][i];
a[j][i]=-c[j]/y;
a[i][j]=a[i][j]/y;
b[j]=a[i][j];
}
a[i][i]=1/y;
j=z[i];
z[i]=z[k];
z[k]=j;
for(k=0; k<n; k++)
if (k!=i)
for (j=0; j<n; j++)
if (j!=i)
a[k][i]-=b[j]*c[k];
}
for (i=0; i<n; i++)
{
k=z[i];
if(k!=1)
{
for (j=0; j<n; j++)
{
w=a[i][j];
a[i][j]=a[k][j];
a[k][j]=w;
}
j=z[i];
z[i]=z[k];
z[k]=j;
d=-d;
}
}
}
void Matrin (int n)
{
int i,j;
for (i=0; i<n; i++)
for (j=0; j<n; j++)
cin>>a[i][j];
}
void Matrout (int n)
{
int i,j;
for (i=0; i<n; i++)
{cout<<"\n";
for (j=0; j<n; j++)
cout<<a[i][j]<<"\t";

}
}
void main ()
{
int n;
float d,e;
cout<<"\n Vvedite razmernost sitemi n=";
cin>>n;
cout<<"\n Vvedite elementi matrici:   \n";
Matrin(n);
cout<<"\n Tochnost vichisleniu e= ";
cin>>e;
InvMatrG(n,e,d);
cout<<"\n Obratnaja matrica;  ";
Matrout (n);
getch ();
}

Не могу понять какие числа надо вводить в переменную "Точность вычисления". Кто разберется подскажите.

Paradoks

Дата 4.5.2006, 20:04 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 5
Регистрация: 11.3.2006

Репутация: нет
Всего: нет

Код


float d,e;
InvMatrG(n,e,d);

Сначало объявляем переменную, а потом передаём её значение в функцию ничего не присвоив?

Количество элементов матрицы должно быть меньше или равно 20, т.е. нужно написать не

Код


int n;
cout<<"\n Vvedite razmernost sitemi n=";    
cin>>n;

а так

Код


int n = 0;
while((n >0) && (n<=20)){
cout<<"\n Vvedite razmernost sitemi n=";    
cin>>n;
}

Вобщем весь код написан не корректно smile

Это сообщение отредактировал(а) Paradoks - 4.5.2006, 20:12

Agentx86

Дата 5.5.2006, 01:17 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 110
Регистрация: 27.4.2006

Репутация: нет
Всего: нет

Цитата(Paradoks @ 4.5.2006, 20:04)

Код


float d,e;
InvMatrG(n,e,d);

Код


int n;
cout<<"\n Vvedite razmernost sitemi n=";    
cin>>n;

а так

Код


int n = 0;
while((n >0) && (n<=20)){
cout<<"\n Vvedite razmernost sitemi n=";    
cin>>n;
}

Вобщем весь код написан не корректно smile

Абсолютно с тобой согласен. Писал не я. Пример переписан один к одному я еще проверил. Переменная "d" не где не используется в void main(), только объявляется. Тем не менее можешь ответить на мой вопрос. Смысл переменной "е"? И какие значения она может принимать.

bsa

Дата 5.5.2006, 09:12 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Модератор
Сообщений: 9185
Регистрация: 6.4.2006
Где: Москва, Россия

Репутация: 63
Всего: 196

Если честно, алгоритм я не очень понял. Но по смыслу e - это "эпселент", т.е. очень малое положительное число, которое ты готов принять за 0.
Т.е.:
Если |A - B| < e, то считается, что A = B

Это проходят еще на 1-м курсе технического ВУЗа.

Это сообщение отредактировал(а) bsa - 5.5.2006, 09:14

--------------------

Правильно заданный вопрос - половина ответа
Отзывы о работе модератора.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Общие вопросы \| Следующая тема »