MATLAB 第四章

v
i 1
n
2 i
。

max { vi } 。
1 ≤i ≤n
设 A 是一个 m ×n 的矩阵，矩阵的 3 种常用范数如下。 1-范数： A 1 max { aij } 。
1 ≤ j ≤n i 1 m
2-范数： A 2 1 ，其中 λ 1 为 A'A 最大特征值。 ∞-范数： A max { aij } 。
【例4.6】先建立5 × 5矩阵A，然后将A的第一行元素乘以1， 第二行乘以2，…，第五行乘以5。 用一个对角矩阵左乘一个矩阵时，相当于用对角阵的第一个 元素乘以该矩阵的第一行，用对角阵的第二个元素乘以该 矩阵的第二行……依此类推，因此，只需按要求构造一个 对角矩阵D，并用D左乘A即可。命令如下： A=[1:5;2:6;3:7;4:8;5:9] D=diag(1:5); D*A %用D左乘A,对A的每行乘以一个指定常数
（2）构造对角矩阵 设V为具有m个元素的向量，diag(V,k)的功能是产生一个 n × n（n = m + k|）对角阵，其第k条对角线的元素即为 向量V的元素。 例如： diag(1:3,-1) ans = 0 0 0 0 1 0 0 0 0 2 0 0 0 0 3 0 省略k时，相当于k为0，其主对角线元素即为向量V的元素。
2．矩阵的秩与迹 （1）矩阵的秩 rank(A) （2）矩阵的迹 矩阵的迹即矩阵的对角线元素之和。 trace(A)。
3．向量和矩阵的范数
设向量 V = (v1 ,v2 ,…,vn )，向量的 3 种常用范数如下。 1-范数： V 2-范数： V ? -范数： V
1
vi 。
i 1
n
2

3．矩阵的转置 所谓转置，即把源矩阵的第一行变成目标矩阵第一列，第二 行变成第二列……依此类推。显然，一个m行n列的矩阵 经过转置运算后，变成一个n行m列的矩阵。MATLAB中， 转置运算符是单撇号(')。

第4章图形处理功能1 内容简介基本内容主要包括：（1）二维图形（2）三维图形（3）图形处理的基本技术2 达到的目标（1）掌握二维图形的绘制。

（2）掌握三维图形的绘制。

（3）掌握图形处理的基本技术3 具体内容３．1 二维图形３．1．1 基本绘图命令（1）当plot函数仅有一个输入变量例4-1y=[5 2 3 8 5]; %y 行矩阵plot(y) %一条线例4-2y=[5 2 3 8 5;2 4 3 1 5;1 1 1 1 1]; %y 矩阵plot(y) %5条线,等于矩阵的列数（2）当plot函数有两个输人变量例4-3x=0:0.01*pi:pi;y=sin(x).*cos(x);plot(x,y)例4-3x=0:0.01*pi:pi;y=[sin(x);cos(x); sin(x).*cos(x)];plot(x,y)例4-4x1=0:0.01*pi:pi;x2=pi:0.01*pi:2*pi;x=[x1' x2'];y=[sin(x1') cos(x2')];plot(x,y)例4-5x1=1:5;x2=6:10;y1=x1;y2=2*x2;plot([x1;x2],[y1;y2])%plot([x1' x2'],[y1' y2'])（3）当plot函数有三个输入变量时MATLAB语言中提供的对曲线的线型、颜色以及标识的控制符如表4．l所示。

例4-6 绘制带有显示属性设置的二维图形。

x=0.5*pi: 0.1*pi:2*pi;y=sin(x);z=cos(x);plot (x, y, '--ko', x, z, '-. r*')３．1．2 特殊的二维图形函数（1）特殊坐标系的二维图形函数（ａ）对数坐标例4-7 绘制X坐标为对数坐标的二维图形。

x=0.5*pi: 0.1*pi:2*pi;y=sin(x);semilogx (x, y, '-ro')（ｂ）极坐标例4-8绘制极坐标下的二维图形。

第4章M A T L A B程序流程控制-习题答案第4章MATLAB程序流程控制习题4一、选择题1．下列关于脚本文件和函数文件的描述中不正确的是（）。

AA．函数文件可以在命令行窗口直接运行B．去掉函数文件第一行的定义行可转变成脚本文件C．脚本文件可以调用函数文件D．函数文件中的第一行必须以function开始2．下列程序的输出结果是（）。

Dy=10;ify==10y=20;elseify>0y=30enddisp(y)A．1B．30C．10D．203．有以下语句：a=eye(5);forn=a(2:end,:)for循环的循环次数是（）。

CA．3B．4C．5D．104．设有程序段k=10;whilekk=k-1end则下面描述中正确的是（）。

AA．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段：x=reshape(1:12,3,4);m=0;n=0; fork=1:4ifx(:,k)<=6m=m+1; elsen=n+1; end end则m 和n 的值分别是（）。

CA ．66B ．21C ．22D ．126．调用函数时，如果函数文件名与函数名不一致用（）。

A A ．函数文件名B ．函数名 C ．函数文件名或函数名均可D ．@函数名 7．如果有函数声明行为“f unction[x,y,z]=f1(a,b,c)”，则下述函数调用格式中错误的是 （）。

B A ．x=f1(a,b,c)B ．[x,y,z,w]=f1(a,b,c) C ．[x,b,z]=f1(a,y,c)D ．[a,b]=f1(x,y,z) 8．执行语句“f n=@(x)10*x;”，则fn 是（）。

A A ．匿名函数B ．函数句柄C ．字符串D ．普通函数 9．执行下列语句后，变量A 的值是（）。

D >>f=@(x,y)log(exp(x+y)); >>A=f(22,3); A ．22,3B ．22C ．3D ．25 10．程序调试时用于设置断点的函数是（）。

4.1.1 M文件的分类
M文件是由若干 Matlab 命令组合在一起构成的，它可 以完成某些操作，也可以实现某种算法
事实上，Matlab 提供的内部函数以及各种工具箱，都是利用 Matlab 语言编写的 M文件 用户也可以结合自己的工作需要，开发自己的程序或工具箱
M文件根据调用方式的不同可以分为两类： Script file：命令文件/脚本文件 Function file：函数文件
例2 输入x，y的值，并将它们的值互换后输出（swap.m）。 x=input('Input x please.'); y=input('Input y please.'); z=x; x=y; y=z; disp(x); disp(y); 例3 求一元二次方程ax2 +bx+c=0的根(root.m)。 a=input('a=?'); b=input('b=?'); c=input('c=?'); d=b*b-4*a*c; x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)]; disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);
例如：
s=0; a=[12 13 14;15 16 17;18 19 20;21 22 23] for k=a s=s+k; end disp(s); 该程序的功能是求矩阵各行元素之和，执行结果是： 39 48 57 66
while语句
while expr （条件） statement（循环体语句） end 若expr成立，则执行循环体的内容，执行后 再判断条件是否为真，如果不成立则跳出循环体。