Есть такая штука: Класс DMatrix .

В двух словах: там сделана перегрузка операторов для матричных вычислений, то есть:

DMatrix A, B, C;

<задаются значения матриц A и B>

C = A * B;

Ну и всякие функции типа обращения матриц и решения систем уравнений.


Вопрос: существуют ли подобные классы на C++ с открытым кодом, чтобы их можно было подключить к проекту и использовать в задачах с большим количеством операций с матрицами? Если - да, то где их можно скачать?

White Owl http://g.zeos.in/?q=c matrix library

Не, так не пойдет... Хочу совет типа "Я пользовался таким-то классом и мне понравилось, потому что..."

NekZOleg77777,

Ну вот от нашей конторы либа , сами написали, сами пользуемся. Работа с матрицами. Правда, либа сама для обработки изображений.

Да, забавно, очень специфично. Такая куча функций! А как все это используется? Может быть, код заточен на прошивку в микросхемы?

BagaBaga lapack C++


Ага, в общем, смотреть здесь: netlib.org/clapack .

А для объектно сориентированных, вроде меня, здесь: math.nist.gov/lapack++

Будем читать.

Теперь такой вопрос, к спецам по мат. методам вычислений.
Вот у меня в компьютер поступает поток информации, в каждый тик времени из значений параметров формируется новое уравнение, то есть строка матрицы + свободный член.
А я хочу, чтобы это было дифференциальное уравнение. То есть каждый раз я считаю все нужные мне производные и помещаю их значения в матрицы.

Вопрос в том, как правильно считать дискретные производные. Если - втупую, получается вот что:

X`(i) = (X(i) - X(i-1)) / dt - первая производная

X``(i) = (X`(i) - X`(i-1)) / (2*dt) = (((X(i) - X(i-1)) / dt) - ((X(i-1) - X(i-2)) / dt)) / (2*dt) =
= (X(i) - X(i-2)) / (2*dt*dt) - вторая производная и т.д.,

то есть взятие n-ой производной получается эквивалентно взятию первой производной на менее подробной решетке, то есть с меньшей частотой дискретизации. Что-то я сомневаюсь, что это правильно.

Может, кто подскажет правильный подход, что почитать и т.д.?