2
голосов
3ответов
1167 просмотров

finding matrix through optimisation

I am looking for algorithm to solve the following problem : I have two sets of vectors, and I want to find the matrix that best approximate the transformation from the input vectors to the output vectors. vectors are 3x1, so matrix is 3x3. This is the general problem. My particular problem is ...

160
голосов
9ответов
121992 просмотров

"Cloning" row or column vectors

Sometimes it is useful to "clone" a row or column vector to a matrix. By cloning I mean converting a row vector such as [1,2,3] Into a matrix [[1,2,3] [1,2,3] [1,2,3] ] or a column vector such as [1 2 3 ] into [[1,1,1] [2,2,2] [3,3,3] ] In matlab or octave this is done pretty ea...

2
голосов
1ответов
677 просмотров

Rotation about the Y-Axis

Doing a couple rotations in Matlab, one which is rotation about the y-axis, however online I've found two different answers: here and here. Which is correct, if both how does one get to the other?

9
голосов
4ответов
17840 просмотров

Solving a matrix in MATLAB?

How does one solve the (non-trivial) solution Ax = 0 for x in MATLAB? A = matrix x = matrix trying to solve for I've tried solve('A * x = 0', 'x') but I only get 0 for an answer.

6
голосов
6ответов
4093 просмотров

C++ libraries for Image Segmentation

I am going to do a project in Data Mining related to image clustering (in C++) .I am looking for a powerful library which is helpful in image processing, linear algebra and 3d graphics. Any thoughts? Thanks.

16
голосов
3ответов
7607 просмотров

Why is Matlab's inv slow and inaccurate?

