matlab
优化模型算法
线性齐次方程组的通解
rank 求秩
% 定义矩阵 A
A = [1 2 3; 4 5 6; 7 8 9];
% 计算矩阵 A 的秩
matrix_rank = rank(A);
disp(['矩阵的秩为:', num2str(matrix_rank)]);
null 求通解
A=[1 2 2 1;2 1 -2 -2;1 -1 -4 -3];
format rat %指定有理式格式输出
B=null(A,'r') %求解空间的有理基
syms k1 k2
X=k1*B(:,1)+k2*B(:,2) %写出方程组的通解
A = [1 1 1 1; 2 3 4 -5; 3 4 5 -6];
format rat %指定有理式格式输出
B=null(A,'r') %求解空间的有理基
syms k1
X=k1*B(:,1) %写出方程组的通解
%结果展示
B =
1
-2
1
0
X =
k1
-2*k1
k1
0