学习线性代数的第一个学期,学了行列式,矩阵,向量,向量组........
这不是要考试了吗,突然在模拟考试的时候发现一个神仙软件 ——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