中南大学matlab习题复习课0401222
中南大学matlab习题复习课
2020/3/3
11
2、利用非线性方程组的符号法求解如下方程组：
x1 0.7sin x1 0.2cos x2 0
x2 0.7 cos x1 0.2sin x2 0
在（0.5，0.5）附近的解。
（a）
a) x1=0.5265 x2=0.5079
b) x1=0.5079 x2=0.5265
ans =
x0=[0.5 0.5] >> fsolve('fc',x0)
0.5265 0.5079
2020/3/3
10
九、方程符号求解
1 求微分方程 y(0)=0的解
x
d2 y dx2
3
dy dx
x2
，y(1)=0，
（c）
a)x^4/5 - x^3/4 b)x^3/3 - x^2/3 c) x^4/3 - x^3/3 d) x^4/5- x^3/4
>> [x,y]=dsolve('Dx=y,Dy=x','t')
2020/3/3
9
非线性方程的符号注解
语法： fsolve(‘fun’, x0) （1）fun为所要求解的函数名，常以M文件形式给出； （2）X0为初始向量或矩阵
Fc.m function y=fc(x) y(1)=x(1)-0.7*sin(x(1))0.2*cos(x(2)); y(2)=x(2)0.7*cos(x(1))+0.2*sin(x(2)); y=[y(1) y(2)]
2020/3/3
8
【例
】 求微分方程
x
d2 y dx2
3
Matlab复习题
M a t l a b复习题(总16页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--一、单项选择题1. x=[1 2 3 4];y=[1 1 1 1];z=x*y，则z等于（）A.1234123412341234⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦B.出错息C.[]1234 D.[]11112. 有一组实验数据x，y，从理论上讲它们应是线性关系，正确的拟合命令应是（）A. p=polyfit(x,y,n)B. p=polyval(x,y,1)C. polyval(x,y,1)D. p=polyfit(x,y,1)3. 用户可以在MATLAB命令窗口用命令编辑功能键对输入的命令进行编辑，编辑键Backspace的功能是（）A.删除光标左侧字符B.删除光标右侧字符C.删除当前行D.光标左移4. format long; pi/4的运行结果是（）A. ans＝B. ans＝C. ans＝D. ans＝5. 已知A=[1 2 3;4 5 6],则命令size(A)的运行结果是（）= = = =3 2 1 6 2 3 66. MATLAB语言可以实现数学上的许多解析功能，正是由于它具备了（）A. 矩阵运算功能B.数值运算功能C.符号运算功能D.数组运算功能7. 计算微分方程221dyyds⎛⎫+=⎪⎝⎭的正确命令是（）A. y=dsolve(‘D2y+y^2=1’,’s’)B. y=dsolve(‘(Dy)^2+y^2=1’,’s’)C. y=dsolve(‘(Dy)^2+y^2=1’)D. y=dsolve(‘D2y+y^2=1’)8. 下列单窗口多曲线子图绘制中的子图分割命令使用正确的时（）A. subplot(2,3,0)B. subplot(1,2,4)C. subplot(225)D. subplot(236)9.三维曲面绘图函数为（）A. surfB. meshC. plot3D. plot使用三维向量[R G B]来表示一种颜色，则黑色为（）A. [1 0 1]B. [1 1 1]C. [0 0 1]D. [0 0 0]11.已知a=[1 2 1;2 2 1;2 1 2];b=[1;2;3];c=a\b,则c等于（）A. []1.00000.33330.6667- B. 错误C.121111/22/31/32/3⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D.1.00000.33330.6667⎡⎤⎢⎥-⎢⎥⎢⎥⎣⎦12.已知a=[1 2 3;4 5 6;7 8 9];b=trace(a)，则b 等于（ ）13.已知a=[1 2 3];c=a^2，则c 等于（ ）A. [1,32,729]B.[1 4 9]C.错误D.[4 10 8]14.已知A=[0 1;-2 -3],则计算矩阵指数函数AT e ，T ＝的命令是（ ）*A) *A) (A) (A)15.已知线性代数方程组1231101011210130014x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦，求解命令为（ ） A. A=[1 1 0;0 1 1;1 0 1;0 0 1];b=[1;2;3;4];X=A/bB. A=[1 1 0;0 1 1;1 0 1;0 0 1];b=[1;2;3;4];X=inv(A)*bC. A=[1 1 0;0 1 1;1 0 1;0 0 1];b=[1;2;3;4];X=inv(A ’*A)*A ’*bD. A=[1 1 0;0 1 1;1 0 1;0 0 1];b=[1;2;3;4];X=pinv(A)*b16.已知()1sin 33x y t e x -=，用数值法计算()40Y y t dt =⎰的正确语句是（ ） A. Y=quad('1./3*exp(-x).*sin(3*x)', '0', '4')B. f=inline('1./3*exp(-x).*sin(3*x)');Y=quad(f,0,4)C. Y=quad(1./3*exp(-x).*sin(3*x),0,4)D. f=inline('1./3*exp(-x).*sin(3*x)');Y=quad('f',0,4)17.能够启动SIMULINK 仿真环境子窗口，展示出SIMULINK 的功能模块组的命令是（ ）A. simulink3B. simulinkC. simulink1D. simulink2为一关于x 的多项式，要求当x=5时多项式的值，正确的MATLAB 命令是（ ）A. polyval(p,5)B. polyfit(p,5)C. value(p,5)D. sum(p,5)21.有一矩阵abc=[1 2 3;2 3 4;3 4 5]，命令mean(abc)后的结果是（ ）A. 提示有错B. 3C. []234D. 234⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦22.在MATLAB 的若干通用操作指令中，清除内存中的变量的是（ ）A. clcB. claC. clfD. clear23.用户可以在MATLAB 命令窗口用命令编辑功能键对输入的命令进行编辑，编辑键Delete 的功能是（ ）A. 删除当前行B. 删除光标左侧字符C. 删除光标右侧字符D. 删除当前命令short e; pi/4 的运行结果是（ ）A. ans ＝B. ans ＝C. ans ＝D. ans ＝25.下面哪个变量是正无穷大变量（ ）A. InfB. NaNC. realmaxD. realmin26.下面命令中能正确定义一个符号变量的是（ ）A. syms w tB. sym w tF=sin(w*t); F=sin(w*t);C. syms(‘F=sin(w*t)’)D. sym(F=sin(w*t))27.计算微分方程2220d y a y dt+=的正确命令是（ ） A. y=dsolve(‘(Dy)^2=－a^2*y’,’s’)B. y=dsolve(‘D2y=－a^2*y’)C. y=dsolve(‘D2y=－a^2*y’,’s’)D. y=dsolve(‘(Dy)^2=－a^2*y’)28.已知t=0:pi/50:10pi,则作三维柱面螺旋线的正确命令是（ ）A. plot(sin(t),cos(t),t)B. plot2(sin(t),cos(t),t)C. plot3(sin(t),cos(t),t)D. plot4(sin(t),cos(t),t)使用三维向量[R G B]来表示一种颜色，则白色为（ ）A. [0 0 0]B. [1 1 1]C. [0 0 1]D. [1 0 1]30.已知a=[1 2 1;2 2 1;2 1 2];b=[3 3 1;3 2 1;1 1 3];c=a+b,则c 等于（ ）A.452542325⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B.[]452542325C. []523245254D.错误31.已知a=[1 2 3;4 5 6;7 8 9];b=fliplr(a),则b 等于（ ）（左右翻转）A. 369258147⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B. 321654987⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C. 789456123⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D.错误32.已知a=[1 2;3 4];b=[5 6;7 8];c=a.*b ，则c 等于（ ）A.5122132⎡⎤⎢⎥⎣⎦ B.19224350⎡⎤⎢⎥⎣⎦C. 错误D.[]5122132 33.已知a=1234⎡⎤⎢⎥⎣⎦;c=a^2，则c 等于（ ） A. 1468⎡⎤⎢⎥⎣⎦B.14916⎡⎤⎢⎥⎣⎦C. 错误D. 7101522⎡⎤⎢⎥⎣⎦ 35. 已知a=[7 3 9 1 0 8];b=sort(a)，则b 等于（ ）A. [0 1 3 7 8 9]B. [9 8 7 3 1 0]C.D. 错误36. 已知A=[3 2 5;7 6 8;9 0 1];b=max(A)，则b 等于（ ）A. [9 6 8]B. [3 0 1]C. 589⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D. 260⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦37.已知A=[0 2 3;1 3 0];B=[1 0 3;1 5 0];C=A.+B ，则C 等于（ ）A.126280⎡⎤⎢⎥⎣⎦B.001101⎡⎤⎢⎥⎣⎦C.错误D.[]12628038. 已知123456x ⎡⎤=⎢⎥⎣⎦，length(x)等于（ ） A. 2 3 B. 3 C. 2 D. 640.有实验数据：t 为时间，L 为函数值序列，用MATLAB 画出L 随t 的变化关系，实验点用“*”表示，正确的命令是（ ）A. plot(t,L,”*”)B. plot(L,t,’*’)C. subplot(t,L,’*’)D. plot(t,L,” *”)41.在图线中标注“图例”的命令是（ ）A. title()B. legend()C. gtext()D. label()42.用户可以在MATLAB 命令窗口用命令编辑功能键对输入的命令进行编辑，编辑键Esc 的功能是（ ）A.删除当前行B.删除光标右面的字符C.删除光标左面的字符D.取消当前命令44.下面哪个变量是非数变量（ ）A. realmaxB. InfC. NaND. realmin46.已知t=0:2*pi/90:2*pi; y=t.*exp(i*t),则绘制复数矢量图的命令是（ ） (y,t) (t,y); (t+i*y) (y)使用三维向量[R G B]来表示一种颜色，则[0 0 1]为（ ）A. 红色B. 蓝色C. 绿色D. 黄色48.已知A=[0 2 3;1 3 0];B=[1 0 3;1 5 0];C=A+B ，则C 等于（ ）A. 126280⎡⎤⎢⎥⎣⎦B. 001101⎡⎤⎢⎥⎣⎦C. 错误D. []12628049.已知a=[1 2 3;4 5 6;7 8 9];b=flipud(a),则b 等于（ ）A. 321654987⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B. 369258147⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C. 789456123⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D. 错误50.已知A=[1 2 3;4 5 0;7 8 9];B=[1 0 3;1 5 0;0 1 2];C=A.*B ，则C 等于（ ）A.22651007911⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B.10942500818⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C.313992512154939⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D.[]22651000911语言中求反正弦的函数是（ ）A. sin()B. asin()C. arcsin()D. asinh()52.已知线性代数方程组1231282313x x x ⎡⎤⎡⎤⎡⎤⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎣⎦，求解命令为（ ）A. A=[1 2;2 3];b=[8;13];X=A/bB. A=[1 2;2 3];b=[8;13];X=inv(A)*bC. A=[1 2;2 3];b=[8;13];X=inv(A’*A)*A’*bD. A=[1 2;2 3];b=[8;13];X=pinv(A)*b53.已知a=[7 3 9 1 0 8];b=sum(a)，则b 等于（ ）A. 28B. 0C. [7 21 189 189 0 0]D. [7 10 19 20 20 28]54.已知x=1:5;y=x.^2;dy=diff(y)，则dy 等于（ ）A. 错误B. [2 2 2]C. 0D. [3 5 7 9]55.已知A=[1,2,6;4,5,10;7,8,9];B=A(1:2,:),则B 等于（ ）A. 610⎡⎤⎢⎥⎣⎦B. 1641079⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C. 1264510⎡⎤⎢⎥⎣⎦D. 1264510789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦ 为一组实验数据，mean(x)的作用是（ ）A.求各列最大值B.求各列数据之和C.求各列中间值D.求各列平均值60.用户可以在MATLAB 命令窗口用命令编辑功能键对输入的命令进行编辑，编辑键Home的功能是（ ）A.光标回首行B.光标回行首C.光标至行尾D.光标至尾行61.创建一个一维行矩阵向量的命令是（ ）=[a b c] =[1;2;3] =[a;b;c] = [1 2 3]62.已知Z ＝zeros （1，2），则Z 等于（ ）A.[]00B.00⎡⎤⎢⎥⎣⎦C.[]01D.10⎡⎤⎢⎥⎣⎦ 64.已知x ＝2＋3i ，则real(x)等于（ ）为色图设定命令，则color([1 0 0])表示设定颜色为（ ）A.红色B.绿色C.蓝色D.黑色语言除了具备基本绘图功能之外，还具有很多其它绘图函数，下面哪个是阶 梯图函数（ ）使用三维向量[R G B]来表示一种颜色，则[0 1 0]为（ ）A.红色B.蓝色C.绿色D.黄色68.已知A=[1 2 3;4 5 0;7 8 9];B=[1 0 3;1 5 0;0 1 2];C=A+B ，则C 等于（ ）A.22651007911⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B.10942500818⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C.313992512154939⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D.[]22651000911 69.已知a=[1 2 3;4 5 6;7 8 9];b=rot90(a),则b 等于（ ）A.321654987⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B.369258147⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C.789456123⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D.错误70.已知a=[1 2;3 4];b=[5 6;7 8];c=a*b ，则c 等于（ ）A.5122132⎡⎤⎢⎥⎣⎦B.19224350⎡⎤⎢⎥⎣⎦C.错误D.[]5122132 语言中求反余弦的函数是（ ）() () () ()73.已知A=[3 2 5;7 6 8;9 0 1];b=mean(A)，则b 等于（ ）A.[ ]B.[7 2 5]C. 3.33337.00003.3333⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D.371⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦ 74.已知x=1:5;y=x.^2;dy=diff(y ，2)，则dy 等于（ ）A.[3 5 7 9]B.[2 2 2] D.错误75.已知A=[1,2,6;4,5,10;7,8,9];B=A(1:2,3),则B 等于（ ）A.610⎡⎤⎢⎥⎣⎦B.1641079⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C. 1264510⎡⎤⎢⎥⎣⎦D. 1264510789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦ 为一组实验数据，median(x)的作用是（ ）A.求各列最大值B.求各列数据之和C.求各列中间值D.求各列平均值为一组实验数据，max(x)的作用是（ ）A ．求各列数据之和B ．求各列最大值C ．求各列中间值D ．求各列平均值82.下列哪个命令是创建一个矩阵向量的正确命令（ ）=’1 2;3 4’ =(1 2;3 4)={1 2;3 4} = [1,2;3,4]83.已知O ＝ones(2,1),则O 等于（ ）A.[]11B.11⎡⎤⎢⎥⎣⎦C.[]01D.10⎡⎤⎢⎥⎣⎦ 84.计算二重不定积分xy xe dxdy -⎰⎰的正确命令是（ ）A. int(int('x*exp(-x*y)','x'),'y')B. int(int(sym(x*exp(-x*y)),'x'),'y')C. int(int(sym(‘x*e xp(-x*y)’),x),y)D. int(int('x*exp(-x*y)',x),y)85.已知x ＝2＋3i ，则imag(x)等于（ ）88. 已知a=[1 2 1;2 2 1;2 1 2];b=[3 3 1;3 2 1;1 1 3];c=a-b,则c等于（）A.210100101--⎡⎤⎢⎥-⎢⎥⎢⎥-⎣⎦B.[]210100101----C. []523245254 D.错误89.已知a=[1 2 3];b=[4 5 6];c=a.\b，则c等于（）A. []579 B.[]0.250.40.5 C.错误 D. []4 2.52语言中求反正切的函数是（）() () () ()91.已知A=[1 2;2 3];P=poly(A)，则P等于（）A.[]141-- B. []282-- C. []141-- D. []141-语言中三次样条插值函数为（）95.已知A=[1,2,6;4,5,10;7,8,9]; A(:,2)=[ ], 则A等于（）A.1641079⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B.1264510⎡⎤⎢⎥⎣⎦C. 610⎡⎤⎢⎥⎣⎦D.1264510789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦100.详细显示工作区变量的命令是（）A. whoB. whosC. whichD. whose 101.下列哪个命令是创建一个矩阵向量的正确命令（）=[1,2;3,4＋5i] =[1 2;3 4+5a]=[1 2a;3 4+5j] =[1 2;3 4+b]102.已知E＝eye（2），则E等于（）A. []01 B.0110⎡⎤⎢⎥⎣⎦C. 1001⎡⎤⎢⎥⎣⎦D.1⎡⎤⎢⎥⎣⎦104.已知x=0:2*pi/90:2*pi; y1=sin(x);y2=cos(x),在同一图上画出y1和y2的命令是（）A. plot(y1,y2)B. plot(x,[y1,y2])C. plot(x,y1);hold on;plot(x,y2)D. plot(x,y1);plot(x,y2)106.已知函数sin(2x)的函数值取值范围是－1～1,自变量x的取值范围是0～6，则绘制该函数曲线的命令是（）A. plot(‘sin(2*x)’,0,6,-1,1)B. fplot(‘sin(2*x)’,0,6,-1,1)C. plot(‘sin(2*x)’,[0,6,-1,1],’r’)D. fplot(‘sin(2*x)’,[0,6,-1,1],’r’)语言中表示暖色色图的色图函数为（）A. hotB. coolC. hsvD. gray语言中表示开平方的函数是（）A. exp()B. sqrt()C. abs()D. log10()111.已知a=[7 3 9 1 0 8];b=max(a)，则b等于（）112.已知x,y为样本向量，xi为已知自变量值，则正确的一维线性插值是（）A. yi=interp1(x’,y’,xi) B. yi=table1(x’,y’,xi)C. yi=spline(x’,y’,xi)D. yi=interp1([x’,y’],xi)114. 在MATLAB的若干通用操作指令中，清除图形窗口的是（）A. clearB. claC. clfD. clc取值从-2π~2π，用MATLAB画sin2x的函数图像正确的命令是（）A．x=-2*π::2*π;plot(x,sin2(x))B．x=-2*pi::2*pi;plot(x,sin(x).^2)C．x=-2*pi:2*pi;plot(x,sinx^2)D．x=-2*::2*;plot(x,(sinx)^2)117.有一矩阵abc=[1 2 3 ;2 3 4 ;3 4 5]，命令mean(abc)后的结果是（）A. 3B. []234 C.提示有错118.用户可以在MATLAB命令窗口用命令编辑功能键对输入的命令进行编辑，编辑键←的功能是（）A. 光标回行首B. 光标右移C. 光标左移D. 光标回行尾119.显示工作区所有变量名称的命令是（）A. whoB. whosC. whichD. whose120.下面把符号表达式赋值给变量f1的正确命令是（）=”exp(-2*x)*sin(x/5)” =exp(-2*x)*sin(x/5)=[exp(-2*x)*sin(x/5)] D. f1=’exp(-2*x)*sin(x/5)’121.初等矩阵函数rand(m,n)所创建的m×n阶随机矩阵的元素的取值范围是（）A. 0～+∞B. -1～1 ～0 D. 0～1中拉普拉斯变换函数是（）123.用plot函数绘制多条曲线时，曲线颜色默认变化值顺序为（）A.红、绿、篮、亮篮、粉红、黄B.篮、红、绿、亮篮、粉红、黄C.篮、绿、红、亮篮、粉红、黄D.红、篮、绿、亮篮、粉红、黄124.绘制三维球面坐标值计算的函数为（）A. grid C. sphere125.在编写绘图程序时，添加图例的函数是（）A. titleB. legendC. xlableD. ylable语言中表示线性灰度色图的色图函数为（）A. hotB. coolC. hsvD. gray127. 已知A=[1 2 3;4 5 0;7 8 9];B=[1 0 3;1 5 0;0 1 2];C=A*B，则C等于（）A.22651007911⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦B.10942500818⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦C.313992512154939⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦D.[]22651000911 128. 已知a=[1 0 0;0 0 1;0 1 0];b=rank(a)，则b 等于（ ）129. 已知a=[1 2 3];b=[4 5 6];c=a.^b ，则c 等于（ ）A.[]132729B.[]149C.错误D.[]41018130. 下列函数中能求矩阵指数的函数是（ ）() () () ()131. x 为一组实验数据，max(x)的作用是（ ）A ．求各列最大值B ．求各列数据之和C ．求各列中间值D ．求各列平均值132. 已知x,y 为样本向量，xi 为已知自变量值，则正确的三次样条插值是（ ） =interp1(x’,y’,xi) =table1(x’,y’,xi)=spline(x,y,xi) =interp1([x’,y’],xi)133.非线性方程组的解析解计算很多情况下无法得到，但用MATLAB 的优化工具箱提供的函数可得到数值解，该函数是（ ）() D.无134. 在MATLAB 中可以用矩阵[C1 C2 C3]来表示颜色，其中C2表示（ ）A.红色B.绿色C.蓝色D.黄色135.已知控制系统的开环传递函数()()215210o G s s s s =++，则绘制根轨迹图的命令是（ ）A. nyquist([15],[1 2 10 0])B. rlocus([15],[1 2 10 0])C. nyquist([15],[1 2 10])D. rlocus([15],[1 2 10])取值从-2π~2π，用MATLAB 画sinx 的函数图象正确的命令是（ ）A. x=-2*pi::2*pi ；plot(x,sin(x));B. x=-2*π::2*π；plot(x,sin(x))C. x=-2*pi:2*pi ；plot(x,sinx)D. x=-2*::2*；plot(x,sinx)137.在MATLAB 的若干通用操作指令中，清除命令窗口的所有显示内容的是（ ）A. clearB. claC. clfD. clc138.用户可以在MATLAB 命令窗口用命令编辑功能键对输入的命令进行编辑，编辑键→的功能是（ ）A.光标右移B.光标左移C.光标回行首D.光标回行尾short; pi/4 的命令结果是（ ）A. ans ＝B. ans ＝C. ans ＝D. ans ＝140.命令m1=[ ； ]；的屏幕显示结果是（ ）= = C. 无 D.错误141.已知x=0:2*pi/90:2*pi; y1=sin(x);y2=cos(x),分别在两张图上画出y1和y2的命令是（ ）A. plot(x,y1);figure(2);plot(x,y2)B. plot(y1,y2)C. plot(x,y1);hold on;plot(x,y2)D. plot(x,y1);plot(x,y2)142.在编写绘图程序时，可以用鼠标在图上添加文字的函数是（ ）143.视角函数的应用格式为view(a,e),其中a 表示方位角，e 表示（ ）A.方位角B.俯视角C.仰视角D.当前视角144.已知a=[1 2;2 3];b=det(a)，则b 等于（ ）A.－1145.已知a=[1 2 3];c=a.^2，则c 等于（ ）A.[]132729B.[]149C.错误D.[]41018146.已知线性代数方程组123110*********x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦，求解命令为（ ） A. A=[1 1 0;1 0 1;0 1 1];b=[1;1;1];X=A/bB. A=[1 1 0;1 0 1;0 1 1];b=[1;1;1];X=inv(A)*bC. A=[1 1 0;1 0 1;0 1 1];b=[1;1;1];X=inv(A’*A)*A’*bD. A=[1 1 0;1 0 1;0 1 1];b=[1;1;1];X=pinv(A)*b147.已知a=[7 3 9 1 0 8];b=mean(a)，则b 等于（ ）为一组实验数据，min(x)的作用是（ ）A. 求各列最大值B. 求各列最小值C. 求各列中间D. 求各列平均值 149.可以对有理多项式进行部分分式化的函数是（ ）A. fsolve C. residue150.在MATLAB 中可以用矩阵[C1 C2 C3]来表示颜色，其中C1表示( )A.红色B.绿色C.蓝色D.黄色151.已知控制系统的开环传递函数()()215210o G s s s s =++，则绘制奈奎斯特图的命令是（ ）A. rlocus([15],[1 2 10])B. rlocus([15],[1 2 10 0])C. nyquist([15],[1 2 10])D. nyquist([15],[1 2 10 0])152.下列随机创建一个3阶2输入2输出系统的状态空间模型的语句是（ ）A. [a,b,c,d]=ord(2,2,3)B. [a,b,c,d]=rmodel(2,2,3)C. [a,b,c,d]=ord(3,2,2)D. [a,b,c,d]=rmodel(3,2,2)二、填空题文件分为独立m文件和两种。
Matlab编程与应用习题和一些参考答案
Matlab编程与应用习题和一些参考答案Matlab 上机实验一、二3.求下列联立方程的解⎪⎪⎩⎪⎪⎨⎧=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\b4．设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。
>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];>> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];>> C1=A*B'>> C2=A'*B>> C3=A.*B>> inv(C1)>> inv(C2)>> inv(C3)5．设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。
>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。
并求该矩阵全体数的平均值和均方差。
（mean var ）a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)10.利用帮助查找limit 函数的用法，并自己编写，验证几个函数极限的例子。
实验一 MATLAB软件应用复习含具体答案
，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证
。 2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品 的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该9种 商品的总收入和总利润。 表1.1
货号
1
2
3. 用两种方法在同一个坐标下作出 这四条曲线的图形，并要求用两种方法在图上加各种标注。 4．用subplot分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题， 1）概率曲线 ； 2）四叶玫瑰线 ； 3）叶形线
4）曳物线
。 5．作出下列曲面的3维图形， 1） ； 2）环面
（2） 新建一个M文件输入命令如下： v=0:0.01:2*pi; u=v; [U,V]=meshgrid(u,v); X=(1+cos(U)).*cos(V); Y=(1+cos(U)).*sin(V); Z=sin(U); mesh(X,Y,Z) 保存并运行，运行结果如下：
T6 新建一个M文件输入命令如下： function shuixianhuashu for m=100:999; a=fix(m/100); b=rem(fix(m/10),10); c=rem(m,10); if m==a.^3+b.^3+c.^3 disp(m) end end 保存，再在命令窗口输入shuixianhuashu结果如下： >> shuixianhuashu 153
探究实验 自由发挥：自己提出问题，实验探索，广泛联想，发现规律，大胆猜想。比如函数cos(1/x)在 x=0附近的振荡现象，有无规律可寻？
T1
新建一个M文件输入命令如下： E=eye(3,3); R=rand(3,2); O=zeros(2,3); S=diag([3,7]); A=[E,R;O,S]; AA=A*A; B=[E,R+R*S;O,S*S]; if B==AA disp('命题成立'); else disp('命题不成立'); end 保存并运行，结果如下： 命题成立 由运行结果可知，得以验证。
中南大学材料院matlab操作题集答案
的491操作题集操作题1 《Matlab与科学计算》P28 例2.1要求计算水在温度为0℃，20℃，40℃，60℃，80℃时的黏度，已知水的黏度随/(1+at+bt2)温度的变化公式为μ=μ为0℃水的黏度，值为1.758*10-3。
其中μ要求分别以长短格式显示计算结果。
操作题2 《Matlab与科学计算》P31 例2.2数值数组【1：5】和字符串转换示例。
操作题3《Matlab与科学计算》P44 例2.9比较用左除和右除法分别求解恰定方程的解。
操作题4 《Matlab与科学计算》 P48 例2.14计算矩阵magic(3)的指数，并比较不同函数的结果。
指数运算expm(a)操作题5 《Matlab与科学计算》 P50 例2.18计算下面矩阵[-149 -50 -154；537 180 546；-27 -9 -25]的特征值条件数。
特征值函数condeig(A)操作题6《Matlab与科学计算》P62 例2.29对矩阵a = 1 1 1 11 2 3 41 3 6 101 4 10 20 进行三角抽取使之分别变成0 0 0 0 1 1 1 01 0 0 0 及 123 41 3 0 0 1 3 6 101 4 10 0 1 4 10 20操作题7《Matlab与科学计算》P62 例2.30用3种方法建立向量v 为[1 2 6 20]多项式的伴随矩阵。
操作题8《Matlab与科学计算》P64 例2.31已知a=[2 1 -3 -1；3 1 0 7；-1 2 4 -2；1 0 -1 5]；求该矩阵的3次方及数组的3次方。
操作题9《Matlab与科学计算》P66 例2.32 数组逻辑运算演示。
已知a=[1:3；4：6；7：9]，b=[0 1 0；1 0 1；0 0 1]， x=5, y=ones （3）*5,求x<=a，ab=a&b，～b操作题10《Matlab与科学计算》P71 例2.41对多项式p =[2 -5 6 -1 9]求根，计算其在x=3的值，进行微分计算并写出其表达式。
Matlab编程与应用习题和一些参考答案
n=rand(1,10) for j=1:10
for i=j+1:10 if n(j)<n(i) a=n(j); n(j)=n(i); n(i)=a; end
end end disp('重新排列后：') n
9.建立一个字符串向量，然后对该向量做如下处理： (1) 取第 1～5 个字符组成的子字符串。 (2) 将字符串倒过来重新排列。 (3) 统计字符串中小写字母的个数。
>> a='abcdef123456'; >> subch=a(1:5) subch = abcde >> revch=a(end:-1:1) revch = 654321fedcba>> k=find(a>='a'&a<='z'); >> a(k)=a(k)-('a'-'A'); >> char(a)
>> z=0:0.01:10;
>> x=z.*sin(3*z);
>> y=z.*cos(3*z);
>> plot3(x,y,z)
12.设 z x e2 (x2 y2 ) ，画出定义域 x=[-2,2],y=[-2,2]内的曲面图。
>> [x,y]=meshgrid([-2:0.1:2]); >> z=x.^2.*exp(-x.^2-y.^2); >> mesh(x,y,z) 13.设 z=0.05x-0.05y+0.1; 画出 z 的曲面(平面)图。 >> [X,Y]=meshgrid(1:10:200,1:10:200); >> Z=0.05*X-0.05*Y+0.1; >> surf(X,Y,Z)
MATLAB复习提纲和练习题
信号的表示和运算 向量法和符号表达式 移位，翻转，累加，求和，卷积，相关 系统模型及转换 传递函数模型，零极点增益模型，状态空间模型
LTI系统的时域分析（连续、离散）
单位冲激响应、阶跃响应、对激励信号的零状态响应 单位脉冲响应、对激励信号的零状态响应 LTI系统的频域分析（连续、离散） 频响函数（幅度响应、相位响应）
36. 绘制二维曲线的函数是 函数是 ;
；绘制离散序列图的
37. 可以使用
38. 使用函数
命令为图形添加网格;
可以调整坐标轴范围；
39. 使用函数
和
添加坐标标签；
；
40. 为图形添加标题的函数是
41. 函数plot(x,y,’option’)中的参数option可以为曲线设置 、 和 ； 42. 为了使两个plot的图形在同一个坐标显示，可以使用 命令进行图形保持；
；
56. 对于离散LTI系统，可以通过函数 impz( )求其单位脉冲响 应；可以通过函数 filter( ) 求对激励信号的零状态响应
57.求系统y( n) 0.25 y( n 1) 0.5 y( n 2) x( n) x( n 1) 的频率特性，绘出其幅 度特性和相位特性
32.用round函数四舍五入对数组[2.48 6.39 3.93 8.52]取整， 结果为 . (A) [2 6 3 8] (B) [2 6 4 8] (C) [2 6 4 9] (D) [3 7 4 9]
33. 已知a=2:2:8, b=2:5，下面的运算表达式中，出错的为 (A) a'*b (B) a .*b (C) a*b (D) a-b
53. 对于连续系统，可以通过函数 impulse( ) 和 step( ) 求其 单位冲激响应和单位阶跃响应；可以通过函数 lsim( ) 求对 激励信号的零状态响应
中南大学matlab课后习题
一行元素，并给出出错信息。
clear
a=rand(5,6);
n=input('请输入一个数:');
if
n>5
n=5
a(n,:) disp('输入超出范围'); else a(n,:) end 请输入一个数:4 ans =
请输入一个数:7 n=
5
ans =
输入超出范围 7.已知 (1)当 f(n)=n+10ln(n^2+5)时，y 的值是多少。 function a=n2a(n) a=n+10*log(n*n+5) y=n2a(40)/(n2a(30)+n2a(20)) a=
23 -10
12 37 26 -15 73 7 a^2-b+eye(3)
ans = -18 -217 17 22 533 109 21 867 526
（2）a*b
ans = 14 14 16 -10 51 21 125 328 180
a.*b
ans = -8 15 4 0 35 24 -9 122 0
表一
特殊角的正弦与正切值表
α（度） 0
sinα
0
tanα
0
15
30
N
1
4
9
16
√N 1
2
3
4
SumA = C=
2，按要求对指定函数进行插值和拟合。
（1） 按表一所示用三次样条方法插值计算 0~90°内整数点的正弦值和 0~75°内整数点的正切
值，然后用 5 次多项式拟合方法计算相同的函数值，并将两种计算结果进行比较。
（2） 按表二所示用 3 次多项式方法插值计算 1~100 之间整数的平方根。
MATLAB系统仿真试题复习资料,有答案及MATLAB课程设计实验体会
PART I 教材复习要点第一章MATLAB入门测试1.11.MATLAB的命令窗口的作用是什么？编辑/调试窗口？图象窗口？MATLAB命令窗口是MATLAB起动时第一个看到的窗口，用户可以在命令窗口提示符">>"后面输入命令，这些命令会被立即执行。
命令窗口也可以用来执行M文件。
编辑/调试窗口是用来新建，修改或调试M文件的。
图像窗口用来显示MATLAB的图形输出。
2.列出几种不同的得到MATLAB帮助的方法。
（1）在命令窗口中输入help <command_name>，本命令将会在命令窗口中显示关于些命令的有信息。
（2）在命令窗口中输入lookfor <keyword>，本命令将会在命令窗口中显示所有的在第一注释行中包含该关键字的命令和函数。
（3）通过在命令窗口输入helpwin或helpdesk启动帮助浏览器，或者是在启动板中选择"Help"。
帮助浏览器包含了基于超文本的MATLAB所有特性的描述，HTML或PDF格式的在线手册，这是MATLAB 最全面的帮助资源。
3.什么是工作区？在同一工作区内，你如何决定它里面存储了什么？工作区是命令、M文件或函数执行时被MATLAB使用的变量或数组的收集器，所有命令都在命令窗口（所有的脚本文件也是从命令窗口执行）共享公共工作区，因此它们也共享所有变量，工作区的内容可以通过whos命令来查看，或者通过工作区浏览器来图形化地查看。
4.你怎样清空MATLAB 工作区内的内容？要清除工作区的内容，只需在命令窗口中输入clear或clear variables即可。
第二章MATLAB基础测试2.11.数组,矩阵,向量有什么区别?数组是在内存中被组织成行和列的数据集合，只有一个名称，数据要通过在数组名后面圆括号里加上表示数据所在行和列的数字来访问。
术语"向量"通常用来描述只有一维的数组，而"矩阵"通常用来描述二维或更多维的数组。
Matlab期末复习题(答案)
M a t l a b期末复习题(答案)-CAL-FENGHAI.-(YICAI)-Company One1一、 选择题（有10个题，共20分）1．有一个4×4矩阵A,表示取得矩阵A 的第一行和第二行的第二列至第四列所有元素的命令是（ A ）(A) A(1:2,2:4) (B) A(2:4,1:2) (C) A(1:2:4) (D) A(4:2:1)2．将一个图形窗口分成3×4个子区并表示第一个子区的命令是（ C ）A ． subplot(12,1)B ．subplot(1,12)C ．subplot(3,4,1)D ．subplot(1,3,4)3．利用Matlab 的符号计算功能计算微分方程初值问题⎪⎩⎪⎨⎧=+-=2)1(42y x x y dx dy 在区间 [1，2]上的解析解时所用到的函数是（ D ）(A) ode23 (B) ode45 (C) solve (D) dsolve4．以下有关Matlab 子函数的性质中不正确的是（ D ）(A) 子函数只能被处于同一个文件的主函数或其他子函数调用；(B) 在M 函数文件中，任何指令通过“名字”对函数进行调用时，子函数的优先级仅次于内装函数；(C) 在M 函数文件内，主函数的位置不可改变，但子函数的排列次序可以任意改变；(D) help, lookfor 等帮助指令都能提供关于子函数的任何帮助信息。
5．以下语句中生成按钮控件的语句是（ B ）(A) h_1=uimenu(gcf,'Label','&Blue');(B) h_1=uicontrol(gcf,'style','push', 'string','grid off','callback','grid off');(C) h_1=uicontrol(gcf,'style','text', 'horizontal','left','string',{'输入'});(D) h_1=axes('unit','normalized','position',[0,0,1,1],'visible','off');6．以下有关Matlab 函数及其句柄的说法中正确的是（ A,C ）(A) 在该函数的目录下可定义该函数的句柄，句柄定义后，可用feval 函数和句柄在任何目录下调用该函数；(B) 函数不在Matlab 的当前工作目录下不可执行，但可以用feval 函数和句柄调用该函数；(C) 函数不在Matlab 的当前工作目录下时，可以定义函数的句柄；(D) 函数不在Matlab 的当前工作目录下时，不可以定义函数的句柄。
MATLAB期末复习资料（终极版）
MATLAB期末复习资料（终极版）一、填空题：(共20分)1、产生4阶全0方阵的命令为 zeros(4)；产生3阶全1方阵的命令为ones(3).2、标点符号；可以使命令行不显示运算结果; % 用来表示该行为注释行。
3、P,Q分别是两个多项式的系数向量，求P对应的多项式的积分(对应的常数为K)，使用的命令是polyint(P,K)；求P/Q的解，商和余数分别保存在k和r中，使用的命令是[k,r]= deconv (P,Q)4、在C盘上创建目录my_dir的命令是mkdir(…c:\?,?my-dir?)；使c:\my_dir成为当前目录的命令是cd c:\my -dir。
5、已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出各指令运行的结果。
A&B的运行结果ans= 0 1 11 1 0A./B的运行结果ans= 0 2.2500 2.00001.0000 0.6000 NaN6、要判断else是否是MATLAB关键字的命令是exist User Name；exist else的运行结果非零，说明else是（填“是”或“不是”) MATLAB的关键字。
7、A=rand(2,5);b=size(A);c=length(A);则b和c的值分别为[ 2 ,5 ] 和 5 。
8、为了使两个plot的图形在同一个坐标显示，可以使用hold on 命令进行图形保持；可以使用grid on 命令为图形添加网格。
9、MATLAB的工作空间中有三个变量v1,v2,v3，写出把它们保存到文件mydata.mat中的指令save mydata.mat；写出把mydata.mat文件中的变量读取到MATLAB工作空间内的指令load mydata.mat。
10、用if语句判断80≥x>60，在MATLAB中if语言后的判断应写为 if x<= 80 & x>6011、写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值。
2024 - Matlab - 复习资料
一、简述下列吩咐在MATLAB中的作用1.clf 清除图对象2.clear 清除工作空间内的全部变量3.clc 清除当前屏幕上显示的全部内容，但不清除工作空间中的数据4.ceil 沿+∞方向取整。
5.factor 符号计算的因式分解。
6.pow2 2的幂。
7.digits(d) 设置当前的可变算术精度的位数为整数d 位。
8.vpa(pi,100) 设置当前pi的小数位数为1009.box on 打开框状坐标轴开10.grid off网格关11.logspace 创建对数刻度分度等间隔行向量12.linspace 创建等间隔行向量13.cat 串接成高维数组14.sym2poly 符号多项式转变为双精度多项式系数向量15.poly2sym 双精度多项式系数转变为向量符号多项式16.plot3 三维线图17.poly2str 以习惯方式显示多项式18.bar 二维直方图19.pie 二维饼图20.zoom on 打开图形缩放模式21.edit M文件编辑22.whos 对当前工作空间变量的信息进行列表23.figure 生成图形窗口24.cart2sph 直角坐标变为球坐标25.pol2cart 极或柱坐标变为直角坐标26.diff 数值差分、符号微分27.dsolve 符号计算解微分方程28.ezplot3 画三维曲线的简捷指令29.fix 向零取整30.factor 符号计算的因式分解二、填空题1、标点符号; 可以使吩咐行不显示运算结果，% 用来表示该行为注释行。
2、x为0～4pi，步长为0.1pi的向量，运用吩咐x=0:0.1*pi:4*pi 创建。
3、输入矩阵A=，运用全下标方式用A(2,2) 取出元素“-5”，运用单下标方式用A(5) 取出元素“-5”。
4、符号表达式sin(2*a+t)+m中主变量为t 。
5、M脚本文件和M函数文件的主要区分是M脚本文件没有函数定义和M函数文件有函数定义。
matlab课后习题与答案
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：不相等。
(2)A./B和B.\A的值是否相等？答：相等。
(3)A/B和B\A的值是否相等？答：不相等。
(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。
答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素。
答：A(7)=[](3)将矩阵A的每个元素值加30。
答：A=A+30;(4)求矩阵A的大小和维数。
答：size(A);ndims(A);(5)将向量t的0元素用机器零来代替。
答：t(find(t==0))=eps;(6)将含有12个元素的向量x转换成34矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8) 求一个ASCII 码所对应的字符。
答：char(49);4. 下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]5. 已知23100.7780414565532503269.5454 3.14A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦完成下列操作：(1) 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角32⨯子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。
中南大学材料学院MATLAB题库答案3.0版
P.28 ex 2.1 数值变量的运算>> format short e>> clear>> muw0=1.785e-3;>> a=0.03368;>> b=0.000221;>> t=0:20:80;>> muw=muw0./(1+a*t+b*t.^2)muw =1.7850e-003 1.0131e-003 6.6092e-004 4.6772e-004 3.4940e-004p.31 ex 2.2 数值驻足和字符串转换>> a=[1:5];>> b=num2str(a);>> aa =1 2 3 4 5>> bb =1 2 3 4 5>> b*2ans =98 64 64 100 64 64 102 64 64 104 64 64 106P.44 ex 2.9 比较用左除和右除分别求解恰定方程（线性方程组如果方程数等于未知数个数，叫做恰定方程组，如果方程多于未知数，叫做超定方程组，反之称为欠定。
换个角度说，系数矩阵如果是方阵，就是恰定方程组）的解见课本P.48 ex 2.14 计算矩阵的指数并比较不同函数的结果>> b=magic(3);>> expm(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897 >> expmdemo2(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897>> expm1(b)ans =1.0e+003 *2.9800 0.0017 0.40240.0191 0.1474 1.09560.0536 8.1021 0.0064>> expmdemo3(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897P50 ex 2.18 计算矩阵的特征值条件数>> a=rand(3)a =0.9649 0.9572 0.14190.1576 0.4854 0.42180.9706 0.8003 0.9157>> [V,D,s]=condeig(a)V =0.4913 0.6696 0.66960.3158 -0.4476 + 0.2831i -0.4476 - 0.2831i0.8117 -0.2332 - 0.4655i -0.2332 + 0.4655iD =1.8146 0 00 0.2757 + 0.3061i 00 0 0.2757 - 0.3061is =1.17921.23661.2366P62 ex 2.29矩阵的抽取、三角抽取>> a=pascal(4)a =1 1 1 11 2 3 41 3 6 101 4 10 20>> diag(a,-2)ans =14>> v=diag(diag(a))v =1 0 0 00 2 0 00 0 6 00 0 0 20%diag简单来说就是抽取矩阵各对角线上的元素，如上是抽取主对角线以下第二条对角线之元素，其另一功能是建立对角矩阵>> tril(a)ans =1 0 0 01 2 0 01 3 6 01 4 10 20>> triu(a)ans =1 1 1 10 2 3 40 0 6 100 0 0 20%triu&tril用法与diag非常类似，用途是提取下、上三角矩阵P62 ex2.30 建立多项式之伴随矩阵这道题有点凌乱了…….求解释P64 ex 2.31 数组的幂运算>> a=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5];>> a^3ans =32 -28 -101 3499 -12 -151 239-1 49 93 851 -17 -98 139>> a.^3ans =8 1 -27 -127 1 0 343-1 8 64 -81 0 -1 125P66 ex 2.32 数组之逻辑运算>> a=[1:3;4:6;7:9];>> b=[0 1 0;1 0 1;0 0 1];>> x=5;y=ones(3)*5;>> x<=aans =0 0 00 1 11 1 1%此处小于等于，对a中元素和x注意比较，大者为假为0，不大者为1>> ab=a&bab =0 1 01 0 10 0 1%此处为与运算，就是同真才为真（同为非零数）>> ~bans =1 0 10 1 01 1 0%逻辑非运算，即全都非，真变假假变真；还有逻辑或运算，看下面即懂：>> a|bans =1 1 11 1 11 1 1%总结多项式运算的函数：poly:Polynomial with specified roots特征多项式的生成p=poly(a)a 为n阶特征矩阵，所得一般为n阶特征多项式；poly2sym数值2符号；polyval 求多项式的值；roots 求多项式的根；conv 多项式的乘法(向量之卷积)conv(p,d);polyder 多项式微分；polyfit 多项式拟合。
中南大学matlab总复习
主讲：唐建国 中南大学材料科学与工程学院 2013.11
第一章 绪 论
1 科学计算 数值计算的定义 数值计算的研究过程 数值计算的误差分析 数值计算的稳定性与收敛性 2 MATLAB简介 MATLAB的类别 MATLAB的概述 3 课程的要求与学习方法 课程的内容 课程的要求 课程的学习与考核 教材与主要参考书 小结
如果结果为整数，则显示没有小数；如果结果不 是整数，则输出形式有： format (short)：短格式（5位定点数）99.1253 format long：长格式（15位定点数 ） 99.12345678900000 format short e：短格式e方式 9.9123e+001 format long e：长格式e方式 9.912345678900000e+001 format bank：2位十进制 99.12 format hex：十六进制格式
（3）正交分解 [Q,R]=qr(A) 将矩阵A做正交化分解，使得 例： a = Q*R=A，其中Q为正交矩阵（其 9 8 范数为1，指令norm(Q)=1)，R为 6 8 对角化的上三角矩阵。 [q,r]=qr(a) q= -0.8321 -0.5547 -0.5547 0.8321 r= -10.8167 -11.0940 0 2.2188
函数3： ezplot：绘制隐函数图形
1(f, [xmin,xmax,ymin,ymax])
说明: 在区间 xmin<x<xmax和ymin<y<ymax绘制 f(x,y) = 0的图 形。,默认区间-2π<x<2π和-2π<y<2π
MATLAB与科学计算期末复习题题库(第一部分).doc
cos(2arccos arccos-)6MATLAB 与科学计算期末复习题题库(第一部分)一、填空1. MATLAB 的主界面是一个高度集成的工作环境，有四个不同职责分工的窗口，分别 为 ________ 、 _________ 、 _________ 、和 _________ 窗口.2. MATLAB 计算屮，欲求空啤2的值，应在命令窗口命令提示符后输入 ___________ ,欲求1+V53空(0•辺的値，可在命令窗口命令提示符后输入计算 log 5 2475 二 ______________________ , ln(l + e'2) = ___________________ , I 23. MATLAB 软件屮，要输入矩阵，应在命令窗口命令提示符后输入 ______________3 4 4. MATLAB 软件屮，生成2x3阶的全零、全一、单位矩阵时，应分别输入 、 ___ 、 ・5 . MATLAB 软件屮，用于求可逆矩阵A 的行列式和逆矩阵的函数分别为 _______ 、 _________ •6. 语句 A=linspace(2/18/9), B=reshape(A,3/3)的执行结果为 B= ____________ .7. MATLAB 中,用命令polyval 计算多项式f(x) = 4x 3+ 2x 2 + x +1在兀= -1,1, 2时的值, 可在命令窗口中输入 _______ .MATLAB 计算屮，欲用命令polyval 计算f (x) = 4x 3 +6x 2 + 1在自变最X 取小于20的正 整数时的函数值，可在命令窗口屮输入 _______ .(提示：先用help polyval 在命令窗口中在线查询函数polyval 的用法)用roots 求/(兀)=5/ +6兀+ 1的根的语句为 _________ •8. 在MATLAB 操作屮，把变量x, y 定义或申明为符号变量的语旬为 ________ •9. MATLAB 语言中，合法的变量名须以 _______ 开头，后可跟 _____ 、 ______ 、 ______10. MATLAB 计算屮，语句 A=[l 2 3]； b=[-l 2 6]； c 二dot(a,b)运行结果为 c= _________ ・(提示：先用help dot 在命令窗口屮在线查询函数dot 的用法)11. 欲将MATLAB 的数值显示格式设置为近似有理数格式，可在命令提不符后输入 ______12. 欲将MATLAB 的数值显示格式设置为长格式，可在命令提示符后输入 _______ ・13. 欲将MATLAB的数值显示格式设置为短格式，可在命令提示符后输入_______ •14・MATLAB符号计算屮，用于求极限、导数及积分的函数分别为___ 、______ 、___________________________________________________________ ・15. MATLAB软件屮，对函数做泰勒展开式的函数为__________ .16. MATLAB的函数文件的第一个可执行语句必须以关键词_______ 开头.17. 欲退出MATLAB环境，可在命令窗口中输入_______ 或 ______ 命令.(提示：实验一下exit、quit函数)18. MATLAB的命令窗口屮的语句»clear x y的作用为_______________ .19. 建立下列数据序列，并对所有序列元素求和/ 1 1 1(1)—, —9°・・,2 3 100⑵斗,21,2,…,201 + z(3) (2n-l)\n = l,2,---,2020 • MATLAB符号计算中，用于求解代数方程和常微分方程的函数分别为、・21. MATLAB软件屮，语句int(f,a,b)表示____________ ・22. 用MATLAB语言编写的m文件有两类：_______ 文件和 ______ 文件.23. 擦除MATLAB命令窗口的所有内容用__________ 命令，清除内存屮的变量和函数用命令，清除MATLAB当前窗口屮的图形用命令，检杳指定名字的变量或函数文件的存在性用命令24. MATLAB绘图时，，在同一个图形窗口屮欲使前途不被麻图擦除，可在前示绘图语句间加上命令_________ .25. MATLAB软件有一些常用的系统预定义的变量，如无穷大、圆周率龙、不定量分别用表26. 已知A=[l 1;2 4];B=[1 1;3 4];运行 A.*B Z A*B , diag(A.*B)的结果分别为 ________ 、 _________ 、 _________ -27. A=zeros(2/4)/A(:)=l:8;s=[2 3 7];贝lj A(s)= __________ .28. 在循环结构中跳出循环，执行循环后面的代码的命令为 ________ .29. 下面程序执行后array 的值为 __________ .for k=l:10if k>6break;elseArray (k)=k;endend30. A=[1 2 3;3 4 5;5 6 7];删除矩阵A 的第三行的语旬为 __________ •31. A=[1 2 3;3 4 5;5 6刀；删除矩阵A 的第三行的语句为 ________ •32. A=[l 2 3;3 4 5;5 6 7]; C=[A,A+1]; C(2,3)= _____ 、C(2,2:3)= _______ 、C([l ?3],[l ：3])= ______ .33. _______________________________________________________________ 求矩阵A=[l 2 3;2 1 3;1 2 3]的特征值和特征向量，特征值为 ____________________________ ,对应的特征向量 为 ________ .(提示：函数eig)34. ________________________________________________ 用gcdjcm 计算135, 279和468的最大公约数 ________________________________________ 和最小公倍数 ______ 。
Matlab 课后习题集
Matlab程序设计课程实验目录Matlab操作实验（第1章）实验指导思考练习矩阵运算实验（第2章）实验指导思考练习程序设计实验（第3章）实验指导思考练习高级绘图实验（第4章）实验指导思考练习数据分析与多项式计算实验（第5章）实验指导思考练习解方程与最优化问题求解实验（第6章）实验指导思考练习数值运算实验（第7章）实验指导思考练习符号实验（第8章）实验指导思考练习图形句柄实验（第9章）实验指导思考练习图形用户界面设计实验（第10章）实验指导思考练习Simulink动态仿真实验（第11章）实验指导思考练习外部接口技术实验（第12章）实验指导思考练习1、Matlab操作实验一、实验内容1. 第一章实验内容1：path(path,’系统盘名:\matlab文件名\要添加的文件名’)如：path(path,’c:\matlab701\my work’)实验内容-2例1-1：x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));例1-2p=[3,7,9,0,-23];x=roots(p)例1-3quad(‘x.*log(1+x)’,0,1)例1-4：实验内容3Help 函数名如：help inv help plot等实验内容4（1）x=0:pi/10:2*pi;y=sin(x);二、思考练习思考题1启动：有三种常见方法，（1）、在Windows桌面，单机任务栏上的‚开始‛按钮，选择‚所有程序‛菜单项，然后选择MATLAB程序组中的MATLAB （版本号）程序选项，即可启动。
（2）、在MATLAB的安装路径中找到MATLAB系统启动程序matlab.exe，然后运行它。
（3）、利用建立快捷方式功能，将MATLAB系统启动程序以快捷方式的形式放在Windows桌面上，以后只要在桌面上双击该图标即可启动MATLAB。
关闭：常见的也有3种方法，（1）、在MATLAB主窗口file菜单中选择exit MATLAB命令。
Matlab课程复习大纲含答案.doc
计算机仿真技术-Matlab应用复习大纲1、产生一个3行4列几在(0,1)Z间均匀分布的随机矩阵A,其命令为_rand(3,4) __________ ;2、将一个5行6列的A矩阵中的第4列所有元素提取岀来赋值给B,其命令为_B=A(: ,4) ___________ ；3、将一个5行6列的A矩阵中的第5行所有元素提取出来赋值给B,其命令为_B=A(5 ,:) _____________ ；4、将一个5行6列的A矩阵屮的第2行到第4行，第3列和第5列所有元索提取出來赋值个B,其命令为_B=A(2:4,3:2:5) ___________ ；5、用一条命令将一个5行6列的A矩阵中的第2行，第4行，第5行所冇元素提取出来赋值个B,其命令为_B=A ([2,4,5],:) _____________ ；6、用一条命令将一个5行6列的A矩阵中的第2列，第4列，第5列所有元素提取出來赋值给B,具命令为B=A(:,[2,4,5]) ____________ ;7、用一条命令将一个5行6列的A矩阵中的第8个元索提取出來赋值给D,其命令为—B=A(8) _________ ；8、用一条命令将一个5行6列的A矩阵中的第8个到第15个元素提取出来赋值给D,其命令为_B=A(8:15) _________ ；9、产生3x4的均匀分布随机数矩阵R,要求其元素在2到14之间取幣数，实现其结果的命令是」'andint(3,4,[2 14]) _________________ 。
9、产生一个5行6列的单位矩阵C,其命令为_C=eye(5,6) ____________ :10、___________________________________________________________提取矩阵A的对角线元素赋值给B,其命令为_B=diag(A) ____________________________________ ；已知A=[5 4 1;20 1;3 5 1],贝ij diag(A r l)= ________ [2 5] _______ , triu(A,l)=[0 4 1;0 0 1;0 00] ___________ ,tril(A,l)=[5 4 0;2 0 1;3 5 1] __________11、已知矩阵A=[2 3 4;5 1 2;2 0 1],则A(3)=_2_, A(6)=—0_。
MATLAB复习题(含答案)
MATLAB语言复习大纲1、掌握以下命令行编辑常用的控制键的用法：Del，Esc，Backspace，Ctrl+C。
答：Del是删除光标右边的字符Esc是删除当前行的全部内容Backspace是删除光标左边的字符Ctrl+c是中断一个matlab任务2、在Matlab中，行注释是以什么符号开头的?答：注释是以%开头，后面是注释的内容3、合法的变量名的命名规则是什么？答：是以字母开头，后接字母、数字或下划线的字符序列，最多允许有63个字符4、在Matlab中，矩阵行与行之间，同行元素之间以什么符号分隔开来？若a是一个m*n的矩阵，则求解矩阵a中最小的元素语句是什么？min(min(a))答：同一行的个元素之间用空格或逗号分隔，不同行的元素用分号分隔5、建立一个均值为3，方差为2的10*10的正态分布随机矩阵命令是什么？答：X=3+sqrt（2）*randn(10)6、产生单位矩阵，全零矩阵，全一矩阵的函数是什么？答：eye ones zeroseye(3)ans =1 0 00 1 00 0 1>> ones(3)ans =1 1 11 1 11 1 1>> zeros(3)ans =0 0 00 0 00 0 0>>7、MATLAB的矩阵有哪两种存储方式？函数文件与命令文件的区别是什么？答：完全存储方式和稀疏存储方式区别：56页8、命令plot(x)，当x为复数向量时，如何画曲线？掌握plotyy函数的用法，它与plot函数有什么区别？P86答：当x为复数向量时，则分别以向量元素实部和虚部为横、纵坐标绘制一条曲线。
9、最常见绘制网线图函数和最常用绘制曲面图函数分别是什么？如何使用？答：mesh（x,y,a,c）,surf(x,y,z,c)10、MATLAB提供哪些函数进行动画制作？答：111页；getframe moviein（n）movie（m,n）11、掌握view函数的使用方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
a）-1/6 b) Inf
（a）
c) –Inf
d）-1
2018/6/9 32
2
（c）
a）-1 b) 1/2
c) 0
d）1
2018/6/9 33
七、符号求导
2018/6/9
34
1 y=
2x t2 x t sin( x ) e
,求x=3，t=2时的dy/dx值
（a）
a)[2.0,0;-1.97998,20.0855] b)[2.0,0;-1.98898,20.0865]
（b）
a) [21 75 54; 4 19 3; 22 72 72]
b) [0 36 18; 1 15 0; 6 25 63] c) [4 19 3;21 75 54; 22 72 72] d) [6 25 63;0 36 18; 1 15 0]
2018/6/9 23
2、求数组b=[1 2 3;4 5 6;7 8 9]的指数运算
程序分析：默认的自由变量是t，C1、C2为任意常数，程序也可指定自由变量， 结果相同： >> [x,y]=dsolve('Dx=y,Dy=x','t')
2018/6/9 9
非线性方程的符号注解
语法： fsolve(‘fun’, x0)
（1）fun为所要求解的函数名，常以M文件形式给出； （2）X0为初始向量或矩阵
2018/6/9 24
3 求矩阵
的特征值
b) [2;3;4] d) [1;3;5]
（a）
a) [1;2;3] c) [3;4;5]
2018/6/9
25
四、数组抽取
2018/6/9
26
四、数组抽取
1、若A为一二维数组，要删除其第34行的 元素，可使用命令 （c）
a)A(34,:)=0 b) A(:,34)=[] c) A(34,:)=[] d) A(:,34)=0
2018/6/9 30
2、 已知A=magic(5),A(:,3)=zeros(5,1)， 则all(A(:,1)<10)
（a）
a)0 b)1 c) [0 0 1 0 0] d)[1 1 0 1 1]
2018/6/9
31
六、符号求极限
1. 求极限
x ln e x 3 3 2 2 lim x x x 1 x x 1 x x
2018/6/9
24 34 36 30
54 84 54 60
36
八、符号求积
1 符号求积
（c）
a)1/2 b) 1
c) 1/4
d) 1/3
2018/6/9 37
2 符号求积
1 1 x 2 dx
（a）
a)pi
b)2*pi c)3*pi d)1/2*pi
2018/6/9 38
dn y y的n阶导数 n dx dn y 或n dt
%求解微分方程 %求解微分方程组
dy
dy
表示为Dy； 表示为Dny。
（2）微分初始条件'con'应写成'y(a)=b，Dy(c)=d'的格式；当初始条件少于 微分方程数时，在所得解中将出现任意常数符C1、C2……，解中任意常数符的数目 等于所缺少的初始条件数。
2018/6/9
5
2018/6/9
6
九、方程符号求解
2018/6/9
7
符号常微分方程
MATLAB提供了dsolve命令，可以用于对符号常微分方程进行求解。 语法： dsolve('eq', 'con', 'v') dsolve('eq1,eq2…', 'con1,con2…', 'v1,v2…') （1）当y是因变量时，微分方程'eq'的表述规定为： y的一阶导数 dx 或dt
（d）
a) 0.4600 b) 0.8400
1.8601 2.9401
0.9730 0.85608
c) 0.4900
d) 0.4900
1.8601
1.2501
0.9730
0.8560
2018/6/9
21
三、矩阵及数组代数运算
2018/6/9
22
三、矩阵及数组代数运算
1、已知A=[0 9 6;1 3 0;3 5 7];B=[1 4 3;1 5 0;2 5 9], 则A.*B
Fc.m function y=fc(x) y(1)=x(1)-0.7*sin(x(1))0.2*cos(x(2)); y(2)=x(2)0.7*cos(x(1))+0.2*sin(x(2)); y=[y(1) y(2)]
2018/6/9
x0=[0.5 0.5] >> fsolve('fc',x0)
2018/6/9
1 对于圆周率值，若取值为3.1415，则其 有数字为 （b） a、3位 b、4位 c、5位 d、2位
2018/6/9
13
2018/6/9
14
2下面的MATLAB语句中正确的有：
（cd）
a) a=sin(x)
b)Record-1=3+4i
c) a=sin(pi) d) C=1+6j
3符号求积
（b）
a) 3.20808 b) 3.10808 c) 3.30808 d) 3.40808
2018/6/9 39
谢谢！
ans =
0.5265 0.5079
10
九、方程符号求解
1 求微分方程 y(0)=0的解
d2 y dy x 2 3 x2 dx dx
，y(1)=0，
（c）
a)x^4/5 - x^3/4
b)x^3/3 - x^2/3 c) x^4/3 - x^3/3 d) x^4/5- x^3/4
2018/6/9
2018/6/9
8
【例
】 求微分方程
d2 y dy x 2 3 x2 ，y(1)=0，y(0)=0的解。 dx dx
>> y=dsolve('x*D2y3*Dy=x^2','x') %求微分方程的通解 y= C3*x^4 - x^3/3 + C2 >> y=dsolve('x*D2y3*Dy=x^2','y(1)=0,y(0)=0','x') %求微分方程的特解 y= x^4/3 - x^3/3 dx dy y , x 的解。 【例 】 求微分方程组 dt dt >> [x,y]=dsolve('Dx=y,Dy=x') x= C8*cos(t) + C7*sin(t) y= C7*cos(t) - C8*sin(t)
2018/6/9 15
3 符合MATLAB变量命名规则的有
a) var!2015 b)VAR15 c)my-var2015 d)2015_VAR
（b）
2018/6/9
16
4在循环结构中跳出循环，但继续下次循环 的命令为
（c ）
(a) return; (b) break (c) continue (d) keyboad
（a）
a)[192 287;416 585] b)[190 285;406 582]
c)[194 292;422 580]
d)[192 288;420 580]
2018/6/9 20
3 试用一个二次多项式拟合下列数据 x y 0.5 1.75 1.0 2.45 1.5 3.81 2.0 4.80 2.5 8.00 3.0 8.60
11
2、利用非线性方程组的符号法求解如下方程组：
x1 0.7sin x1 0.2cos x2 0 x2 0.7 cos x1 0.2sin x2 0
（a）
在（0.5，0.5）附近的解。 x2=0.5079 x2=0.5265 x2=0.7079 x2=0.6265
12
a) x1=0.5265 b) x1=0.5079 c) x1=0.6265 d) x1=0.7079
c)[2.0,0;-1.96998,20.0655]
d)[2.0,0;-1.98998,20.0855]
2018/6/9
35
七、符号求导
2求函数在指定点的导数值。
x 0 x2 2 x3 6x f x 1 2 x 3x 2 , x 1, 2,3
（a）
a) 6 b) 6 c) 4 d) 4
29
五、数组关系逻辑运算
1、已知A=[0 9 6;1 3 0;3 5 7]; B=[1 4 3;1 5 0;2 5 9], 则A<B
（b）
a)[0 1 1; 1 1 0; 1 1 1] b)[1 0 0; 0 1 0; 0 0 1] c)[1 0 1; 1 0 0; 0 1 1] d)[0 0 1; 0 1 0; 1 0 1]
2018/6/9பைடு நூலகம்
17
二、多项式及其运算
2018/6/9
18
二、多项式及其运算
1 求多项式
x 4 12 x 3 25 x 116
在x=3,8处的值
a)-52 b) -54 c) -52 d) -54
2018/6/9
-1730 -1732 -1732 -1730