Matlab考试题库加答案.doc

一、 辨析题（共2小题，每题 5 分，共计 10 分）二、 问答题（共3小题，每题 5 分，共计 15 分）三、 计算题（共7小题，每题 5 分，共计 35 分）四、 编程题（共4小题，每题 10 分，共计 40 分）填空题（计算部分）1、标点符号; _可以使命令行不显示运算结果，%用来表示该行为注释行。2、x 为0 4pi ，步长为0.1pi 的向量，使用命令_创建。x=0:0.1*pi:4*pi3、输入矩阵A= ，使用全下标方式用A(2,2) 取出元素“-5 ”，使用单下标方式用_取出元素“-5 ”。A(5)4、符号表达式sin(2*a+t)+m 中独立的符号变量为_。t5、M脚本文件和M函数文件的主要区别是M脚本文件没有函数定义和M函数文件有函数定义_。6. 设x是一维数组，x的倒数第3个元素表示为_；设y为二维数组，要删除y的第34行和48列，可使用命令_; _；x(_end-2_) y(34,:)= y(:,48)=7. 将变量x以Ascii文本格式存储到文件fname.txt，应使用命令_ _；save _x8. 在while 表达式, 语句体, End 循环语句中，表达式的值_ _时表示循环条件为真，语句体将被执行，否则跳出该循环语句； 非零9.要从键盘读入一个字符串并赋值给变量x，且给出提示“Who is she?”，应使用命令_；x=input(Who is she?,s)_10设A=和B= 和C=均为m*n 矩阵，且存在于WorkSpace中，要产生矩阵D= ，可用命令_ _，计算 可用命令_;D=(A-C)/B.Cdet(inv(A*B)11. 在MATLAB命令窗口中的“>>”标志为MATLAB的_ 提示符，“”标志为_提示符。命令行输入12.已知A=1 2 3；4 5 0；7 8 9；B=1 0 3；1 5 0；0 1 2；写出下列各指令运行的结果。A+B；A.*B；A=B ；ans= 2,2,6;5,10,0;7,9,11ans= 1,0,9;4,25,0;0,8,18ans= 1,0,1;0,1,1;0,0,013.已知A是矩阵，求A的对角矩阵函数是_，求A的下三角矩阵函数是_ 。diag tril14.MATLAB的程序文件和Simulink模型文件的扩展名分别是_、 。.m.mdl15.MATLAB最基本的绘图函数为_ 。plot()16. A=1,2,3;4,5,6; A(:,1,3)=;A=_2；517. fix(-1.5)=_ _, round(-1.5)=_ _.-1 -218.已知A=1 3;4 6;C=A,eye(size(A);A ,ones(size(A) 则_.1 3 1 0；4 6 0 1；1 3 1 1；4 6 1 1_19. A=1,2,3;4,5,6; A(4)=_, A(3,2)=_5_ 无答案20.A=1,2;3,4;B=1,0;0,1; A*B=_, A.*B=_.1 2；3 4 1 0；0 421. A=1 2 3;6 5 4; reshape(A,3,2)=_.1 5；6 3；2 422. A=1,2;3,1; B=1,0;0,1; A=B=_.0 1；1 023. A=1,2,3;1,4,9;1,7,27; triu(A,-1)=_ _.提取矩阵上三角矩阵的函数1 2 3；1 4 9；0 7 2724. A=4,15,-45,10,6;56,0,17,-45,0; find(A>=10 & A<=20)=_ _. 3 6 725.已知A=0 9 6;1 3 0;B=1 4 3;1 5 0;写出下列各指令运行的结果。A & B 的运行结果_；A ./ B 的运行结果_。ans=0 1 1;1 1 0 ans= 0 2.25 2;1 0.6 NaN26. 产生4阶全0方阵的命令为_；产生3阶全1方阵的命令为。zeros(4) ones(3)27. A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为_ 。2 5和 528. 用if判断语句判断80 x >60,在MATLAB中if语言后的判断应写为_。if x<=80 & x>6029.为了使两个plot的图形在同一个坐标显示，可以使用_命令进行图形保持；可以使用_命令为图形添加网格。hold ongrid on30. MATLAB的工作空间中有三个变量v1, v2, v3，写出把它们保存到文件my_data.mat中的指令_； save my_date31. MATLAB提供了两种运行方式，是 按_和把_文件copy到matlab里，回车。F5*.m32MATLAB的含义为_。矩阵实验室（Matrix Laboratory）33.MATLAB中，都包括哪些工具箱，举例不少于6个：_偏微分方程工具箱 、最优化工具箱 、插值运算工具箱、数理统计工具箱 、信号处理工具箱 、通信工具箱，滤波设计工具箱 。34、在MATLAB中，plot命令用于什么， 。 画曲线35、size命令的功能，确定_。 阵列的长度36、a=1 2 2，求a的转置_。 a=1;2;237、在线形图型格式的设置中，字符r表示什么颜色。 红色 38. 在MATLAB中,矩阵的秩由函数 _来实现。 rank39在MATLAB中,单位矩阵由函数_来实现。 eye40. 在MATLAB中,二维饼图由函数_来实现。 pie41. 在线形图型格式的设置中，字符b表示什么颜色 。 蓝色42. MATLAB表达式3*232的结果是_ 。 19243.条形图由函数_来实现。 bar44. 倘若要是对x进行赋值，从5到25，间隔是0.1, x=5:0.1:25;45. 倘若要是对x进行赋值，从2到20，中间间隔100个点, x=linspace(2,20,100);选择题（计算部分）1. 下列变量中的最大数是() Aeps Brealmin Crealmax D-pi2. 下列表示可以作为MATLAB的变量的是() Aabcd- 2 Bxyz_2# Ch DX_1_a3. 如果x=1: 2 : 8,则x(1)和x(4)分别是() A1，8 B1, 7 C2, 8 D2, 74. 下列变量中比0大的最小数是() Ai Brealmin Crealmax Dpi5. 在MATLAB中下列数值的表示不正确的是（） A+99 B1.3e-5 C2-3*e2 D3-2*pi c d b b b6. 下列表示可以作为MATLAB的变量的是() Aabcd-2 Bxyz_2a Ch D3hour7. 如果x=1:-2:-8,则x(1)和x(5)分别是() A1，-8 B-1,-7 C-1,-8 D1,-78. MATLAB表达式2*232的结果是() A128 B4096 C. D256 9. 下列哪个变量的定义是不合法的（）。(A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz 10. 指出下列错误的指令（）。 (A) syms a b; (B) syms a, b; (C) syms(a', b'); (D) syms(a','b','positive'); b d a A B11.下列哪条指令是求矩阵的行列式的值（）。 (A) inv (B) diag (C) det (D) eig inv是求逆矩阵diag 是矩阵对角元素的提取和创建对角阵sqrt(x) 求x的平方根 abs(x) 求x的绝对值det(a) 求行列式的值eig计算矩阵A的特征值和特征向量的函数clf; 用来清除图形的命令12 .清空 Matlab 工作空间内所有变量的指令是（）。(A) clc (B) cls (C) clear (D) clf 13、下列变量名中是合法的。(A) char_1 ; (B) x*y ; (C) xy ; (D) end 14.以下哪个说法是正确的（）。(A)Matlab 进行数值计算的表达精度与其指令窗口中的数值显示精度相同。 (B)Matlab 指令窗口中显示的数值有效位数不能超过 7 位。 (C)输入二维数值数组时，需要用到逗号和分号，它们可以在中文状态下输入。 (D)历史指令窗口所记录的内容 与 diary 指令所产生“日志”内容是不同的。 15. 产生四维单位矩阵的语句为().A.ones(4) B.eye(4) C.zeros(4) D.rand(4) >> eye(4)ans = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1>> ones(4)ans = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1>> zeros(4) ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> rand(4)ans = 0.8147 0.6324 0.9575 0.9572 0.9058 0.0975 0.9649 0.4854 0.1270 0.2785 0.1576 0.8003 0.9134 0.5469 0.9706 0.1419C C a D b16.求解线性方程组Ax=b,当det(A)0时，方程的解是().A. Ab B.b/A C.bA D.A/b17. 运行如下程序后，输入回车，命令窗口（command windows）显示的结果为()c=input('请输入一个字符','s');if c>='A' & c<='Z'disp(setstr(abs(c)+abs('a')-abs('A');elseif c>='a'& c<='z'disp(setstr(abs(c)- abs('a')+abs('A');elseif c>='0'& c<='9'disp(abs(c)-abs('0');elsedisp(c);endA. 9 B. 8 C. 7 D. 618. 运行下面的代码后，输入：回车，回车，回车，则输出结果为（）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);. x1=num2str(1) x2=num2str(1) B. x1=num2str(-1) x2=num2str(1)C. x1= -1 x2=-1D. x1=1 x2=-119. 清除工作空间（wordspace）的命令是（）. clc; B. clear; C. clf; D.delete; 20. 运行如下程序后，输入回车，命令窗口（command windows）显示的结果为( )x=input('请输入x的值:');if x=10y=cos(x+1)+sqrt(x*x+1);elsey=(3(1/2)*sqrt(x+sqrt(x);endyA. 9 B. 8 C. 7 D. 6a a c b d21. 运行如下程序后, 命令窗口（command windows）显示的结果为( b)s=0;a=12,13,14;15,16,17;18,19,20;for k=as=s+k;enddisp(s');A. 144 B. 39 48 57 C.145 D. 45 48 5122. 运行如下程序后, 命令窗口（command windows）显示的结果为(b )k=0;for n=100:200if rem(n,21)=0 %R=rem(X,Y),求余数函数,X,Y应该为正数k=k+1; continueend break;end kA.105 B. 5 C.4 D.10123. 图形窗口（Figure）显示网格的命令是（b）A. axis on B. grid on C. box on D. hold on24. 已知函数文件如下，则factor(4)=(c )function f=factor(n)if n<=1f=1;elsef=factor(n-1)*n; endA. 4 B. 12 C. 24 D.4825. 运行如下程序后, 命令窗口（command windows）显示的结果为( d)A=13,-56,78; 25,63,-735; 78,25,563; 1,0,-1; y=max(max(A) A. y=564 B.y=9 C.y=-735 D.y=56326. 在图形指定位置加标注命令是（c）A. title(x,y,y=sin(x); B. xlabel(x,y,y=sin(x);C. text(x,y,y=sin(x); D. legend(x,y,y=sin(x); %添加图例的标注,27.下列哪个函数为插值函数（b）A. P=polyfit(X,Y,3) B. Y1=interp1(X,Y,X1,'method')C. Y,I=sort(A,dim) D. R=corrcoef(X)28i=2; a=2i; b=2*i; c=2*sqrt(-1); 程序执行后；a, b, c的值分别是多少？c(A)a=4, b=4, c=2.0000i(B)a=4, b=2.0000i, c=2.0000i(C)a=2.0000i, b=4, c=2.0000i (D) a=2.0000i, b=2.0000i, c=2.0000i29. 求解方程x4-4x3+12x-9 = 0 的所有解a(A)1.0000, 3.0000, 1.7321, -1.7321(B)1.0000, 3.0000, 1.7321i, -1.7321i(C)1.0000i, 3.0000i, 1.7321, -1.7321(D)-3.0000i, 3.0000i, 1.7321, -1.732130、在循环结构中跳出循环，但继续下次循环的命令为 。c(A) return; (B) break ; (C) continue ; (D) keyboard31. 用round函数四舍五入对数组2.48 6.39 3.93 8.52取整，结果为 c(A) 2 6 3 8 (B) 2 6 4 8 (C) 2 6 4 9 (D) 3 7 4 932. 已知a=2:2:8, b=2:5，下面的运算表达式中，出错的为 c(A) a' *b (B) a .*b (C) a*b (D) a-b辨析题：1 MATLAB有几种建立矩阵的方法？各有什么优点？。2 MATLAB提供编程语言，请辨析M文件编程中的命令文件与函数文件两种形式。3 辨析matlab中的拟合与插值有什么区别4. 判断下列一段程序的错误，并进行改正。 x=-2:0.2:2;y=x.2;plot(x,y)str1 = num2str(min(x);% 数字转化为字符串str2 = num2str(max(x);% 数字转化为字符串out = Value of f from 'str1' to 'str2'xlabel(out);out = 'Value of f from ' str1 ' to ' str2;或out = 'Value of f from str1 to str2'5根据下面两个矩阵执行下列的矩阵运算命令，并回答有关的问题（1） A.*B和 A*B将分别给出什么结果，它们是否相同，为什么？（2） 得出A.B、A/B及AB的结果，并分别解释它们的物理意义。问答题：1 MATLAB系统包括哪些窗口，各个窗口功能是什么？2 分析下面程序运行后s1,n1,s2,n2,m的值。s1=0;s2=0;n1=0;n2=0；x=1,-4,-8,3,10,-9,7,-3,10, 8,-5,-2,2,0 ;m=length(x);for i=1:mif x(i)<0s1=s1+x(i);n1=n1+1;else s2=s2+x(i); n2=n2+1;endends1,n1,s2,n2,m 请回答s1,n1,s2,n2,m的结果。3 已建立一个简单的SIMULINK模型，如下图，用信号发生器产生一个仿真图，请对下图的符号和功能进行解释。3 A = 1 2; 3 4; B = 1 2; 1 4; 写出下面几条指令的运行结果： C = A*B D = A+B E = A.*B 答：C =3 107 22D =2 44 8E =1 43 164.mysum=0；i=1；while(i<=100)mysum=mysum+i；i=i+1；endmysum答：mysum =50505.在MATLAB窗口中，用哪两种方法启动Simulink？答：在MATLAB命令窗口中输入“simulink”再回车，或者单击工具栏中的SIMULINK图标6. 请修改下面的程序，让他们没有for循环语句！A=1 2 3; 4 5 6; 7 8 9;r c=size(A); for i=1:1:rfor j=1:1:cif (A(i,j)>8 | A(i,j)<2)A(i,j)=0;endendend解：程序如下：a=1 2 3;4 5 6;7 8 9;b=(a<=8&a>=2);a=a.*b我和妹子那些事7 编写一个函数，使其能够产生如下的分段函数：并调用此函数，绘制x=0,+2范围内的f(x)*f(x+2) 。答：function y=f(x)if x<=2y=0.5*x;else if x>6y=0.5;else y =1.5-0.25*x;endend运行结果 x=2f(x)=1x = 0:0.05:2;y = diag(A2(x)'*A2(x+2);plot(x,y);xlabel('bfx');ylabel('bfy');8（1） 将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。,x （曲线图）；f(x,y)=(-2<=x<=2,-4<=x<=4) （曲面图）。答：subplot(2,2,2);ezplot('(cos(x)(1/2)',-pi/2 pi/2);ylabel('y'); subplot(2,2,3);x=-2:0.5:2;y=-4:1:4;ezsurfc('x2/22+y2/42')(2) A 是一个維度m×n的矩阵. 写一段程序, 算出A中有多少个零元素。答：>> A= input ('请输入一个矩阵')m,n= size(A);sig=0;for i=1:mfor j=1:nif A(i,j)=0sig = sig+1;endendend请输入一个矩阵0 1 2;1 0 2; 0 0 0A =0 1 21 0 20 0 0>>sigsig =59 向量 A=. 写一段程序, 找出A中的最小元素。答：A= input ('请输入一个向量')m,n=size(A)min =A(1,n);for i=1:nif A(1,i)<minmin=A(1,i)endend请输入一个向量1 2 3 -5 2 0A =1 2 3 -5 2 0min =-510 写出执行以下代码后C，D，E的值 A=1,2,3;4:6;7:9; C=A;10,11,12,D=C(1:3,2 3)E=C(2,1 2)C = 1 2 3 4 5 6 7 8 9 10 11 12D = 2 3 5 6 8 9E = 4 511 写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值。x=0,1,0,2,0,3,0,4; for k=1:8if x(k)=0x(k)=k;elsex(k)=2*k+1;endenddisp(x);1 5 3 9 5 13 7 1712 创建符号函数并求解，要求写出步骤和运行结果(1)创建符号函数f=ax2+bx+c(2)求f=0的解syms a x b c;f=a*x2+b*x+c;solve(f)ans = 1/2/a*(-b+(b2-4*a*c)(1/2) 1/2/a*(-b-(b2-4*a*c)(1/2)13 求解以下线性方程组，要求写出程序代码和运行结果。 2x1-3x2+x3+2x4=8 x1+3x2+ x4=6 x1-x2+x3+8x4=1 7x1+x2-2x3+2x4=5解：程序代码：a=2 -3 1 2;1 3 0 1;1 -1 1 8;7 1 -2 2;b=8 6 1 5'ra=rank(a);rb=rank(a b);det(a);xx=ab运行结果：xx = 3.1079 1.3278 7.9502 -1.091314. (1) 写出关系运算符中的等号、不等号、小于号、大于号、小于等于号和大于等于号；(2) 写出逻辑操作中逻辑“与”、逻辑“或”及逻辑“非”的符号； (3) 并用语句表达“如果a等于 b且c等于0就让d 等于3，否则如果a大于b且c=1让d等于 0，其余情况让d等于3”；解：（1） = = < > <= >=（2）& | （3）if a=b&c=0d=3else if a>b&c= =1d=0elsed=3end15. 设A= 和B= 是两个10行10列的矩阵（数组），试说明命令A*B, A/B, AB, A.*(B.A), A./B, A.B和sin(A).*exp(B)的涵义。解：A,B的矩阵乘法，矩阵A乘以B的逆，矩阵B乘以A的逆，B的每个元素的对应A次幂乘以A对应的每个元素A的每个元素除以B对应的每个元素，B的每个元素除以A对应的每个元素，A的每个元素的正弦与B对应的每个元素的指数求积，16. 设A为14*17 数组，B为一个行数大于3的数组，请给出（1）删除A的第4、8、12三列的命令； （2）删除B的倒数第3行的命令；解：(1)A(:,4 8 12)=(2) B(end-2)=编程简题1 利用MATLAB数值运算，求解线性方程组(将程序保存为？.m文件)2 求下列联立方程的解 3x+4y-7z-12w=45x-7y+4z+ 2w=-3x +8z- 5w=9-6x+5y-2z+10w=-8求系数矩阵的秩；求出方程组的解。 解：（1）>> a=3 4 -7 -125 -7 4 2 ;1 0 8 -5;-6 5 -2 10;c=4; -3; 9;-8;b=rank(a)b = 4（2）>> d=acd = -1.4841, -0.6816, 0.5337,-1.2429即： x=-1.4841;y= -0.6816;z= 0.5337;w=-1.24293 编写MATALAB程序，完成下列任务（将程序保存为test04.m文件）：（1）在区间 0,4*pi上均匀地取20个点构成向量 ；（2）分别计算函数y1=sin(t) 与y2=2cos(2t) 在向量 t处的函数值；（3）在同一图形窗口绘制曲线y1=sin(t) 与y2=2cos(2t) ，要求 y1曲线为黑色点画线，y2 曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1 and y2”。程序如下：（1）t=linspace(0,4*pi,20);（2）y1=sin(t);y2=2*cos(2*t);（3）plot(t,y1,'k-.');text(1.2,sin(1.2),'y1leftarrow','FontSize',12);hold on;plot(t,y2,'r-o');text(6,1.5,'y2leftarrow','FontSize',12);title('y1 and y2')4 编写M函数文件，利用for循环或while循环完成计算函数y=sum(n)= 的任务，并利用该函数计算 n=20时的和（将总程序保存为？.m文件）。function sum=text4(n);sum=0;m=1;while m<=n; sum=sum+mm m=m+1;end当n=20时 用上述程序求得的和为： 1.0688e+0265 请产生一个100*5的矩阵，矩阵的每一行都是1 2 3 4 5解：repmat(1:5,100,1)6 已知变量：A=ilovematlab；B=matlab, 请找出：（A） B在A中的位置。 解：strfind(A,B):(strfind(A,B)+length(B)-1)（B） 把B放在A后面，形成C=ilovematlabmatlab 解： C = A,B7编写M文件，炼钢基本上是一个氧化脱碳的过程，钢液中原含碳量多少直接影响到冶炼时间的长短，下表是某平炉的熔钢完毕碳(x)与精炼时间(y)的生产记录。现希望从上表的数据中找出x与y变化规律的经验公式，用多项式进行曲线拟合，并给出相应的曲线。>> x=134 150 180 104 190 163 200;>> y=135 168 200 100 215 175 220;>> f=polyfit(x,y,2)f = -0.0016 1.7666 -68.3091>> plot(x,y,'o',x,polyval(f,x)x与y的关系式为y=-0.0016*x2+1.7666*x-68.3091