Модераторы: Poseidon
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> [C++] Класс Matrix, Помогите доделать задачу на классы  
:(
    Опции темы
Гениальный123Робот
Дата 28.1.2008, 12:41 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


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

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



Это решение принадлежит уважаемому тов. Rockie. Я нашел это решение в одной из тем.

Мне же необходимо, чтобы данная задача находила обратную матрицуопределитель матрицы и транспонирование
При реализации класса матрицы, четыре операции над матрицами: + - * /
реализовать в виде перегрузки операторов, т.е. чтобы можно было записать:

Matrix A, B, C;
C = A + B;

Оч прошу помочь в реализации остальных функций!

P.S. Буду очень благодарен, если прокомментируете основные моменты программы.

Заранее большое спасибо!
Код

#include <iostream>

using namespace std;

class Matrix
{
  public:
  int** m_value;
  unsigned m_width, m_height;

  void Get();
  Matrix( unsigned width = 0, unsigned height = 0);
  Matrix(const Matrix&);
  Matrix& operator=(const Matrix&);
  Matrix& operator+=(const Matrix&);
  friend const Matrix operator+ (const Matrix&, const Matrix&);
  ~Matrix();
};


const Matrix operator+( const Matrix& lhs, const Matrix& rhs ){
  Matrix obj(lhs);
  return obj+=rhs;
}

Matrix& Matrix::operator+=(const Matrix& obj)
{
 Matrix T(m_width,m_height);

  for (unsigned  i = 0; i < m_width; i++)
    for (unsigned  j = 0; j <m_height; j++)
      m_value[i][j] = m_value[i][j] + obj.m_value[i][j];

 return *this;
}


Matrix::Matrix( unsigned width, unsigned height) : m_width(width), m_height(height)
   {
     m_value = new int*[m_width];
     for ( unsigned i = 0; i < m_width; i++ )
      {
         m_value[i] = new int[m_height];
         for ( unsigned j = 0; j< m_height; j++ )
           m_value[i][j]=2;
      }
   }

   
Matrix::~Matrix()
   {
    for ( unsigned i = 0; i < m_width; i++ )
      delete [] this->m_value[i];
    delete [] this->m_value;
   }

void Matrix::Get()
   {   cout<<endl;
       for ( unsigned i = 0; i < m_width; ++i )
       {
          for ( unsigned k = 0; k < m_height; ++k )
           cout<<m_value[i][k]<<" ";
           cout<<'\n';
       }
   }


Matrix& Matrix::operator=(const Matrix& obj)
    {
      if(this == &obj)  return *this;

      for(unsigned i=0;i<this->m_width;i++)
      {  for(unsigned j=0;j<this->m_height;j++)
           this->m_value[i][j] = obj.m_value[i][j];
      }
      return *this;
    }


Matrix::Matrix(const Matrix& obj)
{
    // i?iaa?ea ..

     m_width = obj.m_width;
     m_height = obj.m_height;
     m_value = new int*[m_width];

     for ( unsigned i = 0; i < m_width; i++ )
      {  m_value[i] = new int[m_height];
         for ( unsigned j = 0; j< m_height; j++ )
           m_value[i][j]=obj.m_value[i][j];
      }
     m_width = obj.m_width;
     m_height = obj.m_height;

}



int _tmain()
{
      Matrix A(3,3);
      A.m_value[1][1] = 10;
      A.Get();

      Matrix B=A;
      B+=A;
      B.Get();

      Matrix C(3,3);
      C=A+B+B;
      C.Get();

      cin.get();
      return 0;

}

PM MAIL   Вверх
Sergio02
Дата 30.1.2008, 01:10 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


СИпециалист
*


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

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



Это не подойдёт?

Это сообщение отредактировал(а) Sergio02 - 30.1.2008, 01:10

Присоединённый файл ( Кол-во скачиваний: 71 )
Присоединённый файл  matrix.rar 5,22 Kb
PM MAIL   Вверх
Гениальный123Робот
Дата 1.2.2008, 19:49 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


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

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



Большое спасибо, Sergio02! Ты мой спаситель. Оч полезную штуку скинул!
Но тут еще есть нюансы..
Подобавлял я этот код, но разобраться толком не смог. Ошибки не знаю как там исправить. Помогите плиз отладить.
Добавил 3 те операции, что наверху в начале топика подчеркнуты. Их надо отладить...
Спасибо заранее!  


Код

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

//using namespace std;

class Matrix
{
  public:
  int** m_value;
  unsigned m_width, m_height;

  void Get();
  Matrix( unsigned width = 0, unsigned height = 0);
  Matrix(const Matrix&);
  Matrix& operator=(const Matrix&);
  Matrix& operator+=(const Matrix&);
  friend const Matrix operator+ (const Matrix&, const Matrix&);
  ~Matrix();
};



CMatrix& CMatrix::transpon() {
  if (!width || !height) return *this;
  CMatrix &m = *this;
  CMatrix m2(height, width);
  for (unsigned i= 0; i < width; i++)
    for (unsigned j = i + 1; j < height; j++)
      m2[j][i] = m[i][j];
  m = m2;
  return *this;
};


tableitem CMatrix::determinant() const { 
/* TODO : Обработка ошибки */
  if (width != height) return 0;
  if (width == 0) return 0;

  CMatrix m = *this;
  if (width == 1) return m[0][0];

  if (width == 2)
    return (m[0][0] * m[1][1] - m[1][0] * m[0][1]);

// смена знака определителя. по умолчанию - нет
  bool sign = false;
// определитель
  tableitem det = 1;
  tableitem tmp;
  unsigned x, y;
// цикл по всей главной диагонали
  for(unsigned i = 0; i < height; i++) {
  // если элемент на диагонали равен 0, то ищем не нулевой элемент в матрице
    if (m[i][i] == 0) {
    // если все элементы нулевые, то опр. = 0
      if (!m.search(0, false, y, x, i, i)) return 0;
    // меняем i-ую строку с y-ой
      if (i != y) {
        m.swaprows(i, y);
        sign = !sign;
      }
    // меняем i-ый столбец с x-ым
      if (i != x) {
        m.swapcolumns(i, x);
        sign = !sign;
      }
    // таким образом, в m[i][i], теперь ненулевой элемент.
    }

  // выносим элемент m[i][i] за определитель
    det *= m[i][i];
    tmp = m[i][i];
    for (x = i; x < width; x++) {
      m[i][x] = m[i][x] / tmp;
    }
  // таким образом m[i][i] теперь равен 1
  // зануляем все элементы стоящие под (i, i)-ым,
  // при помощи вычитания с опр. коеффициентом
    for (y = i + 1; y < height; y++) {
      tmp = m[y][i];
      for (x = i; x < width; x++)
        m[y][x] -= (m[i][x]*tmp);
    }
  }
  if (sign) return det*(-1);
  return det;
}; // CMatrix::determinant()



unsigned CMatrix::rang() const {
  if (!width || !height) return 0;

  CMatrix m = *this;
  if (width == 1) return 1;

  tableitem tmp;
  unsigned i;
  unsigned x, y;
// выбираем минимум
  unsigned endi = (width > height)? height : width;
// цикл по всей главной диагонали
  for(i = 0; i < endi; i++) {
  // если элемент на диагонали равен 0, то ищем не нулевой элемент в матрице
    if (m[i][i] == 0) {
    // если все элементы нулевые, прерываем цикл
      if (!m.search(0, false, y, x, i, i)) break;
    // меняем i-ую строку с y-ой
      if (i != y) m.swaprows(i, y);
    // меняем i-ый столбец с x-ым
      if (i != x) m.swapcolumns(i, x);
    // таким образом, в m[i][i], теперь ненулевой элемент.
    }

  // выносим элемент m[i][i]
    tmp = m[i][i];
    for (x = i; x < width; x++) {
      m[i][x] = m[i][x] / tmp;
    }
  // таким образом m[i][i] теперь равен 1
  // зануляем все элементы стоящие под (i, i)-ым и справа от него,
  // при помощи вычитания с опр. коеффициентом
    for (y = i + 1; y < height; y++) {
      tmp = m[y][i];
      for (x = i; x < width; x++)
        m[y][x] -= (m[i][x]*tmp);
    }
    for (x = i + 1; x < width; x++) {
      tmp = m[i][x];
      for (y = i; y < height; y++)
        m[y][x] -= (m[y][i]*tmp);
    }
  }


const Matrix operator+( const Matrix& lhs, const Matrix& rhs ){
  Matrix obj(lhs);
  return obj+=rhs;
}

Matrix& Matrix::operator+=(const Matrix& obj)
{
 Matrix T(m_width,m_height);

  for (unsigned  i = 0; i < m_width; i++)
    for (unsigned  j = 0; j <m_height; j++)
      m_value[i][j] = m_value[i][j] + obj.m_value[i][j];

 return *this;
}


Matrix::Matrix( unsigned width, unsigned height) : m_width(width), m_height(height)
   {
     m_value = new int*[m_width];
     for ( unsigned i = 0; i < m_width; i++ )
      {
         m_value[i] = new int[m_height];
         for ( unsigned j = 0; j< m_height; j++ )
           m_value[i][j]=2;
      }
   }

   
Matrix::~Matrix()
   {
    for ( unsigned i = 0; i < m_width; i++ )
      delete [] this->m_value[i];
    delete [] this->m_value;
   }

void Matrix::Get()
   {   cout<<endl;
       for ( unsigned i = 0; i < m_width; ++i )
       {
      for ( unsigned k = 0; k < m_height; ++k )
           cout<<m_value[i][k]<<" ";
           cout<<'\n';
       }
   }


Matrix& Matrix::operator=(const Matrix& obj)
    {
      if(this == &obj)  return *this;

      for(unsigned i=0;i<this->m_width;i++)
      {  for(unsigned j=0;j<this->m_height;j++)
           this->m_value[i][j] = obj.m_value[i][j];
      }
      return *this;
    }


Matrix::Matrix(const Matrix& obj)
{
    // i?iaa?ea ..

     m_width = obj.m_width;
     m_height = obj.m_height;
     m_value = new int*[m_width];

     for ( unsigned i = 0; i < m_width; i++ )
      {  m_value[i] = new int[m_height];
         for ( unsigned j = 0; j< m_height; j++ )
           m_value[i][j]=obj.m_value[i][j];
      }
     m_width = obj.m_width;
     m_height = obj.m_height;

}



int main()
{
      Matrix A(3,3);
      A.m_value[1][1] = 10;
      A.Get();

      Matrix B=A;
      B+=A;
      B.Get();

      Matrix C(3,3);
      C=A+B+B;
      C.Get();

      cin.get();
      return 0;


}

PM MAIL   Вверх
JackYF
Дата 10.2.2008, 22:43 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


полуавантюрист
****


Профиль
Группа: Участник
Сообщений: 5814
Регистрация: 28.8.2004
Где: страна тысячи озё р

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



Навскидку не понял, почему в одних местах Matrix, в других CMatrix.
Если ты в функции создаёшь новую матрицу, то именно её надо и возвращать, а не *this.

Отладка чужого кода - медвежья помощь. Потрать лучше время, научись пользоваться отладчиком.


--------------------
Пожаловаться на меня как модератора можно здесь.
PM MAIL Jabber   Вверх
Гениальный123Робот
Дата 11.2.2008, 18:28 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


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

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



Тут дело в том, что прога работала до добавления этих 3 функций. Ониподписаны. Я не смог вставить куда надо эти алгоритмы и отладить. В итоге остановился на том, что просто их вставил - транспонир., обратную матрицу и нахождение определителя. Вот их надо доработать в этой проге. Помогите плиз.
PM MAIL   Вверх
Гениальный123Робот
Дата 18.2.2008, 10:50 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


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

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



Так-с, ну тут вопрос решился. Только не могу понять, почему не подключается инклюдом файл "matrix.h". Посмотрите, плиз, прикрепленный файл. Запускать - programm.cpp 

Присоединённый файл ( Кол-во скачиваний: 74 )
Присоединённый файл  program.rar 5,29 Kb
PM MAIL   Вверх
phxOne
Дата 26.2.2008, 11:52 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

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



хм... а matrix.h и matrix.cpp  лежат в той же папке, что и program.cpp? если компилятор просто не находит эти файлы, можно попробовать указать путь к ним в двойных кавычках, типа #include "\path\to\files\matrix.h". и что уточни, что именно выдает компилятор/линкер когда ругается.
PM MAIL   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
Правила форума "Центр помощи"

ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Несоблюдение правил может повлечь за собой самые строгие меры от закрытия/удаления темы до бана пользователя!


  • Название темы должно отражать её суть! (Не следует добавлять туда слова "помогите", "срочно" и т.п.)
  • При создании темы, первым делом в квадратных скобках укажите область, из которой исходит вопрос (язык, дисциплина, диплом). Пример: [C++].
  • В названии темы не нужно указывать происхождение задачи (например "школьная задача", "задача из учебника" и т.п.), не нужно указывать ее сложность ("простая задача", "легкий вопрос" и т.п.). Все это можно писать в тексте самой задачи.
  • Если Вы ошиблись при вводе названия темы, отправьте письмо любому из модераторов раздела (через личные сообщения или report).
  • Для подсветки кода пользуйтесь тегами [code][/code] (выделяйте код и нажимаете на кнопку "Код"). Не забывайте выбирать при этом соответствующий язык.
  • Помните: один топик - один вопрос!
  • В данном разделе запрещено поднимать темы, т.е. при отсутствии ответов на Ваш вопрос добавлять новые ответы к теме, тем самым поднимая тему на верх списка.
  • Если вы хотите, чтобы вашу проблему решили при помощи определенного алгоритма, то не забудьте описать его!
  • Если вопрос решён, то воспользуйтесь ссылкой "Пометить как решённый", которая находится под кнопками создания темы или специальным флажком при ответе.

Более подробно с правилами данного раздела Вы можете ознакомится в этой теме.

Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman

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


 




[ Время генерации скрипта: 0.0853 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.