fixed.qlessQRMatrixSolve

Решите систему линейных уравнений (A 'A) x = B для использования x разложение Q-less QR

Описание

пример

x = fixed.qlessQRMatrixSolve(A, B) решает систему линейных уравнений (A 'A) x = B с помощью разложения QR, не вычисляя значение Q.

Результат этого кода эквивалентен вычислению

[~,R] = qr(A,0);
x = R\(R'\B)
или
x = (A'*A)\B

пример

x = fixed.qlessQRMatrixSolve(A, B, outputType) возвращает решение системы линейных уравнений (A 'A) x = B как переменная с выходным типом, заданным outputType.

x = fixed.qlessQRMatrixSolve(A, B, outputType, forgettingFactor) возвращает решение системы линейных уравнений, с forgettingFactor умноженный на R после каждой строки A обрабатывается.

Примеры

свернуть все

В этом примере показано, как решить систему линейных уравнений (AA)x=b использование разложения QR, явным образом не вычисляя фактор Q разложения QR.

rng('default');
m = 6;
n = 3;
p = 1;
A = randn(m,n);
b = randn(n,p);
x = fixed.qlessQRMatrixSolve(A,b)
x = 3×1

    0.2991
    0.0523
    0.4182

fixed.qlessQRMatrixSolve функция эквивалентна следующему коду, однако fixed.qlessQRMatrixSolve функция более эффективна и поддерживает типы данных с фиксированной точкой.

x = (A'*A)\b
x = 3×1

    0.2991
    0.0523
    0.4182

В этом примере показано, как задать, выходные данные вводят, чтобы решить систему уравнений с данными фиксированной точки.

Задайте данные, представляющие систему уравнений. Задайте матричный A как нулевое среднее значение, нормально распределенная случайная матрица со стандартным отклонением 1.

rng('default');
m = 6;
n = 3;
p = 1;
A0 = randn(m,n);
b0 = randn(n,p);

Задайте типы данных с фиксированной точкой для A и b избегать переполнения во время расчета QR.

T.A = fi([],1,22,16);
T.b = fi([],1,22,16);
A = cast(A0, 'like', T.A)
A = 
    0.5377   -0.4336    0.7254
    1.8339    0.3426   -0.0630
   -2.2589    3.5784    0.7147
    0.8622    2.7694   -0.2050
    0.3188   -1.3499   -0.1241
   -1.3077    3.0349    1.4897

          DataTypeMode: Fixed-point: binary point scaling
            Signedness: Signed
            WordLength: 22
        FractionLength: 16
b = cast(b0, 'like', T.b)
b = 
    1.4090
    1.4172
    0.6715

          DataTypeMode: Fixed-point: binary point scaling
            Signedness: Signed
            WordLength: 22
        FractionLength: 16

Задайте тип выходных данных, чтобы избежать переполнения в задней замене.

T.x = fi([],1,29,12);

Используйте fixed.qlessQRMatrixSolve функция, чтобы вычислить решение, x.

x = fixed.qlessQRMatrixSolve(A,b,T.x)
x = 
    0.2988
    0.0522
    0.4180

          DataTypeMode: Fixed-point: binary point scaling
            Signedness: Signed
            WordLength: 29
        FractionLength: 12

Сравните этот результат с результатом встроенных в MATLAB® операций в с двойной точностью, с плавающей точкой.

x0 = (A0'*A0)\b0
x0 = 3×1

    0.2991
    0.0523
    0.4182

Входные параметры

свернуть все

Матрица коэффициентов в линейной системе уравнений (A 'A) x = B.

Типы данных: single | double | fi
Поддержка комплексного числа: Да

Входной вектор или матрица, представляющая B в линейной системе уравнений (A 'A) x = B.

Типы данных: single | double | fi
Поддержка комплексного числа: Да

Выходные данные вводят в виде a numerictype возразите или числовая переменная. Если outputType задан как numerictype объект, выход, x, будет иметь заданный тип данных. Если outputType задан как числовая переменная, x будет иметь совпадающий тип данных как числовую переменную.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | fi | numerictype

Упущение фактора в виде неотрицательного скаляра между 0 и 1. Фактор упущения определяет, сколько веса прошлые данные даны. forgettingFactor значение умножается на выход разложения QR, R после каждой строки A обрабатывается.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | fi

Выходные аргументы

свернуть все

Решение, возвращенное как вектор или матрица. Если A m-by-n матрица и B m-by-p матрица, затем x n- p матрица.

Расширенные возможности

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