ANTICHAT.XYZ    VIDEO.ANTICHAT.XYZ    НОВЫЕ СООБЩЕНИЯ    ФОРУМ  
Баннер 1   Баннер 2
Antichat снова доступен.
Форум Antichat (Античат) возвращается и снова открыт для пользователей. Здесь обсуждаются безопасность, программирование, технологии и многое другое. Сообщество снова собирается вместе.
Новый адрес: forum.antichat.xyz
Вернуться   Форум АНТИЧАТ > Программирование > С/С++, C#, Delphi, .NET, Asm
   
Ответ
 
Опции темы Поиск в этой теме Опции просмотра

Методы Гаусса и Крамера[c++]
  #1  
Старый 02.11.2007, 17:04
Ky3bMu4
Постоянный
Регистрация: 03.02.2007
Сообщений: 520
Провел на форуме:
1777536

Репутация: 932


Отправить сообщение для Ky3bMu4 с помощью ICQ
По умолчанию Методы Гаусса и Крамера[c++]

Как-то на кряклабе нашёл несколько примеров по решению систем линейных уравнений этим методом. Но 404. . Решил написать сам.
Для привидения матрицы к главному ступенчатому виду использовал метод ведущего элемента(для машины он проще чем сначала получение верхне треугольной матрицы , а потом самого вида.)
Т.к. часто попадаются "плохие" числа и приходится работать с дробями, то написал класс для дробей, перегрузил операторы *, / , + , - . Сделал функции 1, 2 и 3 элементарных преобразований.
Также есть алгоритм оброботки нулевых эллемтов на главной диагонале.Для этого написал отдельно класс для матриц.
Ввод матрицы осуществляется в самой программе. Программа таже показывает пошагово значения матриц.

Приводить весь алгоритм небуду, т.к. он очень большой. Исходники и откомпиленная прога ниже(вместе с Крамером):

http://slil.ru/25072002
http://ifolder.ru/4038710
http://rapidshare.com/files/68056699/Gauss_Kramer.rar.html

Последний раз редактировалось Ky3bMu4; 07.11.2007 в 17:27..
 
Ответить с цитированием

  #2  
Старый 02.11.2007, 20:09
Ni0x
Постоянный
Регистрация: 27.08.2006
Сообщений: 367
Провел на форуме:
2009677

Репутация: 472
По умолчанию

Цитата:
сразу извиняюсь за свой ужасный английский
Ты бы лучше за свой ужасный С++ извинился. Что уж говорить, если ты не знаешь как через cout переменную вывести. Использовать код вида:
Код:
printf("Insert %d element of second matrix:",x);
cin>>temp;
cout<<"\r\n";
это норм по твоему? Но больше всего конечно порадовало подключение windows.h ради функции ExitProcess. Ты только не подумай, что я опять наезжаю на тебя или еще чего.
 
Ответить с цитированием

  #3  
Старый 02.11.2007, 22:16
KEZ
Banned
Регистрация: 18.05.2005
Сообщений: 1,981
Провел на форуме:
1941233

Репутация: 2726


По умолчанию

Ni0x, ты че, это кроссплатформеный Опен-Сурс от Кузьмича.
А вот это вообще будет висеть у меня на стенке в корридоре

Цитата:
void echo () const{
printf("%d/%d",chi,zn);
}
Ещё. Крутой деструктор
Цитата:
/*деструктор */
~drob() {}
)))))

Ну и самый прикол
Цитата:
drob::chi*mno.chi
на всякий случай кодер решил везде написать (по всему тексту прямо) "drob::" чтоб вдруг что....гыгы
 
Ответить с цитированием

  #4  
Старый 02.11.2007, 22:38
Constantine
Members of Antichat - Level 5
Регистрация: 24.11.2006
Сообщений: 927
Провел на форуме:
7192869

Репутация: 3033


По умолчанию

Я лучше столбиком буду делить...
__________________

