Реализуйте аппаратно-эффективное решение матрицы реальных пакетов с использованием QR-декомпозиции без Q

Этот пример показывает, как реализовать аппаратно-эффективное решение вещественного матричного уравнения A 'AX = B с использованием решения Real Burst Matrix Solve с использованием блока Q-less QR Decomposition.

Задайте матричные размерности

Укажите количество строк в матрице А, количество столбцов в матрице А и строк В и количество столбцов в матрице B.

m = 100; % Number of rows in A
n = 10;  % Number of columns in A and rows in B
p = 1;   % Number of columns in B

Сгенерируйте матрицы

В данном примере используйте функцию helper realRandomQlessQRMatrices чтобы сгенерировать случайные матрицы А и B для задачи A 'AX = B. Матрицы сгенерированы таким образом, что элементы массива и B находятся между -1 и + 1, а A является полным рангом.

rng('default')
[A,B] = fixed.example.realRandomQlessQRMatrices(m,n,p);

Выбор типов данных с фиксированной точкой

Используйте функцию helper realQlessQRMatrixSolveFixedpointTypes выбрать типы данных с фиксированной точкой для входных матриц А и B и вывести X так, что существует низкая вероятность переполнения во время расчета. Дополнительные сведения о том, как выбираются типы данных, см. в документе FixedPointMatrixLibraryDatypesExample.pdf в текущей директории.

max_abs_A = 1;  % max(abs(A(:))
max_abs_B = 1;  % max(abs(B(:))
f = 24;         % Fraction length (bits of precision)
T = fixed.example.realQlessQRMatrixSolveFixedpointTypes(m,n,max_abs_A,max_abs_B,f);
A = cast(A,'like',T.A);
B = cast(B,'like',T.B);
OutputType = fixed.extractNumericType(T.X);

Откройте модель

model = 'RealBurstQlessQRMatrixSolveModel';
open_system(model);

Подсистема Обработчика Данных в этой модели принимает действительные матрицы А и B как входы. The ready порт запускает обработчик данных. После отправки истинного validIn сигнал, может быть некоторая задержка перед ready задано значение false. Когда Обработчик данных обнаруживает переднее ребро ready сигнал, блок устанавливает validIn true и отправляет следующую строку A и B. Этот протокол позволяет отправлять данные всякий раз, когда переднее ребро ready обнаруживается сигнал, обеспечивающий обработку всех данных.

Установите переменные в рабочем пространстве модели

Используйте функцию helper setModelWorkspace чтобы добавить переменные, определенные выше, в рабочее пространство модели. Эти переменные соответствуют параметрам блоков для решения Real Burst Matrix Solve с использованием блока Q-less QR Decomposition.

numSamples = 1; % Number of samples
fixed.example.setModelWorkspace(model,'A',A,'B',B,'m',m,'n',n,'p',p,...
    'numSamples',numSamples,'OutputType',OutputType);

Симулируйте модель

out = sim(model);

Построение решения из выходных данных

Решение матрицы Real Burst с использованием блока Q-less QR Decomposition выводит данные по одной строке за раз. Когда выводится строка результата, блок устанавливает validOut к true. Строки X выводятся в том же порядке, в котором они вычисляются, первая последняя строка, поэтому необходимо восстановить данные, чтобы интерпретировать результаты. Чтобы восстановить матрицу X из выхода данных, используйте функцию helper matrixSolveModelOutputToArray.

X = fixed.example.matrixSolveModelOutputToArray(out.X,n,p);

Проверьте точность выхода

Чтобы оценить точность решения Real Burst Matrix Solve с использованием блока Q-less QR Decomposition, вычислите относительную погрешность.

relative_error = norm(double(A'*A*X - B))/norm(double(B)) %#ok<NOPTS>
relative_error =

   5.3565e-04

Для просмотра документации необходимо авторизоваться на сайте