matlab相关习题(16页).doc

-1. 已知A=0 9 6;1 3 0;B=1 4 3;1 5 0;写出下列各指令运行的结果。A & B 的运行结果ans= ；A ./ B 的运行结果ans= 。2. 产生4阶全0方阵的命令为 zeros(4) ；产生3阶全1方阵的命令为 ones(3) 。3. A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为 和 2,5,5 。4. 标点符号 分号 可以使命令行不显示运算结果， 百分号 用来表示该行为注释行。5. 用if判断语句判断80 x >60,在MATLAB中if语言后的判断应写为 if x<=80&x>60 。6. P, Q分别是个多项式的系数矢量，求P对应的多项式的积分（对应的常数项为K），使用的命令是 polyint(P,K) ；求P/Q的解，商和余数分别保存在k和r，使用的命令是 k,r=deconv(P,Q) ；7. 为了使两个plot的图形在同一个坐标显示，可以使用 hold on 命令进行图形保持；可以使用 grid on 命令为图形添加网格。8. MATLAB的工作空间中有三个变量v1, v2, v3，写出把它们保存到文件my_data.mat中的指令 save my_data ；写出把my_data.mat文件中的变量读取到MATLAB工作空间内的指令 load my_data 。二、选择题（每空2分，总共20分）1.下列哪个变量的定义是不合法的 A (A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz2. 下列哪条指令是求矩阵的行列式的值 C (A) inv (B) diag (C) det (D) eig3. 在循环结构中跳出循环，执行循环后面代码的命令为 B (A) return (B) break (C) continue (D) keyboard4. 清空Matlab工作空间内所有变量的指令是 C (A) clc (B) cls (C) clear (D) clf5. 用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 96. 已知a=2:2:8, b=2:5，下面的运算表达式中，出错的为 C (A) a'*b (B) a .*b (C) a*b (D) a-b7. 角度，计算其正弦函数的运算为 D (A) SIN（deg2rad(x)） (B) SIN(x) (C) sin(x) (D) sin(deg2rad(x)8. 下面的程序执行后array的值为 ( A )for k=1:10 if k>6 break; elsearray(k) = k; end end(A) array = 1, 2, 3, 4, 5, 6 (B) array = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10(C) array =6 (D) array =10.9i=2; a=2i; b=2*i; c=2*sqrt(-1); 程序执行后；a, b, c的值分别是多少？(A)a=4, b=4, c=2.0000i (C)(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.0000i10. 求解方程x4-4x3+12x-9 = 0 的所有解(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.7321三、写出程序的执行结果或写出给定要求的指令（总共35分）1. 写出执行以下代码后C，D，E的值 (6分)A=1,2,3;4:6;7:9; C=A;10,11,12,D=C(1:3,2 3)E=C(2,1 2)2. 写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 (5分)x=0,1,0,2,0,3,0,4; for k=1:8 if x(k)=0 x(k)=k; else x(k)=2*k+1; endenddisp(x);3. 创建符号函数并求解，要求写出步骤和运行结果（7分）(1)创建符号函数f=ax2+bx+c(2)求f=0的解4. 求解以下线性方程组，要求写出程序代码和运行结果（5分）5绘制函数曲线，要求写出程序代码(12分)(1)在区间0:2均匀的取50个点，构成向量(2)在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求y1曲线为红色点划线，标记点为圆圈；y2为蓝色虚线，标记点为星号四、使用MATLAB语言进行编程（15分）打印出所有的水仙花数。所谓“水仙花数”，是指一个三位数，其各位数字立方之和等于该数本身。一、填空题1、ans = 0 1 1 1 1 0ans = 0 2.2500 2.00001.0000 0.6000 NaN2、 zeros(4)ones(3)3、 2 5 5 4、 ；%5、 if x<= 80 & x>606、 polyint(P,K)k,r= deconv (P,Q)7、 hold ongrid on8、 save my_dataLoad my_data二、选择题1、 A2、 C3、 B4、 C5、 C6、 C7、 D8、 A9、 C10、 A三、写出程序的执行结果或写出给定要求的指令1、 c = 0 + 2.0000iD = 2 3 5 6 8 9E = 4 52、1 5 3 9 5 13 7 173、>> syms a b c x>> s = a*x2 + b*x + c;>> solve(s) ans = 1/2/a*(-b+(b2-4*a*c)(1/2) 1/2/a*(-b-(b2-4*a*c)(1/2)4、>> syms a b c x>> syms x1 x2 x3 x4>> clear>> syms x1 x2 x3 x4>> f = 2*x1 + 3*x2 + x3 + 2*x4 -8;>> g = x1 + 3*x2 + x4 - 6;>> h = x1 - x2 + x3 + 8*x4 - 1;>> i = 7*x1 + x2 - 2*x3 + 2*x4 -5;>> x1,x2,x3,x4 = solve(f,g,h,i) x1 =13/17x2 =320/187x3 =212/187x4 =19/1875、t = linspace(0, 2*pi, 50);plot(t,y1,'r-o' ,t,y2,'b-*');四、使用MATLAB语言进行编程for k = 100:999a = fix(k/100);b = rem(fix(k/10),10);c = rem(k,10);if a.3 + b.3 + c.3 = kfprintf( “%u,tt”, k ) endendMatlab 测试二一、填空题1、MATLAB常用操作界面包括 命令窗口 、工作空间窗口、 命令历史窗口 、当前目录窗口 、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。2、MATLAB中Inf或inf表示 无穷大 、NaN或nan表示 不是一个数 、nargout表示 函数输出变量数目 。3、工作空间浏览器主要用于内存变量的 查阅 、 保持 和 编辑 。4、MATLAB实现将全下标转换为单下标的指令为 sub2ind 、据单下标换算出全下标的指令为 ind2sub 。5、MATLAB中clf用于 清除图形窗口 、clc用于 清除指令窗口中显示内容 、clear用于 清除MATLAB工作空间中保存的变量 。二、简答题（每题5分，共20分） 1、简述MATLAB历史指令窗的主要作用。答：历史指令窗记录着用户在MATLAB指令窗中所输入过的所有指令。历史记录包括：每次开启MATLAB的时间，每次开启MATLAB后在指令窗中运行过的所有指令。应用功能有单行或多行指令的复制和运行、生成M文件等。2、简述MATLAB函数的基本结构。答：典型M函数文件的结构：函数申明行、H1行、在线帮助文本区、编写和修改记录、函数体。3、简述绘制二维图形的一般步骤。绘制二维图形的一般步骤为：曲线数据准备、选定图形窗及子图位置、调用二维曲线绘图指令、设置轴的范围、坐标分格线、图形注释、图形的精细操作。三、阅读程序并回答问题（每题4分，共28分）1、写出下列指令运行结果。A=zeros(2,4); A(:)=1:8; s=2 3 5;A(s)Sa=10 20 30'A(s)=Saans = 2 3 5Sa = 10 20 30A = 1 20 30 7 10 4 6 8 2、写出下列指令运行结果。A=reshape(1:16,2,8) reshape(A,4,4) s=1 3 6 8 9 11 14 16;A(s)=0A = 1 3 5 7 9 11 13 15 2 4 6 8 10 12 14 16ans = 1 5 9 13 2 6 10 14 3 7 11 15 4 8 12 16A = 0 0 5 7 0 0 13 15 2 4 0 0 10 12 0 03、写出下列指令运行结果。A=1,2;3,4; B=-1,-2;2,1;S=3;A.*Bans = -1 -4 6 4A*Bans = 3 0 5 -2S.*Aans = 3 6 9 12S*Bans = -3 -6 6 34、下面的函数主要完成什么功能？ function f=factor(n)if n<=1 f=1;else f=factor(n-1)*n; end利用函数的递归调用求n!5、写出下列指令运行结果。ch=ABc123d4e56Fg9;subch=ch(1:5) revch=ch(end:-1:1) k=find(ch>=a&ch<=z); ch(k)=ch(k)-(a-A); char(ch) subch =ABc12revch =9gF65e4d321cBAans =ABC123D4E56FG9ans = 46、写出下列指令运行结果。A(1,1)='this is cell'A1,2=1 2 3;4 5 6;A2,1=1+2*i;A2,2=A1,21+(A1,21(1,1)+A1,21(2,2);celldisp(A)A1,1 = this is cellA2,1 = 1.0000 + 2.0000iA1,21 = 1 2 3 4 5 6A2,2 = 7 8 9 10 11 127、下面的程序完成功能是什么？t=0:pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b')xlabel(bfit t);ylabel(bfit y);grid on;绘制图形如下：四、编程题（32分）1、用命令来创建GUI并绘制方程y=ax2 +bx+c图形，需要显示绘图结果的坐标系窗口，还能输入a、b、c的值和x的取值范围，同时用命令设置所以对象的属性。2、在同一图上分别用红色实线和绿色虚线绘制y1=sin(x)和y2=cos(x)在区间0，4*pi的曲线，并用星号*标出两条曲线的交点以及建立图例。（5分）x=0:0.01:4*pi;y1=sin(x);y2=cos(x);z=find(abs(y1-y2)<0.007);plot(x,y1,'r.-','Linewidth',2,'Markersize',1);axis(-1.5,15,-1.5,1.5);hold on plot (x,y2,'g.:','Linewidth',2,'Markersize',1);n=size(z);plot(x(z),y1(z),'k*','Markersize',15);3、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)： price<200 没有折扣 200price<500 3%折扣 500price<1000 5%折扣 1000price<2500 8%折扣 2500price<5000 10%折扣5000price 14%折扣输入所售商品的价格，求其实际销售价格。（5分）price=input('please input your price:');a=0;if price<200 a=1;elseif price<500 a=2;elseif price<1000 a=3;elseif price<2500 a=4;elseif price<5000 a=5;else a=6;endswitch a case 0 disp('0 inputs given'); case 1 b=price; case 2 b=0.97*price; case 3 b=0.95*price; case 4 b=0.92*price; case 5 b=0.90*price; case 6 b=0.86*price; otherwise disp('unexpected input!');endfprintf('the actaul price is %f',b);4、二阶微分方程x"+0.2x'+0.4x=0.2u (t), 其中u(t)是单位阶跃函数，试建立系统模型并仿真。第 16 页-