Дети индиго - это бездари, не надо песен! В пять лет едва говорить начинают, мы в этом возрасте стихи наизусть читали!

з.ы http://www.youtube.com/watch?v=sNsQe0KByRY Я ПлакалЪ
 
Ответить с цитированием

  #5  
Старый 02.11.2007, 23:31
alextoun
Постоянный
Регистрация: 07.05.2006
Сообщений: 732
Провел на форуме:
7910701

Репутация: 811


Отправить сообщение для alextoun с помощью ICQ
По умолчанию

не удивительно, шо у меня ошибки при компеляци))
кстате я метод гаус на 5 сдал!!! имея 3 по математике в школе по жизни
 
Ответить с цитированием

  #6  
Старый 02.11.2007, 23:43
.Slip
Leaders of Antichat - Level 4
Регистрация: 16.01.2006
Сообщений: 1,966
Провел на форуме:
21768337

Репутация: 3486


По умолчанию

Идея норм, сурс ужасен=\
__________________
..
 
Ответить с цитированием

  #7  
Старый 03.11.2007, 16:42
Ky3bMu4
Постоянный
Регистрация: 03.02.2007
Сообщений: 520
Провел на форуме:
1777536

Репутация: 932


Отправить сообщение для Ky3bMu4 с помощью ICQ
По умолчанию

Цитата:
не удивительно, шо у меня ошибки при компеляци))
Ну и какие-же это ошибки?

Цитата:
Что уж говорить, если ты не знаешь как через cout переменную вывести.
Так:
Код:
printf("Insert %d element of second matrix:",x);
Намного удобнее и короче, чем:
Код:
cout<<"Insert  "<<x<< " element of second matrix:";
Неправдали?
 
Ответить с цитированием

  #8  
Старый 07.11.2007, 17:17
Ky3bMu4
Постоянный
Регистрация: 03.02.2007
Сообщений: 520
Провел на форуме:
1777536

Репутация: 932


Отправить сообщение для Ky3bMu4 с помощью ICQ
По умолчанию

Большое спасибо за критику. Учёл.
Улучшил Гаусса, добавил обработку нулевых элементов на главной диагонале.

Вот, ещё написал метод Крамера с рекурсивным вычислением определителя.
Код:
#include <tchar.h>
#include <iostream>

using namespace std;

/*функция возведения в определённую степень числа -1*/
int mypow(int x){ 
	int y=1;
	for(int z=1;z<=x;z++){
  y*=-1;
	}
	return y;
}

/*класс матрица*/
class matrix{
public:

/*Конструктор: все элементы 0, размерность 0 */
 matrix(){
  size=0;
  
  for(int x=0;x<=100;x++)
  	for(int y=0;y<=100;y++)
    m[x][y];
 }
/*Рекурсивная функция вычисления определителя */
 int GetDet(){
	int result=0;
  if(size==1){ //Если у матрицы размерность 2(математическая, С`шная 1), то вычисляем определитель
return m[0][0]*m[1][1]-m[0][1]*m[1][0];
  }
  else{
/*Вычисляем определитель путём вычёркивания строк и столбцов. Берём элементы из первой строки. */   
   for(int x=0;x<=size;x++){
result+= mypow(x+2)*m[0][x]*(GetNew(x)).GetDet(); //рекурсия...
   }
  return result;
  }
	}

 /*Функия получения новой матрицы из старой, путём вычёркивания первой(нулевой) строки и х`овой колонки*/
 matrix GetNew(int st){
  matrix result;
  int z=0;
  result.size=size-1;

  for(int x=1;x<=size;x++){
   for(int y=0;y<=size;y++){
  if(y!=st){
   result.m[x-1][z]=m[x][y];
   z++;
  }
  
   }
   z=0;
  }
  return result;
 }

/*Функция вывода на экран квадратной матрицы*/
  void echo(void){
   for(int x=0;x<=size;x++){
    for(int y=0; y<=size;y++){
     cout<<m[x][y]<<" ";
    }
    cout<<"\r\n";
   }
  }
/*Функция вывода на экран матрицы-столбца */
  void echoEx(void){
    for(int x=0;x<=size;x++){
     cout<<m[x][0]<<"\r\n";
    }
   }

/*Функция установки размерности матрицы*/
   inline void SetSize(int r){
    size=r;
   }

/*Функция установки значения элемента*/
   inline void SetValue(int value, int x , int y){
    m[x][y]=value;
   }

