• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

Matlab - 线性方程组求解

原作者: [db:作者] 来自: [db:来源] 收藏 邀请


常用函数:
det 计算矩阵的行列式的值
inv 求矩阵的逆阵
rank 求矩阵的秩
[V D]=eig(A) 求矩阵A的特征值和特征向量
poly 求矩阵的特征多项式
rref 用初等变换将矩阵化成行阶梯形
null(A,’r’) 给出齐次线性方程组Ax=0 的基础解系
fliplr 矩阵左右翻转
flipud 矩阵上下翻转
trace 求矩阵的迹
diag 取得矩阵对角线元素


例子:
1、矩阵函数的应用
A=[3 -4 0; -1 5 2; 4 1 -6]
det (A) %求矩阵的行列式的值
rank (A) %求矩阵的秩
inv (A) %求逆矩阵


2、求解线性方程组Ax=B
A=[3 -4 0; -1 5 2; 4 1 -6];
B=[5;5;16];
解法1 利用矩阵除法: X=A\B
解法2 利用求逆矩阵函数 inv:X1=inv(A)*B
比较:解法1比解法2更简便,
解法1 的算法优于解法2 ,
解法1可用于一般矩阵,而解法2只能用于非奇异的方阵
因此,只需运用解法1 .


3 求线性方程组的通解


3.1 利用除法\和null函数
A=[1 1 -1 -1;2 -5 3 2; 7 -7 3 1];
B=[5;-4;7];
format rat
x1=A\B %求得非齐次方程组Ax=B的一个特解x1
Y=null(A,\'r\') %求得齐次方程组Ax=0 的基础解系Y
输出结果:
x1 =

3
2
0
0

Y =

2/7 3/7
5/7 4/7
1 0
0 1
则方程组Ax=B的通解为: x=x1+k1*Y(:,1)+k2*Y(:,2)

3.2 利用rref函数
format rat
A=[1 1 -1 -1;2 -5 3 2;7 -7 3 1];
B=[5; -4; 7];
%用初等行变换将增广矩阵 [A B] 化成最简行阶梯形T
T=rref([A B])
输出结果:
T =

1 0 -2/7 -3/7 3
0 1 -5/7 -4/7 2
0 0 0 0 0

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
在Delphi中使用C++对象(两种方法,但都要改造C++提供的DLL)发布时间:2022-07-18
下一篇:
DelphiXE2之FireMonkey入门(8)-TImage发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap