MATLAB线性方程组求解

最新推荐文章于 2024-07-29 00:30:00 发布

小游园

最新推荐文章于 2024-07-29 00:30:00 发布

阅读量1.2k

点赞数 2

分类专栏：智能钻完井油田大数据 Matlab在石油工程中应用

本文链接：https://blog.csdn.net/s0302017/article/details/103549159

版权

本文详细介绍了如何使用MATLAB高效地解决线性方程组问题，包括直接方法和迭代方法，适用于科学计算和工程领域的初学者。

摘要由CSDN通过智能技术生成

%% 直接法求线性方程组
% A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];%系数
% b=[13,-9,6,0]';
% x=A\b

%% 矩阵的分解求线性方程组
%UL分解法
% A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];%系数
% b=[13,-9,6,0]';
% [L,U]=lu(A);
% x=U\(L\b)

%QR分解法
% A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];%系数
% b=[13,-9,6,0]';
% [Q,R]=qr(A);
% x=R\(Q\b)

%QR分解法第二种格式
% A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];%系数
% b=[13,-9,6,0]';
% [Q,R,E]=qr(A);
% x=E*(R\(Q\b))

%% cholesky分解法
% A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];%系数
% b=[13,-9,6,0]';
% R=chol(A)  %可以判段矩阵是否为正定
%% Jacobi迭代解法(雅克比)
% A=[10,-1,0;-1,10,-2;0,-2,10];
% b=[9,7,6]';
% [x,n]=jacobi(A,b,[0,0,0]',1.0e-6)
% function[y,n]=jacobi(A,b,x0,eps)
% if nargin == 3
%     eps=1.0e-6;
%  
% else if nargin<3
%         error
%         return
% end
%     D=diag(diag(A));
%     L=-tril(A,-1);
%     U=-triu(A,1);
%     B=D\(L&#