Octave使用矩阵法解线性方程如何实现X=R格式的结果输出
Octave线性方程组求解结果配对输出实现方案
直接修改原有代码如下即可实现变量与结果一一对应展示的效果:
适配任意元方程组的完整代码
A = input('Enter the coeffecient matrix A: '); B = input('Enter the constant matrix B (column matrix form): '); disp('输入的系数矩阵A:'); disp(A); disp('输入的常数矩阵B:'); disp(B); % 若固定求解3元方程组,可将下一行替换为 X = {'x'; 'y'; 'z'}; var_count = size(A, 2); X = strcat('x', cellstr(num2str((1:var_count)'))); R = inv(A) * B; disp('======= 求解结果 ======='); disp('Matrix(X) = Matrix(R) 对应关系:'); for i = 1:var_count fprintf(' %s = %.6f\n', X{i}, R(i)); end % 可选:矩阵形式并排展示 disp('======= 矩阵格式并排展示 ======='); result_mat = [X, num2cell(R)]; disp(result_mat);
修改说明
- 变量名支持动态生成,自动匹配系数矩阵的列数(即未知数个数),不需要手动调整变量列表
- 用循环搭配
fprintf逐行输出变量名和对应计算结果,可通过修改%.6f调整数值保留的小数位数 - 额外提供的元胞矩阵并排展示效果,完全符合你需要的矩阵对应展示逻辑
内容的提问来源于stack exchange,提问作者Sreek
相关产品推荐
相关产品推荐

