ch2 Matlab矩阵及生成及操作.pptVIP

* * * * * * * * * * * * * * * Matlab矩阵的生成与操作 安徽工业大学数理学院 侯为根 wghou@ahut.edu.cn 【例1-1】表示矩阵 一、MATLAB矩阵的生成 1、直接输入法 将矩阵的元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔。 A=[1,2,3;4,5,6;7,8,0] A=[1,2,3;4 5,6;7,8 0] 【例1-2】 试输入复数矩阵 需要避免的语句 B=[1+9*i,2+8*i,3+7*j; 4+6*j 5+5*i,6+4*i;7+3*i,8+2*i,i] B=[1 +9*i,2+8*i,3+7*j; 4+6*j 5+5*i,6+4*i;7+3*i,8+2*i,i] 2、利用M文件建立矩阵 对于比较大且比较复杂的矩阵，可以为它专门建立一个M文件。 【例1-3】利用M文件建立mymat矩阵。 (1) 启动有关编辑程序或MATLAB文本编辑器，并输入待建矩阵. (2) 把输入的内容以纯文本方式存盘(设文件名为mymat.m)。 (3) 运行该M文件，就会自动建立一个名为mymat的矩阵，可供以后使用。 A=[4 10 1 6 2;8 2 9 4 7; 7 5 7 1 5;0 3 4 5 4;23 13 13 0 3] 3、利用MATLAB函数建立矩阵 几个产生特殊矩阵的函数：zeros、ones、 eye、rand、randn、pascal、 magic(n)、 vander(V)、 hilb(n) 、 toeplitz(x,y)、compan(P) 。 这几个函数的调用格式相似，下面以产生零矩阵的zeros函数为例进行说明。其调用格式是： zeros(m) 产生m×m零矩阵 zeros(m,n) 产生m×n零矩阵。 zeros(size(A)) 产生与矩阵A同样大小的零矩阵 相关的函数有：length(A)给出行数和列数中的较大者，即length(A)=max(size(A))；ndims(A)给出A的维数。 指令 说 明 zeros(m,n) 产生阶为 m×n, 元素全为0的矩阵 ones(m,n) 产生阶为 m×n, 元素全为1的矩阵 eye(n) 产生阶为 n×n的单位阵 pascal(m,n) 产生阶为 m×n 的 Pascal 矩阵 vander(m,n) 产生阶为 m×n 的 Vandermonde 矩阵 hilb(n) 产生阶为 n×n 的 Hilbert 矩阵 rand(m,n) 产生 [0,1]均匀分布的随机数矩阵,其阶为 m×n randn(m,n) 产生μ=0,? =1的正态分布随机数矩阵,其阶为 m×n magic(n) 产生阶为n×n的魔方阵 diag(v) 产生以向量v为对角元素的对角阵 【例1-4】 分别建立3×3、3×2和与矩阵A同样大小的零矩阵。 (1) 建立一个3×3零矩阵：zeros(3) (2) 建立一个3×2零矩阵：zeros(3,2) (3) 建立与矩阵A同样大小零矩阵：zeros(size(A)) 此外，常用的函数还有reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m×n的二维矩阵。 【例1-5】 将4阶魔方阵形变为2行8列的矩阵 reshaps(magic(4),2,8) hilb(n)指令可以产生n×n的 Hilbert 矩阵 Hilbert矩阵 Hilbert 矩阵的特性: 当矩阵变大时，其矩阵会接近奇异（即矩阵的行列式会接近于0），Hilbert矩阵常被用来评估各种逆矩阵计算方法的稳定性。 【例1-6】 计算6阶hilbert矩阵的行列式 det(hilb(6)) x1 = rand(10000, 1); x2 = randn(10000, 1); subplot(2,1,1); hist(x1, 40); title(均匀分布); subplot(2,1,2); hist(x2, 40); title(高斯分布); set(findobj(gcf, type, patch),EdgeColor, w); % 改边线为白色 【例1-7】产生10000个均匀均匀与正态分布的随机数 产生随机数矩阵指令rand 和 randn 4、建立大矩阵 大矩阵可由方括号中的小矩阵建立起来。 A=[1 2 3;4 5 6;7 8 0]; 5、冒号表达式 冒号表达式的一般格式：v=s1:s2:s3 还可以用linspace函数产生行向量：linspace(a,b,n) linspace(a,b,n)与a:(b-a)/(n-1):b