【发布时间】:2023-03-23 18:59:01
【问题描述】:
所以我有 AX = B,其中 A 是一个 n x n 矩阵,B = [b1(t); b2(t); b3(t); ... ; bn(t)]。所以 X 将是一个 n x 1 矩阵,每个条目都是 t 的函数。
如何使用高斯消除、LU 分解或任何其他方法找到 X?
【问题讨论】:
-
我的矩阵代数很陈旧,但
X不只是A\B吗? (即逆(A)*B)
所以我有 AX = B,其中 A 是一个 n x n 矩阵,B = [b1(t); b2(t); b3(t); ... ; bn(t)]。所以 X 将是一个 n x 1 矩阵,每个条目都是 t 的函数。
如何使用高斯消除、LU 分解或任何其他方法找到 X?
【问题讨论】:
X 不只是A\B 吗? (即逆(A)*B)
LU 分解是 Matlab 提供的一个内置函数。见here.
【讨论】:
LU 分解:
[L, U, P] = lu(A);
X = (U \ (L \ (P * B)));
你也可以这样做
X = A\B;
利用 A 的任何潜在特殊结构,可以加快计算速度。
【讨论】: