资源描述:
《上机教学一线性代数中的基本运算》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、《线性代数》—上机教学(一)理学院ScienceCollege西南石油大学SouthwestPetroleumUniversity一《线性代数》中的基本运算1上机目的:一、突出教学的重点,减少教学的难点;二、掌握常用计算方法和处理问题的方法.2上机内容:向量与矩阵运算上机软件:Matlab一、向量与矩阵运算的环境1、可以在Matlab指令窗口中运算指令窗工作空间历史纪录/当前路径32、可以在Matlab中建立M文件进行运算M文件的建立:File-New-M-file(或点击工具栏上的白页)4(1)直接输入矩阵用方
2、括号“[]”括起例:>>A=[123;456;789]矩阵同一行中的元素之间用空格或逗号分隔矩阵行与行之间用分号分开直接输入法中,分号可以用回车代替例:二、向量与矩阵运算1、矩阵的生成5(2)由向量生成(4)由函数生成(3)通过编写m文件生成例:>>x=[1,2,3];y=[2,3,4];>>A=[x,y],B=[x;y]例:>>C=magic(3)(魔方矩阵)6C=816357492常见矩阵生成函数zeros(m,n)生成一个m行n列的零矩阵,m=n时可简写为zeros(n)ones(m,n)生成一个m行n列的
3、元素全为1的矩阵,m=n时可写为ones(n)eye(m,n)生成一个主对角线全为1的m行n列矩阵,m=n时可简写为eye(n),即为n维单位矩阵diag(X)若X是矩阵,则diag(X)为X的主对角线向量若X是向量,diag(X)产生以X为主对角线的对角矩阵tril(A)提取一个矩阵的下三角部分triu(A)提取一个矩阵的上三角部分rand(m,n)产生0~1间均匀分布的随机矩阵m=n时简写为rand(n)randn(m,n)产生均值为0,方差为1的标准正态分布随机矩阵m=n时简写为randn(n)其它特殊矩阵
4、生成函数:magic、hilb、pascal72、向量的生成(1)直接输入:a=[1,2,3,4](2)用冒号运算符生成a=[1:4]==>a=[1,2,3,4]b=[0:pi/3:pi]==>b=[0,1.0472,2.0944,3.1416]c=[6:-2:0]==>c=[6,4,2,0]例:(3)从矩阵中抽取行或列8矩阵元素可以是任何数值表达式例:>>x=[-1.3,sqrt(3),(1+2+3)*4/5]3、矩阵元素赋值矩阵元素的单独赋值例:>>x(5)=abs(-3)注:Matlab自动将向量x的长度扩
5、展到5,并将未赋值部分置零.9大矩阵可以把小矩阵作为其元素例:>>A=[123;456;789]B=[A;101112]在原矩阵的下方加一行如何在原矩阵的右边添加一列?3、矩阵元素赋值10单个元素的引用例:>>A(2,3)4、矩阵元素的引用多个元素的引用:冒号的特殊用法利用小括弧和元素所在的位置(下标)x(i):向量x中的第i个元素A(i,j):矩阵A中的第i行,第j列元素a:b:c产生一个由等差序列组成的向量;a是首项,b是公差,c确定最后一项;若b=1,则b可以省略.例:>>x=1:2:5>>y=1:2:6例
6、:>>x=2:1:5>>y=2:511多个元素的引用(续)例:>>A=[123;456;789];x=[123456];>>x(1:3)>>A(3,1:3)>>A(1:2,2:3)可利用冒号提取矩阵的整行或整列.例:>>A(1,:)>>A(:,1:2)>>A(:,:)A(i:j,m:n)表示由矩阵A的第i到第j行和第m到第n列交叉线上的元素组成的子矩阵.4、矩阵元素的引用125、冒号的作用提取矩阵的部分元素:冒号运算符A(:)A的所有元素A(:,:)二维矩阵A的所有元素A(:,k)A的第k列,A(k,:)A的第k
7、行A(k:m)A的第k到第m个元素A(:,k:m)A的第k到第m列组成的子矩阵A(:)与A(:,:)的区别?如何得到由A的第一、三行和第一、二列组成的子矩阵?自己动手136、矩阵的转置矩阵的转置与共轭转置’共轭转置.’转置,矩阵元素不取共轭例:>>A=[12;2i3i]>>B=A'>>C=A.'点与单引号之间不能有空格!在英文状态下输入’或.’147、矩阵的大小查看矩阵的大小:sizesize(A)列出矩阵A的行数和列数size(A,1)返回矩阵A的行数size(A,2)返回矩阵A的列数例:>>A=[123;45
8、6]>>size(A)>>size(A,1)>>size(A,2)length(x)返回向量X的长度length(A)等价于max(size(A))158、矩阵基本运算矩阵的加减:对应分量进行运算要求参与加减运算的矩阵具有相同的维数例:>>clear>>A=[123;456];B=[321;654]>>C=A+B>>D=A-B矩阵的普通乘法要求参与运算的矩阵满足线性代数中