MATLAB数学建模 线性方程式与线性系统
2020-08-25 17:47
986 查看
线性方程
本次内容涉及线性代数,视频中大部分在讲解线性代数的知识,只稍微提及了几个matlab来实现的指令。
学了现代之后再来看一遍(逃~
求解线性方程
将线性方程组用矩阵 Ax=b
表示,则可通过求解矩阵来解方程:
高斯消去法rref()
R = rref(A)使用 Gauss-Jordan 消元法和部分主元消元法返回A的简化行阶梯形。
对增广矩阵
[A b]使用
rref()则可以求解
Ax=b对应的线性方程组
LU因子化
[L,U,P] = lu(A)将满矩阵或稀疏矩阵 A 分解为一个上三角矩阵 U 和一个经过置换的下三角矩阵 L,使得 A = L*U;返回一个置换矩阵 P,并满足 A = P’*L*U。
通过执行 LU 分解,然后使用因子来简化问题,对线性方程组求解。
一些矩阵分解函数
qr():正交三角分解 ldl():Hermitian 不定矩阵的分块 LDL 分解 ilu():不完全 LU 分解 chol():Cholesky 分解 gsvd():广义奇异值分解 svd():奇异值分解
高效mldivide()、\
以上两种方法在对于一般的线性方程组的求解其实并不友好,过于繁琐。实际上,更加高效的方式是使用A\b
(或者mldivide(A,b)
)可直接求得方程组的根 向量x
。
克莱默法则
求解矩阵方程 Ax=b
,x等于A的逆矩阵乘以b,即 x=A-1b
通过
inv(A)对矩阵A求逆,然后直接计算即可:
x = inv(A)*b
需要注意,矩阵A的逆矩阵可能不存在
线性系统
特征值和特征向量eig()
e = eig(A)返回一个列向量,其中包含方阵 A 的特征值;
[V,D] = eig(A)返回特征值的对角矩阵 D 和矩阵 V,其列是对应的右特征向量,使得 AV = VD。
模糊综合评价中利用
eig()求解最大特征根和权向量:
矩阵指数expm()
Y = expm(X)计算 X 的矩阵指数。如果 X 有一组完整的特征向量 V 和对应特征值 D,[V,D] = eig(X),则expm(X) = V*diag(exp(diag(D)))/V;
对于逐个元素的指数运算,使用
exp()
习题
syms R1 R2 R3 R4 R5 V1 V2; A=[R1 0 0 R4 0;0 R2 0 -R4 R5;0 0 -R3 0 R5;1 -1 0 -1 0;0 1 -1 0 -1]; b=[V1;0;V2;0;0]; x=A\b
以上内容为个人笔记,部分图片来源于郭老师课件或课程截图。
笔记汇总:MATLAB基础教程
课程视频:https://www.bilibili.com/video/BV1DA411Y7bN
相关文章推荐
- [数学建模]线性规划与matlab,lingo解法
- [数学建模]线性规划与matlab解法
- 【数学建模集训系列】公交查询系统的matlab实现-运行截图和结果
- matlab 灰色系统预测 GM(1,1) 数学建模
- 【数学建模集训系列】公交查询系统的matlab实现-加入地铁线T1
- 【数学建模集训系列】公交查询系统的matlab实现-加入地铁线T2
- 【数学建模集训系列】公交查询系统的matlab实现-公交与地铁换乘
- 【数学建模集训系列】公交查询系统的matlab实现-站点和站点之间直达矩阵
- 【数学建模集训系列】公交查询系统的matlab实现-问题重述
- 【数学建模集训系列】公交查询系统的matlab实现-只含公交的查询
- 【数学建模集训系列】公交查询系统的matlab实现-公交站点和线路对应矩阵
- 【数学建模集训系列】公交查询系统的matlab实现-价格计算
- 数学建模常用Matlab/Lingo/c代码总结系列——Matlab图形绘制函数汇总
- 数学建模常用的分析法及其MATLAB实现
- 数学建模入门-matlab实现偏微分方程数值解
- 经济金融领域简单数学建模和分析:MATLAB成本曲线方程和销售收入直线方程
- 数学建模专栏 | 第十二篇:MATLAB CUMCM真题求解实例三:机理建模型
- 数学建模与MATLAB-11
- 数学建模神器——Matlab下载安装教程
- Matlab 数学建模方法(七): 连续模型求解(目标跟踪问题实例)