西北工业大学Matlab习题1

1、一、辨析题(共 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 中

2、独立的符号变量 为。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，

3、且给出提示“ Who is she? ”，应使用命令 ；x=input( Wh o is she? , s )_10设 A=错误!未找到引用源。 和 B=错误!未找到引用源。 和 C=错误!未找到 引用源。 均为 m*n 矩阵，且存在于 WorkSpace 中，要产生矩阵 D=错误 !未找，可用命令到引用源。 ，可用命令计算 可用命令错误!未找到引用源。D=(A-C)/B.C det(inv(A *B)11. 在 MATLAB 命令窗口中的“ ”标志为 MATLAB 的 提示符，“标志为 提示符。命令行 输入12. 已知 A=1 2 3 ；4 5 0 ；7 8 9 ；B=1 0 3 ；1 5

4、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 的下三角矩阵函数是 。diagtril14. 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 -21

5、8. 已 知 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)=_

6、.-1 表示为主对角线下一条对角线一下的数全变为 0 提取矩阵上三角矩阵的函数1 2 3；1 4 9； 0 7 2724. A=4,15,-45,10,6;56,0,17,-45,0; find(A=10 & A60, 在 MATLAB中 if 语言后的判断应写为if x60命令进行图形29. 为了使两个 pl ot 的图形在同一个坐标显示， 可以使用 保持；可以使用 命令为图形添加网格。hold on grid on30. MATLAB 的工作空间中有三个变量 my_data.mat中的指令31. MATLAB 提供了两种运行方式，是 matlab 里，回车。save my_date按和把

7、文件 copy 到v1, v2, v3，写出把它们保存到文件34、35、36、37、38. 在 MATLAB 中,矩阵的秩由函数 39在 MATLAB 中 ,单位矩阵由函数40. 在 MATLAB 中,二维饼图由函数41. 在线形图型格式的设置中，字符 b 表示什么颜色42. MATLAB 表达式 3*232 的结果是 。43.条形图由函数 来实现。44. 倘若要是对 x进行赋值，从 5到 25，间隔是 0.1,45. 倘 若 要 是 对 x 进 行 赋 值 ，画曲线阵列的长度a=1;2;2红色rankeyepie蓝色192barx=5:0.1:25;中 间 间 隔 100 个 点 ,1. 下

8、列变量中的最大数是 (C)AepsBrealmin C realmax D -piF5 *.m32MATLAB 的含义为矩阵实验室( Matrix Laboratory )33.MATLAB 中，都包括哪些工具箱，举例不少于 6 个： 偏微分方程工具箱 、最优化工具箱 、插值运算工具箱、 数理统计工具箱 、 信号处理工具箱 、通信工具箱，滤波设计工具箱 在 MATLAB 中， plot 命令用于什么， size 命令的功能，确定 。a=1 2 2，求 a 的转置 在线形图型格式的设置中，字符 r 表示什么颜色。_来实现。_来实现。 来实现。从 2 到 20 ，x=linspace(2,20,1

9、00);选择题 (计算部分)2. 下列表示可以作为 MATLAB 的变量的是 (D)A abcd- 2 B xyz_2# Ch DX_1_a3. 如果 x=1: 2 : 8,则 x(1)和 x(4)分别是 (B)A1，8 B1, 7 C2, 8D2, 74. 下列变量中比 0 大的最小数是 (B)A i B realmin C realmax D pi5. 在 MATLAB 中下列数值的表示不正确的是( B)A +99 B1.3e-5 C2-3*e2 D3-2*pi c d b b b6. 下列表示可以作为 MATLAB 的变量的是 ( B)Aabcd-2 Bxyz_2a C h D3hour

10、7. 如果 x=1:-2:-8,则 x(1)和 x(5)分别是 (D)A 1， -8 B -1,-7C -1,-8D 1,-78. MATLAB 表达式 2*232 的结果是 (A)A128 B 4096 C. 262144 D2569. 下列哪个变量的定义是不合法的( A )(C) abcdef (D) x3yz(Db) )s; yms( a,b,positive); b d a A B(A) abcd-3 (B) xyz_310. 指出下列错误的指令( B)。(A) syms a b; (B) syms a, b; (C) syms( a,11. 下列哪条指令是求矩阵的行列式的值( C)。

11、(A) inv(B) diag (C) det (D) eiginv 是求逆矩阵 diag 是矩阵对角元素的提取和创建对角阵 sqrt(x) 求 x 的平方根 abs(x) 求 x 的绝对值 det (a) 求行列式的值 eig 计算矩阵 A的特征值和特征向量的函数 clf; 用来清除图形的命令12 .清空 Matlab 工作空间内所有变量的指令是( C)。(A) clc (B) cls (C) clear (D) clf13、下列变量名中 A是合法的。(A) char_1 ; (B) x*y ; (C) xy ; (D) end14. 以下哪个说法是正确的()。(A)Matlab 进行数值计

12、算的表达精度与其指令窗口中的数值显示精度相同(B) Matlab 指令窗口中显示的数值有效位数不能超过 7 位。(C) 输入二维数值数组时，需要用到逗号和分号，它们可以在中文状态下输入(D) 历史指令窗口所记录的内容 与 diary 指令所产生 “日志 ”内容是不同的15. 产生四维单位矩阵的语句为 ().A. ones(4) B.eye(4) C.zeros(4) D.rand(4) eye(4)ans = ones(4)1000ans =0100111100101111000111111111 zeros(4) rand(4)ans =ans =00000.81470.63240.9575

13、0.957200000.90580.09750.96490.485400000.12700.27850.15760.800300000.91340.54690.97060.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=a& c=0& c=9 disp(abs(c)-abs(0);else disp(c);endA. 9 B. 8C. 7

14、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. cle

15、ar; C. clf; D.delete;20. 运行如下程序后，输入回车，命令窗口( command windows )显示的结 果为 ( ) x=input( 请输入 x 的值 :);if x=10 y=cos(x+1)+sqrt(x*x+1); else y=(3(1/2)*sqrt(x+sqrt(x);end y A. 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;end disp(s); A.

16、144B. 39 48 57C.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; continue end break; end kD.101A.105B. 5C.423. 图形窗口( Figure)显示网格的命令是( b)A. axis on B. grid onC. box onD. hold on24. 已知函数文件如下，则 factor(4)=(c ) function f

17、=factor(n) if n=1 f=1;else f=factor(n-1)*n; endA. 4B. 12C. 24D.4825. 运行如下程序后 , 命令窗口( command windows)显示的结果为 ( d)A=13,-56,78; 25,63,-735; 78,25,563; 1,0,-1; y=max(max(A)A. y=564B.y=9 C.y=-735 D.y=56326. 在图形指定位置加标注命令是( c )A. title(x,y, y=sin(x) B. x)l;a bel(x,y, y=sin(x) );C. text(x,y, y=sin(x) D.) l;

18、e gend(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,

19、 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

20、 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);% 数字转

21、化为字符串 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 fromstr1 to str2 ;5 根据下面两个矩阵41220123A 124578 ,B4562078136789执行下列的矩阵运算命令，并回答有关的问题（1） A.*B 和 A*B 将分别给出什么结果，它们是否相同，为什么？（2） 得出 A.B、A/B 及 AB 的结果，并分别解释它们的物理意义问答题：1

22、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;elses2=s2+x(i);n2=n2+1;endends1,n1,s2,n2,m请回答 s1,n1,s2,n2,m 的结果。3 已建立一个简单的 SIMULINK 模型，如下图，用信号发生器产生一个仿真图，请对下图 的符号和功能进行解释。3 A = 1 2; 3 4

23、; B = 1 2; 1 4; 写出下面几条指令的运行结果： C = A*BD = A+BE = A.*B答： C =3 107 22 D =2 44 8E =1 43 164. mysum=0； i=1 ； while(i8 | A(i,j)2)A(i,j)=0;endendend解：程序如下： a=1 2 3;4 5 6;7 8 9; b=(a=2);a=a.*b 7 编写一个函数，使其能够产生如下的分段函数： 错误! 未找到引用源 并调用此函数，绘制 x=0,+2 范围内的 f(x)*f(x+2) 。答： function y=f(x) if x6y=0.5;else y =1.5-0.

24、25*x;endend 运行结果 x=2 f(x)=1 x = 0:0.05:2;y = diag(A2(x)*A2(x+2); plot(x,y);xlabel(bfx);ylabel(bfy);0.50.450.40.350.3y 0.250.20.150.10.051.6 1.8 200 0.2 0.4 0.6 0.8 1 1.2 1.48（ 1） 将一个屏幕分 4 幅，选择合适的步长在右上幅与左下幅绘制出下列曲线图）； f（x,y）=函数的图形。错误!未找到引用源。 ,x 错误!未找到引用源错误! 未找到引用源(-2=x=2,-4=x A= input (请输入一个矩阵 )0(cos(

25、x)1/2y10.5-1 0 1xx 2/2 2 +y 2/42y-5-5xm,n= size(A);sig=0;for i=1:m for j=1:n if A(i,j)=0 sig = sig+1; end end end请输入一个矩阵 0 1 2;1 0 2; 0 0 0 A =0 1 21 0 20 0 0sigsig =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

26、 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 34 5 67 8 910 11 12D =2 35 68 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

27、) 创建符号函数 f=ax 2+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=8x1+3x2+ x4=6x1-x2+x3+8x4=17x1+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运行结

28、果：xx =3.10791.32787.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 ab&c= =1d=0elsed=3end15. 设 A=错误!未找到引用源。 和B=错误!未找到引用源。 是两个 10行 10 列的矩

29、阵(数组)，试说明命令 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-