线性代数的数学实验.ppt

线性代数的数学实验1.1.MATLAB特点:1.数值计算和符号计算功能MATLAB的数值计算功能包括：矩阵运算、多项式和有理分式运算、数据统计分析、数值积分、优化处理等。符号计算将得到问题的解析解。2.MATLAB语言MATLAB除了命令行的交互式操作以外，还可以程序方式工作。使用MATLAB可以很容易地实现C或FORTRAN语言的几乎全部功能，包括Windows图形用户界面的设计。MATLAB特点:3.图形功能MATLAB提供了两个层次的图形命令：一种是对图形句柄进行的低级图形命令，另一种是建立在低级图形命令之上的高级图形命令。利用MATLAB的高级图形命令可以轻而易举地绘制二维、三维乃至四维图形，并可进行图形和坐标的标识、视角和光照设计、色彩精细控制等等。4.应用工具箱基本部分和各种可选的工具箱。基本部分中有数百个内部函数。其工具箱分为两大类：功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。学科性工具箱专业性比较强，如控制系统工具箱、信号处理工具箱、神经网络工具箱、最优化工具箱、金融工具箱等，用户可以直接利用这些工具箱进行相关领域的科学研究。常用基本命令:矩阵运算生成向量[1,2,3]生成行向量a:delta:b生成行向量linspace生成矩阵[1,2,3;4,5,6;7,8,9]取第i行第j列a(i,j)或a((j-1)*r+i)取第i行所有列a(i,:)取第i行到第j行所有列a(i:j,:)求矩阵行数列数[r,c]=size(M)生成零矩阵zeros(r,c)生成全1矩阵ones(r,c)生成单位矩阵eye(n)生成以某向量为主对角线的矩阵diag(vec)取矩阵主对角线元素生成列向量diag(M)常用基本命令:矩阵运算求矩阵的秩rank(M)矩阵转置M.’(转置)M’(共轭转置)将矩阵按列优先顺序重新排列reshape(M,r,c)求行列式det(M)求特征多项式poly(M)求矩阵的幂M^n求矩阵每个元素求幂M.^n求矩阵特征值eig(M)求矩阵特征向量(Eigenvectors[M])求矩阵的特征值与特征向量[vec,val]=eig(M)求矩阵的逆inv(M)求线性方程组Ax=b的解(LinearSolve[A,b])求线性方程组Ax=0的解空间的基(NullSpace[M])将矩阵化为行最简(RowReduce[M])1.2、利用MATLAB求行列式二、利用MATLAB进行矩阵运算1、矩阵的加、减2、数与矩阵相乘3、矩阵与矩阵相乘4、求矩阵的逆5、矩阵相除6、矩阵的秩三、利用MATLAB求方程组的解习题习题1已知求：习题2求解方程组4.1求特征值和特征向量如果A是n×n的矩阵，n个λ值满足Ax=λx，则λ为A的特征值，x为A的特征向量。计算特征值使用eig(A)，以列向量形式返回特征值。如果A是实对称矩阵，特征值为实数，如果A不对称，特征值常为复数。利用矩阵函数[V,D]=eig(A)可求得矩阵A的特征值和特征向量.例1求A的特征值各特征向量:程序如下:A=[-1,1,0;-4,3,0;1,0,2]A=-110-430102[V,D]=eig(A)V=00.40820.408200.81650.81651.0000-0.4082-0.4082D=200010001例2将矩阵对角化.程序如下:A=[-1,3,-1;-3,5,-1;-3,3,1]A=-13-1-35-1-331[V,D]=eig(A)V=-0.57740.4264-0.4942-0.57740.6396-0.6763-0.57740.6396-0.5462D=1.00000002.0000