数学建模matlab 课件 9w、10、矩阵.pptVIP

[V,D]=eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365 四、行列式的值 MATLAB提供的内部函数det用来计算矩阵的行列式的值。设矩阵A为一方阵(必须是方阵)，求矩阵A的行列式值的格式为：det(A)。注意：本函数同样能计算通过构造出的稀疏矩阵的行列式的值。关于如何构造稀疏矩阵，将在本章最后一节介绍。 【例6】利用随机函数产生一个三阶方阵A，然后计算方阵之行列式的值。 A=rand(3) A = 0.9501 0.4860 0.4565 0.2311 0.8913 0.0185 0.6068 0.7621 0.8214 det(A) ans = 0.4289 inv —— 矩阵求逆 det —— 行列式的值 eig —— 矩阵的特征值 diag —— 对角矩阵 ’ —— 矩阵转置 sqrt —— 矩阵开方 注： 矩阵的其它运算 矩阵的一些特殊操作 矩阵的变维 a=[1:12];b=reshape(a,3,4) c=zeros(3,4);c(:)=a(:) 矩阵的变向 rot90:旋转; fliplr:上翻; flipud:下翻 矩阵的抽取 diag:抽取主对角线;tril: 抽取主下三角; triu:抽取主上三角 习题 求A的转置矩阵、行列式、特征值特征向量、逆矩阵 2.给出一个5阶矩阵A，并在A上使用命令diag， triu， tril相关的所有格式 到此为止 谢谢? 矩阵及其运算 0、一般矩阵的输入 直接输入需遵循以下基本规则： 整个矩阵应以“ [ ]”为首尾，即整个输入矩阵必须包含在方括号中； 矩阵中，行与行之间必须用分号“ ；”或 Enter 键（ 按 Enter 键）符分隔； 每行中的元素用逗号“ ，”或空格分隔； A=[1,2,3,4;5,6,7,8;9,10,11,12;13,14,15,16] ?A =? 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 一、 特殊矩阵的实现 1.零矩阵:所有元素值为零的矩阵称为零矩阵。零矩阵可以用zeros函数实现。zeros是MATLAB内部函数，使用格式如下： zeros(m)：产生m? m阶零矩阵； zeros(m,n)：产生m? n阶零矩阵，当m=n时等同于zeros(m)； zeros(size(A))：产生与矩阵A同样大小的零矩阵。 一、 特殊矩阵的实现 常见的特殊矩阵有零矩阵、幺矩阵、单位矩阵、三角形矩阵等，这类特殊矩阵在线性代数中具有通用性；还有一类特殊矩阵在专门学科中有用，如有名的希尔伯特(Hilbert)矩阵、范德蒙(Vandermonde) 矩阵等。 2.幺矩阵:所有元素值为1的矩阵称为幺矩阵。幺矩阵可以用ones函数实现。它的调用格式与zeros函数一样。 【例1】 试用ones分别建立3?2阶幺矩阵、和与前例矩阵A同样大小的幺矩阵。 用ones(3,2) 建立一个3? 2阶幺阵： ones(3,2) % 一个3?2阶幺阵 ans =1 1 1 1 1 1 一、 特殊矩阵的实现 3.单位矩阵:主对角线的元素值为1、其余元素值为0的矩阵称为单位矩阵。它可以用MATLAB内部函数eye建立，使用格式与zeros相同。 4.数量矩阵:主对角线的元素值为一常数d、其余元素值为0的矩阵称为数量矩阵。显然，当d=1时，即为单位矩阵，故数量矩阵可以用eye(m)*d或eye(m,n)*d建立。 一、 特殊矩阵的实现 5.对角阵:对角线的元素值为常数、其余元素值为0的矩阵称为对角阵。我们可以通过MATLAB内部函数diag，利用一个向量构成对角阵；或从矩阵中提取某对角线构成一个向量。使用 格式为diag(V)和diag(V,k)两种。 6.用一个向量V构成一个对角阵 设V为具有m个元素的向量，diag(V)将产生一个m?m阶对角阵，其主对角线的元素值