I read at a few places (in the doc and in this blog post : http://blogs.mathworks.com/loren/2007/05/16/purpose-of-inv/ ) that the use of inv in Matlab is not recommended because it is slow and inaccurate. I am trying to find the reason of this inaccuracy. As of now, Google did not give m interes...

1
голосов
1ответов
1984 просмотров

Sparse linear algebra solvers for C#

I'm working on an experimental implementation of Goldenthal et.al's inextensible cloth algorithm in C#. First I used Math.NET Iridium to assemble and solve the matrices, but quickly replaced this with dnAnalytics since the latter allows me to reuse matrices, almost eliminating further memory al...

246
голосов
11ответов
152121 просмотров

What are the most widely used C++ vector/matrix math/linear algebra libraries, and their cost and benefit tradeoffs?

It seems that many projects slowly come upon a need to do matrix math, and fall into the trap of first building some vector classes and slowly adding in functionality until they get caught building a half-assed custom linear algebra library, and depending on it. I'd like to avoid that while not...

0
голосов
4ответов
7678 просмотров

How do I create an identity matrix of arbitrary size in Java?

Is there a utility to create an identity matrix of specified size in Java?

7
голосов
7ответов
3265 просмотров

Programming a smooth change of thrust from current velocity vector to a target vector

TL;dr: "I am not sure how to calculate a smooth transition of thrust between one vector and another." I am programming a simple game where an enemy chases after the player in an open space (no walls). I was calculating the enemy's x & y velocities independently, accelerating them if they wer...

0
голосов
1ответов
232 просмотров

Algorithm that takes 2 'similar' matrices and 'aligns' one to another

First of all, the title is very bad, due to my lack of a concise vocabulary. I'll try to describe what I'm doing and then ask my question again. Background Info Let's say I have 2 matrices of size n x m, where n is the number of experimental observation vectors, each of length m (the time serie...

3
голосов
8ответов
2838 просмотров

Matrix classes in c++

I'm doing some linear algebra math, and was looking for some really lightweight and simple to use matrix class that could handle different dimensions: 2x2, 2x1, 3x1 and 1x2 basically. I presume such class could be implemented with templates and using some specialization in some cases, for perform...

2
голосов
2ответов
1544 просмотров

Compiling PARDISO linear solver test case with GCC

I am trying to compile a linear system solver using PARDISO. The test case (pardiso_sym.c) also downloaded from the same website above. I have the following files inside the directory: [[email protected] my-pardiso]$ ls -lh total 1.3M -rw-r--r-- 1 gv hgc0746 1.3M Aug 7 11:59 libpardiso_GNU_IA64.so -...

9
голосов
6ответов
9678 просмотров

C++ Memory Efficient Solution for Ax=b Linear Algebra System

I am using Numeric Library Bindings for Boost UBlas to solve a simple linear system. The following works fine, except it is limited to handling matrices A(m x m) for relatively small 'm'. In practice I have a much larger matrix with dimension m= 10^6 (up to 10^7). Is there existing C++ approach ...

0
голосов
3ответов
4567 просмотров

Compiling C++ Code With Boost's Numeric Binding Library to Solve Ax=b Linear System

I am using Numeric Library Bindings for Boost UBlas to solve a simple linear system: #include<boost/numeric/ublas/matrix.hpp> #include<boost/numeric/ublas/io.hpp> #include<boost/numeric/bindings/traits/ublas_matrix.hpp> #include<boost/numeric/bindings/lapack/gesv.hpp> #in...

9
голосов
6ответов
18665 просмотров

Boost's Linear Algebra Solution for y=Ax

Does boost have one? Where A, y and x is a matrix (sparse and can be very large) and vectors respectively. Either y or x can be unknown. I can't seem to find it here: http://www.boost.org/doc/libs/1_39_0/libs/numeric/ublas/doc/index.htm

3
голосов
1ответов
799 просмотров

SVG линейная шкала градиента и вопрос о переводе

У меня следующий радиальный градиент: <radialGradient inkscape:collect="always" xlink:href="#linearGradient2454" id="radialGradient2460" cx="4022.8572" cy="5451.2656" fx="4022.8572" fy="5451.2656" r="3559.865" gradientTransform="matrix(-0.1071067,-0.1166362,0.137...

19
голосов
9ответов
6496 просмотров

Как линейная алгебра используется в алгоритмах?

Некоторые из моих коллег отметили, что «линейная алгебра» очень важна при изучении алгоритмов. Я изучил множество алгоритмов и прошел несколько курсов линейной алгебры, но не вижу связи. Так как же линейная алгебра используется в алгоритмах? Например, что интересного можно сделать с помощью м...

14
голосов
7ответов
11050 просмотров

ublas против библиотеки шаблонов матриц (MTL4)

Я пишу программу для гиперболических уравнений в частных производных на C ++. Практически все обозначения векторные и матричные. Кроме того, мне нужен решатель линейной алгебры. И да, размеры вектора и матрицы могут значительно варьироваться (скажем, от 1000 до размеров, которые могут быть решен...

4
голосов
4ответов
1862 просмотров

как сгенерировать псевдослучайную положительно определенную матрицу с ограничениями на недиагональные элементы?

Пользователь хочет установить уникальную, нетривиальную верхнюю / нижнюю границу корреляции между каждой парой переменных в матрице var / covar. Например: мне нужна матрица дисперсии, в которой все переменные имеют 0,9> | rho (x_i, x_j) |> 0,6, rho (x_i, x_j) - корреляция между переменн...

3
голосов
3ответов
11776 просмотров

как сгенерировать псевдослучайную положительно определенную матрицу с ограничениями на недиагональные элементы?

Возможный дубликат: как сгенерировать псевдослучайный положительно определенная матрица с ограничениями на недиагональные элементы? Пользователь хочет установить уникальную, нетривиальную верхнюю / нижнюю границу корреляции между каждой парой переменных в матрице var / covar. ...

10
голосов
4ответов
5332 просмотров

Функция mrdivide в MATLAB: что она делает и как я могу сделать это в Python?

У меня есть эта строка кода MATLAB: a/b Я использую эти входные данные: a = [1,2,3,4,5,6,7,8,9,1,2,3,4,5,6,7,8,9] b = ones(25, 18) Это результат (матрица 1x25): [5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] Что делает MATLAB? Я пытаюсь воспроизвести это поведение в Python...

7
голосов
5ответов
9482 просмотров

Деление массива - перевод из MATLAB в Python

У меня есть эта строка кода в MATLAB, написанная кем-то другим: c=a.'/b Мне нужно перевести это на Python. a, b и c - все массивы. Размеры, которые я сейчас использую для тестирования кода: a: 18 x 1, b: 25x18, что дает мне c с размерами 1x25. Массивы не квадратные, но я бы не х...

2
голосов
1ответов
1301 просмотров

Есть ли в LAPACK / BLAS или где-то еще подпрограмма Fortran для расчета разложения LDL?

Как сказано в заголовке, мне нужно сформировать разложение холецкого ЛПНП для моей положительно определенной матрицы A (как и нормальный холецкий, но есть матрицы, расположенные на одной диагонали от L, а D - диагональная матрица). Я нашел только одну функцию в Lapack, которая делает это, но в н...

3
голосов
3ответов
262 просмотров

Изменение уравнения с плавающей запятой

Я хотел бы знать, есть ли способ повысить точность вычисления наклона. (Это появилось несколько месяцев назад здесь ) . Вроде по изменению: float get_slope(float dXa, float dXb, float dYa, float dYb) { return (dXa - dXb)/(dYa - dYb); } в float get_slope(float dXa, float dXb, fl...

7
голосов
3ответов
2252 просмотров

Отмена проецирования точки на экране обратно в изометрически спроецированный мир

Я выполняю 3D-симуляцию за кулисами во время рендеринга мира в моем 2D-изометрическом движке. Я никогда раньше не делал изометрический движок, и моя матричная математика в целом ржавая, поэтому у меня проблемы. У меня есть матрица проекции, которая в простейшем виде такова: 0.7 0.35...

2
голосов
3ответов
1909 просмотров

основная справка по линейной алгебре

Это то, что мне нужно сделать - У меня есть это уравнение - Ax = y Где A - рациональная матрица размера m * n (m <= n), а x и y - векторы правильный размер. Я знаю A и y, я не знаю, чему равен x. я также знайте, что не существует x, где Ax равно y. Я хочу найти вектор x 'такой, чтобы...

21
голосов
7ответов
9302 просмотров

Алгоритм поиска ближайших точек?

Учитывая набор из нескольких миллионов точек с координатами x, y, какой алгоритм выбрать для быстрого поиска 1000 самых близких точек из местоположения? «Быстро» здесь означает около 100 мсек на домашнем компьютере. Грубая сила означает выполнение миллионов умножений и их сортировку. Хотя даж...

3
голосов
5ответов
1076 просмотров

линейная регрессия высшего порядка

У меня матричная система: A x B = C A - это a на n, а B - это n на b. И A, и B неизвестны, но у меня есть частичная информация о C (у меня есть некоторые значения, но не все), а n выбрано достаточно маленьким, чтобы система ожидается, что он будет чрезмерно ограничен. Не обязательно, чтобы...

7
голосов
2ответов
7854 просмотров

поэлементные операции с матрицами boost c ++ ublas и векторными типами

Я хотел бы выполнять поэлементные функции для типов матрицы и векторов повышения, например возьмите логарифм каждого элемента, возведите в степень каждый элемент, примените специальные функции, такие как гамма и дигамма и т. д. (аналогично тому, как в Matlab эти функции применяются к матрицам и ...