Слияние двух матриц ... в LISP - PullRequest
2 голосов
/ 27 мая 2010
(defun merge-matrix (matrix-1 matrix-2)
    (if (not (or (eql (matrix-rows matrix-1) (matrix-rows matrix-2)) (null matrix-1) (null matrix-2))) (error "Invalid dimensions."))
    (cond
        ((null matrix-1) (copy-tree matrix-2))
        ((null matrix-2) (copy-tree matrix-1))
        (t (let ((result (copy-tree matrix-1)))
                 (dotimes (i (matrix-rows matrix-1))
                     (setf (nth i result) (nconc (nth i result) (nth i matrix-2))))
                 result))))

(матрица слияния '((3 1) (1 3))' ((4 2) (1 1)))

*** - EVAL: переменная NULL не имеет значения

Я получаю такую ​​ошибку, как я могу решить проблему, спасибо

Ответы [ 2 ]

1 голос
/ 02 июля 2011

Код ОП работает у меня. Тем не менее, я чувствовал желание улучшить его и Я реализовал ту же идею (но немного более мощный).

Семантика такая же, как у vertcat Матлаба. Функция добавляет все аргументы в одну большую матрицу.

Обратите внимание, что из-за объявлений мой код должен быть очень эффективным.

(deftype mat ()
  "Non-square matrices. Last index is columns, i.e. row-major order."
  `(simple-array single-float 2))

(defun are-all-elements-typep (type ls)
  (reduce #'(lambda (b x) (and b (typep x type)))
      ls))

(defun are-all-matrix-heights-equalp (ls)
  (let ((first-height (array-dimension (first ls) 0)))
    (reduce #'(lambda (b x) (and b
                (= first-height
                   (array-dimension x 0))))
       ls)))

(defun vertcat (&rest rest)
  (declare (type cons rest))
  (unless (are-all-elements-typep 'mat rest)
    (break "At least one of the arguments isn't a matrix."))
  (unless (are-all-matrix-heights-equalp rest)
    (break "All Matrices must have the same number of rows."))
  (let* ((height (array-dimension (first rest) 0))
     (widths (mapcar #'(lambda (mat) (array-dimension mat 1)) rest))
     (result (make-array (list height
                   (reduce #'+ widths))
                 :element-type 'single-float))
     (current-width 0))
    (dotimes (m (length rest))
      (let ((e (elt rest m)))
    (destructuring-bind (y x) (array-dimensions e)
     (dotimes (j y)
       (dotimes (i x)
         (setf (aref result j (+ current-width i))
           (aref e j i))))
     (incf current-width (elt widths m)))))
    (the mat result)))

#+nil
(let ((a (make-array '(2 3)
             :initial-contents '((1s0 2s0 3s0)
                     (2s0 4s0 5s0))
             :element-type 'single-float))
      (b (make-array '(2 2)
             :initial-contents '((6s0 7s0)
                     (9s0 8s0))
             :element-type 'single-float)))
  (vertcat a b a))
;=> #2A ((1.0 2.0 3.0 6.0 7.0 1.0 2.0 3.0) (2.0 4.0 5.0 9.0 8.0 2.0 4.0 5.0))
0 голосов
/ 29 января 2011

Сообщение об ошибке, которое вы получаете, говорит о том, что lisp пытается обработать один из ваших вызовов null как переменную. Мне удалось воспроизвести это поведение, определив matrix-rows, как это сделал Фрэнк Шиарар, и удалив, например, скобки вокруг ((null matrix-1) (copy-tree matrix-2)) s-выражения. Я бы посоветовал вам проверить свои круглые скобки, либо вручную, либо используя что-то вроде SLIME , что дало мне предупреждение при попытке компилировать функцию.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...