第4章 MATLAB程序设计

第4章 MATLAB程序流程控制习题4一、选择题1．下列关于脚本文件和函数文件的描述中不正确的是（)。

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

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

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

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

CA．6 6 B．2 1 C．2 2 D．1 26．调用函数时,如果函数文件名与函数名不一致，则使用（）。

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

BA．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．执行语句“fn=＠(x）10*x；”，则fn是（）。

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

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

第四章1.a=input('请输入一个4位数:');while (a<1000|a>9999)a=input('输入错误，请重新输入一个4位数:'); endb=fix(a/1000);c=rem(fix(a/100),10);d=rem(fix(a/10),10);e=rem(a,10);b=b+7;c=c+7;d=d+7;e=e+7;b=rem(b,10);c=rem(c,10);d=rem(c,10);e=rem(e,10);g=b;b=d;d=g;g=c;c=e;e=g;a=1000*b+100*c+10*d+e;disp(['加密后：',num2str(a)])2.逻辑表达式法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');x=0.5:1:5.5;x1=(x>=0.5&x<1.5);x2=(x>=1.5&x<3.5);x3=(x>=3.5&x<=5.5);y1=a.*(x.^2)+b.*x+c;y2=a*(sin(b)^c)+x;y3=log(abs(b+c./x));y=y1.*x1+y1.*x2+y3.*x3; disp(y)if语句法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');for x=0.5:1:5.5if x>=0.5 & x<1.5y=a.*(x.^2)+b.*x+celseif x>=1.5 & x<3.5y=a*(sin(b)^c)+xelseif x>=3.5 & x<5.5y=log(abs(b+c./x))endendswitch语句法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');for x=0.5:1:5.5switch floor(x/0.5)case {1,2}y=a.*(x.^2)+b.*x+c;case {3,4,5,6}y=a*(sin(b)^c)+x;case {7,8,9,10}y=log(abs(b+c./x)); enddisp(y)end3.x=fix(rand(1,20)*89)+10;x1=mean(x);n=find(rem(x,2)==0 & x<x1);disp(['小于平均数的偶数是：',num2str(x(n))]);4.(1)A=input('请输入20个数的一个行向量：');a=A(1);b=A(1);for m=Aif a>=ma=m;elseif b<=mb=m;endenddisp(['最小数是：',num2str(a)])disp(['最大数是：',num2str(b)])(2)A=input('请输入20个数的一个行向量：'); maxval=max(A)minval=min(A)5.s=0;for a=0:63c=2^a;s=s+c;enddisp(['2的0次方到63次方的和是：',num2str(s)])k=0:63n=2.^ks=sum(n)6.（1）sum1=0;for n=1:100x=(-1)^(n+1)*(1/n);sum1=sum1+x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;for n=1:1000x=(-1)^(n+1)*(1/n);sum2=sum2+x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;for n=1:10000x=(-1)^(n+1)*(1/n);sum3=sum3+x;enddisp(['当n取10000时：sum=',num2str(sum3)])（2）sum1=0;n1=0;for n=1:2:100x=(-1)^n1*(1/n);sum1=sum1+x;n1=n1+1;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;n2=0;for n=1:2:1000x=(-1)^n2*(1/n);sum2=sum2+x;n2=n2+1;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;n3=0;for n=1:2:10000x=(-1)^n3*(1/n);sum3=sum3+x;n3=n3+1;enddisp(['当n取10000时：sum=',num2str(sum3)])（3）sum1=0;for n=1:100x=1/(4^n);sum1=sum1+x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;for n=1:1000x=1/(4^n);sum2=sum2+x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;for n=1:10000x=1/(4^n);sum3=sum3+x;enddisp(['当n取10000时：sum=',num2str(sum3)])（4）sum1=1;for n=1:100x=4*n*n/(2*n-1)/(2*n+1);sum1=sum1*x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=1;for n=1:1000x=4*n*n/(2*n-1)/(2*n+1);sum2=sum2*x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=1;for n=1:10000x=4*n*n/(2*n-1)/(2*n+1);sum3=sum3*x;enddisp(['当n取10000时：sum=',num2str(sum3)])7.函数文件function f=fibnacci(n)if n==1 | n==2f=1;elsef=fibnacci(n-1)+fibnacci(n-2);end命令文件：shulie=[];for k=1:nshulie=[shulie fibnacci(k)];endshulie8.function [f1,f2]=juzhenji(x1,x2)f1=x1*x2;f2=x1.*x2;命令文件：clear alla=input('请输入一个矩阵：');b=input('请再输入一个矩阵：(注意：两矩阵要可以相乘)'); [f1,f2]=juzhenji(a,b);disp(f1)disp(f2)9.function sum=qiuhe(n,m)if n<=1sum=0;elsesum=n^m+qiuhe(n-1,m);end命令文件：clear ally=qiuhe(100,1)+qiuhe(50,2)+qiuhe(10,-1);disp(y)10.s=0；a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=afor j=1:4if rem(k(j),2)~=0 s=s+k(j);endendendss =108（2）global xx=1:2:5;y=2:2:6;sub(y);xyfunction fun=sub(z) global xz=3*x;x=x+z;x =4 12 20 y =2 4 6。

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中， 转置运算符是单撇号(')。

函数文件中的第一个可执行的行一定是函数 定义行，否则该文件就是一个脚本文件，函 数定义行如下：
将一个文件定义为函数文件 定义了函数的名字 定义了输入参数和输出参数的个数和顺序
25
4.3.1 函数定义行
函数定义行的格式如下:
function[output arguments]=function_name(input arguments)
27
4.3.2 输入和输出参数
输入和输出参数可以用来将数据输入到函数 中或从函数中输出数据。
输入参数写在函数名字后面的括号内，通常 情况下，至少有一个输入参数，虽然有的时 候可以没有输入参数，如果有多个输入参数 ，则参数之间用逗号隔开
28
4.3.2 输入和输出参数
函数文件中用于执行计算的代码是依据输入参数来 编写的，并假设这些参数已经被赋予了数值
31
4.3.2 输入和输出参数
另外，input命令可以用来交互式的输入数据， disp,fprintf,和plot命令可以用来在屏幕上显示输 出，将输出保存至文件中，在脚本文件中画图。表 4-1是有不同输入和输出参数组合的函数定义行示 例：
函数定义行
function
[mpay,tpay]
loan(amount,rate,years)
8
4.1.2 脚本文件的创建与保存
一旦打开窗口，脚本文件中的命令是逐 行输入的，敲击回车键，MATLAB会自 动编号一个新的行
MATLAB的命令可以在任意一个文本 编辑器或文本处理程序中进行编辑，然 后复制粘贴到编辑/调试窗口，编辑/调 试窗口中的小程序如图4-2所示。
9
4.1.2 脚本文件的创建与保存 脚本文件的前几行是注释（不会被执行