В этом примере показано, как решить простую систему линейных уравнений $\mathit{Ax}\mathit{=}$b, используя QR разложение.

В этом примере определите A как матрицу 5 на 3 с большим числом условий. Для решения системы линейных уравнений, включающих неупорядоченные (большое число условий) неквадратные матрицы, необходимо использовать QR-разложение.

rng default;
A = gallery('randsvd', [5,3], 1000000);
b = [1; 1; 1; 1; 1];
x = fixed.qrMatrixSolve(A,b)

x = 3×1
104 ×

   -2.3777
    7.0686
   -2.2703

Сравните результат fixed.qrMatrixSolve с результатом mldivide или \ функция.

x = A\b

x = 3×1
104 ×

   -2.3777
    7.0686
   -2.2703

В этом примере показан эффект параметра регуляризации при решении сверхопределённой системы. В этом примере количество y измеряется при нескольких различных значениях времени t для получения следующих замечаний.

t = [0 .3 .8 1.1 1.6 2.3]';
y = [.82 .72 .63 .60 .55 .50]';

Моделирование данных с помощью затухающей экспоненциальной функции

$y(t){}_{\mathrm{=}}{}_{\mathrm{c1}}{}^{+}$c2e-t.

Предыдущее уравнение говорит, что вектор y должен быть аппроксимирован линейной комбинацией двух других векторов. Один - постоянный вектор, содержащий все, а другой - вектор с компонентами exp(-t). Неизвестные коэффициенты, $${}_{\mathrm{c1}}$$ и $${}_{\mathrm{c2}}$$, могут быть вычислены путем аппроксимации методом наименьших квадратов, что минимизирует сумму квадратов отклонений данных от модели. Существует шесть уравнений и два неизвестных, представленных матрицей 6 на 2.

E = [ones(size(t)) exp(-t)]

E = 6×2

    1.0000    1.0000
    1.0000    0.7408
    1.0000    0.4493
    1.0000    0.3329
    1.0000    0.2019
    1.0000    0.1003

Используйте fixed.qrMatrixSolve функция для получения решения методом наименьших квадратов.

c = fixed.qrMatrixSolve(E, y)

c = 2×1

    0.4760
    0.3413

Другими словами, наименьшие квадраты, соответствующие данным,

Следующие инструкции оценивают модель с регулярно разнесенными приращениями в t, а затем постройте график результата вместе с исходными данными:

T = (0:0.1:2.5)';
Y = [ones(size(T)) exp(-T)]*c;
plot(T,Y,'-',t,y,'o')

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

figure;
lambda = [0:0.1:0.5];
plot(t,y,'o', 'DisplayName', 'Original Data');
for i = 1:length(lambda)
 c = fixed.qrMatrixSolve(E, y, numerictype('double'), lambda(i));
 Y = [ones(size(T)) exp(-T)]*c;
 hold on
 plot(T,Y,'-', 'DisplayName', ['lambda =', num2str(lambda(i))])
end
legend('Original Data', 'lambda = 0', 'lambda = 0.1', 'lambda = 0.2', 'lambda = 0.3', 'lambda = 0.4', 'lambda = 0.5')