   /*Функция замены колонки на матриуц-столбец*/
   inline void SetColumn(int column , matrix source){    
    for(int x=0;x<=size;x++){
     m[x][column]=source.m[x][0];
    }
   }
 
private:
	int m[100][100]; //сама матрица
	int size; //размерность
	
};
int _tmain(int argc, _TCHAR* argv[])
{
matrix first,second,temp; //наши матрицы
int element,size,maindet; //необходимые переменные
bool c=1; //для цикла
cout<<"The Kramer method for the system of linear equations. By Ky3bMu4 (c) 2007.\r\n\r\n";
while(c){
cout<<"\r\nInsert size of matrix:"; 
/*тут небольшая неровность: в математике матрица начинается с 1 столбцы, а в С - с 0-го   */
cin>>size;
first.SetSize(size);
second.SetSize(size);


/* вводим основную матрицу */
for(int x =0; x<=size; x++){
	for( int y=0; y<=size;y++){
  cout<<"\r\nInsert ["<<x<<";"<<y<<"] element of first matrix:";
  cin>>element;
  first.SetValue(element,x,y);
  }
}

/* вводим матрицу-столбец с числами*/
for(int x=0; x<=size; x++){
	cout<<"\r\nInsert "<<x<<" element of second matrix:";
	cin>>element;
	second.SetValue(element,x,0);
}
maindet = first.GetDet(); //вычисляем основной определитель(на него будем делить)
if(maindet!=0){
for(int x=0;x<=size;x++){ //кол-во элементов(переменных в линейном уравнении) == размерности главной матрицы.
//Заменяем колонку, вычисляем определитель , делим на главный определитель, пробуем сократить и выводим 
temp = first;
temp.SetColumn(x,second); //заменянем колонку
cout<<"\r\nElement "<<x<<" = "<<temp.GetDet()<<"/"<<maindet; // выводим на экран	
}
}else{ //облом, главный определитель == 0, Крамер не прокатит :(((
	cout<<"\r\nSorry, but the main determinat of first matrix is 0.Stop.";
}

// The End
cout<<"\r\n\r\nFinish!";
cout<<"\r\nDo you want to again calculate new matrix(1/0)?:";
cin>>c;
}

	return 0;
}

Ссылка на откомпиленную прогу в первом посте.
 
Ответить с цитированием

  #9  
Старый 07.11.2007, 17:30
ak[id]
Участник форума
Регистрация: 22.06.2007
Сообщений: 201
Провел на форуме:
1693011

Репутация: 248
По умолчанию

Спасибо, мы ща в уневере как раз ща решаем матрицы методом Гауса и Крамера.
 
Ответить с цитированием

  #10  
Старый 07.11.2007, 23:15
GlOFF
Постоянный
Регистрация: 08.05.2006
Сообщений: 816
Провел на форуме:
1845671

Репутация: 1338


По умолчанию

Да бред все это! Используйте метод Гаусаса.
 
Ответить с цитированием
Ответ



Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
Мануал к SQL Injection TooL 3 SQLHACK Авторские статьи 21 09.04.2009 15:34
СИ. Метод №2. TALIB ICQ 9 28.01.2007 19:52
Метод шифрования в Web Wiz Forums ysmat Расшифровка хешей 0 14.08.2006 21:06
Сетевой сканер Nmap. Руководство пользователя foreva Чужие Статьи 1 08.02.2005 16:36



Здесь присутствуют: 1 (пользователей: 0 , гостей: 1)
 


Быстрый переход




ANTICHAT.XYZ