第4讲 MATLAB矩阵分析及绘图

x1=x(k);
%取y1与y2相等点的x坐标
y3=0.2*exp(-0.5*x1).*cos(4*pi*x1); %求y1与y2值相等点的y坐标 plot(x,y1,x,y2,'k:',x1,y3,'bp');
5.3 图形标注与坐标控制
1．图形标注
有关图形标注函数的调用格式为：
title(图形名称) xlabel(x轴说明) ylabel(y轴说明) text(x,y,图形说明)
例：用求逆矩阵的方法解线性方程组
x 2 y 3z 5 x 4 y 9 z 2 x 8 y 27 z 6
4.4 矩阵求值
1．方阵的行列式值 把一个方阵看作一个行列式，并对其按行列式的规则求值， 这个值就称为矩阵所对应的行列式的值。 在MATLAB中，求方阵A所对应的行列式的值的函数是 det(A)。
5.2 设置曲线样式
MATLAB提供了一些绘图选项，用于确定所绘曲线
的线型、颜色和数据点标记符号，它们可以组合使用。
线型 线方式： - 实线 :点线 -. 虚点线 - - 波折线。
线型 点方式： . 圆点 +加号 * 星号 x x形 o 小圆 颜色：b蓝 ；g绿；r红； c青； m紫； y黄； k黑； w白；
>> y=[1 2 3;4 5 6;7 8 9];
>> plot(x,y);
(2) 当x,y是同维矩阵时，则以x,y对应列元素为横、纵
坐标分别绘制曲线，曲线条数等于矩阵的列数。 >> x=[1 2 3;4 5 6;7 8 9]; >> y=[7 8 9;1 2 3;4 5 6]; >> plot(x,y); (3) 对只包含一个输入参数的plot函数，当输入参数是实矩阵 时，则按列绘制每列元素值相对其下标的曲线，曲线条数等 于输入参数矩阵的列数。 >> x=[1 2 3;4 5 6;7 8 9]; >> plot(x,y);
yi a bxi
此外，常用的函数还有reshape(A,m,n)，它在矩阵总元 素保持不变的前提下，将矩阵A重新排成m×n的二维矩阵。
2、用于专门学科的特殊矩阵
魔方矩阵
魔方矩阵有一个有趣的性质，其每行、每列及两条 对角线上的元素和都相等。对于n阶魔方阵，其元素由
1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵
第5章 MATLAB绘图
5.1 二维数据曲线图
5.1.1 绘制单根二维曲线 plot函数的基本调用格式为： plot(x,y)
其中x和y为长度相同的向量，分别用于存储x坐标和y坐标数据。
例5-1 在0≤x≤2区间内，绘制曲线 y=2e-0.5xcos(4πx) 程序如下： x=0:pi/100:2*pi; y=2*exp(-0.5*x).*cos(4*pi*x);
>> A=[1 2 3;4 5 6;7 8 9];
>> triu(A) ans = 1 2 0 5 0 0 >> tril(A,-1) ans = 0 0 0 4 0 0 7 8 0 >> tril(A,1) ans = 1 2 0 4 5 6 7 8 9
3 6 9
>> triu(A,1) ans = 0 2 3 0 0 6 0 0 0
2 1.5
1
0.5
0
-0.5
-1
plot(x,y)
-1.5
-2
0
1
2
3
4
5
6
7
例5-2 绘制曲线。 程序如下： t=0:0.1:2*pi; x=t.*sin(3*t); y=t.*sin(t).*sin(t); plot(x,y);
5
4.5
4
3.5
3
2.5
2
1.5
1
0.5
0 -6
-4
-2
0
2
4
6
plot函数最简单的调用格式是只包含一个输入参数：
纵坐标范围默认在最大最 小值之间。
4．图形保持
hold on/off命令控制是保持原有图形还是刷新原有图
形，不带参数的hold命令在两种状态之间进行切换。
例4-5 采用图形保持，在同一坐标内绘制曲线 y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。 程序如下： x=0:pi/100:2*pi; y1=0.2*exp(-0.5*x).*cos(4*pi*x); plot(x,y1) hold on y2=2*exp(-0.5*x).*cos(pi*x); plot(x,y2); hold off
plot(x)
在这种情况下，当x是实向量时，以该向量元素的下标为 横坐标，元素值为纵坐标画出一条连续曲线，这实际上是绘 制折线图。
5.1.2 绘制多根二维曲线
1．plot函数的输入参数是矩阵形式
(1) 当x是向量，y是有一维与x同维的矩阵时，则绘制出多根不 同颜色的曲线。曲线条数等于y矩阵的另一维数，x被作为这 些曲线共同的横坐标。 >> x=[1 2 3];
例：建立3*3，3*2的和矩阵A同样大小的零矩阵 Zero(3,3) Zero(3,2) Zero(size(A))
例(1)在区间[20,50]内均匀分布的5阶随机矩阵
（2）均值为0.6，方差为0.1的5阶正态分布随机矩阵
任意[a,b]区间上均匀分布的随机数： yi a (b a) xi 均值为a，方差为b^2的随机数：
MATLAB对矩阵A实施上下翻转的函数是flipud(A)。
4.3 矩阵的逆
1．矩阵的逆 对于一个方阵A，如果存在一个与其同阶的方阵B，使得： A· B=B· (I为单位矩阵)则称B为A的逆矩阵，当然，A也是B的 A=I 逆矩阵。 求一个矩阵的逆是一件非常烦琐的工作，容易出错，但在 MATLAB中，求一个矩阵的逆非常容易。求方阵A的逆矩阵可调 用函数inv(A)。 例 求方阵A的逆矩阵，并验证其互逆。
2．矩阵的秩和迹
（1）矩阵的秩 矩阵线性无关的行数与列数称为矩阵的秩。在MATLAB中， 求矩阵秩的函数是rank(A)。 （2）矩阵的迹 矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征 值之和。在MATLAB中，求矩阵的迹的函数是trace(A)。
3．矩阵的特征值与特征向量
在MATLAB中，计算矩阵A的特征值和特征向量的函数 是eig(A)，常用的调用格式有3种：
>> plot(x1,y1,x2,y2,x3,y3)
3．具有两个纵坐标标度的图形
在MATLAB中，如果需要绘制出具有不同纵坐标标 度的两个图形，可以使用plotyy绘图函数。调用格式为： plotyy(x1,y1,x2,y2)
其中x1,y1对应一条曲线，x2,y2对应另一条曲线。横坐标的标 度相同，纵坐标有两个，左纵坐标用于x1,y1数据对，右纵坐
(1) E=eig(A)：求矩阵A的全部特征值，构成向量E。
(2) [V,D]=eig(A)：求矩阵A的全部特征值，构成对角阵D， 并求A的特征向量构成V的列向量。
例 用求特征值的方法解方程。3x5-7x4+5x2+2x-18=0 p=[3,-7,0,5,2,-18]; A=compan(p); %求伴随矩阵A x1=eig(A) %求A的特征值 x2=roots(p) %直接求多项式p的零点
例5-6 在同一坐标内，分别用不同线型和颜色绘制曲线 y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)，标记两曲线交叉点。 程序如下：
x=linspace(0,2*pi,1000); y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); k=find(abs(y1-y2)<1e-2); %查找y1与y2相等点(近似相等)的下标
标用于x2,y2数据对。
例5-4 用不同标度在同一坐标内绘制曲线 y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。
程序如下：
x=0:pi/100:2*pi;
y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); plotyy(x,y1,x,y2);
的函数magic(n)，其功能是生成一个n阶魔方阵。 观察magic(5)，其元
素为1~25，每行每列
及对角线的和均为65 例： 将101~125等25个数填 入一个5行5列的表格中，使 其每行每列及对角线的和均 为565。
4.2 矩阵的结构变换
1、对角阵与三角阵
（1）对角阵 只有对角线上有非0元素的矩阵称为对角矩阵，对角线上的元素相等的 对角矩阵称为数量矩阵，对角线上的元素都为1的对角矩阵称为单位矩阵。 提取矩阵的对角线元素
2．含多个输入参数的plot函数 调用格式为： plot(x1,y1,x2,y2,…,xn,yn)
(1)
当输入参数都为向量时，x1和y1，x2和y2，…，xn和yn 分别组成一组向量对，每一组向量对的长度可以不同。 每一向量对可以绘制出一条曲线，这样可以在同一坐标
内绘制出多条曲线。
>> x1=linspace(0,2*pi,100); >> y1=sin(x1); >> x2=linspace(0,3*pi,70); >> y2=1+sin(x2); >> plot(x1,y1,x2,y2)
(2) 当输入参数有矩阵形式时，配对的x,y按对应列元
素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的 列数。
>> x1=linspace(0,2*pi,100); >> y1=sin(x1); >> x2=linspace(0,3*pi,100); >> y2=1+sin(x2); >> x3=[x1;x2]'; >> y3=[y1+2;y2+2]';