数学实验matlab练习题

1、2015-2016数学实验练习题一、选择题1 .清除Matlab工作空间（wordspace）变量的命令是（B）A.clcB.clearC.clf2 .清除当前屏幕上显示的所有内容，但不清除工作空间中的数据的命令是（A）A.clcB.clearC.clf3 .用来清除图形的命令（C）A.clcB.clearC.clf4 .在MATLABi序中，使命令行不显示运算结果的符号是（A）A.;B.%C.#D.&5 .在MATLABS序中，可以将某行表示为注释行的符号是（B）A.;B.%C.#D.&6 .在循环结构中跳出循环，执行循环后面代码的命令为（B）A.returnB.breakC

2、.continueD.Keyboard7 .在循环结构中跳出循环，但继续下次循环的命令为（C）A.returnB.breakC.continueD.Keyboard8 .MATLAB中用于声明全局变量的关键字是（C）A.infB.symsC.globalD.function9 .用户可以通过下面哪项获得指令的使用说明（A）A.helpB.loadC.demoD.lookfor10 .在MATLA命令窗口中键入命令S=zoros（3）;可生成一个三行三列的零矩阵,如果省略了变量名S,MATLAB!现计算结果将用下面的哪一变量名做缺省变量名（A）A.ans;B.pi;C.NaN;D.Eps.11

3、.9/0的结果是（B）A.NAN;B.Inf;C.eps;D.012. 在MATLABP程序或语句的执行结果都可以用不同格式显示，将数据结果显示为分数形式，用下面哪一条命令语句（D）A.formatlong;B.formatlonge;C.formatbank;D.fromatrat13. 下列MATLA命令中是构造1行3列的（-1,1）均匀分布随机矩阵的命令的是D)A.randn(1,3)；B.rand(1,3)；C.ones(3)；D.以上都不对14. 产生四维元素都为1矩阵的语句为(A)A.ones(4)B.eye(4)C.zeros(4)D.rand(4)15. 用round函数对数组

