系统工具与环境（理工科用）.ppt

例如 显示灰度图像月亮(moon.tif)。 I=imread(moon.tif); imshow(I); 这个例子中，使用imread函数将图像文件读入到矩阵I中，原始图像大小可以使用size(I)来求得，然后通过调用imshow显示图像，如图11.12所示。 例如 显示RGB真彩图辣椒(peppers.png) I=imread(pappers.png); imshow(I); 结果如图11.13所示。 例如 将图像辣椒保存为灰度图像，并显示。 I=imread(peppers.png ); I=rgb2gray(I); imwrite(I, mypepper.tif ); imshow(I); 本例中，函数rgb2gray的作用是将矩阵I所表示的图像转换为灰度图像，在MATLAB中有许多图像类型转换的函数，读者可进一步参考相关资料。 2. 随机矩阵和魔方矩阵 使用函数rand()可以生成一个0~1之间的随机数，而使用rand(m,n)则可以生成一个m×n阶随机矩阵，使用rand(m)生成的是m×m阶随机方阵 魔方矩阵是指生成的方阵行、列、对角线之和相等的矩阵，通过函数magic(m)来实现。 例如 生成2×3阶随机方阵以及3阶魔方矩阵。 rand(2,3) ans = 0.7655 0.1869 0.4456 0.7952 0.4898 0.6463 magic(3) ans = 8 1 6 3 5 7 4 9 2 3. 对角矩阵 对角矩阵的生成是通过调用函数diag来实现的，其基本格式为 (1) diag(p)，生成以向量p的元素构成的对角矩阵，矩阵大小由p的元素个数决定。 (2) diag(p, k)，将向量p的元素分布在对角线偏上或偏下的斜列上，k0，代表上部，k0代表下部，k=0时与(1)相同，生成的矩阵阶数=p的元素个数+k。 3、矩阵行列式、秩、逆和条件数 在MATLAB中，通过调用函数det、rank、inv、cond即可分别求的矩阵的行列式的值、秩、矩阵的逆和条件数。 4. 线性方程组求解 考虑线性方程组 ，很显然 ，很显然，利用矩阵的逆可以直接求解线性方程组，这也称为直接法，也可以直接利用之前介绍的反除，即 。 5. 矩阵的三角分解 根据线性代数理论，如果矩阵A是非奇异矩阵，则它一定可以分解为一个非奇异的单位下三角矩阵和非奇异上三角矩阵，使得 成立，但由于需要涉及到选主元来交换某些行，因此还需要增加一个置换矩阵P来构成完整的三角分解 ，在MATLAB中，通过调用函数lu来实现上述过程。 6. 矩阵特征多项式和特征值 求矩阵特征多项式和特征值是线性代数中经常遇到的问题，在MATLAB中通过调用函数eig来实现，其调用格式为 (1) v=eig(A)，求矩阵A的特征值，返回形式为由特征值组成的列向量。 (2) [b,v]=eig(A)，返回结果b是由特征向量组成的矩阵，而v为以特征值为对角线的对角矩阵。 方程求根 1. 函数零点 罗尔(Role)零点定理证明了对于 ，且 ，则必存在区间 内的一点 ，使得 ， 就是函数的零点。在MATLAB中使用函数fzero很容易求得零点，使用方式为 (1) fzero(myfun, x0)，表示求函数fun在x0附近的零点，其中fun为函数句柄。 (2) fzero(myfun, [a b])，表示求函数fun在区间[a,b]上的零点。 二维图形 在 MATLAB 中，plot是最基本的二维绘图函数，其调用格式有 (1) plot(Y)：若Y为实向量，则以该向量元素的下标为横坐标，以Y的各元素值为纵坐标，绘制二维曲线；若Y为复数向量，则等效于plot(real(Y),imag(Y))；若Y为实矩阵，则按列绘制每列元素值相对其下标的二维曲线，曲线的条数等于Y的列数；若Y为复数矩阵，则按列分别以元素实部和虚部为横、纵坐标绘制多条二维曲线。