线性代数怎么过?先教你用好matlab!!!

学习线性代数的第一个学期,学了行列式,矩阵,向量,向量组........

这不是要考试了吗,突然在模拟考试的时候发现一个神仙软件 ——MATLAB,正好这学期在学,好上手,来整理一下,考试的时候,可以用上,哈哈哈,我是天才!!

一,求行列式的值

A = [1 -2 4; -5 2 0; 1 0 3]

d = det(A)

二,行列式初等变换

A = magic(3)

RA = rref(A)

三,解齐次方程组

当齐次线性方程AX=0,rank(A)=r<n时,该方程有无穷多个解。

A=[1 1 1 1 -3 -1 1;1 0 0 0 1 1 0;-2 0 0 -1 0 -1 -2];
r=rank(A);
y=null(A, r )

四,解非齐次方程组

非齐次线性方程组的求解

五,计算矩阵的逆(2个方法)

b=inv(a)
b=a^-1

六,计算矩阵的秩

 a=[1 4 5;3 6 8]

 rank(a)

七,矩阵化成最简型(找出向量组的极大无关组)

 a=[1 4 5;3 6 8]

 rref(a)

八,计算齐次方程组的基础解系

求齐次线性方程组的基础解系的命令:

B=null(A)   B的列向量是AX=0的规范正交的基础解系

B=null(A,'r')  B的列向量是AX=0的有理数形式的基础解系

九,计算矩阵的特征值和特征向量

a=[1 2 3;2 4 5;7 8 9]

[x,y]=eig(a)

y的对角线的数为特征值

十,将二次型标准化

假设题目是这样:f(x1,x2,x3)=x1^2+2x2^2+3x3^2+2x1x3+3x2x3

在Matlab中,我们运用函数eig求出二次型的矩阵A的特征值D和特征向量矩阵P,所求的矩阵D即为系数矩阵A的标准形,

矩阵P即为二次型的变换矩阵。

syms y1 y2 y3
A=[1 0 1;   0 2 3/2;   1 3/2 3];
[P,D]= eig(A)
y=[y1;y2;y3];
x=P*y%所求的正交变换
f=[y1 y2 y3]*D*y
x=vpa(x,5)
f=vpa(f,5)

结果:

x = 0.72551*y1 + 0.64255*y2 + 0.24651*y3

      0.45326*y1 - 0.71565*y2 + 0.5314*y3

      0.2738*y2 - 0.51787*y1 + 0.81046*y3

即 f =0.28619*y1^2 + 1.4261*y2^2 + 4.2877*y3^2

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值