4、取整，结果为(C)A.2638B.2648C.2649D.374916. y=dsolve('Dy=1/(1+xA2)-2*yA2','y(0)=0','x');ezplot(y)的功能是(A)A.求微分方程特解并绘图；B.解代数方程；C.求定积分；D.求微分方程通解.17. MATLAB命令roots(1,0,0,-1)的功能是(D)3A.产生向量1,0,0,1；B.求方程x+1=0的根；C.求多项式x3-1的值；D.求方程x3-1=0的根。18 .在MATLA晞令窗口中，键入命令A=1,2,3;4,5,6;7,8,0;A(2,:)*A(:,2

5、);将出现的结果是(B)A.15;B.81;C.36;D.30.19 .在MATLAB命令窗口中键入命令A=142;312;615;则det(A(1:2,2:3).*A(2:3,2:3)的结果是(C)A.-143;B.60;C.36;D.-19程序执行后a, b, c 的值分别是20 i=2;a=2i;b=2*i;c=2*sqrt(-1);(C)A.a=4,b=4,c=B.a=4,b=,c=C.a=,b=4,c=D.a=,b=,c=21.在MATLABT令窗口中键入命令A=14;35;B=23;38;则A*B的结果是(A)21214352121435A.;B.;C.;D.94021492149

6、940A=1 4;3 5; B=2 3;3 8;则 A*B 的结果22 .在MATLABT令窗口中键入命令是（B）A.2 12; B.9 4014352121435;C.;D.2149214994023 .下面有关MATLA变量名和函数名的说法错误的是（D）A.变量名的第一个字符必须是一个英文字母B.变量名可由英文字母、数字和下划线混合组成C.变量名不得包含空格和标点，但可以有下连字符D.变量名和函数名对于英文的大小使用没有区别24 .正确表达命题A或B都大于C的逻辑表达式应该用下面哪一行()A.A>C;B.B>C;C.A>C&B>C;D.A>C|B>

7、;C。25 .输入方阵A的数据，在MATLABP计算A的逆矩阵用下面哪一命令(D)A.det(A)B.eig(A)C.diag(A)D.inv(A)26 .输入方阵A的数据，在MATLABP计算A的特征值用下面哪一命令(D)A.inv(A)B.diag(A)C.det(A)D.eig(A)27 .输入方阵A的数据，在MATLAB计算A的行列式的信用下面哪一命令(C)A.inv(A)B.diag(A)C.det(A)D.eig(A)28 .输入方阵A的数据，在MATLA呻计算A的秩用下面哪一命令(C)A.inv(A)B.diag(A)C.rank(A)D.eig(A)29 .命令A=magic(

8、3)创建3阶魔方矩阵，求A的特征值绝对值最小用(A)A.min(abs(eig(A);B.min(eig(abs(A);C.min(eig(A);D.min(abs(A);30 .下列哪个变量的定义是不合法的(A)A.abcd-3B.xyz_3C.abcdefD.x3yz31 .下列变量名中合法的是(A)A.char_1B.x*yC.xyD.end32 .在MATLA呻下列数值的表示不正确的是(B)A.+99B.C.2-3*eA2D.3-2*pi33 .图形窗口(Figure)显示网格的命令是(B)A.axisonB.gridonC.boxonD.holdon34 .为了使两个plot的图形在

9、同一个坐标显示，可以使用哪个命令进行图形保持(D)A.axisonB.gridonC.boxonD.holdon35 .假设在当前工作目录中有一名为exal的M文件，其中有些错误。现在为了修改它，用键盘命令将这一文件打开编辑，使用下面的命令不对的是()A.editexal；B.edit;C.typeexal；D.openexal36 .如果x=1:-3:-18,则x(3)和x(7)分别是(D)A-6,-14B-5,-14C-8,-17D-5,-1737 .已知a=2:2:8,b=2:5，下面的运算表达式中，出错的为(C)A.a'*bB.a.*bC.a*bD.a-b38 .设A=243;

10、531;367，则sum(A)，length(A)和size(A)的结果(D)A.101311933B.9916333C.9916933D.10131133339 .输入矩阵A=132;37-5;569;使用全下标方式用A(3,2)取出元素“6”，使用单下标方式用哪个命令可以取得(D)A.A(4)B.A(7)C.A(8)D.A(6)540 .MATLA跄令A=rand(5,5);创建A(aj)55,求max|aj|用(B)ji1A.max(sum(abs(A)B.max(sum(abs(A)C.max(sum(A)D.sum(max(A)41 .MATLAB命令x=1,2,4,5,9;mean

11、(x)的计算结果是(B)A.4B.C.D.2142 .运行如下程序后，命令窗口(commandwindows显示的结果为(D)A=13,-56,78;25,63,-735;78,25,563;1,0,-1;y=max(max(A)A.y=564=9=-735=56343 .运行如下程序后，命令窗口(commandwindows显示的结果为(A)A=13,0,-5;8,-6,2;15,3,7;y=max(A)A.y=1537=13815=1587=133744 .命令factor()用于分解因式，symsx;f=4*xA3+9*xA2-30*x;factor(diff(f)的结果是(B)A.(x

12、-1)*(2*x-5)B.6*(x-1)*(2*x+5)C.6*(x+1)*(2*x+5)D.(x+1)*(2*x-5)45 .MATLAB1达式2*2A3A2的结果是(A)A.128B.4096C.262144D.25646 .已知函数文件如下，则factor(4)=(C)functionf=factor(n)ifn<=1f=1;elsef=factor(n-1)*n;endA.4B.12C.2447 .下面的程序执行后array的值为(A)fork=1:10ifk>6break;elsearray(k)=k;endendA.array=1,2,3,4,5,6B.array=1,

13、2,3,4,5,6,7,8,9,10C.array=6D.array=10.48. 运行如下程序后,命令窗口(commandwindows显示的结果为(B)s=0;a=12,13,14;15,16,17;18,19,20;fork=as=s+k;enddisp(s');A.144B.394857D.45485149. 十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令k=mod(2008,12)+1的结果是(D)A.k指向第二动物牛；B.k指向第三动物虎；C.k指向第四动物兔；D.k指向第五动物龙。50. 在图形指定位置加标注命令是(C)A.title(x,y,y=sin(x);B.xlab

14、el(x,y,y=sin(x);C. text(x,y, 加图例的标注,y=sin(x) ); D. legend(x,y,y=sin(x) ); %添51. 下列哪个函数为插值函数(B)A.P=polyfit(X,Y,3)B.Y1=interp1(X,Y,X1,'method')C.Y,I=sort(A,dim)D.R=corrcoef(X)2、 填空题1. 计算多项式乘法(x22x2)(x25x4)命令是_conv(1,2,1,1,5,4)_;计算多项式除法(3x313x26x8)/(x4)的命令是_k,r=deconv(3,13,6,8,1,4);2. 要从键盘读入一个字

15、符串并赋值给变量x,且给出提示“Whoisshe?”，应使用命令_x=input(whoisshe?,s)；3. 在MATLABr令窗口中的“>>”标志为MATLAB命令行提示符，“I”标志为输入_提示符。4. 求矩阵A的对角矩阵函数是_diag，求A的下三角矩阵函数是_tril。5. MATLA呻作三维曲面网格图的绘图函数为mesh(z)。6. 在MATLA瑜令窗口中键入命令A=1,2,3;4,5,6;A(：,1,3)二口；结果是_2;5。7. 在MATLAB命令窗口中键入命令A=13;46;C=A,eye(size(A);A,ones(size(A);结果是_1310;4601

16、;1311;4611_。8. 在MATLA晞令窗口中键入命令A=123;654;reshape(A,3,2);结果是_15;63;24。9. 生成一个从2到20,中间间隔100个点的数组，并将此数组赋值给x的命令是_x=2:200_。10. P,Q分别是个多项式的系数矢量，求P对应的多项式的积分(对应的常数项为K),使用的命令是;求P/Q的解，商和余数分别保存在k和r,使用的命令是k,r=deconv(p,q)_。3、 程序填空1 .对于任意正实数C,取初始近似值x0>0(例如取x0=2),利用迭代公式xn0.5xn1C/xn1，n1,2,可计算数列。现在程序功能为从键盘输入数据C,计算

17、数歹Uxn的前八项。C=input('inputC:=');x0=2;x=;fork=1:8x1=x=x;x1;x0=end2.反正切函数的Tylor级数展开式arctanx( 1)n2n 1-，取x = 1,可得用2n 1于计算的级数一1-143o当级数的通项绝对值小于给定的误差界errors时，则结束计算。functionz=adds(errors)n=1;an=1;f=1;S=whilean>errorsn=n+2;an=1/n;f=-f;S=endz=S;3 .角谷猜想：任给一个不等于1的正整数n,如果它是偶数则将它除以2,如果它是奇数则将它乘以3并加上1。如此重

18、复操作，最后这个数总变为1。程序如n=input('inputn:=')；disp(n)whilen=1ifmod(n,2)=0elsen=endend4 .编写MATLA毓序，在同一图形窗口中绘制y1和y2两条曲线，并对图形坐标轴进行控制。程序如下：x=-3二3;y1=2*x+5;y2=x.A2-3*x;plot(x,y1);%绘制曲线y1;%呆持曲线y1plot(x,y2);%绘制曲线y2ml=max(y1,y2);m2=min(y1,y2);%fflaxis调制坐标轴，横坐标在-33之间，且完全显示图形;物口入图形注释，横坐标为“时间”，纵坐标为“幅度”5 .用Altke

19、n迭代法求方程x33x10在附近的根，要求精确到4位小数，并给出迭代次数。symsxfxgx;gx=(xA3+1)/3;fx=xA3-3*x+1;；ffx=subs(fx,'x',x);while；u=subs(gx,'x，,x);v=subs(gx,'x，,u);x=v-(v-u)A2/(v-2*u+x);k=k+1;end6 .抛物线的方程为yax2cx,输入c后，绘制该曲线，寻找曲线的最小值和对应的x值。程序如下：clc;clear;closeall;c=input('c=');a=1;x=-2:2;%对变量xnum赋为数组x的个数y=ze

20、ros(1,xnum);%按照函数形式对y赋值figure;%绘制曲线xlabel('x');ylabel('y');ymin,Imin=min(y);%找到最小函数值对应的x值四程序阅读理解1. 程序如下，说明该程序的功能，并写出该程序所用的数学函数symsxf=1/(5+2*cos(x)ezplot(f)df=diff(f,x,1)ezplot(df)2. 下面程序的功能是什么？>x=0:pi/50:4*pi;y1=sin(x);y2=cos(x);plot(x,y1,'-r',x,y2,':b')xlabel('

21、;x');ylabel('y');gridon;3. 解释下面程序的功能，并写出该程序所求解的数学问题。symsxyy=dsolve('Dy=x+y','y(0)=0','x')ezplot(y)4. 下面的函数主要完成什么功能？functionf=factor(n)ifn<=1f=1;elsef=factor(n-1)*n;end,S=3,写出下列指令运行结果。5. 给定A=1,2;3,4，B=-1,-2;2,1A.*BA*B，S.*A和S*B6. 写出下列指令运行结果。A(1,1)='thisiscell

22、'A1,2=123;456;A2,1=1+2*i;A2,2=A1,21+(A1,21(1,1)+A1,21(2,2);celldisp(A)7. 写出下列指令运行结果。a=1127;416;14819;b=(a<=10&a>=5);a=a.*b五.程序设计1.limx 0x sin x2.cosx,求 y(10)3.4.xe sin2xdx4x.dx254x25.6.求数项级数I的和。n 1n(n 1)1sinex求7.1n(n1)(n2)8.因式分解：x45x35x25x69.求代数方程组2axbyc0关于x,y的解。010 .计算多项式4x412x314x25x

23、9的零点、微分、积分以及在x=7处的值11 .用符号函数法求解方程at2+b*t+c=0。12 .求方程20x5+16x4+2x3+x+6=0的所有解13 .将(x-6)(x-3)(x-8)展开为系数多项式的形式。14 .计算多项式乘法(x2+2x+2)(x2+5x+4)。15 .计算多项式除法(3x3+13x2+6x+8)/(x+4)。16.设 f(x,y) = 4 sin (x3y)，求17 .求f(x)(-的二阶导数并化简。x118 .求隐函数F(x,y)xy1siny所确定白导数dy2dx19 .求椭圆的参数方程xaC0st所确定白导数dy。ybsintdxcos2x1d2ydx220

24、.求解微分方程y|x0dydx|x21.分别采用符号法，梯形法，抛物线法及矩形法求定积分22.4矩阵a7364,计算a的行列式和逆矩阵923.求矩阵Aaiia21a12的行列式值、逆和特征根。a2224.用三种不同的方法求解线性方程组2x1xiXiX23x22x24X25x3x37x3x46x42x46x489u的解，并计算系50数矩阵的行列式和秩。25.对于AXB,如果A372628,求解X,A的行列式和秩。26.5x1解线性方程组x1x12x3x23x3x2；4x4x3x45x41231并求系数矩阵的行列式。27.画下列函数的图形:(1)yusintucostt_4202x(2)yzsint(3cost(3sinucosu)cosu)28 .符号函数绘图法绘制函数x=sin(3t)cos(t),y=sin(3t)sin(t)的图形，t的变化范围为0,。29 .绘制极坐标下sin(3*t)*cos(t)