第二章 矩阵和数组
matlab第二章矩阵运算基础
南京信息工程大学
4
例2.1 创建矩阵
>>x=[1 2 3;4 5 6;7 8 9] >>x=[1 2 3 456 7 8 9] >>x=[a b c;e f g;u v w] >>x=[1 2 3;4 5 6]; y=[2 3 4;5 6 7] >>Q=x*y >>a=2;b=3 >>x=a*b
2010-12-29
2010-12-29 南京信息工程大学 6
2.1 矩阵的创建
2、 赋值语句 MATLAB赋值语句有两种格式:
变量=表达式(或数) 表达式
2010-12-29
南京信息工程大学
7
【例2.2】 x=[1,2,3;4,5,6;7,8,9] 与[1,2,3;4,5,6;7,8,9]。
5 + cos 47
【例2.3】计算
2010-12-29
南京信息工程大学
25
§2.2 矩阵和数组的算术运算 六、点运算
C=A.*B C=A.\B
C=A./B C=A.^B
2010-12-29
南京信息工程大学
26
§2.2 矩阵和数组的算术运算 七、幂运算
C=A^B C=A.^B
2010-12-29
南京信息工程大学
27
例2.12 例2.13 例2.14 例2.15
find(x)
检查x是 否全为1
南京信息工程大学 42
2010-12-29
例2.20 建立矩阵A,然后找出大于4的元素位置 (1)建立A >>A=[4 -6 5 -54 0 6 56 0 67 -45 0] (2)找出大于4的元素位置 >>find(A>4)
矩阵论第二章
(2)
则 0 是 经单位化得到的单位向量。 定理1: [cauchy—schwarz不等式]对于内积 空间中任意向量 , ,有 ( , )
(3)
并且, 等号成立的 , 线性相关。
9°(三角不等式)对 向量 , ,有
定义4:设 V 是数域 F上的线性空间, 如果在V 上还定义了一种叫内积的运算:对于V 中任意 向量 , 都有 F 中唯一的数 x 与之对应, 记为
, x, 并且这种内积运算还具有如下性质:
对于任意的 , , V
1) , ,
及任意的 k F
有:
2) k , k , 4) 当 0时, , 0
3) , , ,
此时称 V 为一个内积空间。
n C 对于复数域上的线性空间 , 若规定向量 例1:
a1 , a2 ,, an
1 1 , 2
( 2 , 1 ) 1 2 , [设 2 k1 2 , ( 1 , 1 )
( 2 , 1 ) k ( 1 , 1 )
因 ( 2 , 1 ) 0
],
3
( 3 , 1 ) ( , ) 1 3 2 2 3,…, ( 1 , 1 ) ( 2 , 2 )
定理3: 欧氏空间在一组基下的度量矩阵都 是正定矩阵。
, 证明:设 V 是 n 维欧氏空间,
1 2
,, n
是 V 的一
A 是该基下的度量矩阵。 的一组基, 为证明实对
称矩阵 A 正定, 只须证明实二次型 x
1 1 2 2 n n
T
Ax 正定,
MATLAB基础教程 第2章 数组、矩阵及其运算
写出MATLAB表达式。 解:根据MATLAB的书写规则,以上MATLAB表达式为: (1)y=1/(a*log(1-x-1)+C1) (2)f=2*log(t)*exp(t)*sqrt(pi) (3)z=sin(abs(x)+abs(y))/sqrt(cos(abs(x+y))) (4)F=z/(z-exp(T*log(8)))
命令:X(3:-1:1)
命令:X(find(X>0.5)) 命令:X([1 2 3 4 4 3 2 1])
第二章 数组、矩阵及其运算
2.1 数组(矩阵)的创建和寻访
2. 二维数组的创建和寻访
例2-3 综合练习。将教材P.31~P.44的实例按顺序在MATLAB的 command窗口中练习一遍,观察并体会其输出结果。 (注意变量的大小写要和教材上的严格一致。)
A./B
B.\A
A的元素被B的对应元素相除
(与上相同)
第二章 数组、矩阵及其运算
2.3 数组、矩阵的其他运算
1. 乘方开方运算
数组的乘方运算与power函数 格式:c=a.^k或c=power(a,k) 例如: >> g=[1 2 3;4 5 6] >>g.^2 矩阵的乘方运算与mpower函数 格式:C=A^P或C=mpower(A,P) 注意:A必须为方阵
第二章 数组、矩阵及其运算
2.2 数组、矩阵的运算
3. 矩阵的加法、减法
运算规则是:若A和B矩阵的维数相同,则可以执行矩阵的加减运算, A和B矩阵的相应元素相加减。如果维数不相同,则MATLAB将给出
出错信息。
第二章 数组、矩阵及其运算
2.2 数组、矩阵的运算
3. 矩阵的乘法
MATLAB教学 最新第二章 矩阵与数组2-4
把D的逆阵右乘以B,记作/D,称之为右除.
2.5.3 基本数组运算 1,数组转置 数组转置的操作符是在矩阵转置操作符前加符号".".(实数情 况下等价) 例:数组转置操作
2,数组幂 数组幂运算符 (单个符号自身运算)就是在矩阵运算符前加上符 号".".
3.数组乘法
2.5.4 基本数学函数 在MATLAB中部分函数可以用来进行基本的 数学运算,有三角函数,指数运算函数,复数 运算函数等. 注意:这些函数的参数可以是矩阵,向量或者 多维数组,函数在处理参数时,都是按照数组 运算运算的规则来进行的. 函数数目较多,不一一列出,后面用到时再 作说明. 2.5.5 矩阵(数组)操作函数
例2-5 使用logspace函数创建向量.
上面创建的都是行向量,即创建的都 是一行n列的二维数组.如果需要创建 列向量,即n行一列的数组,则需要使 用分号作为元素与元素之间的间隔或 者直接使用转置运算符" ' ".
2.3 创建矩阵 在编程语言中,矩阵和二维数组一般指的是同一 个概念,在M语言中,矩阵的元素可以为任意的 MATLAB数据类型的数值或者对象.创建矩阵的方 法也有多种,不仅可以直接输入元素,还可以使用 MATLAB MATLAB的数组编辑器编辑矩阵的元素. 2.3.1直接输入法 直接输入矩阵元素创建矩阵的方法适合创建元素较 少的矩阵. 例2-7 用直接输入矩阵元素的方法创建矩阵.
length获取向量长度若输入参数为矩阵或多维数组则返回各个维尺寸的最大值ndims获取矩阵或多维数组的维数numel获取矩阵或数组的元素个数disp显示矩阵或者字符串的内容cat合并不同的矩阵或者数组reshape保持矩阵元素的个数不变修改矩阵的行数和列数repmat复制矩阵元素并扩展矩阵fliplr交换矩阵左右对称位置上的元素flipud交换矩阵上下对称位置上的元素flipdim按照指定的方向翻转交换矩阵元素find获取矩阵或数组中非零元素的索引55例
第二章 MATLAB数值计算
数值
复数 z=a+b*i或z=a+b*j z=a+bi或z=a+bj(当b为标量时) z=r*exp(i*theta) 得出一个复数的实部、虚部、幅值和相角。 a=real(z) %计算实部 b=imag(z) %计算虚部 r=abs(z) %计算幅值 theta=angle(z) %计算相角
变量
矩阵和数组运算
diag(X)产生X矩阵的对角阵 [l,u]=lu(X)方阵分解为一个准下三角方阵和一个上三 角方阵的乘积。l为准下三角阵,必须交换两行才能成 为真的下三角阵。 [q,r]=qr(X) m×n阶矩阵X分解为一个正交方阵q和一个 与X同阶的上三角矩阵r的乘积。方阵q的边长为矩阵X的 n和m中较小者,且其行列式的值为1。 [u,s,v]=svd(X) m×n阶矩阵X分解为三个矩阵的乘积, 其中u,v为n×n阶和m×m阶正交方阵,s为m×n阶的对角 阵,对角线上的元素就是矩阵X的奇异值,其长度为n和 m中的较小者。
第二章
MATLAB数值计算
2.1 变量和数值
数据类型
数据类型包括:数值型、字符串型、元胞型、结构 型等 数值型=双精度型、单精度型和整数类 整数类=无符号类(uint8、uint16、uint32、uint64)和 符号类 (int8、int16、int32、int64)
数值
数据的表达方式 可以用带小数点的形式直接表示 用科学计数法 数值的表示范围是10-309~10309。 以下都是合法的数据表示: -2、5.67、2.56e-56(表示2.56×10^-56)、 4.68e204(表示4.68×10^204)
E1=tril(A) E2=tril(A,1) D=triu(A) E=triu(A,-1)
第二章 MATLAB语言的使用与程序设计
命令历史窗口:显示已执行过的命令。在窗口的某一命令上单击鼠标 右键,会弹出菜单,对所选命令进行操作。
当前路径窗口:提供了当前路径文件的操作
演示
MATLAB的搜索路径
搜索路径是一系列文件路径的组合。当程序和命令执行 时, MATLAB 在搜索路径中查找程序或命令运行所需的函数文 件。 MATLAB 在执行搜索时按照规定的顺序。如:在命令窗口 中输入example,MATLAB将按下面的步骤来处理: 1.检查example是不是一个变量,如果是,则返回变量的值;
本章重点:
MATLAB工作环境掌握 主要文件类型及常用命令
矩阵、变量、表达式、常用函数
MATLAB语言的基本语句结构及程序调试方法
一、MATLAB系统简介
MATLAB的主要组成部分
1.MATLAB语言体系:MATLAB 语言是一种以矩阵运算为基础的高级 语言,具有条件控制、函数调用、数据结构、输入输出及面向对象等 程序语言特征,可以进行程序设计。
6 )对矩阵的特殊操作: rot90(a) 将 a 矩阵旋转 90 度、 fliplr(a) 将 a 矩阵的列反序、 flipud(a) 将 a 矩阵的行反序、diag(a) 将向量 a 构 成对角阵( 元素放在主对角线上 )---a 为向量、triu(a) 提取矩阵的上 三角部分、reshape改变矩阵的阶数,按列的顺序重排。
逻辑运算符: 在MATLAB中,逻辑运算符有3种。 & 逻辑与。当运算双方对应元素都为非零时; 结果为1,否则,结果为0。
| 逻辑或。当运算双方对应元素有一个为非零 时;结果为1,否则,结果为0。
~ 结果为0。 逻辑非。当元素的值为 0 时,结果为 1 ,否则,
例: a=[1 0 3;0 –1 6] , b=[-1 0 0;0 5 0.3] ,计算两矩
线性代数知识点总结(高中)
线性代数知识点总结第一章 行列式二三阶行列式N 阶行列式:行列式中所有不同行、不同列的n 个元素的乘积的和 n nn nj j j j j j j j j nij a a a a ...)1(21212121)..(∑-=τ(奇偶)排列、逆序数、对换行列式的性质:①行列式行列互换,其值不变。
(转置行列式T D D =) ②行列式中某两行(列)互换,行列式变号。
推论:若行列式中某两行(列)对应元素相等,则行列式等于零。
③常数k 乘以行列式的某一行(列),等于k 乘以此行列式。
推论:若行列式中两行(列)成比例,则行列式值为零; 推论:行列式中某一行(列)元素全为零,行列式为零。
④行列式具有分行(列)可加性⑤将行列式某一行(列)的k 倍加到另一行(列)上,值不变 行列式依行(列)展开:余子式ij M 、代数余子式ij j i ij M A +-=)1(定理:行列式中某一行的元素与另一行元素对应余子式乘积之和为零。
克莱姆法则:非齐次线性方程组 :当系数行列式0≠D 时,有唯一解:)21(n j DD x j j ⋯⋯==、齐次线性方程组 :当系数行列式01≠=D 时,则只有零解逆否:若方程组存在非零解,则D 等于零 特殊行列式:①转置行列式:332313322212312111333231232221131211a a a a a a a a a a a a a a a a a a → ②对称行列式:ji ij a a =③反对称行列式:ji ij a a -= 奇数阶的反对称行列式值为零④三线性行列式:333122211312110a a a a a a a 方法:用221a k 把21a 化为零,。
化为三角形行列式⑤上(下)三角形行列式: 行列式运算常用方法(主要)行列式定义法(二三阶或零元素多的) 化零法(比例)化三角形行列式法、降阶法、升阶法、归纳法、第二章 矩阵矩阵的概念:n m A *(零矩阵、负矩阵、行矩阵、列矩阵、n 阶方阵、相等矩阵)矩阵的运算:加法(同型矩阵)---------交换、结合律 数乘n m ij ka kA *)(=---------分配、结合律乘法nm lkj ik n l kj l m ik b a b a B A *1**)()(*)(*∑==注意什么时候有意义一般AB=BA ,不满足消去律;由AB=0,不能得A=0或B=0转置A A T T =)( T T T B A B A +=+)( T T kA kA =)( T T T A B AB =)((反序定理) 方幂:2121k k k k A A A += 2121)(k k k k A A +=几种特殊的矩阵:对角矩阵:若AB 都是N 阶对角阵,k 是数,则kA 、A+B 、 AB 都是n 阶对角阵数量矩阵:相当于一个数(若……)单位矩阵、上(下)三角形矩阵(若……) 对称矩阵 反对称矩阵阶梯型矩阵:每一非零行左数第一个非零元素所在列的下方 都是0分块矩阵:加法,数乘,乘法:类似,转置:每块转置并且每个子块也要转置注:把分出来的小块矩阵看成是元素逆矩阵:设A 是N 阶方阵,若存在N 阶矩阵B 的AB=BA=I 则称A 是可逆的, B A =-1(非奇异矩阵、奇异矩阵|A|=0、伴随矩阵) 初等变换1、交换两行(列)2.、非零k 乘某一行(列)3、将某行(列)的K 倍加到另一行(列)初等变换不改变矩阵的可逆性 初等矩阵都可逆初等矩阵:单位矩阵经过一次初等变换得到的(对换阵 倍乘阵 倍加阵)等价标准形矩阵⎪⎪⎭⎫⎝⎛=O O O I D r r矩阵的秩r(A):满秩矩阵 降秩矩阵 若A 可逆,则满秩若A 是非奇异矩阵,则r (AB )=r (B ) 初等变换不改变矩阵的秩求法:1定义2转化为标准式或阶梯形矩阵与行列式的联系与区别:都是数表;行列式行数列数一样,矩阵不一样;行列式最终是一个数,只要值相等,就相等,矩阵是一个数表,对应元素相等才相等;矩阵n ij n ij a k ka )()(=,行列式n ij nn ij a k ka =逆矩阵注:①AB=BA=I 则A 与B 一定是方阵 ②BA=AB=I 则A 与B 一定互逆; ③不是所有的方阵都存在逆矩阵;④若A 可逆,则其逆矩阵是唯一的。
数组和矩阵(数据结构)
数组和矩阵1、(2分)【单选题】某串的长度小于一个常数,则采用( )存储方式最节省空间A、链式B、顺序C、堆结构D、无法确定参考答案:B解析:串的顺序和链式存储结构2、(2分)【单选题】与线性表相比,串的插入和删除操作的特点是( )。
A、通常以串整体作为操作对象B、需要更多的辅助空间C、算法的时间复杂度较高D、涉及移动的元素更多参考答案:A解析:串的基本运算3、(2分)【单选题】在稀疏矩阵的三元组表示法中,每个三元组表示( )。
A、矩阵中非零元素的值B、矩阵中数据元素的行号和列号C、矩阵中数据元素的行号、列号和值D、矩阵中非零数据元素的行号、列号和值参考答案:D解析:二维数组的存储结构及求址方法4、(2分)【单选题】已知二维数组A8X10,按行存储时,元素a12的地址为1000,每个元素占2个字节,则元素a00的地址为( )A、972B、974C、976D、978参考答案:C解析:二维数组的存储结构及求址方法5、(2分)【单选题】数组通常具有的两种基本操作是( )A、建立和删除B、索引和修改C、查找和修改D、查找和索引参考答案:C解析:二维数组的存储结构及求址方法6、(2分)【单选题】在长度为n的字符串S的第i个位置插入另外一个字符串,i的合法值应该是( )。
A、i>0B、i≤nC、1≤i≤nD、1≤i≤n+1参考答案:D解析:串的基本运算7、(2分)【单选题】两个字符串相等的条件是( )。
A、两串的长度相等B、两串包含的字符相同C、两串的长度相等,并且两串包含的字符相同D、两串的长度相等,并且对应位置上的字符相同参考答案:D解析:串的基本运算8、(2分)【单选题】设有串s=“software”,则其子串的数目是( )。
A、36B、37C、8D、9参考答案:B解析:串的基本运算9、(2分)【单选题】广义表A=((x,(a,b)),((x,(a,b)),y),y),则运算head(head(tail(A)))为( )A、xB、(a,b)C、(x,(a,b))D、A参考答案:C解析:广义表的概念10、(2分)【单选题】串的模式匹配是指( )。
4 矩阵和数组(2):运算与操作
一、矩阵变形
让reshape函数内部自动计算所需的行数
将向量转换成矩阵 (按列排放)
Reshape要求变形前后的元素总数目保持不变
>> A=[1:3;4:6;7:9] A= 1 2 3 4 5 6 7 8 9
>> rot90(A) ans = 3 6 9 2 5 8 1 4 7
逆时针旋转90度
2
3
4
5
8
9
>> B = sort(A,'descend') %降序 B= 9 8 5 4 3 2 0 -7 -10 >> B = sort(A, 'ascend' ) %升序 B= -10 -7 0 2 3 4 5 8
9
>> A = [9 0 -7 5 3 8 -10 4 2] A= 9 0 -7 5 3 8 -10 1 2 3 4 5 6 7
Description Circularly shift matrix contents. Sort array elements in ascending or descending order. Sort rows in ascending order. Determine if matrix elements are in sorted order.
基于多列排序: 前面一列若有相等的数值, 后面一列继续按照数值大小排序
>> B = sortrows(A,[1 7]) B= 76 79 91 0 19 76 44 40 35 60 76 61 93 81 27 95 82 17 5 19 95 1 73 89 20 95 45 92 41 13 >> C = sortrows(A,[1 -7]) C= 76 61 93 81 27 76 44 40 35 60 76 79 91 0 19 95 45 92 41 13 95 1 73 89 20 95 82 17 5 19
数组和矩阵的运算需要掌握运算法则
本节要求
1.掌握数组和矩阵的数值计算,尤其是数 组的“点运算” 2.掌握数组的关系和逻辑运算 3.掌握指令find的使用 4.了解notebook文档的编辑
本课件由飞华健康网/pifu/cc/编辑
>>A=[1 2 3; 4 5 6],B=[4,5,6;1,2,3]; >>C=zeros(2); %生成2阶全0方阵 >>c1=A+B %加法运算 >>c2=A-C %减法相乘 >>c3=A-2 %与标量之间的加减运算
结果如下:
c1 = 5 7
7 7
9 7
c3 = -1 2
0 3
1 4
??? Error using ==> minus Matrix dimensions must agree.
A =
>> A(~B)=0
13 8 12 1
A= 0 2 5 11 0 7 0 0 1 0 0 0 3 13 0 0 0 0 0 0
16 5 9 4
2 11 7 14
3 10 6 15
>> B=isprime(A)
B = 0 1 0 0 1 1 1 0 1 0 0 0
注:isprime是用来检测数值是否为质数。
log10 log2 exp pow2
例1:分析语句a=2+2==4的执行结果。
分析:单个等号表示赋值,后面的双等号表示关系 运算,所以a的值为1.
例2:分析语句a=‘fate’;b=‘cake’;result=a==b 的执行结果。
分析:应用关系运算应该逐个比较字符是否相等。 执行结果如下: result = 0 1 0 1
结果如下:
线性代数第二章,矩阵及其运算
a1n b1
a2n
b2
L L
amn bm
§2 矩阵的运算
一、加法
设 A (ai j )mn , B (bi j )mn 都是m n 矩阵,则加法定义为
a11 b11
A
B
a21
b21
L
a12 b12 L a22 b22 L
LL
am1 bm1 am2 bm2 L
显然,
AB B A
a22
L
L L L
am1 am2 L
a1n
a11 a21 L
a2n
,记
AT
a12
a22
L
L
L L L
amn
a1n an2 L
则称
AT
A
是
的转置矩阵。
am1
am 2
L
amn
显然,
① ( AT )T A ,② ( A B)T AT BT ,③( A)T AT ,④( AB)T BT AT
2. 即使 Amn , Bnm ,则Amn Bnm 是m 阶方阵,而Bnm Amn 是n 阶方阵;
3. 如 果 A , B
都 是n
阶
方
阵
,
例
如
2
A
1
4
2
,
B
2
3
4
6
,则
16
AB
8
32 16
,而BA
0 0
0
0
;
AB BA
综上所述,一般
(即矩阵乘法不满足交换率)。
但是下列性质显然成立:
三、乘法
乘法运算比较复杂,首先看一个例子
设变量t1, t2 到变量 x1, x2 , x3 的线性变换为
《计量经济学》第二章知识
第二章 数学基础 (Mathematics)第一节 矩阵(Matrix)及其二次型(Quadratic Forms)第二节 分布函数(Distribution Function),数学期望(Expectation)及方差(Variance) 第三节 数理统计(Mathematical Statistics ) 第一节 矩阵及其二次型(Matrix and its Quadratic Forms)2.1 矩阵的基本概念与运算 一个m ×n 矩阵可表示为:v a a a a a aa a a a A mn m m n n ij ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡== 212222111211][矩阵的加法较为简单,若C=A +B ,c ij =a ij +b ij但矩阵的乘法的定义比较特殊,若A 是一个m ×n 1的矩阵,B 是一个n 1×n 的矩阵,则C =AB 是一个m ×n 的矩阵,而且∑==nk kj ikij b ac 1,一般来讲,AB ≠BA ,但如下运算是成立的:● 结合律(Associative Law ) (AB )C =A (BC ) ● 分配律(Distributive Law ) A (B +C )=AB +AC 问题:(A+B)2=A 2+2AB+B 2是否成立?向量(Vector )是一个有序的数组,既可以按行,也可以按列排列。
行向量(row ve ctor)是只有一行的向量,列向量(column vector)只有一列的向量。
如果α是一个标量,则αA =[αa ij ]。
矩阵A 的转置矩阵(transpose matrix)记为A ',是通过把A 的行向量变成相应的列向量而得到。
显然(A ')′=A ,而且(A +B )′=A '+B ',● 乘积的转置(Transpose of a production ) A B AB ''=')(,A B C ABC '''=')(。
线性代数第二章 矩阵
(1)
其中 x1, x2 , , xn 是 n 个未知数,m 是方程的个数,
ai(j i 1, 2, , m,j 1, 2, , n)称为线性方程组的系
数,b1, b2, , bm 称为线性方程组的常数项.
由 n 个数 c1, c2 , , cn组成的有序数组 (c1, c2 , , cn ) 称为方程组(1)的解, 是指当 x1, x2 , , xn 分别用 c1, c2 , , cn 替换后,(1)的每个等式都变成了恒
例1 解线性方程组
x1 2x2 2x1 3x2
x3
x4 x4
2, 3,
x1 x2 x3 2x4 3.
我们就可以只考虑方程组的系数和常数项组成的 一个矩形数阵(后面我们称这种矩形数阵为矩阵), 对于方程组(1),其对应的矩形数阵为
a11 a12 a21 a22
等式. 方程组(1)的解的全体组成一个集合,这个集合
称为方程组(1)的解集合. 求解方程组实质上就是找到方程组的所有解,即求
出它的解集合. 把具有相同解集合的两个方程组称为同解的方程组.
定义1 对线性方程组(1)进行如下三种变形,称 为线性方程组的初等变换:
1)用一个非零数 k 乘以某一个方程; 2)用任意数 k 乘以一个方程加到另外一个方程上; 3)交换两个方程的位置.
1.矩阵的加法
第二章 矩阵
第一节 矩阵的基本概念
一、矩阵的引入
所谓具有 m 个方程 n 个未知数的线性方程组的 一般形式是指
a11x1 a12 x2
a21x1
a22 x2
am1x1 am2 x2
第2章数组、矩阵及其运算
,求A B。
2.2.4 集合运算
1.两个集合的交集
命令格式:
c = intersect(a,b) %返回向量a、b的 公共部分,即c= a∩b。 c = intersect(A,B,'rows') %A、B为相同列数 的矩阵,返回元素相同的行。 [c,ia,ib] = intersect(a,b) %c为a、b的公共 元素,ia表示公共元素在a中的位置,ib表示公 %共元素在b中位置。
2.子数组的赋值
>>x(3) = 0 新赋值为0。 x = 0.9501 0.2311 >>x([1 4])=[1 1] 个元素都赋值为1。 x = 1.0000 0.2311 %把上例中的第三个元素重
0 0.4860 0.8913 %把当前x数组的第一、四
0
1.0000
0.8913
2.1.3 二维数组(矩阵)的创建
a 1n
定义为
;数组乘
方A.^P,表示A的每个元素的P次乘方。
2.3 矩阵的关系运算
2.3.1 矩阵的比较关系运算 2.3.2 矩阵的逻辑关系运算
2.3.1 矩阵的比较关系运算
矩阵的比较关系是针对于两个矩阵对应元 素的,所以在使用关系运算时,首先应该保证 两个矩阵的维数一致或其中一个矩阵为标量。 关系运算是对两个矩阵的对应运算进行比 较。 若关系满足,则将结果矩阵中该位置元素 置为1,否则置0。 MATLAB的各种比较关系运算如表2-2所示。
【例2-20】 两集合的并集示例。
6.取集合的单值元素
命令格式:
b = unique (a) 构成的向量。 b = unique (A,'rows') 成的矩阵。 [b,i,j] = unique (…) 向量(矩阵)中的位置。 %取集合a的不重复元素
第二章 矩阵运算
据全下标换算出单下标。 据单下标换算出全下标。
广西大学电气工程学院
Matlab程序设计
“逻辑1”标识
所谓“逻辑1”标识 法是:通过与A同样大小的 逻辑数组L中“逻辑值1”所在的位置,指出A中 元素的位置。 L的元素或是0或是1,它是“逻辑数组(Logical Array)”。是一种特殊的数据类型。 例
(3)定数线性采样法
在设定的“总点数”下,均匀采样生成一维“行”数组。 – 通用格式:x=linspace(a,b,n) [说明]该指令生成(1 ×n)数组,其作用与x=a:(b-a)/(n-1):b同。
–
广西大学电气工程学院
Matlab程序设计
一维数组的子数组寻访和赋值
例1:子数组的寻访(Address)。 例2:子数组的赋值(Assign)。
其显示结果是一样的。
广西大学电气工程学院
Matlab程序设计
特殊矩阵(P42)
非数NaN “空”数组(空阵) 全 0阵 单位阵 全 1阵 随机阵 其他特殊矩阵(表2-6)
广西大学电气工程学院
Matlab程序设计
非数NaN
按IEEE规定,0/0,∞/∞,0×∞等运算都会产生非数 (Not a Number)。该非数在MATLAB中用NaN或 nan记述。 根据IEEE数学规范,NaN具有以下性质:
A.^n
A.^p p.^A A+B A-B
A^n
A^p p^A A+B A-B
Matlab程序设计
数组运算和矩阵运算指令对照汇总(二)
数 指 令 组 运 含
对应元素相乘 A的元素被B的对应元素除 (一定与上相同) 以自然数e为底,分别以A的 元素为指数,求幂 对A的各元素求对数 对A的各元素求平方根 求A各个元素的函数值。f(.) 表示为上节所列各函数 A、B阵对应元素间的关系 运算。#代表关系运算符
MATLAB教程及实训
4 7
5 8
6 ,则a(:,end)是指
9
A. 所有元素 B. 第一行元素 C. 第三列元素 D. 第三行元素
答案: C
3.数组的赋值
数组的赋值包括全下标方式、单下标方式 和全元素方式。
• 全下标方式:a(i,j,k…)=b,给a数组的部分 元素赋值,则b数组的行列数必须等于a数组 的行列数。
表示范围 0~28 -1 0~216 -1 0~232 -1 0~264 -1 2-7~27 -1 2-15~215 -1 2-31~231 -1 2-63~263 -1
字节数 1 2 4 8 1 2 4 8
类型转换函数 uint8() uint16() uint32() uint64() int8() int16() int32() int64()
第2章 MATLAB基本运算
2.1 数据类型 2.2 矩阵和数组的算术运算 2 .3 字符串 2.4 日期和时间 2.5 结构体和元胞数组 2.6 多维数组 2.7 关系运算和逻辑运算 2.8 数组的信息获取 2.9 多项式
2.1数据类型
MATLAB 定义了15种基本的数据类型,包括整型、浮 点型、字符型和逻辑型等,用户甚至可以定义自己的数 据类型。
2. 字符串与数值的转换 abs将字符串转换为ASCII码数值 str2num将字符串转换为数值 str2double将元胞字符串数组转换为数值
matlab课件--第2讲-数组和矩阵
Matlab 软件实习
矩阵的创建
1) 矩阵的直接定义
键入:A=[1 2 3;4 5 6]
输出:A= 123 456
或键入:A=[1 2 3 4 5 6]
第9页,共66页。
Matlab 软件实习
>> A=[1 2 3;4 5 6;7 8 9]
A= 123 456
789
第10页,共66页。
Matlab 软件实习
第15页,共66页。
Matlab 软件实习
c. 三角矩阵
命令
triu(A) triu(A,k) tril(A) tril(A,k)
运行结果
生成一个和A维数相同的上三角矩阵。该矩 阵主对角线及以上元素取自A中相应元素。 其余元素为0。
生成一个和A维数相同的上三角矩阵。该矩 阵第k条对角线及以上元素取自A中相应元 素。其余元素为0。
max(v) min (v) sum(v) mean(v) sort(v)
求最大值
求最小值
求和
求平均值
按升序排列
第5页,共66页。
Matlab 软件实习
矩阵的创建
标量的创建
直接输入:
>>x = 7 x=
7
第6页,共66页。
Matlab 软件实习
行、列向量的创建
1、逐个元素输入法
❖ >> x=[2 pi/2 3+5i]
第20页,共66页。
Matlab 软件实习
用于专门学科的特殊矩阵
(1) 魔方矩阵 魔方矩阵有一个有趣的性质,其每行、
每列及两条对角线上的元素和都相等。对于 n阶魔方阵,其元素由1,2,3,…,n2共n2个整数 组成。MATLAB提供了求魔方矩阵的函数 magic(n),其功能是生成一个n阶魔方阵。
自考线性代数经管类要点考点
线性代数〔经管类〕考点逐个击破第一章 行列式〔一〕行列式的定义行列式是指一个由假设干个数排列成同样的行数与列数后所得到的一个式子,它实质上表示把这些数按一定的规则进展运算,其结果为一个确定的数.1.二阶行列式由4个数)2,1,(=j i a ij 得到以下式子:11122122a a a a 称为一个二阶行列式,其运算规则为2.三阶行列式由9个数)3,2,1,(=j i a ij 得到以下式子:333231232221131211a a a a a a a a a称为一个三阶行列式,它如何进展运算呢.教材上有类似于二阶行列式的所谓对角线法,我们采用递归法,为此先要定义行列式中元素的余子式及代数余子式的概念.3.余子式及代数余子式设有三阶行列式 3332312322211312113a a a a a a a a a D =对任何一个元素ij a ,我们划去它所在的第i 行及第j 列,剩下的元素按原先次序组成一个二阶行列式,称它为元素ij a 的余子式,记成ij M例如 3332232211a a a a M =,3332131221a a a a M =,2322131231a a a a M =再记 ij ji ij M A +-=)1( ,称ij A 为元素ij a 的代数余子式.例如 1111M A =,2121M A -=,3131M A = 则 ,三阶行列式3D 定义为我们把它称为3D 按第一列的展开式,经常简写成∑∑=+=-==3111131113)1(i i i i i i i M a A aD4.n 阶行列式一阶行列式 11111a a D ==n 阶行列式 1121211111212222111211n n nnn n n nn A a A a A a a a a a a a a a a D +++==其中(,1,2,,)ij A i j n =为元素ij a 的代数余子式.5.特殊行列式上三角行列式111212*********n n nn nn a a a a a a a a a =下三角行列式1122112212000nn n n nn a a a a a a a a a =21对角行列式1122112200000nn nna a a a a a =〔二〕行列式的性质性质1 行列式和它的转置行列式相等,即TD D =性质2 用数k 乘行列式D 中*一行〔列〕的所有元素所得到的行列式等于kD ,也就是说,行列式可以按行和列提出公因数.性质3 互换行列式的任意两行〔列〕,行列式的值改变符号.推论1 如果行列式中有*两行〔列〕一样,则此行列式的值等于零.推论2 如果行列式中*两行〔列〕的对应元素成比例,则此行列式的值等于零. 性质4 行列式可以按行〔列〕拆开.性质5 把行列式D 的*一行〔列〕的所有元素都乘以同一个数以后加到另一行〔列〕的对应元素上去,所得的行列式仍为D.定理1〔行列式展开定理〕n 阶行列式nija D =等于它的任意一行〔列〕的各元素与其对应的代数余子式的乘积的和,即),,2,1(2211n i A a A a A a D in in i i i i =+++=或),,2,1(2211n j A a A a A a D nj nj j j j j =+++=前一式称为D 按第i 行的展开式,后一式称为D 按第j 列的展开式.本定理说明,行列式可以按其任意一行或按其任意一列展开来求出它的值.定理2 n 阶行列式nij a D =的任意一行〔列〕各元素与另一行〔列〕对应元素的代数余子式的乘积之和等于零.即)(02211k i A a A a A a kn in k i k i ≠=+++ 或)(02211s j A a A a A a ns nj s j s j ≠=+++〔三〕行列式的计算行列式的计算主要采用以下两种根本方法:〔1〕利用行列式性质,把原行列式化为上三角〔或下三角〕行列式再求值,此时要注意的是,在互换两行或两列时,必须在新的行列式的前面乘上〔-1〕,在按行或按列提取公因子k 时,必须在新的行列式前面乘上k.〔2〕把原行列式按选定的*一行或*一列展开,把行列式的阶数降低,再求出它的值,通常是利用性质在*一行或*一列中产生很多个“0〞元素,再按这一行或这一列展开:例1 计算行列式 52072325121314124-=D解:观察到第二列第四行的元素为0,而且第二列第一行的元素是112=a ,利用这个元素可以把这一列其它两个非零元素化为0,然后按第二列展开.例2 计算行列式 ab b b b a b b b b a b bb b a D =4解:方法1 这个行列式的元素含有文字,在计算它的值时,切忌用文字作字母,因为文字可能取0值.要注意观察其特点,这个行列式的特点是它的每一行元素之和均为b a 3+〔我们把它称为行和一样行列式〕,我们可以先把后三列都加到第一列上去,提出第一列的公因子b a 3+,再将后三行都减去第一行: 方法2 观察到这个行列式每一行元素中有多个b ,我们采用“加边法〞来计算,即是构造一个与4D 有一样值的五阶行列式:这样得到一个“箭形〞行列式,如果b a =,则原行列式的值为零,故不妨假设b a ≠,即0≠-b a ,把后四列的ba -1倍加到第一列上,可以把第一列的〔-1〕化为零. 例3 三阶德蒙德行列式 ))()((1112313122322213213x x x x x x x x x x x x V ---== 〔四〕克拉默法则定理1〔克拉默法则〕设含有n 个方程的n 元线性方程组为如果其系数行列式0≠=nija D ,则方程组必有唯一解:n j DD x j j ,,2,1, ==其中j D 是把D 中第j 列换成常数项n b b b ,,,21 后得到的行列式. 把这个法则应用于齐次线性方程组,则有定理2 设有含n 个方程的n 元齐次线性方程组如果其系数行列式0≠D ,则该方程组只有零解:021====n x x x换句话说,假设齐次线性方程组有非零解,则必有0=D ,在教材第二章中,将要证明,n 个方程的n 元齐次线性方程组有非零解的充分必要条件是系数行列式等于零.第二章 矩阵〔一〕矩阵的定义1.矩阵的概念由n m ⨯个数),,2,1;,,2,1(n j m i a ij ==排成的一个m 行n 列的数表 称为一个m 行n 列矩阵或n m ⨯矩阵当n m =时,称()nn ija A ⨯=为n 阶矩阵或n 阶方阵元素全为零的矩阵称为零矩阵,用n m O ⨯或O 表示2.3个常用的特殊方阵:①n 阶对角矩阵是指形如 ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=nn a a a A0000002211的矩阵 ②n 阶单位方阵是指形如 ⎪⎪⎪⎪⎪⎭⎫⎝⎛=100010001 n E 的矩阵③n 阶三角矩阵是指形如 ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛nn n n nn n n a a a a a a a a a a a a2122211122211211000,000的矩阵 3.矩阵与行列式的差异矩阵仅是一个数表,而n 阶行列式的最后结果为一个数,因而矩阵与行列式是两个完全不同的概念,只有一阶方阵是一个数,而且行列式记号“*〞与矩阵记号“()*〞也不同,不能用错.〔二〕矩阵的运算1.矩阵的同型与相等设有矩阵n m ij a A ⨯=)(,λ⨯=k ij b B )(,假设k m =,λ=n ,则说A 与B 是同型矩阵.假设A 与B 同型,且对应元素相等,即ij ij b a =,则称矩阵A 与B 相等,记为B A =因而只有当两个矩阵从型号到元素全一样的矩阵,才能说相等.2.矩阵的加、减法设n m ij a A ⨯=)(,n m ij b B ⨯=)(是两个同型矩阵则规定注意:只有A 与B 为同型矩阵,它们才可以相加或相减.由于矩阵的相加表达为元素的相加,因而与普通数的加法运算有一样的运算律.3.数乘运算设n m ij a A ⨯=)(,k 为任一个数,则规定n m ij ka kA ⨯=)(故数k 与矩阵A 的乘积就是A 中所有元素都乘以k ,要注意数k 与行列式D 的乘积,只是用k 乘行列式中*一行或*一列,这两种数乘截然不同.矩阵的数乘运算具有普通数的乘法所具有的运算律.4.乘法运算设k m ij a A ⨯=)(,n k ij b B ⨯=)(,则规定n m ij c AB ⨯=)(其中kj ik j i j i ij b a b a b a c +++= 2211),,2,1;,,2,1(n j m i ==由此定义可知,只有当左矩阵A 的列数与右矩阵B 的行数相等时,AB 才有意义,而且矩阵AB 的行数为A 的行数,AB 的列数为B 的列数,而矩阵AB 中的元素是由左矩阵A 中*一行元素与右矩阵B 中*一列元素对应相乘再相加而得到.故矩阵乘法与普通数的乘法有所不同,一般地: ①不满足交换律,即BA AB ≠②在0=AB 时,不能推出0=A 或0=B ,因而也不满足消去律.特别,假设矩阵A 与B 满足BA AB =,则称A 与B 可交换,此时A 与B 必为同阶方阵. 矩阵乘法满足结合律,分配律及与数乘的结合律.5.方阵的乘幂与多项式方阵设A 为n 阶方阵,则规定m A AAA =m 个特别E A =0又假设1110()m m m m f x a x a x a x a --=++++,则规定称)(A f 为A 的方阵多项式,它也是一个n 阶方阵6.矩阵的转置设A 为一个n m ⨯矩阵,把A 中行与列互换,得到一个m n ⨯矩阵,称为A 的转置矩阵,记为TA ,转置运算满足以下运算律:A A T =T )(,T T TB A B A +=+)(,T T kA kA =)(,T T T A B AB =)(由转置运算给出对称矩阵,反对称矩阵的定义设A 为一个n 阶方阵,假设A 满足A A T =,则称A 为对称矩阵,假设A 满足A A T -=,则称A 为反对称矩阵.7.方阵的行列式矩阵与行列式是两个完全不同的概念,但对于n 阶方阵,有方阵的行列式的概念.设)(ij a A =为一个n 阶方阵,则由A 中元素构成一个n 阶行列式nij a ,称为方阵A 的行列式,记为A方阵的行列式具有以下性质:设A ,B 为n 阶方阵,k 为数,则①A A T =; ②A k kA n = ③B A AB ⋅=〔三〕方阵的逆矩阵1.可逆矩阵的概念与性质设A 为一个n 阶方阵,假设存在另一个n 阶方阵B ,使满足E BA AB ==,则把B 称为A 的逆矩阵,且说A 为一个可逆矩阵,意指A 是一个可以存在逆矩阵的矩阵,把A 的逆矩阵B 记为1-A ,从而A 与1-A 首先必可交换,且乘积为单位方阵E.逆矩阵具有以下性质:设A ,B 为同阶可逆矩阵,0≠k 为常数,则①1-A 是可逆矩阵,且A A =--11)(; ②AB 是可逆矩阵,且111)(---=A B AB ;③kA 是可逆矩阵,且111)(--=A kkA ④TA 是可逆矩阵,且T TA A )()(11--=⑤可逆矩阵可从矩阵等式的同侧消去,即设P 为可逆矩阵,则B A PB PA =⇔=B A BP AP =⇔=2.伴随矩阵设)(ij a A =为一个n 阶方阵,ij A 为A 的行列式n ij a A =中元素ij a 的代数余子式,则矩阵⎪⎪⎪⎪⎪⎭⎫ ⎝⎛nn nn n n A A A A A A A A A212221212111称为A 的伴随矩阵,记为*A 〔务必注意*A 中元素排列的特点〕伴随矩阵必满足1*-=n AA 〔n 为A 的阶数〕3.n 阶阵可逆的条件与逆矩阵的求法定理:n 阶方阵A 可逆⇔0≠A ,且*11A AA=- 推论:设A ,B 均为n 阶方阵,且满足E AB =,则A ,B 都可逆,且B A =-1,A B =-1例1 设⎪⎪⎭⎫⎝⎛=d c b a A〔1〕求A 的伴随矩阵*A〔2〕a ,b ,c ,d 满足什么条件时,A 可逆.此时求1-A解:〔1〕对二阶方阵A ,求*A 的口诀为“主交换,次变号〞即⎪⎪⎭⎫⎝⎛--=a c b d A *〔2〕由bc ad dc b a A -==,故当0≠-bc ad 时,即0≠A ,A 为可逆矩阵此时⎪⎪⎭⎫⎝⎛---==-a c b d bc ad A A A 11*1〔四〕分块矩阵1.分块矩阵的概念与运算对于行数和列数较高的矩阵,为了表示方便和运算简洁,常用一些贯穿于矩阵的横线和纵线把矩阵分割成假设干小块,每个小块叫做矩阵的子块,以子块为元素的形式上的矩阵叫做分块矩阵.在作分块矩阵的运算时,加、减法,数乘及转置是完全类似的,特别在乘法时,要注意到应使左矩阵A 的列分块方式与右矩阵B 的行分块方式一致,然后把子块当作元素来对待,相乘时A 的各子块分别左乘B 的对应的子块.2.准对角矩阵的逆矩阵形如 ⎪⎪⎪⎪⎪⎭⎫⎝⎛r A A A 21的分块矩阵称为准对角矩阵,其中r A A A ,,,21 均为方阵空白处都是零块.假设r A A A ,,,21 都是可逆矩阵,则这个准对角矩阵也可逆,并且〔五〕矩阵的初等变换与初等方阵1.初等变换对一个矩阵A 施行以下三种类型的变换,称为矩阵的初等行〔列〕变换,统称为初等变换,〔1〕交换A 的*两行〔列〕;〔2〕用一个非零数k 乘A 的*一行〔列〕;〔3〕把A 中*一行〔列〕的k 倍加到另一行〔列〕上.注意:矩阵的初等变换与行列式计算有本质区别,行列式计算是求值过程,用等号连接,而对矩阵施行初等变换是变换过程用“→〞连接前后矩阵.初等变换是矩阵理论中一个常用的运算,而且最常见的是利用矩阵的初等行变换把矩阵化成阶梯形矩阵,以至于化为行简化的阶梯形矩阵.2.初等方阵由单位方阵E 经过一次初等变换得到的矩阵称为初等方阵.由于初等变换有三种类型,相应的有三种类型的初等方阵,依次记为ij P ,)(k D i 和)(k T ij ,容易证明,初等方阵都是可逆矩阵,且它们的逆矩阵还是同一类的初等方阵.3.初等变换与初等方阵的关系设A 为任一个矩阵,当在A 的左边乘一个初等方阵的乘积相当于对A 作同类型的初等行变换;在A 的右边乘一个初等方阵的乘积相当于对A 作同类型的初等列变换.4.矩阵的等价与等价标准形假设矩阵A 经过假设干次初等变换变为B ,则称A 与B 等价,记为B A ≅ 对任一个n m ⨯矩阵A ,必与分块矩阵⎪⎪⎭⎫⎝⎛O O O E r 等价,称这个分块矩阵为A 的等价标准形.即对任一个n m ⨯矩阵A ,必存在n 阶可逆矩阵P 及n 阶可逆矩阵Q ,使得 ⎪⎪⎭⎫⎝⎛=O O O E PAQ r5.用初等行变换求可逆矩阵的逆矩阵设A 为任一个n 阶可逆矩阵,构造n n 2⨯矩阵〔A ,E 〕 然后 ),(),(1-→A E E A注意:这里的初等变换必须是初等行变换.例2 求⎪⎪⎪⎭⎫⎝⎛----=421412311A 的逆矩阵解:则 ⎪⎪⎪⎭⎫⎝⎛----=-1132141241A例3 求解矩阵方程解:令⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛----=213411,421412311B A ,则矩阵方程为B AX =,这里A 即为例2中矩阵,是可逆的,在矩阵方程两边左乘1-A ,得也能用初等行变换法,不用求出1A -,而直接求B A 1-则 ⎪⎪⎪⎭⎫ ⎝⎛==-2052031B A X〔六〕矩阵的秩1.秩的定义设A 为n m ⨯矩阵,把A 中非零子式的最高阶数称为A 的秩,记为秩)(A 或)(A r零矩阵的秩为0,因而{}n m A ,m in )(0≤≤秩,对n 阶方阵A ,假设秩n A =)(,称A 为满秩矩阵,否则称为降秩矩阵.2. 秩的求法由于阶梯形矩阵的秩就是矩阵中非零行的行数,又矩阵初等变换不改变矩阵的秩.对任一个矩阵A ,只要用初等行变换把A 化成阶梯形矩阵T ,则秩(A)=秩(T)=T 中非零行的行数.3.与满秩矩阵等价的条件n 阶方阵A 满秩⇔A 可逆,即存在B ,使E BA AB ==⇔A 非奇异,即0≠A⇔A 的等价标准形为E⇔A 可以表示为有限个初等方阵的乘积 ⇔齐次线性方程组0=AX 只有零解⇔对任意非零列向量b ,非齐次线性方程组b AX =有唯一解 ⇔A 的行〔列〕向量组线性无关 ⇔A 的行〔列〕向量组为n R 的一个基⇔任意n 维行〔列〕向量均可以表示为A 的行〔列〕向量组的线性组合,且表示法唯一. ⇔A 的特征值均不为零⇔A A T 为正定矩阵.〔七〕线性方程组的消元法.对任一个线性方程组⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++mn mn m m n n n n b x a x a x a b x a x a x a b x a x a x a 22112222212*********可以表示成矩阵形式b AX =,其中n m ij a A ⨯=)(为系数矩阵,Tm b b b b ),,,(21 =为常数列矩阵,T n x x x X ),,,(21 =为未知元列矩阵.从而线性方程组b AX =与增广矩阵),(b A A =一一对应.对于给定的线性方程组,可利用矩阵的初等行变换,把它的增广矩阵化成简化阶梯形矩阵,从而得到易于求解的同解线性方程组,然后求出方程组的解.第三章 向量空间〔一〕n 维向量的定义与向量组的线性组合1.n 维向量的定义与向量的线性运算由n 个数组成的一个有序数组称为一个n 维向量,假设用一行表示,称为n 维行向量,即n ⨯1矩阵,假设用一列表示,称为n 维列向量,即1⨯n 矩阵与矩阵线性运算类似,有向量的线性运算及运算律.2.向量的线性组合设m ααα,,,21 是一组n 维向量,m k k k ,,,21 是一组常数,则称 为m ααα,,,21 的一个线性组合,常数m k k k ,,,21 称为组合系数.假设一个向量β可以表示成则称β是m ααα,,,21 的线性组合,或称β可用m ααα,,,21 线性表出.3.矩阵的行、列向量组设A 为一个n m ⨯矩阵,假设把A 按列分块,可得一个m 维列向量组称之为A 的列向量组. 假设把A 按行分块,可得一个n 维行向量组称之为A 的行向量组.4.线性表示的判断及表出系数的求法.向量β能用m ααα,,,21 线性表出的充要条件是线性方程组βααα=+++m m x x x 2211有解,且每一个解就是一个组合系数.例1 问T )5,1,1(-=β能否表示成T )3,2,1(1=α,T)4,1,0(2=α,T )6,3,2(3=α的线性组合.解:设线性方程组为 βααα=++332211x x x对方程组的增广矩阵作初等行变换: 则方程组有唯一解1,2,1321-===x x x所以β可以唯一地表示成321,,ααα的线性组合,且3212αααβ-+=〔二〕向量组的线性相关与线性无关1.线性相关性概念设m ααα,,,21 是m 个n 维向量,如果存在m 个不全为零的数m k k k ,,,21 ,使得02211=+++m m k k k ααα ,则称向量组m ααα,,,21 线性相关,称m k k k ,,,21 为相关系数.否则,称向量m ααα,,,21 线性无关.由定义可知,m ααα,,,21 线性无关就是指向量等式02211=+++m m k k k ααα 当且仅当021====m k k k 时成立.特别 单个向量α线性相关⇔0=α;单个向量α线性无关⇔0≠α2.求相关系数的方法设m ααα,,,21 为m 个n 维列向量,则m ααα,,,21 线性相关⇔m 元齐次线性方程组02211=+++m m x x x ααα 有非零解,且每一个非零解就是一个相关系数⇔矩阵),,,(21m A ααα =的秩小于m例2 设向量组123(2,1,7),(1,4,11),(3,6,3)T T Tααα=-==-,试讨论其线性相关性.解:考虑方程组0332211=++αααx x x其系数矩阵 ⎪⎪⎪⎭⎫ ⎝⎛-→⎪⎪⎪⎭⎫ ⎝⎛--==0001102013117641312),,(321αααA于是,秩32)(<=A ,所以向量组线性相关,与方程组同解的方程组为 令13=x ,得一个非零解为1,1,2321==-=x x x 则02321=++-ααα3.线性相关性的假设干根本定理定理1 n 维向量组m ααα,,,21 线性相关⇔至少有一个向量是其余向量的线性组合.即m ααα,,,21 线性无关⇔任一个向量都不能表示为其余向量的线性组合.定理2 如果向量组m ααα,,,21 线性无关,又m αααβ,,,,21 线性相关,则β可以用m ααα,,,21 线性表出,且表示法是唯一的.定理3 假设向量组中有局部组线性相关,则整体组也必相关,或者整体无关,局部必无关. 定理4 无关组的接长向量组必无关.〔三〕向量组的极大无关组和向量组的秩1.向量组等价的概念假设向量组S 可以由向量组R 线性表出,向量组R 也可以由向量组S 线性表出,则称这两个向量组等价.2.向量组的极大无关组设T 为一个向量组,假设存在T 的一个局部组S ,它是线性无关的,且T 中任一个向量都能由S 线性表示,则称局部向量组S 为T 的一个极大无关组.显然,线性无关向量组的极大无关组就是其本身.对于线性相关的向量组,一般地,它的极大无关组不是唯一的,但有以下性质:定理1 向量组T 与它的任一个极大无关组等价,因而T 的任意两个极大无关组等价. 定理2 向量组T 的任意两个极大无关组所含向量的个数一样.3.向量组的秩与矩阵的秩的关系把向量组T 的任意一个极大无关组中的所含向量的个数称为向量组T 的秩.把矩阵A 的行向量组的秩,称为A 的行秩,把A 的列向量组的秩称为A 的列秩. 定理:对任一个矩阵A ,A 的列秩=A 的行秩=秩〔A 〕此定理说明,对于给定的向量组,可以按照列构造一个矩阵A ,然后用矩阵的初等行变换法来求出向量组的秩和极大无关组.例3 求出以下向量组的秩和一个极大无关组,并将其余向量用极大无关组线性表出:解:把所有的行向量都转置成列向量,构造一个54⨯矩阵,再用初等行变换把它化成简化阶梯形矩阵()B A TT T T T =⎪⎪⎪⎪⎪⎭⎫⎝⎛--→⎪⎪⎪⎪⎪⎭⎫⎝⎛------==1000001100010100000133697446224112122111,,,,54321ααααα易见B 的秩为4,A 的秩为4,从而秩{}4,,,,54321=ααααα,而且B 中主元位于第一、二、三、五列,则相应地5321,,,αααα为向量组的一个极大无关组,而且324ααα--=〔四〕向量空间1.向量空间及其子空间的定义定义1 n 维实列向量全体〔或实行向量全体〕构成的集合称为实n 维向量空间,记作nR定义2 设V 是n 维向量构成的非空集合,假设V 对于向量的线性运算封闭,则称集合V 是nR 的子空间,也称为向量空间.2. 向量空间的基与维数设V 为一个向量空间,它首先是一个向量组,把该向量组的任意一个极大无关组称为向量空间V 的一个基,把向量组的秩称为向量空间的维数.显然,n 维向量空间nR 的维数为n ,且nR 中任意n 个线性无关的向量都是nR 的一个基.3. 向量在*个基下的坐标设r ααα,,,21 是向量空间V 的一个基,则V 中任一个向量α都可以用r ααα,,,21 唯一地线性表出,由r 个表出系数组成的r 维列向量称为向量α在此基下的坐标.第四章 线性方程组(一) 线性方程组关于解的结论定理1 设b AX =为n 元非齐次线性方程组,则它有解的充要条件是)(),(A r b A r = 定理2 当n 元非齐次线性方程组b AX =有解时,即r A r b A r ==)(),(时,则〔1〕b AX =有唯一解⇔n r =; 〔2〕b AX =有无穷多解⇔n r <.定理3 n 元齐次线性方程组0=AX 有非零解的充要条件是n r A r <=)( 推论1 设A 为n 阶方阵,则n 元齐次线性方程组0=AX 有非零解⇔0=A 推论2 设A 为n m ⨯矩阵,且n m <,则n 元齐次线性方程组必有非零解〔二〕齐次线性方程组解的性质与解空间首先对任一个线性方程组,我们把它的任一个解用一个列向量表示,称为该方程组的解向量,也简称为方程组的解.考虑由齐次线性方程组0=AX 的解的全体所组成的向量集合显然V 是非空的,因为V 中有零向量,即零解,而且容易证明V 对向量的加法运算及数乘运算封闭,即解向量的和仍为解,解向量的倍数仍为解,于是V 成为n 维列向量空间nR 的一个子空间,我们称V 为方程组0=AX 的解空间〔三〕齐次线性方程组的根底解系与通解把n 元齐次线性方程组0=AX 的解空间的任一个基,称为该齐次线性方程组的一个根底解系.当n 元齐次线性方程组0=AX 有非零解时,即n r A r <=)(时,就一定存在根底解系,且根底解系中所含有线性无关解向量的个数为r n -求根底解系与通解的方法是:对方程组0=AX 先由消元法,求出一般解,再把一般解写成向量形式,即为方程组的通解,从中也能求出一个根底解系.例1 求⎪⎩⎪⎨⎧=-++=+-+=+-+0022*********43214321x x x x x x x x x x x x 的通解解:对系数矩阵A ,作初等行变换化成简化阶梯形矩阵:42)(<=A r ,有非零解,取43,x x 为自由未知量,可得一般解为⎪⎪⎩⎪⎪⎨⎧==+-=-=4433432431,54,43x x x x x x x x x x写成向量形式,令13k x =,24k x =为任意常数,则通解为⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=1054014321k k X 可见,⎪⎪⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=1054,014321ξξ为方程组的一个根底解系. 〔四〕非齐次线性方程组1.非齐次线性方程组与它对应的齐次线性方程组〔即导出组〕的解之间的关系设b AX =为一个n 元非齐次线性方程组,0=AX 为它的导出组,则它们的解之间有以下性质: 性质1 如果21,ηη是b AX =的解,则21ηηξ-=是0=AX 的解性质2 如果η是b AX =的解,ξ是0=AX 的解,则ηξ+是b AX =的解 由这两个性质,可以得到b AX =的解的构造定理:定理 设A 是n m ⨯矩阵,且r A r b A r ==)(),(,则方程组b AX =的通解为其中*η为b AX =的任一个解〔称为特解〕,r n -ξξξ,,,21 为导出组0=AX 的一个根底解系.2.求非齐次线性方程组的通解的方法对非齐次线性方程组b AX =,由消元法求出其一般解,再把一般解改写为向量形式,就得到方程组的通解.例2 当参数a ,b 为何值时,线性方程组⎪⎪⎩⎪⎪⎨⎧-=+++=--+-=++=+++1232)3(122043214324324321ax x x x b x x a x x x x x x x x有唯一解.有无穷多解.无解.在有无穷多解时,求出通解.解:对方程组的增广矩阵施行初等行变换,把它化成阶梯形矩阵:当1≠a 时,4)(),(==A r b A r ,有唯一解; 当1,1≠=b a 时,3),(=b A r ,2)(=A r ,无解; 当1,1-==b a 时,2)(),(==A r b A r ,有无穷多解.此时,方程组的一般解为 ⎪⎪⎩⎪⎪⎨⎧==--=++-=44334324312211x x x x x x x x x x令2413,k x k x ==为任意常数,故一般解为向量形式,得方程组通解为。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
repmat
fliplr flipud flipdim
复制矩阵元素并扩展矩阵
交换矩阵左右对称位置上的元素 交换矩阵上下对称位置上的元素 获取指定的方向翻转交换矩阵元素
find
获取矩阵或者数组中非零元素的索引
2.4基本数学函数
函数的主要类别
三角函数 指数运算函数 复数运算函数
圆整和求余函数
2.2向量分析
2.2.1向量概念 从编程语言的角度上看,向量其实就是一维数组 从数学的角度上看,向量就是1×N或者N×1的矩
阵,即行向量或列向量 b1,1
b2,1
B= b3,1
∶ ∶
和B=[b1,1 b1,2 b1,3 ······ b1,n]
bn,1
2.2.2 创建向量
1. 在命令窗口逐个输入元素 例2-1:>>X=[ 1 3 pi 3+5i ]; Y=[ 1; 3; pi; 3+5i ] 2. 利用冒号运算符创建向量 X=J:INC:K J为向量的第一个元素,K为向量的最后一个元素, INC为向量元素递增的步长 J、INC、K之间必须用“:”间隔 若忽略INC,则默认的递增步长为1 INC可以为正数,也可以为负数 例2-2:>>X=1:10
>> A(end,:) ans = 5 10 15 20 25 >>I=[1 3 5];J=[2 4]; >>A(I,J) >>A([1 3 5], [2 4]) ans = 6 16 8 18 10 20
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
函数在处理参数时,是按照数组运算的规则进
行的
基本数学函数(续)
三角函数 函数
sin
sinh asin asinh cos cosh acos acosh tan
说明
正弦函数
双曲正弦函数 反正弦函数 反双曲正弦函数 余弦函数 双曲余弦函数 反余弦函数 反双曲余弦函数 正切函数
函数
tanh
atan atan2 atanh sec sech asec asech csc
>> C = A*B C = 14 32 23 >> D = [ -1; 0; 2 ]; >> F = pi*D F = -3.1416 0 6.2832
基本矩阵运算(续)
3. 矩阵除运算及线性方程组的解 在线性代数中没有矩阵的除运算,只有矩阵逆的运算,在MATLAB 中有两种矩阵除运算。 A/B — 矩阵右除,相当于 Ainv(B) A\B — 矩阵左除,相当于 inv(A)B 因此,x = A\B 是线性方程组Ax=B的解。
求方阵的行列式
求矩阵的秩 求矩阵的特征向量和特征值 对矩阵进行奇异值分解 求矩阵的范数
基本矩阵运算
A.*B 矩阵A与矩阵B对应元素相乘 A./B 矩阵A与矩阵B对应元素做除法 A.^B B的每一个元素作为A对应元素的幂 次
基本矩阵运算
1. 矩阵加、减运算 (A+B、A-B) 规则: 相加、减的两矩阵必须有相同的行和列,两矩阵 对应元素相加减。 MATLAB允许参与运算的两矩阵之一是标量,标量 与矩阵的所有元素分别进行加减操作。 例:A=[1 2 3;4 5 6] B=[3 4 5;7 8 9] C=3 A+B=[4 6 8;11 13 15] A+C=[4 5 6;7 8 9]
说明
双曲正切函数
反正切函数 四象限反正切函数 反双曲正切函数 正割函数 双曲正割函数 反正割函数 双曲反正割函数 余割函数
函数
csch
acsc acsch cot coth acot acoth
说明
双曲余割函数
反余割函数 反双曲余割函数 余切函数 双曲余切函数 反余切函数 反双曲余切函数
基本数学函数(续)
指数运算函数 函数
exp log log10 log2 pow2 指数函数 自然对数函数 常用对数函数 以2为底的对数函数 2的幂函数
说明
函数
realpow reallog realsqrt sqrt nextpow2
说明
实数幂运算函数 实数自然对数函数 实数平方根函数 平方根函数 求大于输入参数的第一个2的幂
基本矩阵运算(续)
4.矩阵乘方 A^n —— A自乘n次幂
方阵 >1的整数
例 >>a = [ 1, 2, 3 ; 4, 5, 6 ; 7, 8, 9 ]; >>a^2 ans = 30 36 42 66 81 96 102 126 150
2.3.4矩阵的数据操作
max(A) min(A) mean(A) median(A) std(A) sum(A) prod(A) cumsum cumprod sort sortrows 求矩阵A每列的最大元素 求矩阵A每列的最小元素 求矩阵A列元素的平均值 求矩阵A列元素的中值 求矩阵A元素的标准差 求矩阵A各列元素的和 求矩阵A各列元素的积 求列元素的累计和 求列元素的累计积 按升序对元素进行排序 按升序排列矩阵各行
使用索引访问矩阵元素的方法
矩阵元素的访问
A(i,j) A(I,J) A(i,:) A(:,j)
说 明
访问矩阵A的第i行第j列上的元素,其中i和j为标量 访问由向量I和J指定的矩阵A中的元素 访问矩阵A中第i行的所有元素 访问矩阵A中第j列的所有元素
A( : )
A( l ) A(L)
访问矩阵A中的所有元素,将矩阵看成一个向量
矩阵元素的访问 例 A=
A(1,2) A(5)
4
1 2
10 1
5 6 9
6
13
2
17
8 7
3
2 5
7
9
10
4 1
5
14
7 5
8
18
7
11
15
19
A(1:4,5) A(:,5) A(:,end) A(17:20)'
A(2:4,2:3) 4 A([2 3 4],[2 3])
0
3
8
4
12
16
20
矩阵元素的访问
2.3.1创建矩阵
例: z=zeros(3,4) x=ones(4) y=rand(5) A=fix(100*rand(5)) e1=eye(5) m1=magic(5)
2.3.2矩阵的操作
1.直接修改 在命令行窗口中,可用↑键找到所要修改的矩阵, 用←键移动到要修改的矩阵元素上即可修改。 在数组编辑器中,可用↑、↓、←、→键找到所要 修改的矩阵元素进行修改。 2.指令修改:用A(,)= 来修改。 例: >>A = [ 1 2 0; 3 0 5; 7 8 9 ] A=1 2 0 3 0 5 7 8 9 >>A (3, 3) = 0 A=1 2 0 3 0 5 7 8 0
2.3.1创建矩阵
函 数 zeros 说 明 产生元素全为0的矩阵
ones
eye rand randn diag tril triu pascal magic
产生元素全为1的矩阵
产生单位矩阵 产生均匀分布的随机数矩阵,数值范围(0,1) 产生均值为0,方差为1的正态分布随机数矩阵 获取矩阵的对角线元素,也可生成对角矩阵 产生下三角矩阵 产生上三角矩阵 产生帕斯卡矩阵 产生幻方阵
b.^a
2.3 矩阵分析
2.3.1创建矩阵
1.通过元素列表输入 输入方式: A=[a11 a12 a13; a21 a22 a23; a31 a32 a33] 说明:元素用空格或逗号隔开,换行用分号分割 例:A=[1 2 3; 4 5 6; 7 8 9] B=[1:3; 4 5 6; 7 8 9] X = [ 2 pi/2 ; sqrt(3) 3+5i ] 2.通过函数产生矩阵
B+C=[6 7 8;10 11 12]
基本矩阵运算(续)
2.矩阵乘运算 A*B:A矩阵的列数必须等于B矩阵的行数。 s*A 或 A*s:标量可与任何矩阵相乘,标量s分别与 矩阵A每个元素相乘。 例:>> A = [ 1 2 3; 4 5 6; 7 8 0 ]; B=[ 1; 2; 3 ];
1.用冒号法输入向量x=(-3 -2 -1 0 1 2 3); 2.y1=x(abs(x)>1) 3.y2=x([1 1 1 1]) 4.x(abs(x)>1)=[]
2.2.4向量运算
向量a加向量b a+b 向量a减向量b a-b 数λ乘以向量a λ*a 向量a与b的数量积 dot(a,b) 三维向量a与b的向量积 cross(a,b) 向量a的模 norm(a) 向量a与向量b对应元素相乘 a.*b 向量a与向量b对应元素相除 a./b 向量a的元素是向量b对应元素的方幂 向量a的每个元素的k次幂 a.^k
调用函数:x=rand(n, m) 说明:n为行数,m为列数,随机数在0~1之间 例:创建6个元素的行向量a,6个元素的列向量b a=rand(1,6) b=rand(6,1) 创建随机整数向量的方法: 调用格式:x=fix(rand(1,n)*30) 例:C=fix(rand(1,5)*30)
2.2.3向量元素的操作
例:求解方程组 3x1 + x2 - x3 = 3.6 x1 + 2x2 + 4x3 = 2.1 -x1 + 4x2 + 5x3 = -1.4 >> A = [ 3 1 -1 ; 1 2 4 ; -1 4 5 ]; >> B = [ 3.6 ; 2.1 ; -1.4 ]; >> x = A\B x= 1.4818 -0.4606 0.3848