MATLAB考试题.docx
《MATLAB考试题.docx》由会员分享,可在线阅读,更多相关《MATLAB考试题.docx(16页珍藏版)》请在冰豆网上搜索。
MATLAB考试题
MATLAB考试题
一、填空题
1、在MATLAB中,清除所有变量的命令为_clear___________。
2、在MATLAB中,保存和载入数据分别使用命令___save_________、____load________。
3、矩阵
的条件数记为___________,定义为___________。
4、在MATLAB中,对图形窗口进行灵活分割的命令为_suplot___________,对函数自适应采样的绘图函数为___________,绘制三维曲面的函数为_surf___________、____________。
5、在MATLAB中,求矩阵特征值及对应的特征向量的命令为__[X,Lamda]=eig(A)_________。
6、
,
=[334],A.*B=[223]。
7、在MATLAB中,zeros(1,3)=,ones(1,2)=。
8、给定
个点可以唯一确定一个
阶多项式,在MATLAB中,利用命令_________可以确定多项式的系数,对多项式求导的命令为____________。
9、在MATLAB中,使用命令____linspace________可做线性插值。
10、在MATLAB中,输入命令:
sum([215;985]),输出结果:
ans=。
11、在MATLAB中,画二维图形的常用函数为_plot()___________,画三维图形常用函数为__plot3()________、surf()___________、_sphere()______________。
12、在MATLAB中,求矩阵
的逆矩阵可用命令_inv(A)_________。
13、
,
=6,A.^2=。
14、在MATLAB中,可以用命令polyval计算多项式的值,用命令roots求多项式的根。
15、将区间
三等分,步长为
,各个节点的函数值为
,则3/8辛普森法的公式为。
1/3辛普森法的公式为。
16、在MATLAB中,数据输出时可以用命令__format_______改变数据输出格式。
17、在MATLAB的工具箱中,求定积分的函数为____________、____________、____________。
18、在MATLAB中,函数文件的基本结构为:
_____function_______输出形参表=函数名(输入形参表)
注释说明部分
函数体语句
19、在MATLAB中,求向量或矩阵范数的函数为。
20、在MATLAB中,eye(3,3)=。
二、选择题
1、在MATLAB中,列出当前工作空间变量名的函数是
(A)who(s)(B)clc(C)path(D)draw
2、在MATLAB中,图形窗口色图的设置和改变需调用函数
(A)color(B)colormap(C)plot(D)change
3、在MATLAB中,求齐次线性方程组
的基础解系可直接调用函数
(A)numl(B)solve(C)null(D)polyfit
4、在MATLAB中,不能画三维图形的函数为
(A)text(B)mesh(C)surf(D)plot3
5、MATLAB中M文件的分类为
(A)脚本文件、函数文件(B)仅有脚本文件
(C)仅有函数文件(D)函数文件、命令文件
6、MATLAB中生成平面网格坐标矩阵的函数为
(A)meshgrid(B)gridmesh(C)subplot(D)sub
7、给定
个点可唯一确定一个
次多项式,需调用的命令
(A)poly(x,y,n)(B)polyfit(x,y,n)(C)polyfit(n)(D)poly(n)
8、MATLAB中‘legend’函数的功能为
(A)给图形加标题(B)在指定位置给图形加注释
(C)给图形变色(D)给图形添加图例
9、已知
,则在MATLAB中调用命令A(3)的结果为
(A)3(B)2(C)无结果(D)7
10、MATLAB中调用命令[L,U,P]=lu(A)后,矩阵A满足
(A)
(B)
(C)
(D)
11、MATLAB中调用命令v=eye(2,1)后的结果为
(A)v=[1;0](B)v=[1,0](C)v=[0;1](D)v=[1,0;0,1]
12、MATLAB中计算积分
,则下列程序中的空格处应填的命令是
symt ;f=3/((t-1)^2*(t-2));F=_____________,y=double(F)
(A)(f,t,2,6)(B)qad(f,t,2,6)(C)int(t,2,6)(D)int(f,t,2,6)
13、如果在MATLAB的命令窗口输入命令a=input('a=?
'),出来的结果为
(A)a=?
(B)没有任何显示(C)a=(D)’a=?
’
14、’try’语句为MATLAB所特有的试探性语句,下面的程序MATLAB能执行吗,执行的命令为A=[1,2,3;4,5,6];B=[7,8,9;10,11,12];tryC=A*B;catchC=A.*B;end
(A)不执行(B)执行的是’catch’后的命令
(C)执行,但没有结果(D)执行的是’try’后的命令
15、MATLAB中与linspace(3,7,11)等价的命令为
(A)3:
7(B)3:
0.4:
7(C)3:
4/11:
7(D)3:
7/10:
7
16、对于下面的程序,在MATLAB中运行后f的形式为
a=0;b=3*pi;n=1000;h=(b-a)/n;x=a:
h:
b;f=exp(-0.5*x).*sin(x+pi/6)
(A)数(B)矩阵(C)向量(D)没有结果
17、已知
,
,MATLAB中若想由x生成矩阵
可调用命令
(A)A=[x,x.^2,x.^3](B)A=[x,x^2,x^3](C)A=[x;x^2;x^3](D)A=[x;x.^2;x.^3]
18、已知数据点为x=-6:
6,y=sin(x),x0=-6 :
0.25 :
6,过
、
做分段三次样条插值,并求该多项式再x0处的函数值,在MATLAB中应使用的命令是
(A)interp1(x,y,x0,’pchip’)(B)interp1(x,y,x0,pchip)
(C)interp1(x,y,x0,’spline’)(D)interp1(x,y,x0,’cubic’)
19、下面命令中不能求解线性方程组
的命令是
(A)A/b(B)A\b(C)inv(A)*b(D)A^(-1)*b
20、MATLAB中命令exp
(1)的结果为
(A)0(B)1(C)e(D)inf
21、MATLAB中命令sum([412;036])的结果为
(A)16(B)448(C)79(D)程序出错,没有结果
22、MATLAB中生成希尔伯特矩阵的函数为
(A)invhilb(n)(B)company(n)(C)eye(n)(D)hilb(n)
23、MATLAB中给变量w赋空矩阵的语句为
(A)w=zeros()(B)w=[](C)w=zeros(0,0)(D)clearw
24、在MATLAB中,如需查询一个函数的基本用法,命令为
(A)help(B)version(C)demo(D)path
25、A=[121;223;500],则MATLAB中求A的逆矩阵的命令为
(A)det(A)(B)uv(A)(C)A*(-1)(D)inv(A)
26、b=[5,23,7,2,3,6,0],要把b向量的元素按从小到大的顺序排列,MATLAB中用的命令为
(A)sort(B)put(C)fprintf(D)mesh
27、已知A=[40],MATLAB中C=[A,ones(size(A))]为
(A)4011(B)4010(C)4000(D)4040
28、在MATLAB中,列出当前工作空间变量名的函数是
(A)path(B)clc(C)who(s)(D)draw
29、计算矩阵、向量范数的命令为
(A)length(B)norm(C)size(D)anorm
30、MATLAB对图形窗口进行灵活分割的命令为
(A)plot(B)subplot(C)aplot(D)ploval
31、MATLAB中生成一个2行4列的单位矩阵的命令为
(A)eye(2,4)(B)eye(2,2)(C)eye
(2)(D)zeros
32、已知x=[4,2,],求以x为根的多项式的语句是
(A)x=[4,2,];polval(x)(B)x=[4,2,];poly(x)
(C)p=roots(x);poly(p)(D)x=[4,2,];roots(x)
33、在MATLAB中,清除所有变量的命令为
(A)clear(B)clc(C)delete(D)holdoff
34、命令floor(3.4)及ceil(3.4)的值为多少
(A)34(B)44(C)33(D)43
35、在MATLAB中,打开示例程序的命令为
(A)demo(B)help(C)clear(D)plot
36、MATLAB中保持图形窗口不关闭的命令为
(A)alter(B)holdon(C)holdoff(D)cd
37、在MATLAB中,强行终止程序运行的命令为
(A)ctrl+v(B)ctrl+shift(C)ctrl+c(D)shift+c
38、
,
,则A.*B=
(A)运算不执行,因有错(B)683(C)574(D)663
39、双纵坐标函数为
(A)plot(B)yplotty(C)plot3(D)plotty
40、已知数据点为x=[1,1.4,1.6,1.8,2.0]',y=[3.75,3.99,2.07,1.52,4.22]',拟合一条过这些点的3次曲线,应使用的命令是
(A)polyfit(x,y)(B)polyfit
(2)(C)polyfit(x,y,1)(D)polyfit(x,y,3)
41、MATLAB中表示无穷的符号是
(A)
(B)int(C)inf(D)+in
42、MATLAB中给图形加网格线的命令为
(A)gridoff(B)holdon(C)boxon(D)gridon
43、MATLAB中程序A=[123;-121;352];B=[346;589;453];C=(B>1)&(B<7),运行后结果为
(A)C=[1,1,1;1,0,0;1,1,1](B)C=[1,0,0;1,0,0;1,0,1]
(C)无结果(D)C=[0,1,0;1,0,0;1,1,1]
44、MATLAB中生成一个从0到20的等差为0.5的数列应使用的命令为
(A)x=linspace(0,20)(B)x=0:
20(C)x=linspace(0,20,41)(D)x=linspace(0,20,0.5)
45、积分区间被等分为奇数个小区间时,合理使用1/3辛普森法和3/8辛普森法的是
(A)只用3/8辛普森法
(B)1/3辛普森法和3/8辛普森法都不能使用,用cotes方法
(C)前三个区间用3/8辛普森法,剩下的区间用1/3辛普森法(D)只用1/3辛普森法
46、MATLAB中c=[1234],要得到c的转置向量所用命令为
(A)b=c’(B)b=[1,2,3,4](C)b=inv(c)(D)b==c
47、MATLAB中计算积分
,则下列程序中的空格处应填的命令是
symt ;f=5/((t-1)*(t-2)*(t-3));F=_____________,y=double(F)
(A)(f,t,4,5)(B)qad(f,t,4,5)(C)int(t,4,5)(D)int(f,t,4,5)
48、MATLAB表示圆周率的符号是
(A)
(B)int(C)pi(D)没有符号可表示
49、MATLAB中a=[1,3,2],则b=a.^2的结果为
(A)b=[1,6,4](B)b=[1,9,4](C)b=[2,6,4](D)执行出错,无结果
50、求解
在x=1附近的根的错误的语句为
(A)fzero('x^3-3*x+1=0',1)(B)fzero('x^3-3*x+1',1)
(C)fzero('x^3-3*x+1=0',[1,2])(D)P=[1,-3,1];roots(p)
三、阅读下述MATLAB脚本程序
A=[3579;0241;6448];
B=ones(3,3);
b=A([12],2:
3);
B(4,4)=5;
A([1,2],:
)=[];
C=[A;B(:
1:
3)'];
D=reshape(C,2,8);
写出b、A、B、C、D的结果。
四、已知多项式
和
,编写脚本程序计算:
1)多项式a(x)+b(x)
2)多项式a(x)*b(x)。
3)多项式a(x)的根。
五、用MATLAB编写函数,计算
,函数的输入参数为x,输出参数为f,函数名为fun。
六、编程求
的值,要求:
(1)使用循环语句编程,但不能使用数组运算符和sum;
(2)使用数组运算符和sum编程,不能使用循环结构。
七、已知超定方程组
写出用Matlab中语言求解上述线性方程组的程序。
八、Hilbert矩阵是非常著名的病态矩阵,写出其定义,并编出Matlab程序求
到
的Hilbert矩阵的条件数、行列式。
九、数值求解定积分时,当积分区间被等分为奇数个小区间时,如何合理使用1/3辛普森法和3/8辛普森法。
十、写出将矩阵
进行LU分解,并求出该矩阵的行列式的语句
x=-5:
5;y=1./(1+x.^2);
t=-5:
0.05:
5;
y1=1./(1+t.^2);
n=length(t);
fori=1:
n
z=t(i);s=0;
for
Lk=1;
u=x(k);
for
if
Lk=Lk*(z-x(j))/(u-x(j));
end
end
;
end
y2(i)=s;
end
plot(x,y,'ko',t,y1,t,y2,'r')
十一、写出
个插值节点
的Lagrange插值多项式的公式,并将下面的Lagrange插值多项式的Matlab实现程序补充完整。
(每空3分,共15分)
Lagrange插值多项式公式:
十二、已知
的表达式
,当
时,求
的值,请补充完整。
n=100;i=1:
n;
f=;y=
十三、
是数据点的横坐标数组,
是拟合数据点
,
的多项式,
是一个矩阵,它的第
行即为
的幂系数,请补充完整。
functionp=shape_pw(x)
np=length(x);
forj=1:
np
y=zeros(1,np);y(j)=1;
p(j,:
)=;
end
十四、下面的程序是解线性方程组的一个函数文件,请补充完整。
function[x,y]=line_solution(A,b)
[m,n]=size(A);y=[];
ifnorm(b)>0
ifrank(A)==rank([A,b])
ifrank(A)==
x=A\b;
else
disp('原方程组有有无穷个解,其齐次方程组的基础解系为y,特解为x');
y=null(A,'r');x=;end
else
disp('方程组的最小二乘法解是:
');
x=A\b;end
else
ifrank(A)==n
x=;
else
disp('方程组有无穷个解,基础解系为x');
x=null(A,'r');end
end
return
十五、下面是求解非线性方程的牛顿法的函数文件,公式为
,
请补充完整。
functionx=Newt_n(f_name,x0)
x=x0;xb=x-999;n=0;h=0.01;
whileabs(x-xb)>0.0001
n=;xb=x;
ifn>300break;end
y=feval(f_name,x);
yd=(feval(f_name,x+h)-y)/;
x=;end
fprintf('\nFinalanswer=%12.6e\n',x);
十六、编写函数文件求半径为
的圆的周长和面积。
十七、求解下列微分方程
设著名的Lorenz方程可表示为
其中,
。
若令其初值为
(1)写出用desolve命令求解该方程的语句,并指出能否用desolve求出解析解
(2)写出用ode45求解该微分方程组的解并分别绘出状态变量的时间响应曲线和相空间三维图的程序